Difference between revisions of "LB RunGroupC Vertex"

From New IAC Wiki
Jump to navigation Jump to search
Line 291: Line 291:
 
|| GenPart.x || GenPart.y || GenPart.z || Event.x || Event.y || Event.z  
 
|| GenPart.x || GenPart.y || GenPart.z || Event.x || Event.y || Event.z  
 
|-
 
|-
−
|| [[ File:GenPart.x 270DegRot.png|200px]] || [[ File:GenPart.y 270DegRot.png|200px]] || [[File:GenPart.z 270DegRot.png|200px]] || [[File:Event.x 270DegRot.png|200px]] || [[File:Event.y 270DegRot.png|200px]] || [[File:Event.z 270DegRot.png|200px]]
+
|| || || || || ||
 
|-
 
|-
 
|}<br>
 
|}<br>

Revision as of 19:06, 8 August 2016

Vertex Reconstruction studies

Inclusive electron GEMC 2.4 & Coatjava 2.4

LB_GEMC2.4CJ2.4_Vertex



Summary

Set all histogram ranges to -6,6 cm  and bin sizes of 0.1 cm

Point target X,Y,Z=0 cm

No Rastering
Vx Difference NoRaster VxDifference sol.png
Vy Difference NoRaster VyDifference sol.png
Vz Difference NoRaster VzDifference sol.png
X Resolution (cm) 0.2362 +/- 0.0008
Y Resolution(cm) 0.09809 +/- 0.00033
Z Resolution (cm) 0.05435 +/- 0.00017


These histograms raise some questions. To begin I created a LUND file with 6 GeV incident electrons at 25 degrees in theta and 0 degrees in phi. All of the vertex positions were set to 0. I then ran GEMC 2.4 using the command line ~/src/CLAS/GEMC/source/gemc -USE_GUI=0 -INPUT_GEN_FILE="LUND,No_Raster.LUND" -N=75000 eg12_sol_No_Raster.gcard

which created an output file called eg12_sol_75k_No_Raster.ev. After that the reconstruction command line used was

~/src/CLAS/coatjava-2.4/bin/clas12-reconstruction -i eg12_sol_75k_No_Raster.ev-config GEOM::new=true -config MAG::torus=-1.0 -config MAG::solenoid=1.0 -o eg12_sol_75k_No_Raster_rec.evio -s DCHB:DCTB:EC:FTOF:EB -config DATA::mc=true -config DCTB::useRaster=true

After the reconstruction a root file was created using

~/src/CLAS/evio2root/bin/evio2root eg12_sol_75k_No_Raster_rec.evio eg12_sol_75k_No_Raster_rec.root 75000

A plot was then created to show the X Vertex and Y Vertex Differences

NoRaster Vx YDifference.png

The first question is why is there a shift in the Y Difference. What would make X more centered around 0? The next question is when rastering begins, why do the resolutions in the X Vertices remain unchanged while the resolutions in the Y Vertices change?

Point in 2-D but extended target in 1-D

Two vertex variables are fixed at zero and the third vertex variable is altered

2cm X Raster, Y=0,Z=0 2 cm Y Raster, X=0,Z=0 2 cm Z Raster ,X=0,Y=0 6 cm Z Raster ,X=0,Y=0
Vx Difference 1cm XRaster VxDifference sol 75k.png 1cm YRaster VxDifference sol.png 2cm ZRaster VxDifference sol.png 3cm ZRaster VxDifference sol.png
Vy Difference 1cm XRaster VyDifference sol.png YRaster Eventy.png 2cm ZRaster VyDifference sol.png 3cm ZRaster VyDifference sol.png
Vz Difference 1cm XRaster VzDifference sol.png 1cm YRaster VzDifference sol.png 2cm ZRaster VzDifference.png 3cm ZRaster VZDifference sol 75k.png
X Resolution (cm) 0.2658 +/- 0.0010 0.2531 +/- 0.0009 0.2544 +/- 0.0009 0.2539 +/- 0.0009
Y Resolution(cm) 0.1728 +/- 0.0007 0.2368 +/- 0.0007 0.2644 +/- 0.0007 0.6191 +/- 0.0013
Z Resolution (cm) 0.1792 +/- 0.0006 0.1998 +/- 0.0004 0.4318 +/- 0.0009 1.105 +/- 0.002


