Difference between revisions of "Run in GEMC"
Jump to navigation
Jump to search

| Line 77: | Line 77: | ||
</pre> | </pre> | ||
| + | |||
| + | Using cat to combine the reconstructed particles found using Analysis.groovy | ||
| + | |||
| + | For 0T: | ||
| + | <pre> | ||
| + | 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 | ||
| + | </pre> | ||
| + | |||
| + | For 5T: | ||
| + | <pre> | ||
| + | 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 | ||
| + | </pre> | ||
| + | <pre> | ||
=5.5GeV incident electron= | =5.5GeV incident electron= | ||
<center>[[File:Mol_E_Theta_Lab_Spread_half.png]]</center> | <center>[[File:Mol_E_Theta_Lab_Spread_half.png]]</center> | ||
Revision as of 20:32, 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
