DV Running GEMC

gcards

For 0T Solenoid Field: eg12.gcard

<gcard>

        <detector name="../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/dc/dc"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/ec/ec"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/ctof/ctof"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/ftof/ftof"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/htcc/htcc"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/pcal/pcal"            factory="TEXT" variation="original"/>
        <option name="SCALE_FIELD" value="clas12-torus-big, -1.0"/>
        <option name="OUTPUT" value="evio,eg12.ev"/>

</gcard>

For 5T Solenoid Field: eg12_sol.gcard

 
<gcard>

        <detector name="../../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/dc/dc"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/ec/ec"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/ctof/ctof"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/ftof/ftof"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/htcc/htcc"            factory="TEXT" variation="original"/>
        <detector name="../../../../../clas12/pcal/pcal"            factory="TEXT" variation="original"/>
        <option name="SCALE_FIELD" value="clas12-torus-big, -1.0"/>
        <option name="HALL_FIELD"  value="clas12-solenoid"/>
        <option name="SCALE_FIELD" value="clas12-solenoid, 1.0"/>
        <option name="OUTPUT" value="evio,eg12.ev"/>

</gcard>

The Solenoid can vary between -1.2 to 1.2 on its field strength, where 1 corresponds to 5T and -1 to -5T.

Running GEMC

Setup

Pre-commands

cd
cd ~/src/CLAS/GEMC
tcsh
source setup
cd experiments/eg12/MolrBckGrd/DV

3 directories are created in this area

No_Solenoid  Solenoid  Solenoid_half

Within No_Solenoid:

trial1  trial2  trial3  trial4  trial5

Within Solenoid are additional directories for positive/negative 5T:

Solenoid_negative  Solenoid_positive

Within Solenoid_half are additional directories for positive/negative 2.5T:

Solenoid_negative  Solenoid_positive

The subdirectories for positive and negative have 5 trials within as well similar to No_Solenoid

Running

The gcards are distributed from within the positive 5T, trial1 subdirector to the remaining solenoid field strength matching directories via:

echo ../Solenoid_positive/trial* | xargs -n 1 cp eg12_sol.gcard

This is repeated for all remaining field stregths.

The 5 different LUND files are distributed in a similar fashion:

echo *Sol*/trial1 | xargs -n 1 cp MollerScattering_set0.LUND
echo Sol*/Sol*/trial1 | xargs -n 1 cp MollerScattering_set0.LUND

From within each separate subdirectory for each test's conditions; the commands are for No_Solenoid/test1:

gemc -Hall_Material="Vacuum" -INPUT_GEN_FILE="LUND, MollerScattering_set0.LUND" -N=75000 eg12.gcard 
~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -i eg12.ev -config DCHB::torus=1.0 -config DCHB::solenoid=0.0 -config DCTB::kalman=true -o eg12_rec.ev  -s DCHB:DCTB:EC:FTOF:EB
~/src/CLAS/coatjava-1.0/bin/run-groovy Analysis.groovy eg12_rec.evio

---

Links

Back