Difference between revisions of "Run in GEMC"
Jump to navigation
Jump to search
Line 95: | Line 95: | ||
<center>[[File:Mol_E_Theta_Lab_Spread_half.png]]</center> | <center>[[File:Mol_E_Theta_Lab_Spread_half.png]]</center> | ||
+ | |||
+ | For 0T: | ||
+ | <pre> | ||
+ | 16544 0.10778422619732184 5.519947956962607 9.999999999974188 0.19995739 25.982777 96.17941 | ||
+ | </pre> | ||
+ | |||
+ | |||
+ | For 5T: | ||
+ | <pre> | ||
+ | 354 0.09640233076149667 5.8423283009263045 9.999999999976643 0.09051078 49.163692 -85.56053 | ||
+ | 5526 0.03948624565299794 9.164813676343726 9.999999999971514 0.12387912 52.14526 50.051956 | ||
+ | 21892 0.02741437437384611 11.000629446532491 9.99999999997668 0.26457193 81.77688 179.14293 | ||
+ | </pre> |
Revision as of 20:34, 20 June 2016
11GeV incident electron
For an incident electron at 11GeV, the Moller electron in the range of 5 to 40 degrees for the scattering angle Theta, has energies lower than 200MeV.
Spliting the LUND file into 5 parts of 25000 events, or 75000 lines
split -d -l 75000 -a 2 LUND_Spread.LUND LUND_Spread_
Naming these files LUND files
prename 's/(LUND_Spread_\d{2})/$1.LUND/' LUND_Spread_*
Placing each of these files into its own directory, within a directory named LUND_Spread
find . -name "*.LUND" -exec sh -c 'mkdir "${1%.*}" ; mv "$1" "${1%.*}" ' _ {} \;
Defining the gcards
<gcard> <!-- minumal detectors for HB tracking --> <detector name="../../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/> <detector name="../../../../../../clas12/dc/dc" factory="TEXT" variation="ccdb"/> <detector name="../../../../../../clas12/ftof/ftof" factory="TEXT" variation="java"/> <option name="SCALE_FIELD" value="clas12-torus-big, -1"/> <option name="OUTPUT" value="evio, Molr.evio"/> <!-- Solenoid apparatus and field --> <!--detector name="../../../../../../clas12/magnets/solenoid" factory="TEXT" variation="original"/--> <!--option name="HALL_FIELD" value="clas12-solenoid"/--> <!-- other CLAS detectors --> <detector name="../../../../../../clas12/ec/ec" factory="TEXT" variation="original"/> <detector name="../../../../../../clas12/ctof/ctof" factory="TEXT" variation="original"/> <detector name="../../../../../../clas12/htcc/htcc" factory="TEXT" variation="original"/> <detector name="../../../../../../clas12/pcal/pcal" factory="TEXT" variation="javageom"/> <!--detector name="../../../../../../clas12/micromegas/micromegas" factory="TEXT" variation="original"/--> <!-- Beamline material --> <!--detector name="../../../../../../clas12/beamline/beamline" factory="TEXT" variation="ft"/--> </gcard>
Using a file to store the needed commands
/src/CLAS/GEMC/source/gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,MolTest.LUND" -N=100 eg12.gcard Solenoid off ~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -i ../files/MolrBckGrd/Molr.evio -config GEOM::new=true -config MAG::torus=1.0 -config MAG::solenoid=0.0 -o ../files/MolrBckGrd/Molr_rec.evio -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true Solenoid ON ~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -i ../files/MolrBckGrd/Molr.evio -config GEOM::new=true -config MAG::torus=-1.0 -config MAG::solenoid=1.0 -o ../files/MolrBckGrd/Molr_rec.evio -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true -config DCTB::useRaster=true ~/src/CLAS/coatjava-2.4/bin/evio-dump -i ../files/MolrBckGrd/Molr_rec.0.evio
Using cat to combine the reconstructed particles found using Analysis.groovy
For 0T:
571 0.12103204139108564 5.232292493596922 9.99999999997547 0.2397734 64.37257 163.2439 23025 0.02376743729734356 11.825408412314742 9.999999999978474 0.29301766 89.180336 -164.40204
For 5T:
24726 0.08576327895104396 6.223977074710134 9.99999999997547 0.17311224 20.298353 39.036583 19970 0.026540097007028158 11.193535798440614 9.999999999969427 0.15106007 49.16122 14.2969885
5.5GeV incident electron
For 0T:
16544 0.10778422619732184 5.519947956962607 9.999999999974188 0.19995739 25.982777 96.17941
For 5T:
354 0.09640233076149667 5.8423283009263045 9.999999999976643 0.09051078 49.163692 -85.56053 5526 0.03948624565299794 9.164813676343726 9.999999999971514 0.12387912 52.14526 50.051956 21892 0.02741437437384611 11.000629446532491 9.99999999997668 0.26457193 81.77688 179.14293