New IAC Wiki - User contributions [en]

TS RunGroupC

2016-04-27T09:44:12Z

Stoktyle: /* Sweet spot */

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Detector Images
|-
|[[File:Re1eetSpotOn.png | 300 px]] ||[[File:RecSwetOff.png | 300 px]]
|-
| No Events|| Solenoid On with a proton and electron event||
|}

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

1-D Phi and Theta histograms for Solenoid Off

Overlay of Solenoid off and On for Theta and Phi.

Difference between Generated and Reconstructed Theta and Phi.

The generated plots are going to be centered at one point. For Theta that is 25 degrees. For Phi that is 0 degrees.
There are 323,066 generated events for no solenoid.
There are 573,362 generated events for solenoid on.
{| class="wikitable" align="center" border=1
|+ Reconstructed 1D No Solenoid
|-
|[[File:RecPhiOff1D.png | 300 px]] ||[[File:RecThetaOff1D.png | 300 px]]
|-
|Rec Phi Off|| Rec Theta Off||
|}

{| class="wikitable" align="center" border=1
|+ Overlay
|-
|[[File:OnOverOff1D.png | 300 px]] ||[[File:ThetaOnOverOff1D.png | 300 px]]
|-
|Phi On Over Off|| Theta On Over Off||
|}

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Reconstructed with solenoid|| Reconstructed with no solenoid|| solenoid off divided by 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 spans 25 +/- 15

[[File:TS_NoSolGcard_022916.txt]]

[[File:TS_SolGcard_022916.txt]]
gcard for Solenoid

=== Solenoid Off===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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"

CLAS12Sim->Draw("evt.Theta >> Gen(20,15,35)");
CLAS12Sim->Draw("evt.Theta >> Rec(20,15,35)");

Gen->Draw();
Rec->Draw("same");

===scp to laptop===

scp stoktyle@minerve2.cose.isu.edu:src/CLAS/GEMC/experiments/eg12/lh2/Sim.root ./

