12 Dec
December 2008 posts
2008.12.08 Run 8 EEMC QA
Ilya Selyuzhenkov December 08, 2008
Data sets:
- pp2006 - STAR 2006 pp longitudinal data (~ 3.164 pb^1)
Trigger: eemc-http-mb-L2gamma [id:137641]
- pp2008 - STAR 2008 pp data
Trigger: etot-mb-l2 [id:7]
Days: 53-70; ~0.5M triggered events (1/3 of available statistics)
Detector subsystems involved in analysis:
- TPC (vertex, jets, charge particle veto)
- Endcap EMC (triggering, photon candidate reconstruction)
- Barrel EMC (away side jet reconstruction)
Gamma-jet analysis cuts:
- Select only di-jet events
- cos(phi_gamma - phi_jet) < -0.8 : gamma-jet opposite in phi
- R_EM^jet < 0.9 : neutral energy fraction cut for the away side jet
- N_ch=0 : no charge tracks associated with a gamma candidate
- N_bTow = 0 : no barrel towers associated with a gamma candidate (gamma in the endcap)
- N_(5-strip cluster)^u > 2 : minimum number of strips in EEMC SMD u-plane cluster around peak
- N_(5-strip cluster)^v > 2 : minimum number of strips in EEMC SMD v-plane cluster around peak
- gamma-algo fail : my algorithm failed to match tower with SMD uv-intersection, etc...
- Tow:SMD match : SMD uv-intersection has a tower which is not in a 3x3 cluster
- R_{3x3cluster}: Energy in 3x3 cluster of EEMC tower to the total jet energy (not applied here)
Figure 1: EEMC x vs. y position of photon candidate for 2008 data sample
Problem with pre-shower layer in Sector 10 can been seen in the upper left corner
Figure 2: EEMC x vs. y position of photon candidate for 2006 data sample
Figure 3: Average < E_pre1 * E_pre2 > for 3x3 cluster around high tower
vs run number for sectors 9, 10 and 11
Note, zero pre-shower energy for sector 10 (black points) for days 61, 62, 64, and 67.
All di-jet events for pp2008 data are shown (no gamma-jet cuts)
Figure 3a: Same as Fig.3, zoom into day 61
Figure 3b: Same as Fig.3, zoom into day 62
Figure 3c: Same as Fig.3, zoom into day 64
Figure 3d: Same as Fig.3, zoom into day 67
Figure 4: EEMC x vs. y position of photon candidate for 2008 data sample
Same as Fig. 1, but excluding days: 61, 62, 64, and 67
Conclusion on QA:
No problem with pp2008 data have been found,
except that for some runs (mostly on days 61, 62, 64, and 67)
EEMC pre-shower layer for sector 10 was off.
Comparison between 2006 and 2008 data
Figure 5: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 6: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 5):
Figure 7: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 5, no scaling):
2008.12.09 pp Run 8 vs. Run 6 SMD shower shapes
Ilya Selyuzhenkov December 09, 2008
Data sets:
- pp2006 - STAR 2006 pp longitudinal data (~ 3.164 pb^1)
Trigger: eemc-http-mb-L2gamma [id:137641]
- pp2008 - STAR 2008 pp data
Trigger: etot-mb-l2 [id:7]
Days: 53-70; ~0.5M triggered events (1/3 of available statistics)
2008.12.09 pp Run 8 vs. Run 6 SMD shower shapes
Ilya Selyuzhenkov December 09, 2008
Data sets:
- pp2006 - STAR 2006 pp longitudinal data (~ 3.164 pb^1)
Trigger: eemc-http-mb-L2gamma [id:137641]
- pp2008 - STAR 2008 pp data
Trigger: etot-mb-l2 [id:7]
Days: 53-70; ~0.5M triggered events (1/3 of available statistics)
Gamma-jet analysis cuts:
- Select only di-jet events
- cos(phi_gamma - phi_jet) < -0.8 : gamma-jet opposite in phi
- R_EM^jet < 0.9 : neutral energy fraction cut for the away side jet
- N_ch=0 : no charge tracks associated with a gamma candidate
- N_bTow = 0 : no barrel towers associated with a gamma candidate (gamma in the endcap)
- N_(5-strip cluster)^u > 2 : minimum number of strips in EEMC SMD u-plane cluster around peak
- N_(5-strip cluster)^v > 2 : minimum number of strips in EEMC SMD v-plane cluster around peak
- gamma-algo fail : my algorithm failed to match tower with SMD uv-intersection, etc...
- Tow:SMD match : SMD uv-intersection has a tower which is not in a 3x3 cluster
- R_{3x3cluster}: Energy in 3x3 cluster of EEMC tower to the total jet energy (not applied here)
Comparison between 2006 and 2008 data
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
2008.12.09 pp Run 8 vs. Run 6 SMD shower shapes
Ilya Selyuzhenkov December 09, 2008
Data sets:
- pp2006 - STAR 2006 pp longitudinal data (~ 3.164 pb^1)
Trigger: eemc-http-mb-L2gamma [id:137641]
- pp2008 - STAR 2008 pp data
Trigger: etot-mb-l2 [id:7]
Days: 53-70; ~0.5M triggered events (1/3 of available statistics)
Gamma-jet analysis cuts:
- Select only di-jet events
- cos(phi_gamma - phi_jet) < -0.8 : gamma-jet opposite in phi
- R_EM^jet < 0.9 : neutral energy fraction cut for the away side jet
- N_ch=0 : no charge tracks associated with a gamma candidate
- N_bTow = 0 : no barrel towers associated with a gamma candidate (gamma in the endcap)
- N_(5-strip cluster)^u > 2 : minimum number of strips in EEMC SMD u-plane cluster around peak
- N_(5-strip cluster)^v > 2 : minimum number of strips in EEMC SMD v-plane cluster around peak
- gamma-algo fail : my algorithm failed to match tower with SMD uv-intersection, etc...
- Tow:SMD match : SMD uv-intersection has a tower which is not in a 3x3 cluster
- R_{3x3cluster}: Energy in 3x3 cluster of EEMC tower to the total jet energy (not applied here)
Comparison between 2006 and 2008 data
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
2008.12.09 pp Run 8 vs. Run 6 SMD shower shapes
Ilya Selyuzhenkov December 09, 2008
Data sets:
- pp2006 - STAR 2006 pp longitudinal data (~ 3.164 pb^1)
Trigger: eemc-http-mb-L2gamma [id:137641]
- pp2008 - STAR 2008 pp data
Trigger: etot-mb-l2 [id:7]
Days: 53-70; ~0.5M triggered events (1/3 of available statistics)
Gamma-jet analysis cuts:
- Select only di-jet events
- cos(phi_gamma - phi_jet) < -0.8 : gamma-jet opposite in phi
- R_EM^jet < 0.9 : neutral energy fraction cut for the away side jet
- N_ch=0 : no charge tracks associated with a gamma candidate
- N_bTow = 0 : no barrel towers associated with a gamma candidate (gamma in the endcap)
- N_(5-strip cluster)^u > 2 : minimum number of strips in EEMC SMD u-plane cluster around peak
- N_(5-strip cluster)^v > 2 : minimum number of strips in EEMC SMD v-plane cluster around peak
- gamma-algo fail : my algorithm failed to match tower with SMD uv-intersection, etc...
- Tow:SMD match : SMD uv-intersection has a tower which is not in a 3x3 cluster
- R_{3x3cluster}: Energy in 3x3 cluster of EEMC tower to the total jet energy (not applied here)
Comparison between 2006 and 2008 data
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
2008.12.09 pp Run 8 vs. Run 6 SMD shower shapes
Ilya Selyuzhenkov December 09, 2008
Data sets:
- pp2006 - STAR 2006 pp longitudinal data (~ 3.164 pb^1)
Trigger: eemc-http-mb-L2gamma [id:137641]
- pp2008 - STAR 2008 pp data
Trigger: etot-mb-l2 [id:7]
Days: 53-70; ~0.5M triggered events (1/3 of available statistics)
Gamma-jet analysis cuts:
- Select only di-jet events
- cos(phi_gamma - phi_jet) < -0.8 : gamma-jet opposite in phi
- R_EM^jet < 0.9 : neutral energy fraction cut for the away side jet
- N_ch=0 : no charge tracks associated with a gamma candidate
- N_bTow = 0 : no barrel towers associated with a gamma candidate (gamma in the endcap)
- N_(5-strip cluster)^u > 2 : minimum number of strips in EEMC SMD u-plane cluster around peak
- N_(5-strip cluster)^v > 2 : minimum number of strips in EEMC SMD v-plane cluster around peak
- gamma-algo fail : my algorithm failed to match tower with SMD uv-intersection, etc...
- Tow:SMD match : SMD uv-intersection has a tower which is not in a 3x3 cluster
- R_{3x3cluster}: Energy in 3x3 cluster of EEMC tower to the total jet energy (not applied here)
Comparison between 2006 and 2008 data
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
Conclusion:
2008.12.09 pp Run 8 vs. Run 6 shower shapes
Ilya Selyuzhenkov December 09, 2008
Data sets:
- pp2006 - STAR 2006 pp longitudinal data (~ 3.164 pb^1)
Trigger: eemc-http-mb-L2gamma [id:137641]
- pp2008 - STAR 2008 pp data
Trigger: etot-mb-l2 [id:7]
Days: 53-70; ~0.5M triggered events (1/3 of available statistics)
Gamma-jet analysis cuts:
- Select only di-jet events
- cos(phi_gamma - phi_jet) < -0.8 : gamma-jet opposite in phi
- R_EM^jet < 0.9 : neutral energy fraction cut for the away side jet
- N_ch=0 : no charge tracks associated with a gamma candidate
- N_bTow = 0 : no barrel towers associated with a gamma candidate (gamma in the endcap)
- N_(5-strip cluster)^u > 2 : minimum number of strips in EEMC SMD u-plane cluster around peak
- N_(5-strip cluster)^v > 2 : minimum number of strips in EEMC SMD v-plane cluster around peak
- gamma-algo fail : my algorithm failed to match tower with SMD uv-intersection, etc...
- Tow:SMD match : SMD uv-intersection has a tower which is not in a 3x3 cluster
- R_{3x3cluster}: Energy in 3x3 cluster of EEMC tower to the total jet energy (not applied here)
Comparison between 2006 and 2008 data
Figure 1: Vertex z distribution:
All gamma-jet cuts applied, plus pt_gamma>7 and pt_jet > 5 GeV (exlcuding days 61, 62, 64, and 67)
Results are shown for pp2008 data sample (black), vs. pp2006 data (red).
pp2008 data scaled to the same total number of candidates as in pp2006 data.
Figure 2: Shower shapes within +/- 30 strips from high strip (same cuts as in Fig. 1):
Figure 3: Shower shapes within +/- 5 strips from high strip
(same cuts as in Fig. 1, no scaling):
2008.12.11 Run 8 EEMC QA (presentation for Spin PWG)
Ilya Selyuzhenkov December 11, 2008
Run 8 QA with EEMC gamma-jet candidates
Presentation in pdf or open office file format
2008.12.16 Effect of L2gamma trigger in simulations
Ilya Selyuzhenkov December 16, 2008
Data sets
- pp2006 - STAR 2006 pp longitudinal data (~ 3.164 pb^1)
Trigger: eemc-http-mb-L2gamma [id:137641]
- gamma-jet[gamma-filtered] - data-driven Prompt Photon [p6410EemcGammaFilter] events.
Partonic pt range 2-25 GeV.
- QCD jets[gamma-filtered] - data-driven QCD [p6410EemcGammaFilter] events.
Partonic pt range 2-25 GeV.
Cuts applied
Gamma-jet isolation cuts except 3x3/r=0.7 energy isolation cut
Data driven shower shape replacement maker fix
Figure 1: (reproducing old results with dd-maker fix)
transverse momentum and vertex z distributions
before (with ideal gains/pedestals) and
after (with realistic gains/pedestal tables) dd-maker fix are in a good agreement.
For details on "dd-maker problem", read these hyper news threads:
emc2:2905, emc2:2900, and phana:294
Now we can run L2gamma trigger emulation and Eemc SMD ddMaker
in the same analysis chain.
l2-gamma trigger effect in simulation
Figure 2: Vertex z distribution with and without trigger condition in simulations
(emulated trigger: eemc-http-mb-L2gamma [id:137641]).
Solid red/green symbols show results with l2gamma condition applied,
while red/green lines show results for the same analysis cuts but without trigger condition.
Note, good agreement between MC QCD jets with trigger condition on (green solid squared)
and pp2006 data (black solid circles) for pre-shower1>0 case.
Figure 3: pt distribution with/without trigger condition in simulations.
Same color coding as in Fig. 2
Figure 4: Same as Fig. 3 just on a log scale
One can clearly see large trigger effect when applied for QCD jet events,
and a little effect for direct gammas.
Figure 5: gamma candidate pt QCD (right) and prompt photon (left) Monte-Carlo:
no (upper) with (lower) L2e-gamma trigger condition
No photon pt and no jet pt cuts
Figure 6: gamma candidate pt for QCD Monte-Carlo: no L2e-gamma trigger condition
No photon pt and no jet pt cuts
Figure 7: gamma candidate pt for QCD Monte-Carlo: L2e-gamma trigger condition applied (id:137641)
No photon pt and no jet pt cuts
2008.12.19 Parton pt distribution for Pythia QCD and gamma-jet events
Ilya Selyuzhenkov December 19, 2008
Data sets
Cuts applied
Gamma-jet isolation cuts except 3x3/r=0.7 energy isolation cut
Figure 1: Parton pt distibution for gamma-jet candidates from Pythia QCD sample
with various pt and l2gamma trigger conditions
Figure 2: Parton pt distibution for gamma-jet candidates from Pythia prompt photon sample
with various pt and l2gamma trigger conditions