LB RunGroupC Vertex
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
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
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