Difference between revisions of "Run in GEMC"

From New IAC Wiki
Jump to navigation Jump to search
Line 21: Line 21:
  
 
<pre>find . -name "*.LUND" -exec sh -c 'mkdir "${1%.*}" ; mv "$1" "${1%.*}" ' _ {} \;</pre>
 
<pre>find . -name "*.LUND" -exec sh -c 'mkdir "${1%.*}" ; mv "$1" "${1%.*}" ' _ {} \;</pre>
 +
 +
Defining the gcards
 +
<pre>
 +
<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>
 +
</pre>
 +
 +
Using a file to store the needed commands
 +
<pre>
 +
/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
 +
 +
</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:30, 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.

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

5.5GeV incident electron

Mol E Theta Lab Spread half.png