2009.12.01 BEMC 1x1, 2x1, 2x2, 3x3 clustering
Monte-Carlo setup:
- Throwing one photon per event
- Full y2009 STAR geometry configurations with/without LOW_EM option
- Throw particles flat in eta (-1,1), phi (0, 2pi),
with energy: 30GeV with flat (+/-0.5 GeV) spread - bfc.C options:
trs,fss,y2009a,Idst,IAna,l0,tpcI,fcf,ftpc,Tree,logger,ITTF,Sti,MakeEvent,McEvent,
geant,evout,IdTruth,tags,bbcSim,tofsim,emcY2,EEfs,
GeantOut,big,-dstout,fzin,-MiniMcMk,beamLine,clearmem,eemcDB,VFPPVnoCTB - Use fixed (7%) sampling fraction in StEmcSimpleSimulator.cxx
mSF[0] = 1/0.07;
mSF[1] = 0.;
mSF[2] = 0.; - Vertex z=0
- 50K/per particle type
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");
Note: for BEMC ideal pedSigma set to 0, so effectively
there is no effect when I apply 3-sigma threshold above pedestal.
Figure 1: Energy sampling of various cluster in the Barrel EMC
E_reco is the total energy in the BEMC towers from mMuDstMaker->muDst()->muEmcCollection()
eta_thrown - rapidity of the thrown photon from the Geant record
Cuts: |eta| < 0.97 && |eta|>0.01 && total energy s.f. > 0.055
Figure 2: Various cluster energy ratios
Groups:
- Printer-friendly version
- Login or register to post comments