2010.02.11 BFC and Pythia QA: Gain spread=0.05, mean=0.95
QCD and gamma-jet data samples are described here
Resultys without gain spread can be found here
(Note: ignore parton pt=25-35GeV for the gamma-jet sample since all jobs failed)
Gain spread implementation in StEEmcSlowMaker.cxx (private version):
void StEEmcSlowMaker::setTowerGainSpread(Float_t s, Float_t mTowerGainMean)
{
LOG_INFO << "setTowerGainSpread(): gain spread: " << s << "; gain mean value: " << mTowerGainMean << endm;
// initialize tower gain factors to 1
for ( Int_t sec=0;sec<kEEmcNumSectors;sec++ )
for ( Int_t sub=0;sub<kEEmcNumSubSectors;sub++ )
for ( Int_t eta=0;eta<kEEmcNumEtas;eta++ )
{
// mTowerGainFact[sec][sub][eta]=1.0;
Float_t f = -1.0E9;
while ( f <= -1. || f >= 1.0 )
f = gRandom->Gaus(0., s);
mTowerGainFact[sec][sub][eta] = mTowerGainMean + f;
}
}
Pythia filter configuration
StEemcGammaFilter:: running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events
StEemcGammaFilter:: mConeRadius 0.22 mSeedThreshold 2.1 mClusterThreshold 3.25 mEtaLow 0.95 mEtaHigh 2.1
StEemcGammaFilter:: mCalDepth 279.5 mMinPartEnergy 1e-05 mHadronScale 0.4 mFilterMode 0 mPrintLevel 0
BFC filter configuration
StChain:INFO - Init() : Seed energy threshold = 2.8 GeV
StChain:INFO - Init() : Cluster eT threshold = 4.2 GeV
StChain:INFO - Init() : Maximum vertex = +/- 120 cm
StChain:INFO - Init() : Running the TEST mode (accepting all events). Set mFilterMode=1 to actually reject events in BFC
Accept/Reject relative to the total number of Pythia generated events
Figure 1: Fraction of accepted events
Figure 2: False rejection (Y-axis scale is 10^-3)
Accept/Reject relative to the number of triggered events
Figure 3: Fraction of accepted events (relative to triggered events)
- Printer-friendly version
- Login or register to post comments