Hi Bjorn, On Mon, 27 May 2002 08:37:09 -0400 (EDT) Bjorn H Samset <bjornhs@rcf2.rhic.bnl.gov> wrote concerning "AMPT tree": > fY = -999; > if (fE != fPz) fY = 0.5*TMath::Log((fE + fPz)/(fE - fPz)); Why not fY = fE != fPz ? .5 * TMath::Log((fE + fPz) / (fE - fPz)) : -999; > PS: Since I know this one is coming - no, there is unfortunalely no easy > way of turning this data into something BRAG can read. Not unless someone > has a DumpZdat-routine for ROOT... BRAG does not _need_ the input stored in ZEBRA/HBOOK banks. BRAG can read ASCII files with columns, OSC1997 formated files, and a few others. So, if you can write the chain into something like that, BRAG can read it. The format of the ASCII column files are file : num_particles_in_event '\n' particle_list particle_list : /* empty */ | particle particle_list particle : geant_pid px py pz '\n' These can be read by setting USER/CONTROL/CARDS UKIN ! 'myfile.dat' to either 2 or 3, depending on wether you want to set cuts or not (try help USER/CONTROL/CARDS in BRAG). The format of OSC1997 files are file : format1 format2 code event_list format1 : 'OSC1997A' '\n' format2 : 'final_id_p_x' '\n' code : code_name version '(' aproj ',' zproj ')+(' atarg ',' ztarg ')' refframe ',' ebeam ',' ntestpart '\n' event_list : /* empty */ | event event_list event : event_info particle_list event_info : event_number number_of_particles impact_paramter reaction_plane '\n' particle_list: /* empty */ | particle particle_list particle : number particle_id px py pz p0 mass x y z t The full description can be found at [1]. I think the header appropiate for AMPT is: OSC1997A final_id_p_x AMPT 3.14 (197,79)+(197,79) aacm 200 1 # Comments are allowed before the first event header, so you can put # stuff here too :-) # I have no idea what the version number of AMPT is, so I just put # in pi :-) Micheal, perhaps you could convince the AMPT people to write the output files in the format above, or perhaps OSC1999? The later format has room for full event history too, so in a way that could be interresting for them too. You can pick up a small C++ program that converts the AMPT format to OSC1997A in [2]. Compile it with: prompt% g++ ampt2osc1997a.cc -o ampt2osc1997a Use it as a filter or with explicit arguments: prompt% cat ampt_file | ./ampt2osc1997a > osc1997a_file prompt% ./ampt2osc1997a ampt_file > osc1997a_file prompt% ./ampt2osc1997a ampt_file osc1997a_file You can also pick up an OSCAR<->ROOT conversion utility at [3] Yours, ____ | Christian Holm Christensen |_| | ------------------------------------------------------------- | | Address: Sankt Hansgade 23, 1. th. Phone: (+45) 35 35 96 91 _| DK-2200 Copenhagen N Cell: (+45) 24 61 85 91 _| Denmark Office: (+45) 353 25 305 ____| Email: cholm@nbi.dk Web: www.nbi.dk/~cholm | | [1] http://www-cunuke.phys.columbia.edu/people/molnard/mirror-OSCAR/oscar/ [2] rcf:/brahms/u/cholm/ampt2osc1997a.cc [3] http://cholm.home.cern.ch/cholm/root
This archive was generated by hypermail 2b30 : Mon May 27 2002 - 10:53:26 EDT