Difference between revisions of "LB RunGroupC Vertex"

From New IAC Wiki
Jump to navigation Jump to search
Line 142: Line 142:
 
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.   
 
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 line target in X
+
 
 +
 
 +
 
 +
Below are the histograms for the 2cm line target in X
  
 
{|border="5"
 
{|border="5"

Revision as of 18:38, 27 July 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.2731 +/- 0.0010 0.2702 +/- 0.0009 0.2736 +/- 0.0010
Y Resolution (cm) 0.2041 +/- 0.0010 0.2282 +/- 0.0011 ** 0.433 +/- 0.002 0.5584 +/- 0.0018 0.6454 +/- 0.0015
Z Resolution (cm) 0.28 +/-0.00 0.3582 +/- 0.0013 ** 0.7927 +/- 0.0023 1.034 +/- 0.003 1.253 +/-0.003
X Centroid (cm) 0.0139 +/- 0.001 0.0139 +/- 0.00107 ** 0.01217 +/- 0.00107 0.00991 +/- 0.00104 0.009022 +/- 0.00105
Y Centroid (cm) -0.1893 +/- 0.0011 -0.2004 +/- 0.0012 ** -0.2936 +/- 0.0022 -0.3648 +/- 0.0027 -0.3871 +/- 0.0026
Z Centroid (cm) 0.0928 +/- 0.0014 0.06464 +/- 0.0013 ** -0.1582 +/- 0.0040 -0.3016 +/- 0.0053 -0.486 +/- 0.006
Generated Vx (cm) 05cmZTarg GenPartx.png 1cmZTarg GenPartx.png ** 3cmZTarg GenPartx.png 4cmZTarg GenPartx.png 5cmZTarg GenPartx.png
Generated Vy (cm) 0.5ZTarg GenPartY.png 1cmZTarg GenParty.png ** 3cmZTarg GenParty.png 4cmZTarg GenPartY.png 5cmZTarg GenPartY.png
Generated Vz (cm) 05ZTarg GenPartz.png 1cmZTarg GenPartz.png ** 3cmZTarg GenPartz.png 4cmZTarg GenPartz.png 5cmZTarg GenPartz.png
Reconstructed Vx (cm) 0.5ZTarg Eventx.png 1cmZTarg Eventx.png ** 3cmZTarg Eventx.png 4cmZTarg Eventx.png 5cmZTarg Eventx.png
Reconstructed Vy (cm) 0.5ZTarg Eventy.png 1cmZTarg Eventy.png ** 3cmZTarg Eventy.png 4cmZTarg Eventy.png 5cmZTarg EventY.png
Reconstructed Vz (cm) 0.5ZTarg Eventz.png 1cmZTarg Eventz.png ** 3cmZTarg Eventz.png 4cmZTarg Eventz.png 5cmZTarg Eventz.png
Vx Difference (cm) 0.5ZTarg VxDiff.png 1cmZTarg VxDiff.png ** 3cmZTarg VxDiff.png 4cmZTarg VxDiff.png 5cmZTarg VxDiff.png
Vy Difference (cm) 0.5ZTarg VyDiff.png 1cmZTarg VyDiff.png ** 3cmZTarg VyDiff.png 4cmZTarg VyDiff.png 5cmZTarg VyDiff.png
Vz Difference (cm) 0.5ZTarg VzDiff.png 1cmZTarg 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.

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 X

[-1.0 , -0.9] [-0.9 , -0.8] [-0.8 , -0.7] [-0.7 , -0.6] [-0.6 , -0.5] [-0.5 , -0.4] [-0.4 , -0.3]
Resolution (cm) 0.2512 +/- 0.0042 0.3018 +/- 0.0049 0.2873 +/- 0.0049 0.2475 +/- 0.0038 0.2398 +/- 0.0038 0.2203 +/- 0.0035 0.2178 +/- 0.0034
Vx Difference VxDiff --1.0,-0.9-.png VxDiff --0.9,-0.8-.png VxDiff--0.8,-0.7-.png VxDiff--0.7,-0.6-.png VxDiff--0.6,-0.5-.png VxDiff--0.5,-0.4-.png VxDiff--0.4,-0.3-.png
* * * * * * * *
[-0.3 , -0.2] [-0.2 , -0.1] [-0.1 , 0.0] [0.0 , 0.1] [0.1 , 0.2] [0.2 , 0.3] [0.3 , 0.4]
Resolution (cm) 0.227 +/- 0.003 0.225 +/- 0.003 0.2357 +/- 0.0035 0.2256 +/- 0.0035 0.234 +/- 0.003 0.254 +/- 0.004 0.3227 +/- 0.0049
Vx Difference VxDiff--0.3,-0.2-.png VxDiff--0.2,-0.1-.png VxDiff--0.1,0-.png VxDiff-0,0.1-.png VxDiff-0.1,0.2-.png VxDiff-0.2,0.3-.png VxDiff-0.3,0.4-.png


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