bool TopPairMuPlusJetsReferenceSelection2011::passesTriggerSelection(const EventPtr event) const {
unsigned int runNumber(event->runnumber());
if (event->isRealData()) {
if (runNumber >= 160404 && runNumber < 173236)
return event->HLT(HLTriggers::HLT_IsoMu24);
else if (runNumber >= 173236)
return event->HLT(HLTriggers::HLT_IsoMu24_eta2p1);
else
return false;
} else {
//Fall11 MC
return event->HLT(HLTriggers::HLT_IsoMu24);
}
}
bool QCDPFRelIsoEPlusJetsSelection::passesTriggerSelection(const EventPtr event) const {
if (!useNonIsoTrigger_)
return TopPairEPlusJetsReferenceSelection::passesTriggerSelection(event);
else {
unsigned int runNumber(event->runnumber());
if (event->isRealData()) {
if (runNumber >= 160404 && runNumber <= 163869)
return event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_CentralTriJet30);
else if (runNumber > 163869 && runNumber <= 178380)
return event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralJet30);
else if (runNumber > 178380)
return event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralPFJet30);
else
return false;
} else {
//Fall11 MC
return event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralJet30);
}
}
}
bool TopPairEMuReferenceSelection::passesTriggerSelection(const EventPtr event) const {
return event->HLT(HLTriggers::HLT_Mu17_Ele8_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL) || event->HLT(HLTriggers::HLT_Mu8_Ele17_CaloIdT_CaloIsoVL_TrkIdVL_TrkIsoVL);
}
void HLTriggerQCDAnalyser::analyse(const EventPtr event) {
unsigned long run(event->runnumber());
//only do this analysis for runs above 193834 as previous runs don't have all triggers
//after run 194225 the trigger has been changed to HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30_30_20
if (!(run >= 193834 && run <= 194225 && event->isRealData()))
return;
// if (!(run >= 194270 && event->isRealData()))
// return;
if (event->Electrons().size() == 0)
return;
const ElectronPointer mostEnergeticElectron(event->Electrons().front());
bool passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT = event->HLT(
HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30)
|| event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30_30_20);
bool passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT = event->HLT(
HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30);
bool passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT = event->HLT(
HLTriggers::HLT_Ele25_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30);
bool passCaloIdVT_TrkIdT = event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralPFNoPUJet30)
|| event->HLT(HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralPFNoPUJet30_30_20);
if (passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT
&& (!studyExclusiveEvents_
|| (!passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT && !passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT
&& !passCaloIdVT_TrkIdT))) {
int prescale = event->HLTPrescale(HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30);
if (run >= 194270)
prescale = event->HLTPrescale(
HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_TriCentralPFNoPUJet30_30_20);
eleAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->setPrescale(prescale);
metNonIsoRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->setPrescale(prescale);
metAntiIDRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->setPrescale(prescale);
if (passesTriggerAnalysisSelection(event)) {
eleAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyse(event);
eleAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyseElectron(mostEnergeticElectron, event->weight());
}
if (passesNonIsoWithoutBtagAndHLT(event)) {
QCDNonIsoRegionCount_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_++;
metNonIsoRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyse(event, qcdNonIsoSelection_->signalLepton(event));
// metNonIsoRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyseTransverseMass(event,
// qcdNonIsoSelection_->signalLepton(event));
}
if (passesAntiIDWithoutBtagAndHLT(event)) {
QCDAntiIDRegionCount_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_++;
metAntiIDRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyse(event,
qcdAntiIDSelection_->signalLepton(event));
// metAntiIDRegionAnalyser_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_->analyseTransverseMass(event,
// qcdAntiIDSelection_->signalLepton(event));
}
if (passesSignalSelectionWithoutBtagAndHLT(event)) {
TopSignalRegionCount_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_++;
if (topSignalSelection_->hasAtLeastTwoGoodBJets(event))
TopSignal_TwoBtagsRegionCount_CaloIdVT_CaloIsoT_TrkIdT_TrkIsoT_++;
}
}
if (passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT
&& (!studyExclusiveEvents_
|| (!passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT && !passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT
&& !passCaloIdVT_TrkIdT))) {
int prescale = event->HLTPrescale(
(HLTriggers::HLT_Ele25_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30));
eleAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->setPrescale(prescale);
metNonIsoRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->setPrescale(prescale);
metAntiIDRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->setPrescale(prescale);
if (passesTriggerAnalysisSelection(event)) {
eleAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyse(event);
eleAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyseElectron(mostEnergeticElectron, event->weight());
}
if (passesNonIsoWithoutBtagAndHLT(event)) {
QCDNonIsoRegionCount_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_++;
metNonIsoRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyse(event,
qcdNonIsoSelection_->signalLepton(event));
// metNonIsoRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyseTransverseMass(event,
// qcdNonIsoSelection_->signalLepton(event));
}
if (passesAntiIDWithoutBtagAndHLT(event)) {
QCDAntiIDRegionCount_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_++;
metAntiIDRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyse(event,
qcdAntiIDSelection_->signalLepton(event));
// metAntiIDRegionAnalyser_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_->analyseTransverseMass(event,
// qcdAntiIDSelection_->signalLepton(event));
}
if (passesSignalSelectionWithoutBtagAndHLT(event)) {
TopSignalRegionCount_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_++;
if (topSignalSelection_->hasAtLeastTwoGoodBJets(event))
TopSignal_TwoBtagsRegionCount_CaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT_++;
}
}
if (passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT
&& (!studyExclusiveEvents_
|| (passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT && !passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT
&& !passCaloIdVT_TrkIdT))) {
int prescale = event->HLTPrescale(
(HLTriggers::HLT_Ele25_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_TriCentralPFNoPUJet30));
eleAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->setPrescale(prescale);
metNonIsoRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->setPrescale(prescale);
metAntiIDRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->setPrescale(prescale);
if (passesTriggerAnalysisSelection(event)) {
eleAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyse(event);
eleAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyseElectron(mostEnergeticElectron, event->weight());
}
if (passesNonIsoWithoutBtagAndHLT(event)) {
QCDNonIsoRegionCount_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_++;
metNonIsoRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyse(event,
qcdNonIsoSelection_->signalLepton(event));
// metNonIsoRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyseTransverseMass(event,
// qcdNonIsoSelection_->signalLepton(event));
}
if (passesAntiIDWithoutBtagAndHLT(event)) {
QCDAntiIDRegionCount_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_++;
metAntiIDRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyse(event,
qcdAntiIDSelection_->signalLepton(event));
// metAntiIDRegionAnalyser_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_->analyseTransverseMass(event,
// qcdAntiIDSelection_->signalLepton(event));
}
if (passesSignalSelectionWithoutBtagAndHLT(event)) {
TopSignalRegionCount_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_++;
if (topSignalSelection_->hasAtLeastTwoGoodBJets(event))
TopSignal_TwoBtagsRegionCount_CaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT_++;
}
}
if (passCaloIdVT_TrkIdT
&& (!studyExclusiveEvents_
|| (!passCaloIdVT_CaloIsoVL_TrkIdVL_TrkIsoT && !passCaloIdVT_CaloIsoT_TrkIdT_TrkIsoT
&& !passCaloIdVL_CaloIsoT_TrkIdVL_TrkIsoT))) {
int prescale = event->HLTPrescale((HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralPFNoPUJet30));
if (run >= 194270)
prescale = event->HLTPrescale((HLTriggers::HLT_Ele25_CaloIdVT_TrkIdT_TriCentralPFNoPUJet30_30_20));
eleAnalyser_CaloIdVT_TrkIdT_->setPrescale(prescale);
metNonIsoRegionAnalyser_CaloIdVT_TrkIdT_->setPrescale(prescale);
metAntiIDRegionAnalyser_CaloIdVT_TrkIdT_->setPrescale(prescale);
if (passesTriggerAnalysisSelection(event)) {
eleAnalyser_CaloIdVT_TrkIdT_->analyse(event);
eleAnalyser_CaloIdVT_TrkIdT_->analyseElectron(mostEnergeticElectron, event->weight());
}
if (passesNonIsoWithoutBtagAndHLT(event)) {
QCDNonIsoRegionCount_CaloIdVT_TrkIdT_++;
metNonIsoRegionAnalyser_CaloIdVT_TrkIdT_->analyse(event,
qcdNonIsoSelection_->signalLepton(event));
// metNonIsoRegionAnalyser_CaloIdVT_TrkIdT_->analyseTransverseMass(event,
// qcdNonIsoSelection_->signalLepton(event));
}
if (passesAntiIDWithoutBtagAndHLT(event)) {
QCDAntiIDRegionCount_CaloIdVT_TrkIdT_++;
metAntiIDRegionAnalyser_CaloIdVT_TrkIdT_->analyse(event,
qcdAntiIDSelection_->signalLepton(event));
// metAntiIDRegionAnalyser_CaloIdVT_TrkIdT_->analyseTransverseMass(event,
// qcdAntiIDSelection_->signalLepton(event));
}
if (passesSignalSelectionWithoutBtagAndHLT(event)) {
TopSignalRegionCount_CaloIdVT_TrkIdT_++;
if (topSignalSelection_->hasAtLeastTwoGoodBJets(event))
TopSignal_TwoBtagsRegionCount_CaloIdVT_TrkIdT_++;
}
}
}
