2009.09.04 Test of corrected EEMC geometry: SVT, slow-simulator on/off, pre-shower migration

Test of corrected EEMC geometry: shower shapes (bug 1618)

Monte-Carlo setup:

• One particle per event (photons, electrons, pions, and eta-meson)
• Full STAR 2006 geometry (with/without SVT)
  In Kumac file: detp geom y2006g; gexec $STAR_LIB/geometry.so (remove SVT with SVTT_OFF option)
• Throw particles flat in eta (1.08, 2.0), phi (0, 2pi), and pt (6-10 GeV)
• Using A2Emaker to get reconstructed Tower/SMD energy
  (with/without EEMC SlowSimulator in chain)
• Vertex z=0
• ~50K/per particle type
• Non-zero energy: 3 sigma above pedestal

Color coding:

• Photon (single particle MC)
• Electron (single particle MC)
• Neutral pion (single particle MC)
• Eta-meson (single particle MC)
• Eta-meson [pp2006 data] (single photons from eta-meson decay)

Pre-shower bins:

1. Ep1 = 0, Ep2 = 0 (no energy in both EEMC pre-shower layers)
2. Ep1 = 0, Ep2 > 0
3. 0 < Ep1 < 4 MeV
4. 4 < Ep1 < 10 MeV
5. Ep1 > 10 MeV
6. All pre-shower bins combined

Ep1/Ep2 is the energy deposited in the 1st/2nd EEMC pre-shower layer.
For a single particle MC it is a sum over
all pre-shower tiles in the EEMC with energy of 3 sigma above pedestal.
For eta-meson from pp2006 data the sum is over 3x3 tower patch

Shower shapes

Single particle kinematic cuts: pt=7-8GeV, eta=1.2-1.4
Eta-meson shower shapes (blue) taken from Fig. 1 from here of this post
All shapes are normalized to 1 at peak (central strip)

Figure 1: Shower shape sorted by pre-shower conditions.
cAir-Fixed EEMC geometry (NO slow simulator, WITH SVT)
Ratio plot

Figure 2: Shower shape sorted by pre-shower conditions.
cAir-Fixed EEMC geometry (NO slow simulator, NO SVT)
Ratio plot

Figure 3: Shower shape sorted by pre-shower conditions.
cAir-Fixed EEMC geometry (WITH slow simulator, WITH SVT)
Ratio plot

Figure 4: Shower shape sorted by pre-shower conditions.
Old cAir-bug EEMC geometry (NO slow simulator, WITH SVT)
Click here to see the plot

Pre-shower migration with/without SVT

Starting with a fixed (50K events) for each type of particle.
Change in number of counts for a given pre-shower bin
with different detector configuration shows pre-shower migration

Figure 5: Pre-shower migration.
cAir-Fixed EEMC geometry (WITH SVT)

Figure 6: Pre-shower migration.
cAir-Fixed EEMC geometry (WITHOUT SVT)

Sampling fraction with/without Slow-simulator

Figure 7: Sampling fraction (0.05 E_reco / E_thrown).
cAir-Fixed EEMC geometry (WITHOUT Slow-simulator)

Figure 8: Sampling fraction (0.05 E_reco / E_thrown).
cAir-Fixed EEMC geometry (WITH Slow-simulator)
Slow simulator introduce some non-linearity in the sampling fraction

Figure 9: Sampling fraction (0.05 E_reco / E_thrown).
cAir-Fixed EEMC geometry (WITHOUT SVT, WITHOUT Slow-simulator)
Click here to see the plot

Figure 10: Sampling fraction (0.05 E_reco / E_thrown).
Old cAir-bug EEMC geometry (NO slow simulator, WITH SVT)
Click here to see the plot

2009.09.11 Test of corrected EEMC geometry: LOW_EM cuts

Test of corrected EEMC geometry: SVT and LOW_EM cuts

Monte-Carlo setup:

