Difference between revisions of "Run in GEMC"

From New IAC Wiki
Jump to navigation Jump to search
Line 136: Line 136:
 
   | 571
 
   | 571
 
   | 1
 
   | 1
 +
  |
 +
  |
 
   |
 
   |
 
   |
 
   |
Line 141: Line 143:
 
   | 571
 
   | 571
 
   | 1
 
   | 1
 +
  |
 +
  |
 
   |
 
   |
 
   |
 
   |
Line 146: Line 150:
 
   | 571
 
   | 571
 
   | 1
 
   | 1
 +
  |
 +
  |
 
   |
 
   |
 
   |
 
   |
Line 151: Line 157:
 
   | C
 
   | C
 
   | D
 
   | D
 +
  |
 +
  |
 
   |
 
   |
 
   |
 
   |

Revision as of 15:32, 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.

Mol E Theta Lab Spread.png


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:

[math] 571*1=571*3=1713-2=1711 23025*3=69075*3=207225-2=207223 24726*1=24726*3=74178-2=74176 19970*3=59910*3=179730-2=179728 [/math]

Multiplying the line number by 3 accounts for the 3 lines per particle in the LUND file. From this LUND particle we find the weights associated with each particle (Specific Energy and angle Theta) to be:

using an awk command:

awk 'NR == 1711 { print $6 }' LUND_Spread.LUND

Gives the weight:

<math>

1097.227223687062

25951.273801947842

2579.181846387747

14255.222465310777


Caption
Particle number Run number LUND Moller electron number Overall LUND particle number Corrected LUND header line number Weight
571 1
571 1
571 1
C D

5.5GeV incident electron

Mol E Theta Lab Spread half.png

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