LB GEMC2.4CJ2.4 Vertex

From New IAC Wiki
Revision as of 18:22, 16 June 2016 by Foretony (talk | contribs)
Jump to navigation Jump to search

Reconstruction Efficiencies

Simulation commands

Solenoid Off

~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -SPREAD_V="(0.1, 10)cm" -N=100 eg12_nosol_noLUND.gcard

~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -i eg12_nosol_75k_noLUND.ev -config GEOM::new=true -config MAG::torus=-1.0 -config MAG::solenoid=0.0 -o temp.evio -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true

Solenoid On

~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -SPREAD_V="(0.1, 10)cm" -N=100 eg12_sol_noLUND.gcard

~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -i eg12_nosol_75k_noLUND.ev -config GEOM::new=true -config MAG::torus=-1.0 -config MAG::solenoid=1.0 -o temp.evio -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true


Sweet spot

TS_RunGroupC#Sweet_spot Below is the command line I used for GEMC

gemc -USE_GUI=0 -HALL_MATERIAL="Vacuum" , -N=75000 eg12_nosol1.gcard

Here is the 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="javageom"/>
      <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
     <option name="OUTPUT" value="evio,eg12_nosol_1.ev"/>
     <option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>
     <option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

</gcard>

Here is the reconstruction command

~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -s FMT:DCHB:DCTB:FTOF:EC:EB -config DCTB::useRaster=true -config MAG::solenoid=0 -i eg12_nosol_1.ev -o eg12_nosol_1_rec.ev


Below is the command and the gcard used for when the solenoid was on

gemc -USE_GUI=0 -HALL_MATERIAL="Vacuum" , -N=100000 eg12_sol1.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="javageom"/>
       <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
       <option name="HALL_FIELD"  value="clas12-solenoid"/>
       <option name="OUTPUT" value="evio,eg12_sol_1.ev"/>
       <option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>
       <option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

</gcard>

~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -s FMT:DCHB:DCTB:FTOF:EC:EB -config MAG::fields=nominal -config DATA::mc=true -i eg12_sol_1.ev -o eg12_sol_1_rec.evio

Below are the reconstructed events with the solenoid on and off for 1k events

Rec theta 1k.png

Rec Theta Sol 1k.png


Rec phi nosol 1k.png


Rec phi sol 1k.png

Raster in Z

Since the LUND files are not currently operating correctly, I will use a different method to raster in Z. I will use the GEMC command line to reduce the radius of the target to 0 cm and randomly place vertex positions along Z (like a line target). The command line for GEMC is shown below along with the gcard when the solenoid is off:

~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -SPREAD_V="(0.0, 3)cm" -N=100 eg12_nosol_noLUND.gcard


<gcard>

      <detector name="../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
      <detector name="../../clas12/dc/dc"            factory="TEXT" variation="ccdb"/>
      <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="java"/>
      <detector name="../../clas12/htcc/htcc"            factory="TEXT" variation="original"/>
      <detector name="../../clas12/pcal/pcal"            factory="TEXT" variation="javageom"/>
      <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
     <option name="OUTPUT" value="evio,eg12_nosol_75k_noLUND.ev"/>
     <option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>
     <option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

</gcard>

and here is the reconstruction command line ~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -i eg12_nosol_75k_noLUND.ev -config GEOM::new=true -config MAG::torus=-1.0 -config MAG::solenoid=0.0 -o eg12_nosol_75k_noLUND_rec.ev -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true

Below are the differences between the generated and reconstructed angles for 75k events:

PhiDifference75kZraster.png ThetaDifference75kZraster.png

Below are the command lines and gcard used for GEMC and the reconstruction software when the solenoid is on

~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -SPREAD_V="(0.0, 3)cm" -N=75000 eg12_sol_noLUND.gcard


~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -i eg12_sol_75k_noLUND_3cm_Zraster.ev -config GEOM::new=true -config MAG::torus=-1.0 -config MAG::solenoid=1.0 -o eg12_sol_75k_noLUND_3cm_Zraster_rec.evio -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true


<gcard>

      <detector name="../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
      <detector name="../../../clas12/dc/dc"            factory="TEXT" variation="ccdb"/>
      <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="java"/>
      <detector name="../../../clas12/htcc/htcc"            factory="TEXT" variation="original"/>
      <detector name="../../../clas12/pcal/pcal"            factory="TEXT" variation="javageom"/>
      <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
      <option name="HALL_FIELD"  value="clas12-solenoid"/>
      <option name="OUTPUT" value="evio,eg12_sol_75k_noLUND_3cm_Zraster.ev"/>
      <option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>
      <option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

</gcard>


Below are the plots for the differences between the generated and the reconstructed values for theta and phi with the solenoid on for 75k events:


ThetaDifferenceZRasterSol.png PhiDifferenceZrasterSol.png

Raster in X

Raster in Y

Raster in X & Y Zis zero

Below is the command line used for rastering in X and Y while leaving the z-component of the vertex at 0 with the solenoid off:

~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -SPREAD_V="(0.5, 0)cm" -N=75000 eg12_nosol_noLUND.gcard

and here is the reconstruction command line

~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -i eg12_nosol_75k_noLUND_xyraster.ev -config GEOM::new=true -config MAG::torus=-1.0 -config MAG::solenoid=0.0 -o eg12_nosol_75k_noLUND_xyraster_rec.evio -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true

Below are the differences between the reconstructed and the generated angles for rastering 1cm in XY:

ThetaDifferenceXYRaster.png PhiDifferenceXYRaster.png

Below are the command lines and the gcard used for when the solenoid was on:

~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -SPREAD_V="(0.5, 0)cm" -N=75000 eg12_sol_noLUND.gcard

~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -i eg12_sol_75k_noLUND_1cm_XYraster.ev -config GEOM::new=true -config MAG::torus=-1.0 -config MAG::solenoid=1.0 -o eg12_sol_75k_noLUND_1cm_XYraster_rec.evio -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true


<gcard>

      <detector name="../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
      <detector name="../../../../clas12/dc/dc"            factory="TEXT" variation="ccdb"/>
      <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="java"/>
      <detector name="../../../../clas12/htcc/htcc"            factory="TEXT" variation="original"/>
      <detector name="../../../../clas12/pcal/pcal"            factory="TEXT" variation="javageom"/>
      <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
      <option name="HALL_FIELD"  value="clas12-solenoid"/>
      <option name="OUTPUT" value="evio,eg12_sol_75k_noLUND_1cm_XYraster.ev"/>
      <option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>
      <option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

</gcard>



Here are the plots for the differences in the generated and reconstructed values for theta and phi with the solenoid on:

ThetaDifferenceXYrasterSol.png PhiDifferenceXYrasterSol.png

Raster in XYZ

For this simulation I rastered the beam in the XY plane around 1cm, and z was rastered +/- 3cm. Below are the differences in the the generated and reconstructed angles for when the solenoid was off:

ThetaDifferenceXYZRasterNoSol.png PhiDifferenceXYZRasterNoSol.png



and here are the plots for when the solenoid was on:


ThetaDifferenceXYZSol.png PhiDifferenceXYZRasterSol.png


LB_RunGroupC_Vertex#GEMC_2.4_.26_Coatjava_2.4