Difference between revisions of "TS RunGroupC"

From New IAC Wiki
Jump to navigation Jump to search
Line 33: Line 33:
 
   |+ Reconsstruction efficiency
 
   |+ Reconsstruction efficiency
 
|-
 
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]] |
+
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]  
 
|-
 
|-
 
|Generated with no solenoid. || Reconstructed with no solenoid.|| reconstructed divided by generated with no solenoid 64 Bins, 6.4 Bins .
 
|Generated with no solenoid. || Reconstructed with no solenoid.|| reconstructed divided by generated with no solenoid 64 Bins, 6.4 Bins .

Revision as of 07:30, 15 March 2016

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 NoSolGcard 022916.txt


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.


Reconsstruction efficiency
ThetaPhiNoSolGen.png ThetaPhiRecNoSol.png ThetaPhiDivNoSol.png NSEffZoom.png
Generated with no solenoid. Reconstructed with no solenoid. reconstructed divided by generated with no solenoid 64 Bins, 6.4 Bins .

Solenoid On

Reconsstruction efficiency
ThetaPhiGenSol.png ThetaPhiRecSol.png DivThetaPhiSol.png
Generated with solenoid. Reconstructed with solenoid. reconstructed divided by generated with solenoid. 20x Bins 50y 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. NSRecDivRec.png

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 ./

EG12#Tyler_Stokes_EG12_studies