TS RunGroupC
Inclusive electron scattering
Theta and Phi graphs
Solenoid On
Inclusive Track Reconstruction efficiencies
Plot the reconstruction efficiency for an electron as a function of the angles Theta and Phi
I used the following: UBUNTU 14, coat-java-1.0
No Solenoid gcard:
angle spane 25 +/- 15
File:TS SolGcard 022916.txt
gcard for Solenoid
Solenoid Off
Why is efficiency grater than 1?
Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.
NOTE: The reason why the graphs do not extend down to 0 degrees in the new uploads is because the gcard is 25 +- 15. So, 25-15=10 25+15=40 10 and 40 are the boundaries.
 |
 |
 |
|
| Generated with no solenoid. | Reconstructed with no solenoid. | reconstructed divided by generated with no solenoid 64 Bins, 6.4 Bins . |
Solenoid On
 |
 |
 |
|
| Generated with solenoid. | Reconstructed with solenoid. | reconstructed divided by generated with solenoid. 71 Bins, 7.1 Bins |
Ratio Solenoid Off / On
NOTE: The reason the division is at around 5 is because for some reason there were 5x the amount of no solenoid events vs solenoid events. I ran 975,000 events for both no solenoid and solenoid, but a reasonable ~500,000 came out for no solenoid. For some reason only ~100,000 came out for solenoid.
No solenoid reconstructed divided by solenoid reconstructed.
Phi angle shift
Phi angle Solenoid on - off difference
spread from 0 to 5 GeV.
Observations
Why does the Solenoid off events only extend to 25 degrees whereas the solenoid on go to 40.
Vertex Reconstruction
For particle at origin
plot Generated xpos - Reconstructed x position
With Beam rastering
Commands
Simulation commands
tcsh cd ~/src/CLAS/GEMC source setup cd experiments/eg12/lh2 gemc -HALL_MATERIAL="Vacuum" eg12.gcard gemc -USE_GUI=0 -HALL_MATERIAL="Vacuum" -N=10 eg12.gcard
to edit files
emacs -nw filename
without solenoid
~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -config DCHB::torus=1.0 -config DCHB::solenoid=0.0 -config DCTB::kalman=true -s DCHB:DCTB:EC:FTOF:EB -i eg12.ev -o eg12_rec.evio
with solenoid
~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -s FMT:DCHB:DCTB:FTOF:EC:EB -config MAG::fields=nominal -i eg12_5GeVe.evio -o eg12_rec.evio
evio analysis program
~/src/CLAS/coatjava-1.0/bin/run-groovy Analysis.groovy eg12_rec.0.evio
dump file contents
~/src/CLAS/coatjava-1.0/bin/evio-dump -i eg12_rec.0.evio
batch job cheat
nohup gemc -USE_GUI=0 -HALL_MATERIAL="Vacuum" -N=75000 eg12.gcard > /dev/null &
Analysis commands
cat Analysis.groovy Analysis.groovy.old1 Sim.dat > trash
Converting files to ROOT
~/src/CLAS/coatjava-1.0/bin/run-groovy Evio2Text.groovy eg12_rec.0.evio > Sim.dat
edit Sim.dat and remove the first few lines until you have only numbers
minerve2:~/src/CLAS/GEMC/experiments/eg12/lh2> root -l
root [0] .L Asci2Root.C
root [1] Asci2Root("Sim.dat");
Using ROOT to analyze a file
minerve2:~/src/CLAS/GEMC/experiments/eg12/lh2> root -l Sim.root
root [0] Attaching file Sim.root as _file0...
root [1] new TBrowser();
root [2] CLAS12Sim->Draw("evt.Theta");
root [3] CLAS12Sim->Draw("evt.Theta:evt.Phi");
root [4] CLAS12Sim->Draw("evt.Theta:evt.Phi","","colz");
root [5] TH2D *Rec=new TH2D("Rec","Rec",400,-200,200, 50, 0,50);
root [6] CLAS12Sim->Draw("evt.Theta:evt.Phi>>Rec","","colz");
root [7] TH2D *Gen=new TH2D("Gen","Gen",400,-200,200, 50, 0,50);
root [9] CLAS12Sim->Draw("evt.Theta:evt.Phi>>Gen","","colz");
root [10] Gen->Draw("colz")
root [11] Rec->Draw("colz")
root [12] Rec->Divide(Gen)
root [4] TH1D *RecPhi=new TH1D("RecPhi","RecPhi",400, -200,200);
root [5] CLAS12Sim->Draw("evt.Phi>>RecPhi");
go to Generated Events File
root [8] CLAS12Sim->Draw("evt.Phi>>GenPhi"); root [9] RecPhi->Divide(GenPhi);
return to ROOT file with reconstructed events and do
root [6] RecPhi->Divide(GenPhi); root [7] RecPhi->Draw();
TColor*col;
Int_t colors[3];
colors[0]=1;
colors[1]=2;
colors[2]=3;
gStyle->SetPalette(3,colors);
File Type Search:
find . -type f -name "*.txt"
scp to laptop
scp stoktyle@minerve2.cose.isu.edu:src/CLAS/GEMC/experiments/eg12/lh2/Sim.root ./