Difference between revisions of "TS RunGroupC"

From New IAC Wiki
Jump to navigation Jump to search
 
(81 intermediate revisions by 2 users not shown)
Line 1: Line 1:
 +
=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=
 
=Inclusive electron scattering=
  
==Theta and Phi graphs==
+
==Sweet spot==
  
===Solenoid off ===
+
Send electrons only at the center of a single sector.
[[File:TS_11162015_GenTheta.png]][[File:TS_11162015_Theta.png]]
 
  
 +
<option name="BEAM_P"  value="e-, 6.0*GeV, 25.0*deg, 0*deg"/>
  
===Solenoid On===
+
<option name="SPREAD_P" value="4.0*GeV, 0*deg, 0*deg"/>
  
==Inclusive Track Reconstruction efficiencies==
+
 +
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||
 +
|}
  
Plot the reconstruction efficiency for an electron as a function of the angles Theta and Phi
 
  
I used the following:
+
{| class="wikitable" align="center" border=1
UBUNTU 14, coat-java-1.0
+
  |+ Overlay
 +
|-
 +
|[[File:OnOverOff1D.png | 300 px]] ||[[File:ThetaOnOverOff1D.png | 300 px]]
 +
|-
 +
|Phi On Over Off|| Theta On Over Off||
 +
|}
  
No Solenoid gcard:
+
zoom in 24 < theta < 26 and -15 < phi < 15
  
<pre>
+
{| class="wikitable" align="center" border=1
<!-- minumal detectors for HB tracking -->
+
  |+ Reconstruction efficiency
 +
|-
 +
|[[File:RecSweetSpotOn.png | 300 px]] ||[[File:RecSweetSpotOff.png | 300 px]] || [[File:SweetSpotEff.png | 300 px]]
 +
|-
 +
| Reconstructed with no solenoid|| Reconstructed with solenoid|| solenoid off divided by solenoid on||
 +
|}
  
        <detector name="../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
+
zoom in like above and do efficiency plots like below for on/off solenoid
        <detector name="../../../../clas12/dc/dc"            factory="TEXT" variation="original"/>
 
        <detector name="../../../../clas12/ftof/ftof"            factory="TEXT" variation="original"/>
 
        <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
 
        <option name="OUTPUT" value="evio,eg12.ev"/>
 
  
        <!-- Solenoid apparatus and field -->
+
==Inclusive Track Reconstruction efficiencies==
  
 +
Plot the reconstruction efficiency for an electron as a function of the angles Theta and Phi
  
        <!--detector name="../../../../clas12/magnets/solenoid" factory="TEXT" variation="original"/-->
+
I used the following:
          <!--option name="HALL_FIELD"  value="clas12-solenoid"/-->
+
UBUNTU 14, coat-java-1.0
  
<!-- other CLAS detectors -->
+
No Solenoid gcard:
  
c      <detector name="../../../../clas12/ec/ec"            factory="TEXT" variation="original"/>
+
  angle spans 25 +/- 15
      <detector name="../../../../clas12/ctof/ctof"            factory="TEXT" variation="original"/>
 
      <detector name="../../../../clas12/htcc/htcc"            factory="TEXT" variation="original"/>
 
      <detector name="../../../../clas12/pcal/pcal"            factory="TEXT" variation="javageom"/>
 
      <detector name="../../../../clas12/micromegas/micromegas" factory="TEXT" variation="original"/>
 
  
 +
[[File:TS_NoSolGcard_022916.txt]]
  
<!-- Beamline material -->
 
        <!--detector name="../../../../clas12/beamline/beamline" factory="TEXT" variation="ft"/-->
 
  
 +
[[File:TS_SolGcard_022916.txt]]
 +
gcard for Solenoid
  
          <option name="BEAM_P"  value="e-, 4.0*GeV, 20.0*deg, 0*deg"/>
+
=== Solenoid Off===
          <option name="SPREAD_P" value="0.0*GeV, 20*deg, 180*deg"/>
 
          <!-- option name="RANDOM" value="TIME"/ -->
 
</gcard>
 
  
 +
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.
  
  
<!-- minumal detectors for HB tracking -->
+
{| 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.
 +
|}
  
        <detector name="../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
