LB Summer 2017 JLab Vertex Work
Jump to navigation
Jump to search
Simulation Setup
The following simulations were done using GEMC 4a.0.0.It is useful to note that the LUND file vertices are in units of cm, while the reconstructed events are in mm.Below is a sample gcard that was used for many of my simulations
<gcard>
<detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/fc/forwardCarriage" factory="TEXT" variation="original"/> <detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/dc/dc" factory="TEXT" variation="java"/> <detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/ec/ec" factory="TEXT" variation="java"/> <detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/ctof/ctof" factory="TEXT" variation="cad"/>
<detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/ctof/cad/" factory="CAD"/> <detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/ftof/ftof" factory="TEXT" variation="java"/> <detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/htcc/htcc" factory="TEXT" variation="original"/> <detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/pcal/pcal" factory="TEXT" variation="java"/> <detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/magnets/solenoid" factory="TEXT" variation="original"/>
<detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/magnets/torus" factory="TEXT" variation="original"/>
<detector name ="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/targets/target" factory="TEXT" variation="lH2"/>
<detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/bst/bst" factory="TEXT" variation="java"/>
<detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/cnd/cnd" factory="TEXT" variation="original"/>
<detector name="../../../../../../../clas12Tags/4a.0.0/experiments/clas12/micromegas/micromegas" 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_sol_75k_NoRaster_30DegRot.ev"/>
</gcard>
and here is a sample GEMC command line that was used for the simulations. Obviously the name of the LUND file, the number of events, and the output filename will change.
The LUND files were created using WriteLUND.C. The vertex positions are generated by a RNG whose output is between 2 numbers. This will determine if the target is 1D,2D,3D. To check the angles on the LUND files, run a GEMC simulation with only a few events to see where the particles hit within the detector. For angles
theta=
phi=