Extended target -3 < Z < 3 cm

1cm X Raster 1 cm Y Raster 1 cm X & Y Raster
X Resolution (cm) 0.2696 +/- 0.0009 0.2588 +/- 0.0009 0.2734 +/- 0.0009
Y Resolution (cm) 0.6234 +/- 0.0014 0.6851 +/- 0.0014 0.6927 +/- 0.0013
Z Resolution (cm) 1.136 +/- 0.002 1.361 +/- 0.002 1.381 +/- 0.003
Generated Vx (cm) 3cm ZRaster 1cm XRaster GenX sol.png 3cm ZRaster 1cm YRaster GenX sol.png 3cm ZRaster 1cm XYRaster sol GenPartx.png
Generated Vy (cm) 3cm ZRaster 1cm XRaster GenY sol.png 3cm ZRaster 1cm YRaster GenY sol.png 3cm ZRaster 1cm XYRaster GenParty.png
Generated Vz (cm) 3cm ZRaster 1cm XRaster GenZ sol.png 3cm ZRaster 1cm YRaster GenZ sol.png 3cm ZRaster 1cm XYRaster GenPartz sol.png
Reconstructed Vx (cm) 3cm ZRaster 1cm XRaster EventX sol.png 3cm ZRaster 1cm YRaster EventX sol.png 3cm ZRaster 1cm XYRaster Eventx.png
Reconstructed Vy (cm) 3cm ZRaster 1cm XRaster EventY sol.png 3cm ZRaster 1cm YRaster EventY sol.png 3cm ZRaster 1cm XYRaster Eventy.png
Reconstructed Vz (cm) 3cm ZRaster 1cm XRaster EventZ sol.png 3cm ZRaster 1cm YRaster EventZ sol.png 3cm ZRaster 1cm XYRaster eventz.png
Vx Difference (cm) XZRaster VxDifference 1mmbins.png YZRaster VxDifference 1mmbins.png XYZRaster VxDifference 1mmbins.png
Vy Difference (cm) XZRaster VyDifference 1mmbins.png YZRaster VyDifference 1mmbins.png XYZRaster VyDifference 1mmbins.png
Vz Difference (cm) XZRaster VzDifference 1mmbins.png YZRaster VzDifference 1mmbins.png XYZRaster VzDifference 1mmbins.png


Extended target in Z. Vy shift analysis

All Histograms below were created by simply making the target longer along the Z axis. Each case uses vertex points for X and Y ranging from -1cm to 1cm

