bool PseudoTopAnalyser::passesEventSelection( const MCParticlePointer pseudoLepton, const ParticlePointer pseudoNeutrino, const JetCollection pseudoJets, const MCParticleCollection pseudoBs, const ParticleCollection allPseudoLeptons, const ParticlePointer pseudoMET ) { // Event selection taken from here : https://twiki.cern.ch/twiki/bin/view/LHCPhysics/ParticleLevelTopDefinitions unsigned int numberGoodLeptons = 0; unsigned int numberVetoLeptons = 0; ParticlePointer leadingLepton; for ( unsigned int leptonIndex = 0; leptonIndex < allPseudoLeptons.size(); ++ leptonIndex ) { const ParticlePointer lepton = allPseudoLeptons.at(leptonIndex); // Check if this is a good signal type lepton if ( lepton->pt() > minLeptonPt_ && fabs(lepton->eta()) < maxLeptonAbsEta_ ) { ++numberGoodLeptons; if ( leadingLepton == 0 ) leadingLepton = lepton; } // Check if this is a veto lepton if ( lepton->pt() > minVetoLeptonPt_ && fabs(lepton->eta()) < maxVetoLeptonAbsEta_ ) { ++numberVetoLeptons; } } // Neutrino pt sum bool passesNeutrinoSumPt = false; if ( pseudoMET != 0 ) { if ( pseudoMET->pt() > minNeutrinoSumPt_ ) passesNeutrinoSumPt = true; } // W MT bool passesWMT = false; if ( leadingLepton != 0 && pseudoMET != 0 ) { double genMT = sqrt( 2 * leadingLepton->pt() * pseudoMET->pt() * ( 1 - cos(leadingLepton->phi() - pseudoMET->phi() ) ) ); if (genMT > minWMt_) passesWMT = true; } // Jets unsigned int numberGoodJets = 0; unsigned int numberGoodBJets = 0; for ( unsigned int jetIndex = 0; jetIndex < pseudoJets.size(); ++ jetIndex ) { const JetPointer jet = pseudoJets.at(jetIndex); // Check if this is a good jet if ( jet->pt() > minJetPt_ && fabs(jet->eta()) < maxJetAbsEta_ ) { ++numberGoodJets; // Check if this is also a good b jet if ( fabs( jet->partonFlavour() ) == 5 ) { ++numberGoodBJets; } } } if ( numberGoodLeptons == 1 && numberVetoLeptons <= 1 && passesNeutrinoSumPt && passesWMT && numberGoodJets >= minNJets_ && numberGoodBJets >= minNBJets_ ) { return true; } else return false; }
void METAnalyser::analyseTransverseMass(const EventPtr event, const ParticlePointer particle) { histMan_->setCurrentHistogramFolder(histogramFolder_); weight_ = event->weight() * prescale_ * scale_; for (unsigned index = 0; index < METAlgorithm::NUMBER_OF_METALGORITHMS; ++index) { std::string prefix = METAlgorithm::prefixes.at(index); METAlgorithm::value metType = (METAlgorithm::value) index; if (!MET::isAvailableInNTupleVersion(Globals::NTupleVersion, index)) continue; if (MET::isMCOnlyMETType(index) && event->isRealData()) continue; //skip MC only METs for real data const METPointer met(event->MET(metType)); histMan_->setCurrentHistogramFolder(histogramFolder_ + "/" + prefix); double MT = Event::MT(particle, met); double angle = met->angle(particle); double delPhi = met->deltaPhi(particle); histMan_->H1D("Transverse_Mass")->Fill(MT, weight_); histMan_->H1D("Angle_lepton_MET")->Fill(angle, weight_); histMan_->H1D("DeltaPhi_lepton_MET")->Fill(delPhi, weight_); if (met->et() < 20) histMan_->H1D("Transverse_Mass_MET20")->Fill(MT, weight_); histMan_->H2D("MET_vs_leptonPt")->Fill(particle->pt(), met->et(), weight_); histMan_->H1D("MET_plus_leptonPt")->Fill(particle->pt() + met->et(), weight_); histMan_->H2D("MET_vs_leptonEta")->Fill(particle->eta(), met->et(), weight_); } }
double Particle::deltaEta(const ParticlePointer other) const { return eta() - other->eta(); }