[[EG12#Tyler_Stokes_EG12_studies]]

TS RunGroupC

2016-04-27T09:41:44Z

Stoktyle: /* Sweet spot */

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Detector Images
|-
|[[File:Re1eetSpotOn.png | 300 px]] ||[[File:RecSwetOff.png | 300 px]]
|-
| No Events|| Solenoid On with a proton and electron event||
|}

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

1-D Phi and Theta histograms for Solenoid Off

Overlay of Solenoid off and On for Theta and Phi.

Difference between Generated and Reconstructed Theta and Phi.

The generated plots are going to be centered at one point. For Theta that is 25 degrees. For Phi that is 0 degrees.

{| class="wikitable" align="center" border=1
|+ Reconstructed 1D No Solenoid
|-
|[[File:RecPhiOff1D.png | 300 px]] ||[[File:RecThetaOff1D.png | 300 px]]
|-
|Rec Phi Off|| Rec Theta Off||
|}

{| class="wikitable" align="center" border=1
|+ Overlay
|-
|[[File:OnOverOff1D.png | 300 px]] ||[[File:ThetaOnOverOff1D.png | 300 px]]
|-
|Phi On Over Off|| Theta On Over Off||
|}

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Reconstructed with solenoid|| Reconstructed with no solenoid|| solenoid off divided by 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 spans 25 +/- 15

[[File:TS_NoSolGcard_022916.txt]]

[[File:TS_SolGcard_022916.txt]]
gcard for Solenoid

=== Solenoid Off===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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"

CLAS12Sim->Draw("evt.Theta >> Gen(20,15,35)");
CLAS12Sim->Draw("evt.Theta >> Rec(20,15,35)");

Gen->Draw();
Rec->Draw("same");

===scp to laptop===

scp stoktyle@minerve2.cose.isu.edu:src/CLAS/GEMC/experiments/eg12/lh2/Sim.root ./

[[EG12#Tyler_Stokes_EG12_studies]]

File:RecThetaOff1D.png

2016-04-27T09:36:30Z

Stoktyle:

File:RecPhiOff1D.png

2016-04-27T09:35:58Z

Stoktyle:

TS RunGroupC

2016-04-27T09:35:36Z

Stoktyle: /* Sweet spot */

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Detector Images
|-
|[[File:Re1eetSpotOn.png | 300 px]] ||[[File:RecSwetOff.png | 300 px]]
|-
| No Events|| Solenoid On with a proton and electron event||
|}

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

1-D Phi and Theta histograms for Solenoid Off

Overlay of Solenoid off and On for Theta and Phi.

Difference between Generated and Reconstructed Theta and Phi.

{| class="wikitable" align="center" border=1
|+ Reconstructed 1D No Solenoid
|-
|[[File:RecPhiOff1D.png | 300 px]] ||[[File:RecThetaOff1D.png | 300 px]]
|-
|Rec Phi Off|| Rec Theta Off||
|}

{| class="wikitable" align="center" border=1
|+ Overlay
|-
|[[File:OnOverOff1D.png | 300 px]] ||[[File:ThetaOnOverOff1D.png | 300 px]]
|-
|Phi On Over Off|| Theta On Over Off||
|}

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Reconstructed with solenoid|| Reconstructed with no solenoid|| solenoid off divided by 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 spans 25 +/- 15

[[File:TS_NoSolGcard_022916.txt]]

[[File:TS_SolGcard_022916.txt]]
gcard for Solenoid

=== Solenoid Off===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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"

CLAS12Sim->Draw("evt.Theta >> Gen(20,15,35)");
CLAS12Sim->Draw("evt.Theta >> Rec(20,15,35)");

Gen->Draw();
Rec->Draw("same");

===scp to laptop===

scp stoktyle@minerve2.cose.isu.edu:src/CLAS/GEMC/experiments/eg12/lh2/Sim.root ./

[[EG12#Tyler_Stokes_EG12_studies]]

TS RunGroupC

2016-04-27T09:35:06Z

Stoktyle: /* Sweet spot */

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Detector Images
|-
|[[File:Re1eetSpotOn.png | 300 px]] ||[[File:RecSwetOff.png | 300 px]]
|-
| No Events|| Solenoid On with a proton and electron event||
|}

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

1-D Phi and Theta histograms for Solenoid Off

Overlay of Solenoid off and On for Theta and Phi.

Difference between Generated and Reconstructed Theta and Phi.

{| class="wikitable" align="center" border=1
|+ Overlay
|-
|[[File:RecPhiOff1D.png | 300 px]] ||[[File:RecThetaOff1D.png | 300 px]]
|-
|Rec Phi Off|| Rec Theta Off||
|}

{| class="wikitable" align="center" border=1
|+ Overlay
|-
|[[File:OnOverOff1D.png | 300 px]] ||[[File:ThetaOnOverOff1D.png | 300 px]]
|-
|Phi On Over Off|| Theta On Over Off||
|}

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Reconstructed with solenoid|| Reconstructed with no solenoid|| solenoid off divided by 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 spans 25 +/- 15

[[File:TS_NoSolGcard_022916.txt]]

[[File:TS_SolGcard_022916.txt]]
gcard for Solenoid

=== Solenoid Off===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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"

CLAS12Sim->Draw("evt.Theta >> Gen(20,15,35)");
CLAS12Sim->Draw("evt.Theta >> Rec(20,15,35)");

Gen->Draw();
Rec->Draw("same");

===scp to laptop===

scp stoktyle@minerve2.cose.isu.edu:src/CLAS/GEMC/experiments/eg12/lh2/Sim.root ./

[[EG12#Tyler_Stokes_EG12_studies]]

File:ThetaOnOverOff1D.png

2016-04-27T09:30:44Z

Stoktyle:

File:OnOverOff1D.png

2016-04-27T09:30:14Z

Stoktyle:

TS RunGroupC

2016-04-27T09:29:53Z

Stoktyle: /* Sweet spot */

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Detector Images
|-
|[[File:Re1eetSpotOn.png | 300 px]] ||[[File:RecSwetOff.png | 300 px]]
|-
| No Events|| Solenoid On with a proton and electron event||
|}

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

1-D Phi and Theta histograms for Solenoid Off

Overlay of Solenoid off and On for Theta and Phi.

Difference between Generated and Reconstructed Theta and Phi.

{| class="wikitable" align="center" border=1
|+ Overlay
|-
|[[File:OnOverOff1D.png | 300 px]] ||[[File:ThetaOnOverOff1D.png | 300 px]]
|-
|Phi on over off|| Theta on over off||
|}

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Reconstructed with solenoid|| Reconstructed with no solenoid|| solenoid off divided by 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 spans 25 +/- 15

[[File:TS_NoSolGcard_022916.txt]]

[[File:TS_SolGcard_022916.txt]]
gcard for Solenoid

=== Solenoid Off===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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"

CLAS12Sim->Draw("evt.Theta >> Gen(20,15,35)");
CLAS12Sim->Draw("evt.Theta >> Rec(20,15,35)");

Gen->Draw();
Rec->Draw("same");

===scp to laptop===

scp stoktyle@minerve2.cose.isu.edu:src/CLAS/GEMC/experiments/eg12/lh2/Sim.root ./

[[EG12#Tyler_Stokes_EG12_studies]]

TS RunGroupC

2016-04-13T07:52:16Z

Stoktyle: /* Detector Images */

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Detector Images
|-
|[[File:Re1eetSpotOn.png | 300 px]] ||[[File:RecSwetOff.png | 300 px]]
|-
| No Events|| Solenoid On with a proton and electron event||
|}

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Reconstructed with solenoid|| Reconstructed with no solenoid|| solenoid off divided by 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 spans 25 +/- 15

[[File:TS_NoSolGcard_022916.txt]]

[[File:TS_SolGcard_022916.txt]]
gcard for Solenoid

=== Solenoid Off===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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]]

TS RunGroupC

2016-04-13T07:50:55Z

Stoktyle: /* Sweet spot */

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Detector Images
|-
|[[File:Re1eetSpotOn.png | 300 px]] ||[[File:RecSwetOff.png | 300 px]]
|-
| ||
|}

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Reconstructed with solenoid|| Reconstructed with no solenoid|| solenoid off divided by 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 spans 25 +/- 15

[[File:TS_NoSolGcard_022916.txt]]

[[File:TS_SolGcard_022916.txt]]
gcard for Solenoid

=== Solenoid Off===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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]]