0.5cm Z Target 1cm Z Target 2cm Z Target 3cm Z Target 4cm Z Target 5cm Z Target
X Resolution (cm) 0.2705 +/- 0.0010 0.2726 +/- 0.0010 0.2716 +/- 0.0010 0.2731 +/- 0.0010 0.2702 +/- 0.0009 0.2736 +/- 0.0010
Y Resolution (cm) 0.2041 +/- 0.0010 0.2282 +/- 0.0011 0.3152 +/- 0.0014 0.433 +/- 0.002 0.5584 +/- 0.0018 0.6454 +/- 0.0015
Z Resolution (cm) 0.28 +/-0.00 0.3582 +/- 0.0013 0.5659 +/- 0.0018 0.7927 +/- 0.0023 1.034 +/- 0.003 1.253 +/-0.003
X Mean (cm) 0.0139 +/- 0.001 0.0139 +/- 0.00107 0.01436 +/- 0.00104 0.01217 +/- 0.00107 0.00991 +/- 0.00104 0.009022 +/- 0.00105
Y Mean (cm) -0.1893 +/- 0.0011 -0.2004 +/- 0.0012 -0.2379 +/- 0.0016 -0.2936 +/- 0.0022 -0.3648 +/- 0.0027 -0.3871 +/- 0.0026
Z Mean (cm) 0.0928 +/- 0.0014 0.06464 +/- 0.0013 -0.0337 +/- 0.0028 -0.1582 +/- 0.0040 -0.3016 +/- 0.0053 -0.486 +/- 0.006
Generated Vx (cm) 05cmZTarg GenPartx.png 1cmZTarg GenPartx.png 2cmZTarg GenPartx.png 3cmZTarg GenPartx.png 4cmZTarg GenPartx.png 5cmZTarg GenPartx.png
Generated Vy (cm) 0.5ZTarg GenPartY.png 1cmZTarg GenParty.png 2cmZTarg GenPartY.png 3cmZTarg GenParty.png 4cmZTarg GenPartY.png 5cmZTarg GenPartY.png
Generated Vz (cm) 05ZTarg GenPartz.png 1cmZTarg GenPartz.png 2cmZTarg GenPartz.png 3cmZTarg GenPartz.png 4cmZTarg GenPartz.png 5cmZTarg GenPartz.png
Reconstructed Vx (cm) 0.5ZTarg Eventx.png 1cmZTarg Eventx.png 2cmZTarg Eventx.png 3cmZTarg Eventx.png 4cmZTarg Eventx.png 5cmZTarg Eventx.png
Reconstructed Vy (cm) 0.5ZTarg Eventy.png 1cmZTarg Eventy.png 2cmZTarg EventY.png 3cmZTarg Eventy.png 4cmZTarg Eventy.png 5cmZTarg EventY.png
Reconstructed Vz (cm) 0.5ZTarg Eventz.png 1cmZTarg Eventz.png 2cmZTarg Eventz.png 3cmZTarg Eventz.png 4cmZTarg Eventz.png 5cmZTarg Eventz.png
Vx Difference (cm) 0.5ZTarg VxDiff.png 1cmZTarg VxDiff.png 2cmZTarg VxDiff.png 3cmZTarg VxDiff.png 4cmZTarg VxDiff.png 5cmZTarg VxDiff.png
Vy Difference (cm) 0.5ZTarg VyDiff.png 1cmZTarg VyDiff.png 2cmZTarg VyDiff.png 3cmZTarg VyDiff.png 4cmZTarg VyDiff.png 5cmZTarg VyDiff.png
Vz Difference (cm) 0.5ZTarg VzDiff.png 1cmZTarg VzDiff.png 2cmZTarg VzDiff.png 3cmZTarg VzDiff.png 4cmZTarg VzDiff.png 5cmZTarg VzDiff.png


It would seem that increasing the length of the target in the Z direction has an effect on the resolution of the Y vertex position along with how far the center of the Gaussian peak will shift for the Y vertex differences.

Localized Vertex Resolution Investigation of Line Targets

GEMC was ran to create electrons being produced in a line along certain axes. These include three separate 2cm targets centered at (0,0,0) in the X, Y , and Z directions. The final target is a line along the Z axis centered at (0,0,0) that is 6cm long. The generated electrons were shot at 25 degrees in theta and 0 degrees in phi at an energy of 6GeV.

Here is an example of how the cuts were made. Consider the case of a 2cm line target in the Z direction centered at (0,0,0). Closed intervals begin at the end of the target and are 0.1cm long. So the first interval would be [-1.0,-0.9], then the next would be [-0.9, -0.8] and so on until the entire target is covered. The events in the EVENTHB bank were cut so that in any given interval there would be a margin of +/- 1cm to remove any outliers that may cause the histograms to have long tails.



Below are the histograms for the 2cm line target in Y,X=0,Z=0. The electrons were fired with an energy of 6 GeV at 25 degrees in theta at 0 degrees in phi.

GenPart.x GenPart.y GenPart.z Event.x Event.y Event.z
GenPart.x Local Vert.png GenPart.y Local Vert.png GenPart.z Local Vert.png Event.x Local Vert.png Event.y Local Vert.png Event.z Local Vert.png


