const ElectronCollection TopPairEMuReferenceSelection::signalElectrons(const EventPtr event) const {
const MuonCollection muons(goodMuons(event));
const ElectronCollection electrons(goodElectrons(event));
ElectronCollection signalElectrons;
double ptMax = 0;
int storeIndexA = -1;
// cout << "loop for max pt, muons: " << muons.size() << "electrons: " << electrons.size() << endl;
if(electrons.size() >= 1 && muons.size() >= 1){
for (unsigned int indexA = 0; indexA < electrons.size(); ++indexA) {
const ElectronPointer electron(electrons.at(indexA));
for (unsigned int indexB = 0; indexB < muons.size(); ++indexB) {
const MuonPointer muon(muons.at(indexB));
if((electron->charge()*muon->charge() < 0) && ((electron->pt()+muon->pt())>ptMax)){
ptMax = electron->pt()+muon->pt();
storeIndexA = indexA;
}
}
}
if(storeIndexA != -1)
signalElectrons.push_back(electrons.at(storeIndexA));
}
return signalElectrons;
}
const ElectronCollection TopPairEMuReferenceSelection::goodElectrons(const EventPtr event) const {
const ElectronCollection allElectrons(event->Electrons());
ElectronCollection goodIsolatedElectrons;
for (unsigned int index = 0; index < allElectrons.size(); ++index) {
const ElectronPointer electron(allElectrons.at(index));
if (isGoodElectron(electron) && isIsolatedElectron(electron)) {
goodIsolatedElectrons.push_back(electron);
}
}
return goodIsolatedElectrons;
}
bool TopPairEMuReferenceSelection::passesZmassVeto(const EventPtr event) const {
const ElectronCollection electrons(signalElectrons(event));
const MuonCollection muons(signalMuons(event));
ElectronCollection electronColl;
MuonCollection muonColl;
double ptMax = 0;
int storeIndexA = -1;
int storeIndexB = -1;
if(electrons.size() >= 1 && muons.size() >= 1){
for (unsigned int indexA = 0; indexA < electrons.size(); ++indexA) {
const ElectronPointer electron(electrons.at(indexA));
for (unsigned int indexB = 0; indexB < muons.size(); ++indexB) {
const MuonPointer muon(muons.at(indexB));
if((electron->charge() == -muon->charge()) && ((electron->pt()+muon->pt())>ptMax)){
ptMax = electron->pt()+muon->pt();
storeIndexA = indexA;
storeIndexB = indexB;
}
}
if(storeIndexA != storeIndexB){
electronColl.push_back(electrons.at(storeIndexA));
muonColl.push_back(muons.at(storeIndexB));
}
}
}
ElectronPointer electron = electronColl.front();
MuonPointer muon = muonColl.front();
double mass = 0;
ParticlePointer dilepton;
dilepton = ParticlePointer(new Particle(*electron + *muon));
mass = dilepton->mass();
return mass < 76 || mass > 106;
}
void DiElectronAnalyser::analyse(const EventPtr event) {
histMan_->setCurrentHistogramFolder(histogramFolder_);
ElectronCollection electrons = event->Electrons();
weight_ = event->weight() * prescale_ * scale_;
if (electrons.size() == 2) {
ElectronPointer leadingElectron = electrons.front();
ElectronPointer secondElectron = electrons.at(1);
histMan_->H1D_BJetBinned("diElectronMass")->Fill(leadingElectron->invariantMass(secondElectron), weight_);
}
ElectronCollection isolatedElectrons;
for (unsigned int index = 0; index < electrons.size(); ++index) {
const ElectronPointer electron(electrons.at(index));
if (electron->pfRelativeIsolation(0.3) < 0.1)
isolatedElectrons.push_back(electron);
}
if (isolatedElectrons.size() == 2) {
ElectronPointer leadingElectron = isolatedElectrons.front();
ElectronPointer secondElectron = isolatedElectrons.at(1);
histMan_->H1D_BJetBinned("diElectronMass_iso")->Fill(leadingElectron->invariantMass(secondElectron), weight_);
}
}
const LeptonPointer QCDPFRelIsoEPlusJetsSelection::signalLepton(const EventPtr event) const {
// if (!hasExactlyOneIsolatedLepton(event)) {
// cerr << "An error occurred in QCD*Selection in event = " << event->eventnumber();
// cerr << ", run = " << event->runnumber() << ", lumi = " << event->lumiblock() << "!" << endl;
// cerr << "File = " << event->file() << endl;
// cerr
// << "Access exception: No signal lepton could be found. Event doesn't pass 'hasExactlyOneIsolatedLepton' selection"
// << endl;
// throw "Access exception: No signal lepton could be found. Event doesn't pass 'hasExactlyOneIsolatedLepton' selection";
// }
const ElectronCollection allElectrons(event->Electrons());
ElectronCollection goodElectrons;
for (unsigned int index = 0; index < allElectrons.size(); ++index) {
const ElectronPointer electron(allElectrons.at(index));
if (isGoodElectron(electron)) {
goodElectrons.push_back(electron);
}
}
return MostIsolatedElectron(goodElectrons);
}