+
===Solenoid On ===
        <detector name="../../../../clas12/dc/dc"            factory="TEXT" variation="original"/>
 
        <detector name="../../../../clas12/ftof/ftof"            factory="TEXT" variation="original"/>
 
        <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
 
        <option name="OUTPUT" value="evio,eg12.ev"/>
 
  
 +
{| 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.
 +
|}
  
</pre>
+
=== Ratio Solenoid Off / On ===
  
gcard for Solenoid
+
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.
  
<pre>
+
directly divide reconstructed solenoid off by reconstructed solenoid on
<!-- minumal detectors for HB tracking -->
 
  
        <detector name="../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="origi\
 
\
 
nal"/>
 
        <detector name="../../../../clas12/dc/dc"            factory="TEXT" variation="origina\
 
\
 
l"/>
 
        <detector name="../../../../clas12/ftof/ftof"            factory="TEXT" variation="ori\
 
\
 
ginal"/>
 
        <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
 
        <option name="OUTPUT" value="evio,eg12.ev"/>
 
  
        <!-- Solenoid apparatus and field -->
+
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 .||
 +
|}
  
        <!detector name="../../../../clas12/magnets/solenoid" factory="TEXT" variation="original"/>
 
          <option name="HALL_FIELD"  value="clas12-solenoid"/>
 
  
<!-- other CLAS detectors -->
+
directly divide reconstructed solenoid off by reconstructed solenoid on that has been multiplied by (825/150)
  
c      <detector name="../../../../clas12/ec/ec"            factory="TEXT" variation="original"/>
+
===Phi angle shift===
      <detector name="../../../../clas12/ctof/ctof"            factory="TEXT" variation="original"/>
 
      <detector name="../../../../clas12/htcc/htcc"            factory="TEXT" variation="original"/>
 
      <detector name="../../../../clas12/pcal/pcal"            factory="TEXT" variation="javageom"/>
 
      <detector name="../../../../clas12/micromegas/micromegas"  factory="TEXT" variation="original"/>
 
  
 +
Phi angle Solenoid on - off difference
  
<!-- Beamline material -->
+
spread from 0 to 5 GeV.
        <!--detector name="../../../../clas12/beamline/beamline" factory="TEXT" variation="ft"/-->
 
  
 +
=== Observations===
  
          <option name="BEAM_P"  value="e-, 4.0*GeV, 20.0*deg, 0*deg"/>
+
Why does the Solenoid off events only extend to 25 degrees whereas the solenoid on go to 40.
          <option name="SPREAD_P" value="0.0*GeV, 20*deg, 180*deg"/>
 
          <!-- option name="RANDOM" value="TIME"/ -->
 
</gcard>
 
  
 +
=Vertex Reconstruction=
  
 +
== Elastic Electrons from origin==
  
<!-- minumal detectors for HB tracking -->
 
  
        <detector name="../../../../clas12/fc/forwardCarriage" factory="TEXT" variation="original"/>
+
Efficiency plots for Elasticlly scattered electrons when incident electron energy was 11 GeV
        <detector name="../../../../clas12/dc/dc"            factory="TEXT" variation="original"/>
 
        <detector name="../../../../clas12/ftof/ftof"            factory="TEXT" variation="original"/>
 
        <option name="SCALE_FIELD" value="clas12-torus-big, -1"/>
 
        <option name="OUTPUT" value="evio,eg12.ev"/>
 
  
  
 
+
using the file ReconElasticExclusive_FromOrigin.LUND
</pre>
 
 
 
=== Solenoid Off===
 
Generated with no solenoid.
 
[[File:ThetaPhiNoSolGen.png]]
 
Reconstructed with no solenoid.
 
[[File:ThetaPhiRecNoSol.png]]
 
reconstructed divided by generated with no solenoid.
 
[[File:ThetaPhiDivNoSol.png]]
 
  
 
===Solenoid On ===
 
===Solenoid On ===
 +
nohup gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,ReconElasticExclusive_FromOrigin.LUND" -N=75000 shift.gcard > /dev/null &
  
Generated with solenoid.
 
[[File:ThetaPhiGenSol.png]]
 
