Difference between revisions of "Run in GEMC"
| Line 174: | Line 174: | ||
! Corrected LUND header line number | ! Corrected LUND header line number | ||
! Weight | ! Weight | ||
| + | ! Number of Particles | ||
|- | |- | ||
| 16544 | | 16544 | ||
| Line 181: | Line 182: | ||
| 49630 | | 49630 | ||
| 814.898350440404 | | 814.898350440404 | ||
| + | | | ||
|- | |- | ||
| 354 | | 354 | ||
| Line 188: | Line 190: | ||
| 1060 | | 1060 | ||
| 556.656736645754 | | 556.656736645754 | ||
| + | | | ||
|- | |- | ||
| 5526 | | 5526 | ||
| Line 195: | Line 198: | ||
| 66310 | | 66310 | ||
| 937.328333457728 | | 937.328333457728 | ||
| + | | | ||
|- | |- | ||
| 21892 | | 21892 | ||
| Line 202: | Line 206: | ||
| 262702 | | 262702 | ||
| 8343.893471134810 | | 8343.893471134810 | ||
| + | | | ||
|} | |} | ||
Revision as of 16:12, 21 June 2016
Old runs
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, we can output the particle number from the LUND file, Total Energy, scattering angles Theta, and Phi for the generated and reconstructed particles respectively.
For 0T:
571-1 0.12103204139108564 5.232292493596922 9.99999999997547 0.2397734 64.37257 163.2439 23025-3 0.02376743729734356 11.825408412314742 9.999999999978474 0.29301766 89.180336 -164.40204
For 5T:
24726-1 0.08576327895104396 6.223977074710134 9.99999999997547 0.17311224 20.298353 39.036583 19970-3 0.026540097007028158 11.193535798440614 9.999999999969427 0.15106007 49.16122 14.2969885
Converting the particle number per run into the particle number via the LUND file:
Multiplying the line number by 3 accounts for the 3 lines per particle in the LUND file. From this LUND particle we find the differential cross section associated with each particle (Specific Energy and angle Theta).
To find the differential cross section we can use an awk command:
awk 'NR == 1711 { print $6 }' LUND_Spread.LUND
From this, assuming an NH3 target we can find the weight by dividing by 1.19E-6 barns
| Particle number | Run number | LUND Moller electron number | Overall LUND particle number | Corrected LUND header line number | X-Section | Number of Particles |
|---|---|---|---|---|---|---|
| 571 | 1 | 571 | 1713 | 1711 | 1097.227223687062 | 922039684 |
| 23025 | 3 | 69075 | 207225 | 207223 | 25951.273801947842 | 21807793111 |
| 24726 | 1 | 24726 | 74178 | 74176 | 2579.181846387747 | 2167379703 |
| 19970 | 3 | 59910 | 179730 | 179728 | 14255.222465310777 | 11979178542 |
5.5GeV incident electron
For 0T:
16544-1 0.10778422619732184 5.519947956962607 9.999999999974188 0.19995739 25.982777 96.17941
For 5T:
354-1 0.09640233076149667 5.8423283009263045 9.999999999976643 0.09051078 49.163692 -85.56053 5526-4 0.03948624565299794 9.164813676343726 9.999999999971514 0.12387912 52.14526 50.051956 21892-4 0.02741437437384611 11.000629446532491 9.99999999997668 0.26457193 81.77688 179.14293
| Particle number | Run number | LUND Moller electron number | Overall LUND particle number | Corrected LUND header line number | Weight | Number of Particles |
|---|---|---|---|---|---|---|
| 16544 | 1 | 16544 | 49632 | 49630 | 814.898350440404 | |
| 354 | 1 | 354 | 1062 | 1060 | 556.656736645754 | |
| 5526 | 4 | 22104 | 66312 | 66310 | 937.328333457728 | |
| 21892 | 4 | 87568 | 262704 | 262702 | 8343.893471134810 |