Recreation of Plots Using GEMC 2.5

From New IAC Wiki
Revision as of 16:23, 21 October 2016 by Brenleyt (talk | contribs)
Jump to navigation Jump to search

Here are the plots that I am trying to recreate

https://clasweb.jlab.org/wiki/index.php/File:EG12_SweetSpotVertex_ZdiffDist_10-6-16.png

Here is the gcard used

<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/magnets/solenoid" 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_1cm_DiskTarg_30DegRot_GEMC_2.5.ev"/>
       <option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 30*deg"/>
               <option name="SPREAD_V" value="(1.0,0.0)cm"/>

</gcard>

Here is the command line for GEMC:

~/src/CLAS/GEMC/2.0.Linux__Ubuntu16.04-x86_64-gcc5.4.0/gemc/2.5/gemc -USE_GUI=0 -N=75000 eg12.gcard

I used the same coatjava command (no fields) as seen below

~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -config GEOM::new=true -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true -config DCTB::useRaster=true -i eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5.ev -o eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5_rec.evio

and finally the evio2root command line (done on lds2):

~/src/CLAS/evioRoot/bin/evio2root eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5_rec.evio eg12_sol_75k_1cm_DiskTarg_30DegRot_GEMC_2.5_rec.root 75000

For the plot below I required that Event.z is non-zero and that Event.z < 2cm.

VzDiff Z=0 DeltaR=0.2.png