void DiElectronAnalyser::analyse(const EventPtr event, const ElectronCollection electrons) {
histMan_->setCurrentHistogramFolder(histogramFolder_);
weight_ = event->weight() * prescale_ * scale_;
if (electrons.size() == 2) {
ElectronPointer leadingElectron = electrons.front();
ElectronPointer secondElectron = electrons.at(1);
//both
histMan_->H1D_BJetBinned("diLepton_Mass")->Fill(leadingElectron->invariantMass(secondElectron), weight_);
histMan_->H1D_BJetBinned("diLepton_DeltaR")->Fill(leadingElectron->deltaR(secondElectron), weight_);
//lead
histMan_->H1D_BJetBinned("LeadLepton_Pt")->Fill(leadingElectron->pt(), weight_);
histMan_->H1D_BJetBinned("LeadLepton_Eta")->Fill(leadingElectron->eta(), weight_);
histMan_->H1D_BJetBinned("LeadLepton_AbsEta")->Fill(abs(leadingElectron->eta()), weight_);
histMan_->H1D_BJetBinned("LeadLepton_Phi")->Fill(leadingElectron->phi(), weight_);
histMan_->H1D_BJetBinned("LeadLepton_RelIso")->Fill(leadingElectron->pfRelativeIsolationRhoCorrected(), weight_);
histMan_->H1D_BJetBinned("LeadLepton_HadOverEM")->Fill(leadingElectron->HadOverEm(), weight_);
histMan_->H1D_BJetBinned("LeadLepton_dB")->Fill(leadingElectron->d0_wrtBeamSpot(), weight_);
//second
histMan_->H1D_BJetBinned("SecondLepton_Pt")->Fill(secondElectron->pt(), weight_);
histMan_->H1D_BJetBinned("SecondLepton_Eta")->Fill(secondElectron->eta(), weight_);
histMan_->H1D_BJetBinned("SecondLepton_AbsEta")->Fill(abs(secondElectron->eta()), weight_);
histMan_->H1D_BJetBinned("SecondLepton_Phi")->Fill(secondElectron->phi(), weight_);
histMan_->H1D_BJetBinned("SecondLepton_RelIso")->Fill(secondElectron->pfRelativeIsolationRhoCorrected(), weight_);
histMan_->H1D_BJetBinned("SecondLepton_HadOverEM")->Fill(secondElectron->HadOverEm(), weight_);
histMan_->H1D_BJetBinned("SecondLepton_dB")->Fill(secondElectron->d0_wrtBeamSpot(), weight_);
}
}
ElectronPointer TestObjectFactory::goodIsolatedElectron2() {
ElectronPointer electron = TestObjectFactory::goodIsolatedElectron();
ElectronPointer electron2 = TestObjectFactory::goodIsolatedElectron();
FourVector vec = electron->getFourVector();
vec.SetPx(vec.Px() - 20);
vec.SetPy(-vec.Py());
vec.SetPz(-vec.Pz());
electron2->setFourVector(vec);
//Z mass requirement
assert(electron->invariantMass(electron2) > 80);
assert(electron->invariantMass(electron2) < 100);
assert(electron2->relativeIsolation() < 0.1);
return electron2;
}
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_);
}
}
