[ Home | Contents | Search | Post | Reply | Next | Previous | Up ]
From: nikolai.zaitsev@cern.ch
Date: 11/10/97
Time: 10:34:40 AM
Remote Name: 192.16.199.49
Remote User: lhcb
Dear Colleagues,
I comitted an example of Bd-Jpsi+Ks reconstruction into task/task2 directory. It is simplifyed version of reconstruction routines I did during preparation of LOI. I hope it will help for beginners to understand how to work with AX-routines and will give idea about HEAP and STACK.
The program creates J/Psi,Ks,Bd as HEAP objects within large window ~50MeV (actual value is taken from sicb.dat, see sudata.F). Mass distribution is put into histograms 1,2,3. CWN-tuple is also filled with some usefull info (see below).
More complex example with detailed analysis is useless, since now everything is different.
I tryed to give good on-line comments and I believe you may start even without more deep explanations. For more questions, look unique, though out of date, manual written by J.Zweizig.
The example consist of following routines (simple flow chart):
include files are in the same directory.
sudata - init of ffkey cards (variables ievdmp,myprint,window can be taken from sicb.dat) suanal - steering of job \make_evhead - inits, fills CWN-block /evhead/ with header of event \make_evreco - inits, fills CWN-block /evreco/ with general information; makes Primary Vertex and stores its coordinates (true positions as well) \make_bkine - inits, fills CWN-block /bkine/ with information about b-particles: oscillations flag, acceptance flag, kinematics, decay position etc. \bjk - steering of Bd-JPsi+Ks reconstruction \didcay - general routine for reconstruction of particle ID out of ID1 and ID2. Inside performs AX-scan of ID1,ID2. \subundle - routine called by didcay for every combination of ID1,ID2. It is upto you (of course) to decide whether to bundle ID or not. Fills common block /trtrak/ with reconstructed information and HEAP-like references. Didcay routine can be easily extended for recosntruction of 3/4-prong decays (I foresaw such possibility in calling parameters of subundle).
Please, look carefully every routine I give (they are not so many), before to complain. For a moment, I can garantee only technical functioning of reconstruction I give, but no of AX-library !!!
Make_* routines are of general use and I plan to give similar routines filling TRIG/TAG information (as soon as I will update my own). Though anybody is welcome to share the experience.
the executable is made with 'gmake PACK=task PATCH=task2 VERSION=111'
sicb.dat = task/sicb_jks.dat
In principle, would be nice to extend this package as general, so aby suggestions, comments are welcome.
Good luck, Nikolai.
------------------------------------------------------------------ work address: | home address: -------------------------------+---------------------------------- NIKHEF, | Paramaribostraat,12,II Kruislaan, 409, | 1058 AN Amsterdam 1009 DB Amsterdam, | The Netherlands The Netherlands | | tel:+31 20 592 2050, | tel:+31 20 689 60 57 fax:+31 20 592 5155 | ------------------------------------------------------------------ e-mail: Nikolai.Zaitsev@cern.ch or n.zaitsev@nikhef.nl
http://wwwcn1.cern.ch/~nzaitsev