Reconstructed with solenoid.
 
[[File:ThetaPhiRecSol.png]]
 
Reconstructed divided by generated with solenoid.
 
[[File:DivThetaPhiSol.png]]
 
 
=== Ratio Solenoid Off / On ===
 
 
No solenoid reconstructed divided by solenoid reconstructed.
 
[[File:NSRecDivRec.png]]
 
 
===Phi angle shift===
 
  
Phi angle Solenoid on - off difference
+
{| 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. ||
 +
|}
  
spread from 0 to 5 GeV.
+
== Elastics from Dual Target==
  
=== Observations===
+
  Generated and reconstructed separate Plots of  sqrt(Vx^2+Vy^2) -vs- Vz
  
Why does the Solenoid off events only extend to 25 degrees whereas the solenoid on go to 40.
 
  
Insert the two cards and normalized graphs.
 
  
=Vertex Reconstruction=
 
 
== For particle at origin==
 
 
plot Generated xpos - Reconstructed x position
 
  
  
Line 180: Line 171:
 
without solenoid
 
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
+
  ~/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
 
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
+
  ~/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
 
evio analysis program
Line 251: Line 242:
 
root [4] TH1D *RecPhi=new TH1D("RecPhi","RecPhi",400, -200,200);
 
root [4] TH1D *RecPhi=new TH1D("RecPhi","RecPhi",400, -200,200);
 
root [5] CLAS12Sim->Draw("evt.Phi>>RecPhi");
 
root [5] CLAS12Sim->Draw("evt.Phi>>RecPhi");
 +
 +
root [5] CLAS12Sim->Draw("evt.Phi>>RecPhi","10");
  
 
go to Generated Events File
 
go to Generated Events File
Line 261: Line 254:
 
root [6] RecPhi->Divide(GenPhi);
 
root [6] RecPhi->Divide(GenPhi);
 
root [7] RecPhi->Draw();
 
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 to laptop===

Latest revision as of 19:33, 27 April 2016

Detector Images

Detector Images
Re1eetSpotOn.png RecSwetOff.png
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"/>


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.

Reconstructed 1D No Solenoid
RecPhiOff1D.png RecThetaOff1D.png
Rec Phi Off Rec Theta Off


Overlay
OnOverOff1D.png ThetaOnOverOff1D.png
Phi On Over Off Theta On Over Off
zoom in 24 < theta < 26 and -15 < phi < 15
Reconstruction efficiency
RecSweetSpotOn.png RecSweetSpotOff.png SweetSpotEff.png
Reconstructed with no solenoid Reconstructed with solenoid solenoid off divided by solenoid on
zoom in like above and do efficiency plots like below for on/off solenoid

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.


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

Solenoid On

Reconstruction efficiency
ThetaPhiGenSol.png ThetaPhiRecSol.png DivThetaPhiSol.png SolEffZoom.png
Generated with solenoid. Reconstructed with solenoid. reconstructed divided by generated with solenoid. 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.

Ratio Off/On
RecOFfDivRecOnUnzoomed.png RecOffRecOnZoomed.png.png
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 &


Reconstruction efficiency
Gen Elastic.png Recon Elastic.png RecDivGen.png|
Generated with solenoid. Reconstructed with solenoid. reconstructed divided by generated with solenoid. 16 Bins, 4 Bins.

Elastics from Dual Target

 Generated and reconstructed separate Plots of  sqrt(Vx^2+Vy^2) -vs- Vz



With Beam rastering

Commands

Simulation commands

tcsh
cd ~/src/CLAS/GEMC
source setup
cd experiments/eg12/lh2
gemc -HALL_MATERIAL="Vacuum" eg12.gcard
gemc -USE_GUI=0 -HALL_MATERIAL="Vacuum" -N=10 eg12.gcard

to edit files

emacs -nw filename

without solenoid

~/src/CLAS/coatjava-1.0/bin/clas12-reconstruction -config DCHB::torus=1.0 -config DCHB::solenoid=0.0 -config DCTB::kalman=true  -s DCHB:DCTB:EC:FTOF:EB -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");

root [5] CLAS12Sim->Draw("evt.Phi>>RecPhi","10");

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