LB RunGroupC Vertex

From New IAC Wiki
Jump to navigation Jump to search

Vertex Reconstruction studies

Reconstruction Efficiencies

Sweet spot

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

gemc -USE_GUI=0 -HALL_MATERIAL="Vacuum" , -N=100000 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

Electron Vertex reconstruction

I have created histograms of the z-vertex resolutions for both the proton and the electron using elastic events from the origin. I used the gemc command line to output a file called eg12_lb_NoSol.ev

gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND, ElasticEvents_FromOrigin.LUND" -N=75000 eg12.gcard


My gcard looks like this

      <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"/>
     <option name="OUTPUT" value="evio,eg12_lb_NoSol.ev"/>

then I ran the reconstruction program using

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

After that I used the eg12_lb_NoSol_rec.ev to create the histograms

Here are the histograms when the solenoid is off

Z Vertex NoSol.png

Then I ran the GEMC simulation with the solenoid on. My gcard was

<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"/>
      <option name="HALL_FIELD"  value="clas12-solenoid"/>
     <option name="OUTPUT" value="evio, eg12_lb_Sol.ev"/>


</gcard>


Which gives the output file eg12_lb_Sol.ev. I used the GEMC command line

gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND, ElasticEvents_FromOrigin.LUND" -N=75000 eg12_Sol.gcard

Then I made the reconstructed file, eg12_lb_Sol_rec.ev using the command line

~/src/CLAS/coatjava//bin/clas12-reconstruction -s FMT:DCHB:DCTB:FTOF:EC:EB -config MAG::fields=nominal -i eg12_Sol.ev -o eg12_lb_Sol_rec.ev

The histogram is below

Z Vertex Sol.png


Now I will get more simulated events and use the java geometry for the pcal and the LUND file Vertex0.075k_10-40.dat and investigate the Z vertex resolution using different radius cuts.

References

https://clasweb.jlab.org/wiki/index.php/TF_EG12_Vertex#Z_resolution_With_micro-megas

EG12#Vertex_Studies