File:RecSwetOff.png

2016-04-13T07:34:07Z

Stoktyle:

File:Re1eetSpotOn.png

2016-04-13T07:33:35Z

Stoktyle:

TS RunGroupC

2016-04-13T07:29:56Z

Stoktyle:

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Detector Images
|-
|[[File:Re1eetSpotOn.png | 300 px]] ||[[File:RecSwetOff.png | 300 px]]
|-
| ||
|}

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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]]

TS RunGroupC

2016-04-13T07:29:08Z

Stoktyle:

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Detector Images
|-
|[[File:Re1eetSpotOn.png | 300 px]] ||[[File:RecSwetOff.png | 300 px]] || [[File:SweeotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
| ||
|}

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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]]

TS RunGroupC

2016-04-13T07:26:59Z

Stoktyle:

=Detector Images=
{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpon.png | 300 px]] ||[[File:RecSweet.png | 300 px]] ||}
=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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]]

File:SweetSpotEff.png

2016-04-13T07:24:12Z

Stoktyle: uploaded a new version of "File:SweetSpotEff.png"

File:RecSweetSpotOn.png

2016-04-13T07:16:38Z

Stoktyle: uploaded a new version of "File:RecSweetSpotOn.png"

File:SweetSpotEff.png

2016-04-13T07:14:25Z

