Build: 1. set environment to any hydra2 version (only ROOT + Makefile are needed). . /cvmfs/hades.gsi.de/install/5.34.34/hydra2-XXX/defall.sh 2. svn checkout https://subversion.gsi.de/hades/hydra2/trunk/programs/converter/smash smashconverter_make cd smashconverter_make setting up the environment: ROOT + optional PLUTO 1. ROOTSYS has to be set 2. Makefile : if PLUTO should be used (version 1.2) PLUTODIR has to be set make Run: 1. set environment to any hydra2 version (ROOT) . /cvmfs/hades.gsi.de/install/5.34.34/hydra2-XXX/defall.sh 2. version 1.3 USAGE: convert-smash -i infile -o outfile [-h,v,u,a,r,w,c,p,n] [-e nevents num] [-b beam_energy ] [-s seed] [-d databasefile] -i (--input) infile : smash OSCAR13/OSCAR13Extended ascii/Root file format or UrQmd f14 input (Center of mass system). -o (--output) outfile: will be in HGeant ascii evt format. -v (--verbose) : verbose output level (0 default, 1,2,3) -e (--events) num : run num events (default is all events) -h (--help) : print this help -a (--incSpec) : include spectators (no default) -r (--recoRho) : reconstruct rho mesons for external decay (no default, needs OSCAR13 format) -w (--recoOmega) : reconstruct rho,omega mesons for external decay (no default, needs OSCAR13 format) -p (--recoPhi) : reconstruct phi mesons for external decay (no default, needs OSCAR13 format) -c (--secondDecay) : if PLUTO decays are used, known channels daughters will be be decayed too (no default, needs OSCAR13 format, compile #define WITHPLUTO ) -u (--pluto) : required my reco options (no default, needs OSCAR13 format, compile #define WITHPLUTO ) -n (--ntuple) : produce PLUTO Tree (no default) -b (--boost) : boost CM to LAB frame (no default) format : \"EKin/u\" [GeV] -s (--seed) num : run with this random seed (no default) -d (--databasefile) : decay data base file (required when -p is used) HINT: options recoRho,recoOmega,recoPhi,pluto,ntuple can be set in databasefile in the form keyword=value (0 or 1) options set in the databasefile will be overwritten by command line options. FORMAT database file example lepton decays using the reco options: Extension has to be .dat (active=0 decays will be ignored) #------------------------------------------------------------------ doRecoRho=1 doRecoOmega=1 doRecoPhi=1 doPluto=1 doSecondDecay=1 doTree=0 type=pi0dalitz ;channel=pi0:dilepton[e+,e-],g;branching=1.198e-2 ;active=1; type=etadalitz ;channel=eta:dilepton[e+,e-],g;branching=0.6e-2 ;active=1; type=eta ;channel=eta:e+,e- ;active=0; type=rho0 ;channel=rho0:e+,e- ;branching=0.00448e-2;active=1; type=w ;channel=w:e+,e- ;branching=0.00715e-2;active=1; type=wdalitz ;channel=w:dilepton[e+,e-],pi0;branching=0.059e-2 ;active=1; type=phi ;channel=phi:e+,e- ;branching=0.0299e-2 ;active=1; #------------------------------------------------------------------ FORMAT database as string example (line seperator \"#\"): option -d "doRecoRho=1#doRecoOmega=1#doRecoPhi=1#doPluto=1#doTree=0#type=pi0dalitz;channel=pi0:dilepton[e+,e-],g;active=1;" The key word "branching=" for branching ratio is optional, by default branching is assumed to be 1. The branching ratio is used when calculating the weights of a particle : weight = active nchannels for this mother particle * branching. For the above example w+wdalitz = 2 active channels. weight for daughters of w = 2*0.00715e-2. HELP on geninfo, geninfo1, geninfo2 in output to identify sources: Standard smash output: 1. single particles keep Int_t sourceID = pow(10,1*2) * parentID1 sourceID += pow(10,0*2) * parentID2 genInfo1 = sourceID genInfo2 = 1000000 + proc_type_origin genInfo3 = 1000000 + proc_id_origin example proton from elastic scattering : 1414 100001 100001 2. decay daughters from smash rho (0,+,-), w, phi (without reco options) sourceID = motherID * pow(10,1*3) + daughter1ID * pow(10,1*2) + daughter2ID * pow(10,1*1) genInfo2 = 1000000 + proc_type_origin genInfo3 = 1000000 + proc_id_origin example : rho+ -> pi+ p- 420807 100005 100003 3. decay daughters from PLUTO decays sourceID = motherID * pow(10,1*3) + daughter1ID * pow(10,1*2) + daughter2ID * pow(10,1*1) genInfo2 = like PLUTO genInfo3 = counter for each PLUTO decay in event example : pi0 -> dipleton[e+,e-],g 75101 7051(e+,e-) or 7(g) 1 OUTPUT FORMAT from OSCAR13 extended input contains generator infos Event Header : Evt Nr NParticle Collision energy Impact parameter HGeant format flag Particle line : E, px, py, pz, geantID, genInfo1, genInfo2, genInfo3, weight 5 6 4.500000 0.396000 -2 1.311937 0.505859 0.246233 0.724438 13 1313 1000003 1000001 0.000000 2.010489 -0.408474 0.391084 1.924376 8 808 1000005 1000002 0.000000 2.473349 0.123280 -0.655407 2.189259 14 1414 1000005 1000003 0.000000 0.132957 0.003891 -0.009491 0.132561 1 75101 7 1 0.011980 0.027938 -0.013071 0.001796 0.024621 2 75101 7051 1 0.011980 0.417129 -0.211485 0.025785 0.358616 3 75101 7051 1 0.011980 OUTPUT FORMAT from UrQMD f14 input contains no generator infos Event Header : Evt Nr NParticle Collision energy Impact parameter HGeant format flag Particle line : E, px, py, pz, geantID, weight convert-smash expects : smash output should be set to "fixed target" !! see example config_au123au_oscar13_smash_1.8.yaml (smash >=1.8) config_au123au_oscar13_smash_1.5.yaml (smash 1.5) test input files are located in data/0 ./convert-smash -i data/0/particle_lists.oscar -o test.evt version: 1.0 : initial commit to svn 1.1 : added K0L K0_SS handling (50:50) 1.2 : - config option removed - added --boost to boost from CMS to LAB. The boost beta is calculated by the reaction string, for example -b 4.5:1:1 for p+p at 4.5AGeV or 1.23:197:197 for Au+Au at 1.23AGeV - added --recoMesons : reconstruct mesons from oscar13 from their decay products and modify the event - added --pluto : neeeds --recomesons. Decay the mesons using PLUTO into predfined channels. WITHPLUTO in the top of convert-smash.cc has to be defined to compile with PLUTO - added --ntuple : needs WITHPLUTO : write the event in PLUTO Tree format. - added --seed : set a random seed, by default a random seed will be generated by /dev/urandom 1.3 : added decay of particle (PLUTO) , generic reaction build (leptons+hadrons)+ database input file added geninfo, geninfo1, geninfo, weight to output (PLUTO to output) boost option changed : -b 4.5 for p+p at 4.5AGeV or -b 1.23 for Au+Au at 1.23AGeV 1.4 : in case recoRho etc, check event after omega reco for remaining rho+- and remove them(putting org smash particles in place) 1.5 : adding branching rations for weight calculations. added doSecondDecay option 1.6 : -incSpec option (by default spectators are not not propagated to the output)