Difference between revisions of "LB RunGroupC Vertex"

From New IAC Wiki
Jump to navigation Jump to search
Line 3: Line 3:
 
=Electron Vertex reconstruction=
 
=Electron Vertex reconstruction=
  
Here are the plots of differences between the simulated x, y, and z vertex positions and the reconstructed vertex positions for a 2cm target. Each histogram corresponds to a difference in the matching coordinate i.e. Simulated Vx - Reconstructed Vx
+
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
  
Histogram for Vx differences
+
gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND, ElasticEvents_FromOrigin.LUND" -N=75000 eg12.gcard
 
 
[[File:eVxDiffCuts_fnl.jpg]]
 
 
 
Histogram for Vy differences
 
 
 
[[File:eVyDiffCuts_fnl.jpg]]
 
 
 
Histogram for Vz differences
 
 
 
[[File:eVzDiffCuts.jpg]]
 
 
 
Now I will show the histograms for the difference between the simulated and the reconstructed vertex positions for the proton
 
 
 
Histogram for Vx differences
 
 
 
[[File:pVxDiffCuts_fnl.jpg]]
 
 
 
Histogram for the Vy differences
 
 
 
[[File:pVyDiffCuts_fnl.jpg]]
 
 
 
Histogram for the Vz differences
 
 
 
[[File:pVzDiffCuts.jpg]]
 
 
 
Now I will attempt to find the phi shift between the vertex position and the detector. First, I created an ev file called eg12_lb_NoSol.ev with GEMC using the command line
 
gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND, ElasticEvents_FromOrigin.LUND" -N=7000 eg12.gcard
 
  
  
Line 45: Line 18:
 
         <detector name="../../../clas12/pcal/pcal"            factory="TEXT" variation="original"/>
 
         <detector name="../../../clas12/pcal/pcal"            factory="TEXT" variation="original"/>
 
       <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
 
       <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
       <option name="OUTPUT" value="evio,eg12.ev"/>
+
       <option name="OUTPUT" value="evio,eg12_lb_NoSol.ev"/>
  
 
then I ran the reconstruction program using
 
then I ran the reconstruction program using
Line 51: Line 24:
 
~/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
 
~/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
  
Once these two evio files have been created, the investigation into the phi shift can begin to be studied by looking at the events in the ced of coatjava. Below I will show a few of the different events for the simulated and the reconstructed particles as seen in the ced.
+
After that I used the eg12_lb_NoSol_rec.ev to create the histograms
 
 
 
 
Simulated
 
 
 
 
 
[[File:Event2_sim.png]]
 
 
 
Reconstructed
 
 
 
 
 
[[File:Event2_rec.png]]
 
 
 
 
 
Simulated
 
 
 
[[File:Event3_sim.png]]
 
 
 
Reconstructed
 
  
[[File:Event3_rec.png]]
+
Here are the histograms when the solenoid is off
  
 
=References=
 
=References=

Revision as of 20:32, 30 October 2015

Vertex Reconstruction studies

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

References

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

EG12#Vertex_Studies