Stoktyle:

File:RecSweetSpotOff.png

2016-04-13T07:13:52Z

Stoktyle:

File:RecSweetSpotOn.png

2016-04-13T07:13:23Z

Stoktyle:

TS RunGroupC

2016-04-13T07:12:54Z

Stoktyle: /* Sweet spot */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]] ||[[File:Phi.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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 -config DATA::mc=true -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=nomina -config DATA::mc=true -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]]

TS RunGroupC

2016-04-11T10:49:44Z

Stoktyle: /* Solenoid On */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 80 Bins, 20 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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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]]

File:DivThetaPhiSol.png

2016-04-11T10:48:51Z

Stoktyle: uploaded a new version of "File:DivThetaPhiSol.png"

TS RunGroupC

2016-04-11T10:19:19Z

Stoktyle: /* Solenoid Off */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 200 Bins, 40 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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]]

File:ThetaPhiDivNoSol.png

2016-04-11T10:18:24Z

Stoktyle: uploaded a new version of "File:ThetaPhiDivNoSol.png"

File:RecOffRecOnZoomed.png.png

2016-04-06T19:15:25Z

Stoktyle:

File:RecOFfDivRecOnUnzoomed.png

2016-04-06T19:14:48Z

Stoktyle:

TS RunGroupC

2016-04-06T19:14:19Z

Stoktyle: /* Ratio Solenoid Off / On */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 150 Bins, 15.0 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:RecOFfDivRecOnUnzoomed.png | 300 px]] ||[[File:RecOffRecOnZoomed.png.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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]]

TS RunGroupC

2016-04-06T19:10:03Z

Stoktyle: /* Ratio Solenoid Off / On */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 150 Bins, 15.0 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.

{| class="wikitable" align="center" border=1
|+ Ratio Off/On
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]]
|-
|Bins 400 x 50. || Bins 3 x 1 .||
|}

[[File:NSRecDivRec.png]] [[File:RecOffRecOnZoomed.png]] This image has the following bins: 3 , 1

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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]]

TS RunGroupC

2016-04-06T19:06:49Z

Stoktyle: /* Ratio Solenoid Off / On */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 150 Bins, 15.0 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.
[[File:NSRecDivRec.png]] [[File:RecOffRecOnZoomed.png]] This image has the following bins: 3 , 1

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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]]

File:RecOffRecOnZoomed.png

2016-04-06T19:05:32Z

Stoktyle:

TS RunGroupC

2016-04-06T19:05:09Z

Stoktyle: /* Ratio Solenoid Off / On */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 150 Bins, 15.0 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.
[[File:NSRecDivRec.png]] [[File:RecOffRecOnZoomed.png]]

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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]]

File:NSRecDivRec.png

2016-04-06T18:41:46Z

Stoktyle: uploaded a new version of "File:NSRecDivRec.png"

No solenoid reconstructed divided by reconstructed with solenoid.

TS RunGroupC

2016-04-06T18:38:56Z

Stoktyle: /* Solenoid On */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 150 Bins, 15.0 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.
[[File:NSRecDivRec.png]]

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. ||
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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]]

TS RunGroupC

2016-04-06T18:38:13Z

Stoktyle: /* Solenoid On */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 150 Bins, 15.0 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.
[[File:NSRecDivRec.png]]

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]|
|-
|Generated with solenoid. || Reconstructed with solenoid.|| reconstructed divided by generated with solenoid. 16 Bins, 4 Bins. |
|}

==With Beam rastering==

=Commands=
==Simulation commands==

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

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]]

File:RecDivGen.png

2016-04-06T18:36:46Z

Stoktyle: uploaded a new version of "File:RecDivGen.png"

File:Recon Elastic.png

2016-04-06T17:52:24Z

