Difference between revisions of "TS RunGroupC"

From New IAC Wiki
Jump to navigation Jump to search
Line 29: Line 29:
 
|[[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. || Zoomed in version of the efficiency.
+
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 64 Bins, 6.4 Bins. || Zoomed in version of the efficiency.
 
|}
 
|}
  

Revision as of 10:43, 3 April 2016

Inclusive electron scattering

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 spans 25 +/- 15

File:TS NoSolGcard 022916.txt


File:TS SolGcard 022916.txt gcard for Solenoid

Solenoid Off

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.


Reconstruction efficiency
ThetaPhiNoSolGen.png ThetaPhiRecNoSol.png ThetaPhiDivNoSol.png NSEffZoom.png
Generated with no solenoid. 824,979 Events Reconstructed with no solenoid. 609,106 Events reconstructed divided by generated with no solenoid 64 Bins, 6.4 Bins. Zoomed in version of the efficiency.

Solenoid On

Reconstruction efficiency
ThetaPhiGenSol.png ThetaPhiRecSol.png DivThetaPhiSol.png SolEffZoom.png
Generated with solenoid. Reconstructed with solenoid. reconstructed divided by generated with solenoid. 60 Bins, 5.0 Bins. Zoomed in version of the efficiency.

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

Elastic Electrons from origin

Efficiency plots for Elasticlly scattered electrons when incident electron energy was 11 GeV


using the file ReconElasticExclusive_FromOrigin.LUND

Solenoid On

Reconstruction efficiency
Gen Elastic.png Recon Elastic.png RecDivGen.png Efficiency Elastic.png
Generated with solenoid. Reconstructed with solenoid. reconstructed divided by generated with solenoid. 60 Bins, 5.0 Bins. Zoomed in version of the efficiency.

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