• One particle per event (only photons in this post)
• Full STAR 2006 geometry (with/without SVT, LOW_EM cuts)
  In Kumac file: detp geom y2006g; gexec $STAR_LIB/geometry.so (vary SVTT_OFF, LOW_EM)
  LOW_EM cut definition is given at the end of this page
• Throw particles flat in eta (1.2, 1.9), phi (0, 2pi), and pt (6-10 GeV)
• Using A2Emaker to get reconstructed Tower/SMD energy
  (this post: no EEMC SlowSimulator)
• Vertex z=0
• ~50K/per iteration
• Non-zero energy: 3 sigma above pedestal

Color coding:

• SVT, LOW_EM marked in legend as LowEM (single photon MC)
• STV, no-LOW_EM marked in legend as default (single photon MC)
• no-SVT, no-LOW_EM marked in legend as no-SVT (single photon MC)
• photon-jet candidates [pp2006] (used data points from this post)
• photons from eta-meson [pp2006]

Pre-shower bins:

1. Ep1 = Ep2 = 0 (no energy in both EEMC pre-shower layers)
2. Ep1 = 0, Ep2 > 0
3. 0 < Ep1 < 4 MeV
4. 4 < Ep1 < 10 MeV
5. Ep1 > 10 MeV
6. All pre-shower bins combined

Note: Ep1/Ep2 is the energy deposited in the 1st/2nd EEMC pre-shower layer.
For a single photon MC it is a sum over
all pre-shower tiles in the EEMC with energy of 3 sigma above pedestal.
For eta-meson/gamma-jet candidates from pp2006 data the sum is over 3x3 tower patch

Shower shapes

Single particle kinematic cuts: pt=7-8GeV, eta=1.2-1.4
Eta-meson shower shapes (blue) taken from Fig. 1 from here of this post
All shapes are normalized to 1 at peak (central strip)

Figure 1: Shower shape sorted by pre-shower conditions.

Figure 2: Shower shape ratio. All shapes in Fig. 1 are divided by single photon shape
for "SVT+LOW_EM" configuration (black circles in Fig. 1)

Sampling fraction

Figure 3: Sampling fraction (0.05 * E_reco/ E_thrown)

Pre/post-shower energy and migration

Figure 4: Pre-shower1 energy (all tiles)

Figure 5: Pre-shower2 energy (all tiles)

Figure 6: Post-shower energy (all tiles)

Figure 7: Pre-shower bin photon migration

Tower energy profile

Figure 8a: Energy ratio in 2x1 to 3x3 cluster
For the first 4 pre-shower bins total yield in MC is normalized to that of the data
Blue circles indicate photon-jet candidates [pp2006] (points from this post)
Same data on a linear scale

Figure 8b: Energy ratio in 2x1 to 3x3 cluster: 7 < pt < 8 and 1.2 < eta < 1.4

 

Figure 8c: Energy ratio in 2x1 to 3x3 cluster: 7 < pt < 8 and 1.6 < eta < 1.8

Figure 9: Average energy ratio in 2x1 to 3x3 cluster vs. thrown energy

Figure 10: Average energy ratio in 2x1 to 3x3 cluster vs. thrown energy

LOW_EM cut definition

LOW_EM option for the STAR geometry (Low cuts on Electro-Magnetic processes)
is equivalent to the following set of GEANT cuts:

• CUTGAM=0.00001
• CUTELE=0.00001
• BCUTE =0.00001
• BCUTM =0.00001
• DCUTE =0.00001
• DCUTM =0.00001

All these values are for kinetic energy in GeV.

 

Cut meaning and GEANT default values:

• CUTGAM threshold for gamma transport (0.001);
• CUTELE threshold for electron and positron transport (0.001);
• BCUTE threshold for photons produced by electron bremsstrahlung (-1,);
• BCUTM threshold for photons produced by muon bremsstrahlung (-1);
• DCUTE threshold for electrons produced by electron delta-rays (-1);
• DCUTM threshold for electrons produced by muon or hadron delta-rays (-1);

Some details can be found at this link and in the GEANT manual