09 Sep

September 2009 posts

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