Difference between revisions of "LB RunGroupC Vertex"
Line 613: | Line 613: | ||
For the plots below, a 6 GeV electron was fired at angles 25 degrees in theta and 30 degrees in phi. The target has an XY radius of 1cm and a length in the Z direction of 6cm. Setting GenPart.z = 0cm, the resolution was investigated as a function of the target's XY radius along with a table for the values. | For the plots below, a 6 GeV electron was fired at angles 25 degrees in theta and 30 degrees in phi. The target has an XY radius of 1cm and a length in the Z direction of 6cm. Setting GenPart.z = 0cm, the resolution was investigated as a function of the target's XY radius along with a table for the values. | ||
− | [[ File: | + | [[File:Annular Rad vs Width Z=0.png|200px]] |
clas12->Draw("Event.z >> (100,-5,5)","sqrt(GenPart.x*GenPart.x+GenPart.y*GenPart.y)<0.95 && sqrt(GenPart.x*GenPart.x+GenPart.y*GenPart.y)>0.85"); | clas12->Draw("Event.z >> (100,-5,5)","sqrt(GenPart.x*GenPart.x+GenPart.y*GenPart.y)<0.95 && sqrt(GenPart.x*GenPart.x+GenPart.y*GenPart.y)>0.85"); | ||
Line 619: | Line 619: | ||
{|border="5" | {|border="5" | ||
|- | |- | ||
− | ||R (cm) || 0. | + | ||R (cm) || 0.05 +/- 0.05 || 0.15 +/- 0.05 || 0.25 +/- 0.05 || 0.35 +/- 0.05 || 0.45 +/- 0.05 ||0.55 +/- 0.05||0.65 +/- 0.05||0.75 +/- 0.05||0.85 +/- 0.05||0.95 +/- 0.05 |
|- | |- | ||
− | ||Width (cm) ||0. | + | ||Width (cm) || 0.05094 +/- 0.0019 || 0.07812 +/- 0.0018 || 0.3688 +/- 0.0063 || 0.5231 +/- 0.0077 || 0.718 +/- 0.012 || 0.7987 +/- 0.0092 || 0.8709 +/- 0.0085 || 1.038 +/- 0.01 || 1.183 +/- 0.012 || 1.296 +/- 0.013 |
|- | |- | ||
|}<br> | |}<br> |
Revision as of 21:36, 9 September 2016
Vertex Reconstruction studies
Inclusive electron GEMC 2.4 & Coatjava 2.4
Summary
Set all histogram ranges to -6,6 cm and bin sizes of 0.1 cm
Point target X,Y,Z=0 cm
Below are examples of the reconstruction for several line targets as seen in the CED
We should change it to +30,150,210, and 330 degrees so the electron is not going in between two drift chambers but rather into the center, keep 90 and 270.
The histograms below are using a 6 GeV electron fired at an angles of 25 degrees in theta and 0 degrees in phi.
No Rastering | |
---|---|
Vx Difference | ![]() |
Vy Difference | ![]() |
Vz Difference | ![]() |
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
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. A 6 GeV electron was fired at angles of 25 degrees in theta and 0 degrees in phi.
Extended target -3 < Z < 3 cm
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
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 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
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 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
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 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
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 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Below are the histograms for a 6cm target in the Z direction using a 6 GeV electron fired at 25 degrees in theta and 0 degrees in phi.
GenPart.x | GenPart.y | GenPart.z | Event.x | Event.y | Event.z |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
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 |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Below are the histograms for a 6 GeV electron fired at angles 25 degrees in theta and 270 degrees in phi. The target was 6cm in the Z direction.
GenPart.x | GenPart.y | GenPart.z | Event.x | Event.y | Event.z |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Below are the histograms for a 6 GeV electron fired at 25 degrees in theta and 0 degrees in phi. The target is 2cm long on the Y axis.
GenPart.x | GenPart.y | GenPart.z | Event.x | Event.y | Event.z |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
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 30 degrees in phi.
GenPart.x | GenPart.y | GenPart.z | Event.x | Event.y | Event.z |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Below are the histograms for a 2cm Y target with an incident electron energy of 6 GeV. The angles were theta = 25 degrees and phi = 210 degrees
GenPart.x | GenPart.y | GenPart.z | Event.x | Event.y | Event.z |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Localized Vertex Resolution Investigation of 2 Dimensional Targets
The histograms below were created using a 6 GeV electron fired at 25 degrees in theta and 30 degrees in phi. To make the cuts I simply made circles of increasing radius centered at the origin. I used 0.1cm increments to make the 10 histograms for each vertex position.
Below are the Vx and Vy Differences for a 2cm XY target, Z=0.
GenPart.x | GenPart.y | GenPart.z | Event.x | Event.y | Event.z |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Below are the Vy and Vz differences for a 2cm target in Y and Z, X=0. The electrons were fired with an energy of 6 GeV and angles of 25 degrees in theta and 30 degrees in phi.
GenPart.x | GenPart.y | GenPart.z | Event.x | Event.y | Event.z |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Below are the Vx and Vz differences for a 2cm XZ target, Y = 0. The electrons were fired with an energy of 6 GeV at angles theta = 25 degrees and phi = 30 degrees.
GenPart.x | GenPart.y | GenPart.z | Event.x | Event.y | Event.z |
![]() |
![]() |
![]() |
![]() |
![]() |
![]() |
Localized Vertex Resolution of 3 Dimensional targets
For the plots below, a 6 GeV electron was fired at angles 25 degrees in theta and 30 degrees in phi. The target has an XY radius of 1cm and a length in the Z direction of 6cm. Setting GenPart.z = 0cm, the resolution was investigated as a function of the target's XY radius along with a table for the values.
clas12->Draw("Event.z >> (100,-5,5)","sqrt(GenPart.x*GenPart.x+GenPart.y*GenPart.y)<0.95 && sqrt(GenPart.x*GenPart.x+GenPart.y*GenPart.y)>0.85");
R (cm) | 0.05 +/- 0.05 | 0.15 +/- 0.05 | 0.25 +/- 0.05 | 0.35 +/- 0.05 | 0.45 +/- 0.05 | 0.55 +/- 0.05 | 0.65 +/- 0.05 | 0.75 +/- 0.05 | 0.85 +/- 0.05 | 0.95 +/- 0.05 |
Width (cm) | 0.05094 +/- 0.0019 | 0.07812 +/- 0.0018 | 0.3688 +/- 0.0063 | 0.5231 +/- 0.0077 | 0.718 +/- 0.012 | 0.7987 +/- 0.0092 | 0.8709 +/- 0.0085 | 1.038 +/- 0.01 | 1.183 +/- 0.012 | 1.296 +/- 0.013 |
Here is the Radius vs Width plot for when abs(GenPart.z) <= 0.5cm along with a table of values.
R (cm) | 0.1 | 0.2 | 0.3 | 0.4 | 0.5 | 0.6 | 0.7 | 0.8 | 0.9 | 1.0 |
Width (cm) | 0.345 +/- 0.032 | 0.2889 +/- 0.0101 | 0.2704 +/- 0.0053 | 0.2821 +/- 0.0042 | 0.2782 +/- 0.0035 | 0.287 +/- 0.003 | 0.294 +/- 0.003 | 0.294 +/- 0.003 | 0.2967 +/- 0.003 |
Here is the Radius vs Width plot for when abs(GenPart.z) <= 1.0 cm along with a table of values.
R (cm) | 0.1 | 0.2 | 0.3 | 0.4 | 0.4 | 0.5 | 0.6 | 0.7 | 0.8 | 0.9 | 1.0 |
Width (cm) | 0.5467 +/- 0.0351 | 0.4475 +/- 0.0102 | 0.4976 +/- 0.0084 | 0.505 +/- 0.006 | 0.4966 +/- 0.0055 | 0.4884 +/- 0.0045 | 0.502 +/- 0.004 | 0.5033 +/- 0.0038 | 0.4788 +/- 0.0029 | 0.4766 +/- 0.0026 |
For this section I created a target that is 2cm in X and Y, and 6cm in Z. To produce the histograms I made a cut in the radius of a cylinder (R=0.1cm, R=0.2cm, .... , R=1.0cm) and also cut Z (Z=1cm, Z=2cm,...,Z=6cm) so that for each fixed radius there will be 6 histograms corresponding to changes in Z. For example, assume I am using R=0.1cm. There will be 6 histograms for X,Y,and Z corresponding to changing Z from 1cm to 6cm in 1cm increments. I also placed a cut for the reconstructed events using the commands abs(Event.z) <= 4.0 and sqrt(Event.x*Event.x+Event.y*Event.y) <=1.
For the plots below, a 6 GeV electron was fired at angles 25 degrees in theta and 30 degrees in phi. The target has an XY radius of 1cm and a length in the Z direction of 6cm. Setting GenPart.z = 0cm, the resolution was investigated as a function of the target's XY radius along with a table for the values.
Below are the Vx Differences for R=1.0cm as an example.
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