2009.11.03 BEMC sampling fraction: with and without LOW_EM option

Monte-Carlo setup:

  • Throwing one photon per event
  • Full STAR geometry (y2006g) configurations with and without LOW_EM option.
    Note: LOW_EM cuts are listed at the bottom of this page,
    and some related discussion can be found in this phana thread
  • Throw particles flat in eta (-1,1), phi (0, 2pi), and energy (30 +/- 0.5 GeV)
  • Vertex z=0
  • 50K/per particle type

Geometry configurations and notations:

  1. BEMC-noLOW_EM: Full STAR y2006g without LOW_EM option
  2. BEMC-LOW_EM: Full STAR y2006g with LOW_EM option

data base settings (same settings in bfc.C (Jan's trick) and in my MuDst reader):
dbMk->SetFlavor("sim","bemcPed");
dbMk->SetFlavor("Wbose","bemcCalib");
dbMk->SetFlavor("sim","bemcGain");
dbMk->SetFlavor("sim","bemcStatus");

dbMk->SetFlavor("sim","bprsPed");
dbMk->SetFlavor("Wbose","bprsCalib");
dbMk->SetFlavor("sim","bprsGain");
dbMk->SetFlavor("sim","bprsStatus");

dbMk->SetFlavor("sim","bsmdePed");
dbMk->SetFlavor("Wbose","bsmdeCalib");
dbMk->SetFlavor("sim","bsmdeGain");
dbMk->SetFlavor("sim","bsmdeStatus");

dbMk->SetFlavor("sim","bsmdpPed");
dbMk->SetFlavor("Wbose","bsmdpCalib");
dbMk->SetFlavor("sim","bsmdpGain");
dbMk->SetFlavor("sim","bsmdpStatus");

Note: for BEMC ideal pedSigma set to 0, so effectively
there is no effect when I apply 3-sigma threshold above pedestal.

Figure 1: E_reco/E_thrown distribution.
E_reco is the total energy in the BEMC towers from mMuDstMaker->muDst()->muEmcCollection()
E_thrown energy of the thrown photon from tne GEant record
No cut (yet) applied to exclude otliers in the average
Outliers in E_reco/E_thrown

Figure 2: Average E_reco/E_thrown vs. thrown photon eta (left) and phi (right)
Average is taken over a slice in eta or phi (no gaussian fits)

Figure 3: Average E_reco/E_thrown vs. thrown position (eta and phi)
Left: without LOW_EM option; right: with LOW_EM option
No cut applied to exclude otliers