E1039 simulation tutorial 201904 - E1039-Collaboration/e1039-wiki GitHub Wiki

Obtain a FermiLab computing account

All the related libraries are compiled on the spinquestgpvm01.fnal.gov machine. To access those machines, please start with obtaining a FermiLab computing account if you haven't got one.

login, setup and checkout the code

init kerberos key, <usr> below means your user name.

kinit -f <usr>@FNAL.GOV
<type your password>

ssh -Y <usr>@spinquestgpvm01.fnal.gov

cd to your work area

mkdir -p /e906/app/users/$USER/tutorial1904
cd /e906/app/users/$USER/tutorial1904

source the setup macro

source /e906/app/software/osg/software/e1039/this-e1039.sh

check out analysis repository from github

git clone https://github.com/E1039-Collaboration/e1039-analysis.git

How to run

Run the simulation + reconstruction macro

cd e1039-analysis/SimChainDev/
root -l Fun4Sim.C\(100\)

You will see several root output files. 'trk_eval.root' contains evaluation output from the TrkEval module. This will be the input for the final plotting.

Run 'int_run.C' to view the event display:

root -l int_run.C

You should get a event display like this: G4 Event Display - DY

To get a display for mac OS users, please refer to this report from John Haggerty @ BNL. The key step for me is this:

defaults write org.macosforge.xquartz.X11 enable_iglx -bool true

Some simple analysis

Open the root file 'trk_eval.root' generated by the TrkEval module

root -l trk_eval.root

Hit reconstruction eval:

Truth->Draw("pos+driftDistance-truth_pos:detectorID","detectorID<31","box")

hit-eval

Single track momentum reconstruction eval:

Truth->Draw("(pz-gpz)/gpz>>h(100,-0.1,0.1)","")

track-eval

Dimuon lab phi reconstruction eval:

Truth->Draw("dimu_phi-dimu_gphi","dimu_phi<1000")

dimu-eval