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)

Figure 4: False rejection relative to triggered events