Difference between revisions of "LB RunGroupC Vertex"

From New IAC Wiki
Jump to navigation Jump to search
Line 6: Line 6:
  
  
−
=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
 
−
 
−
[[File:rec_theta_1k.png | 200 px]]
 
−
 
−
[[File:Rec_Theta_Sol_1k.png | 200 px]]
 
−
 
−
 
−
[[File:Rec_phi_nosol_1k.png | 200 px]]
 
−
 
−
 
−
[[File:Rec_phi_sol_1k.png | 200 px]]
 
−
 
−
==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:
 
−
 
−
[[File:PhiDifference75kZraster.png|200px]]
 
−
[[File:ThetaDifference75kZraster.png | 200px]]
 
−
 
−
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:
 
−
 
−
 
−
[[File:ThetaDifferenceZRasterSol.png | 200px]]
 
−
[[File:PhiDifferenceZrasterSol.png | 200px]]
 
−
 
−
==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:
 
−
 
−
[[File:ThetaDifferenceXYRaster.png | 200px]]
 
−
[[File:PhiDifferenceXYRaster.png | 200px]]
 
−
 
−
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:
 
−
 
−
[[File:ThetaDifferenceXYrasterSol.png | 200px]]
 
−
[[File:PhiDifferenceXYrasterSol.png | 200px]]
 
−
 
−
== 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:
 
−
 
−
[[File:ThetaDifferenceXYZRasterNoSol.png | 200px]]
 
−
[[File:PhiDifferenceXYZRasterNoSol.png | 200px]]
 
−
 
−
 
−
 
−
 
−
and here are the plots for when the solenoid was on:
 
−
 
−
 
−
[[File:ThetaDifferenceXYZSol.png | 200px]]
 
−
[[File:PhiDifferenceXYZRasterSol.png | 200px]]
 
  
 
==Summary==
 
==Summary==

Revision as of 18:22, 16 June 2016

Vertex Reconstruction studies

GEMC 2.4 & Coatjava 2.4

LB_GEMC2.4CJ2.4_Vertex



Summary

Var X(cm) Var Y(cm) Var Z(cm) X resold(cm) Y resolution (cm) Z resolution (cm)


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

Retrieved from "https://wiki.iac.isu.edu/index.php?title=LB_RunGroupC_Vertex&oldid=106495"