Interval Resolution (cm) Vy Difference *** Interval Resolution (cm) Vy Difference
[-1.0,-0.9] 0.09406 +/- 0.00128 VyDiff--1.0,-0.9-.png *** [0.9,1.0] 0.1001 +/- 0.0016 VyDiff-0.9,1.0-.png
[-0.9,-0.8] 0.0868 +/- 0.0013 VyDiff--0.9,-0.8-.png *** [0.8,0.9] 0.1221 +/- 0.0018 VyDiff-0.8,0.9-.png
[-0.8,-0.7] 0.08961 +/- 0.00128 VyDiff--0.8,-0.7-.png *** [0.7,0.8] 0.1144 +/- 0.0017 VyDiff-0.7,0.8-.png
[-0.7,-0.6] 0.09853 +/- 0.00148 VyDiff--0.7,-0.6-.png *** [0.6,0.7] 0.169 +/- 0.003 VyDiff-0.6,0.7-.png
[-0.6,-0.5] 0.1038 +/- 0.0016 VyDiff--0.6,-0.5-.png *** [0.5,0.6] 0.171 +/- 0.003 VyDiff-0.5,0.6-.png
[-0.5,-0.4] 0.1202 +/- 0.0018 VyDiff--0.5,-0.4-.png *** [0.4,0.5] 0.1105 +/- 0.0016 VyDiff-0.4,0.5-.png
[-0.4,-0.3] 0.1755 +/- 0.0026 VyDiff--0.4,-0.3-.png *** [0.3,0.4] 0.09508 +/- 0.00137 VyDiff-0.3,0.4-.png
[-0.3,-0.2] 0.1847 +/- 0.0028 VyDiff--0.3,-0.2-.png *** [0.2,0.3] 0.09065 +/- 0.0140 VyDiff-0.2,0.3-.png
[-0.2,-0.1] 0.1255 +/- 0.0019 VyDiff--0.2,-0.1-.png *** [0.1,0.2] 0.093002 +/- 0.0013 VyDiff-0.1,0.2-.png
[-0.1,0.0] 0.09663 +/- 0.00160 VyDiff--0.1,0.0-.png *** [0.0,0.1] 0.09008 +/- 0.0012 VyDiff-0.0,0.1-.png


Below are the histograms for a 2cm target along the Y axis, X = 0 Z = 0. The electrons were fired with an energy of 6 GeV at 25 degrees in theta and 180 degrees in phi.

GenPart.x GenPart.y GenPart.z Event.x Event.y Event.z
GenPart.x 180DegRot.png GenPart.y 180DegRot.png GenPart.z 180DegRot.png Event.x 180DegRot.png Event.y 180DegRot.png Event.z 180DegRot.png