Stoktyle: uploaded a new version of "File:Recon Elastic.png"

File:Gen Elastic.png

2016-04-06T17:51:26Z

Stoktyle: uploaded a new version of "File:Gen Elastic.png"

TS RunGroupC

2016-04-05T10:06:23Z

Stoktyle: /* Solenoid On */

=Inclusive electron scattering=

==Sweet spot==

Send electrons only at the center of a single sector.

<option name="BEAM_P" value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>

<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>

==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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 150 Bins, 15.0 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.

directly divide reconstructed solenoid off by reconstructed solenoid on

No solenoid reconstructed divided by solenoid reconstructed.
[[File:NSRecDivRec.png]]

directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)

===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 ===
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]||[[File:Efficiency_Elastic.png | 300 px]]
|-
|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==

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

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]]

TS RunGroupC

2016-04-03T11:21:38Z

Stoktyle: /* Solenoid Off */

=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===

Question for Same Number of Events ran:Why is their a discrepancy between the number of events in solenoid vs no solenoid? I ran 975,000 events for both.

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 150 Bins, 15.0 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.
[[File: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 ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]||[[File:Efficiency_Elastic.png | 300 px]]
|-
|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==

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

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]]

TS RunGroupC

2016-04-03T11:20:39Z

Stoktyle: /* Solenoid Off */

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

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|Generated with no solenoid. 824,979 Events || Reconstructed with no solenoid. 609,106 Events|| reconstructed divided by generated with no solenoid 150 Bins, 15.0 Bins. || Zoomed in version of the efficiency.
|}

===Solenoid On ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.
[[File: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 ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]||[[File:Efficiency_Elastic.png | 300 px]]
|-
|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==

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

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]]

File:NSEffZoom.png

2016-04-03T11:19:30Z

Stoktyle: uploaded a new version of "File:NSEffZoom.png"

File:ThetaPhiDivNoSol.png

2016-04-03T11:18:49Z

Stoktyle: uploaded a new version of "File:ThetaPhiDivNoSol.png"

File:DivThetaPhiSol.png

2016-04-03T11:17:56Z

Stoktyle: uploaded a new version of "File:DivThetaPhiSol.png": Reverted to version as of 23:11, 2 April 2016

File:SolEffZoom.png

2016-04-03T11:17:15Z

Stoktyle: uploaded a new version of "File:SolEffZoom.png": Reverted to version as of 07:31, 15 March 2016

File:SolEffZoom.png

2016-04-03T11:16:38Z

Stoktyle: uploaded a new version of "File:SolEffZoom.png"

File:DivThetaPhiSol.png

2016-04-03T11:12:11Z

Stoktyle: uploaded a new version of "File:DivThetaPhiSol.png"

TS RunGroupC

2016-04-03T10:43:30Z

Stoktyle: /* Solenoid Off */

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

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiNoSolGen.png | 300 px]] ||[[File:ThetaPhiRecNoSol.png | 300 px]] || [[File:ThetaPhiDivNoSol.png | 300 px]] ||[[File:NSEffZoom.png | 300 px]]
|-
|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 ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:ThetaPhiGenSol.png | 300 px]] ||[[File:ThetaPhiRecSol.png | 300 px]] || [[File:DivThetaPhiSol.png | 300 px]]||[[File:SolEffZoom.png | 300 px]]
|-
|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.
[[File: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 ===

{| class="wikitable" align="center" border=1
|+ Reconstruction efficiency
|-
|[[File:Gen_Elastic.png | 300 px]] ||[[File:Recon_Elastic.png | 300 px]] || [[File:RecDivGen.png | 300 px]]||[[File:Efficiency_Elastic.png | 300 px]]
|-
|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==

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

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]]

File:ThetaPhiDivNoSol.png

2016-04-03T10:42:07Z

Stoktyle: uploaded a new version of "File:ThetaPhiDivNoSol.png"