Interval Resolution (cm) Vy Difference *** Interval Resolution (cm) Vy Difference
[-1.0,-0.9] 0.07741 +/- 0.00166 VyDiff--1.0,-0.9- 180Deg.png *** [0.9,1.0] 0.2319 +/- 0.0189 VyDiff-0.9,1.0- 180DegRot.png
[-0.9,-0.8] 0.0976 +/- 0.0017 VyDiff--0.9,-0.8- 180Deg.png *** [0.8,0.9] 0.2147 +/- 0.0185 VyDiff-0.8,0.9- 180DegRot.png
[-0.8,-0.7] 0.09924 +/- 0.00150 VyDiff--0.8,-0.7- 180Deg.png *** [0.7,0.8] 0.2364 +/- 0.0188 VyDiff-0.7,0.8- 180Deg.png
[-0.7,-0.6] 0.1386 +/- 0.0024 VyDiff--0.7,-0.6- 180Deg.png *** [0.6,0.7] 0.2323 +/- 0.0189 VyDiff-0.6,0.7- 180Deg.png
[-0.6,-0.5] 0.1278 +/- 0.0027 VyDiff--0.6,-0.5- 180Deg.png *** [0.5,0.6] 0.06016 +/- 0.00299 VyDiff-0.5,0.6- 180Deg.png
[-0.5,-0.4] 0.09882 +/- 0.00142 VyDiff--0.5,-0.4- 180Deg.png *** [0.4,0.5] 0.3986 +/- 0.00104 VyDiff-0.4,0.5- 180Deg.png
[-0.4,-0.3] 0.09322 +/- 0.0013 VyDiff--0.4,-0.3- 180Deg.png *** [0.3,0.4] 0.3056 +/- 0.0047 VyDiff-0.3,0.4- 180Deg.png
[-0.3,-0.2] 0.09399 +/- 0.00145 VyDiff--0.3,-0.2- 180Deg.png *** [0.2,0.3] 0.1097 +/- 0.0029 VyDiff-0.2,0.3- 180Deg.png
[-0.2,-0.1] 0.08935 +/- 0.00131 VyDiff--0.2,-0.1- 180Deg.png *** [0.1,0.2] 0.1146 +/- 0.0025 VyDiff-0.1,0.2- 180Deg.png
[-0.1,0.0] 0.08913 +/- 0.00130 VyDiff--0.1,0.0- 180Deg.png *** [0.0,0.1] 0.09557 +/- 0.0015 VyDiff-0.0,0.1- 180Deg.png



Below are the Histograms for a 2cm long target in the X Direction with a 6GeV electron fired at 25 degrees in theta and 90 degrees in phi

GenPart.x GenPart.y GenPart.z Event.x Event.y Event.z
GenPart.x 90DegRot.png GenPart.y 90DegRot.png GenPart.z 90DegRot.png Event.x 90DegRot.png Event.y 90DegRot.png Event.z 90DegRot.png


Interval Resolution (cm) Vx Difference *** Interval Resolution (cm) Vx Difference
[-1.0,-0.9] 0.1341 +/- 0.0140 VxDiff--1.0,-0.9-.png *** [0.9,1.0] 0.01572 +/- 0.0002 VxDiff-0.9,1.0-.png
[-0.9,-0.8] 0.1304 +/- 0.0051 VxDiff--0.9,-0.8-.png *** [0.8,0.9] 0.01844 +/- 0.00031 VxDiff-0.8,0.9-.png
[-0.8,-0.7] 0.1182 +/- 0.0036 VxDiff--0.8,-0.7-.png *** [0.7,0.8] 0.01722 +/- 0.0028 VxDiff-0.7,0.8-.png
[-0.7,-0.6] 0.09062 +/- 0.0031 VxDiff--0.7,-0.6-.png *** [0.6,0.7] 0.02386 +/- 0.00088 VxDiff-0.6,0.7-.png
[-0.6,-0.5] 0.06816 +/- 0.00191 VxDiff--0.6,-0.5-.png *** [0.5,0.6] 0.024 +/- 0.001 VxDiff-0.5,0.6-.png
[-0.5,-0.4] n/a VxDiff--0.5,-0.4-.png *** [0.4,0.5] 0.04406 +/- 0.00182 VxDiff-0.4,0.5-.png
[-0.4,-0.3] 0.0635 +/- 0.0013 VxDiff--0.4,-0.3-.png *** [0.3,0.4] 0.04517 +/- 0.00118 VxDiff-0.3,0.4-.png
[-0.3,-0.2] 0.0618 +/- 0.0012 VxDiff--0.3,-0.2-.png *** [0.2,0.3] 0.04722 +/- 0.00107 VxDiff-0.2,0.3-.png
[-0.2,-0.1] 0.07473 +/- 0.00178 VxDiff--0.2,-0.1-.png *** [0.1,0.2] 0.05667 +/- 0.00106 VxDiff-0.1,0.2-.png
[-0.1,0.0] 0.08787 +/- 0.00132 VxDiff--0.1,0.0-.png *** [0.0,0.1] 0.08356 +/- 0.00159 VxDiff-0.0,0.1-.png


Below are the Vx Differences using a 2cm target in X with incident 6 GeV electrons at 25 degrees in theta and 270 degrees in phi


GenPart.x GenPart.y GenPart.z Event.x Event.y Event.z
GenPart.x 270DegRot.png GenPart.y 270DegRot.png GenPart.z 270DegRot.png Event.x 270DegRot.png Event.y 270DegRot.png Event.z 270DegRot.png


Interval Resolution (cm) Vx Difference *** Interval Resolution (cm) Vx Difference
[-1.0,-0.9] 0.01571 +/- 0.0002 VxDiff--1.0,-0.9- 270Deg.png *** [0.9,1.0] 0.1257 +/- 0.0100 VxDiff-0.9,1.0- 270Deg.png
[-0.9,-0.8] 0.0163 +/- 0.0003 VxDiff--0.9,-0.8- 270Deg.png *** [0.8,0.9] 0.112 +/- 0.004 VxDiff-0.8,0.9- 270Deg.png
[-0.8,-0.7] 0.01821 +/- 0.0003 VxDiff--0.8,-0.7- 270Deg.png *** [0.7,0.8] 0.1162 +/- 0.0034 VxDiff-0.7,0.8- 270Deg.png
[-0.7,-0.6] 0.03065 +/- .000163 VxDiff--0.7,-0.6- 270Deg.png *** [0.6,0.7] 0.08958 +/- 0.00316 VxDiff-0.6,0.7- 270Deg.png
[-0.6,-0.5] 0.02109 +/- 0.00217 VxDiff--0.6,-0.5- 270Deg.png *** [0.5,0.6] 0.07347 +/- 0.00210 VxDiff-0.5,0.6- 270Deg.png
[-0.5,-0.4] 0.03772 +/- 0.00133 VxDiff--0.5,-0.4- 270Deg.png *** [0.4,0.5] 0.07271 +/- 0.00159 VxDiff-0.4,0.5- 270Deg.png
[-0.4,-0.3] 0.04817 +/- 0.00125 VxDiff--0.4,-0.3- 270Deg.png *** [0.3,0.4] 0.06196 +/- 0.00146 VxDiff-0.3,0.4- 270Deg.png
[-0.3,-0.2] 0.05154 +/- 0.00113 VxDiff--0.3,-0.2- 270Deg.png *** [0.2,0.3] 0.06087 +/- 0.00114 VxDiff-0.2,0.3- 270Deg.png
[-0.2,-0.1] 0.05677 +/- 0.00106 VxDiff--0.2,-0.1- 270Deg.png *** [0.1,0.2] 0.07298 +/- 0.00149 VxDiff-0.1,0.2- 270Deg.png
[-0.1,0.0] 0.08348 +/- 0.00144 VxDiff--0.1,0.0- 270Deg.png *** [0.0,0.1] 0.09017 +/- 0.00128 VxDiff-0.0,0.1- 270Deg.png


Below are the histograms for a 6cm target in the Z direction using a 6 GeV electron fired at 25 degrees in theta and 90 degrees in phi.


GenPart.x GenPart.y GenPart.z Event.x Event.y Event.z


Interval Resolution (cm) Vx Difference *** Interval Resolution (cm) Vx Difference
[-3.0,-2.5] *** [2.5,3.0]
[-2.5,-2.0] *** [2.0,2.5]
[-2.0,-1.5] *** [1.5,2.0]
[-1.5,-1.0] *** [1.0,1.5]
[-1.0,-0.5] *** [0.5,1.0]
[-0.5,0.0] *** [0.0,0.5]


Elastic electron proton GEMC 2.4 & Coatjava 2.4

References

https://clasweb.jlab.org/wiki/index.php/TF_EG12_Vertex#Z_resolution_With_micro-megas

EG12#Vertex_Studies

Retrieved from "https://wiki.iac.isu.edu/index.php?title=LB_RunGroupC_Vertex&oldid=107739"