#include <RivetAnalysis.h>

Inheritance diagram for Rivet::RivetAnalysis:

Public Member Functions
void	analyze (const Event &event) override

Jets	fatjets () const

void	finalize () override

void	init () override

Jets	jets () const

Particles	leptons () const

Vector3	met () const

Particles	neutrinos () const

Particles	photons () const

	RivetAnalysis (const edm::ParameterSet &pset)

std::string	status () const override

Private Attributes
bool	_doJetClustering

bool	_excludeNeutrinosFromJetClustering

bool	_excludePromptLeptonsFromJetClustering

double	_fatJetConeSize

double	_fatJetMaxEta

double	_fatJetMinPt

Jets	_fatjets

double	_jetConeSize

double	_jetMaxEta

double	_jetMinPt

Jets	_jets

double	_lepConeSize

double	_lepMaxEta

double	_lepMinPt

Particles	_leptons

Vector3	_met

Particles	_neutrinos

double	_particleMaxEta

double	_particleMinPt

double	_phoIsoConeSize

double	_phoMaxEta

double	_phoMaxRelIso

double	_phoMinPt

Particles	_photons

bool	_usePromptFinalStates

Detailed Description

Definition at line 23 of file RivetAnalysis.h.

Constructor & Destructor Documentation

◆ RivetAnalysis()

Rivet::RivetAnalysis::RivetAnalysis ( const edm::ParameterSet & pset )

inline

Definition at line 50 of file RivetAnalysis.h.

         : Analysis("RivetAnalysis"),
           _usePromptFinalStates(pset.getParameter<bool>("usePromptFinalStates")),
           _excludePromptLeptonsFromJetClustering(pset.getParameter<bool>("excludePromptLeptonsFromJetClustering")),
           _excludeNeutrinosFromJetClustering(pset.getParameter<bool>("excludeNeutrinosFromJetClustering")),
           _doJetClustering(pset.getParameter<bool>("doJetClustering")),
  
           _particleMinPt(pset.getParameter<double>("particleMinPt")),
           _particleMaxEta(pset.getParameter<double>("particleMaxEta")),
  
           _lepConeSize(pset.getParameter<double>("lepConeSize")),
           _lepMinPt(pset.getParameter<double>("lepMinPt")),
           _lepMaxEta(pset.getParameter<double>("lepMaxEta")),
  
           _jetConeSize(pset.getParameter<double>("jetConeSize")),
           _jetMinPt(pset.getParameter<double>("jetMinPt")),
           _jetMaxEta(pset.getParameter<double>("jetMaxEta")),
  
           _fatJetConeSize(pset.getParameter<double>("fatJetConeSize")),
           _fatJetMinPt(pset.getParameter<double>("fatJetMinPt")),
           _fatJetMaxEta(pset.getParameter<double>("fatJetMaxEta")),
  
           _phoMinPt(pset.getParameter<double>("phoMinPt")),
           _phoMaxEta(pset.getParameter<double>("phoMaxEta")),
           _phoIsoConeSize(pset.getParameter<double>("phoIsoConeSize")),
           _phoMaxRelIso(pset.getParameter<double>("phoMaxRelIso")) {}

Member Function Documentation

◆ analyze()

void Rivet::RivetAnalysis::analyze ( const Event & event )

inlineoverride

Definition at line 141 of file RivetAnalysis.h.

                                               {
       _jets.clear();
       _fatjets.clear();
       _leptons.clear();
       _photons.clear();
       _neutrinos.clear();
  
       // Get analysis objects from projections
       Cut jet_cut = (Cuts::abseta < _jetMaxEta) and (Cuts::pT > _jetMinPt * GeV);
       Cut fatjet_cut = (Cuts::abseta < _fatJetMaxEta) and (Cuts::pT > _fatJetMinPt * GeV);
  
       _leptons = apply<DressedLeptons>(event, "DressedLeptons").particlesByPt();
  
       // search tau ancestors
       Particles promptleptons = apply<PromptFinalState>(event, "PromptLeptons").particles();
       for (auto& lepton : _leptons) {
         const auto& cl = lepton.constituents().front();  // this has no ancestors anymore :(
         for (auto const& pl : promptleptons) {
           if (cl.momentum() == pl.momentum()) {
             for (auto& p : pl.ancestors()) {
               if (p.abspid() == 15) {
                 p.setMomentum(p.momentum() * 10e-20);
                 lepton.addConstituent(p, false);
               }
             }
             break;
           }
         }
       }
  
       // Photons
       Particles fsparticles = apply<FinalState>(event, "FS").particles();
  
       for (auto& photon : apply<FinalState>(event, "Photons").particlesByPt()) {
         if (photon.pt() < _phoMinPt)
           continue;
  
         if (abs(photon.eta()) > _phoMaxEta)
           continue;
  
         double photonptsum = 0;
  
         for (auto& fsparticle : fsparticles) {
           if (deltaR(fsparticle, photon) == 0)
             continue;
  
           if (deltaR(fsparticle, photon) > _phoIsoConeSize)
             continue;
  
           photonptsum += fsparticle.pt();
         }
  
         if (photonptsum / photon.pt() > _phoMaxRelIso)
           continue;
  
         _photons.push_back(photon);
       }
  
       if (_doJetClustering) {
         _jets = apply<FastJets>(event, "Jets").jetsByPt(jet_cut);
         _fatjets = apply<FastJets>(event, "FatJets").jetsByPt(fatjet_cut);
       }
       _neutrinos = apply<FinalState>(event, "Neutrinos").particlesByPt();
       _met = apply<MissingMomentum>(event, "MET").missingMomentum().p3();
  
       // check for leptonic decays of tag particles
       for (auto& jet : _jets) {
         for (auto& pt : jet.tags()) {
           for (auto& p : pt.children(isLepton)) {
             pt.addConstituent(p, false);
           }
         }
       }
       for (auto& jet : _fatjets) {
         for (auto& pt : jet.tags()) {
           for (auto& p : pt.children(isLepton)) {
             pt.addConstituent(p, false);
           }
         }
       }
     };

References _doJetClustering, _fatJetMaxEta, _fatJetMinPt, _fatjets, _jetMaxEta, _jetMinPt, _jets, _leptons, _met, _neutrinos, _phoIsoConeSize, _phoMaxEta, _phoMaxRelIso, _phoMinPt, _photons, funct::abs(), GetRecoTauVFromDQM_MC_cff::cl, PbPb_ZMuSkimMuonDPG_cff::deltaR, MillePedeFileConverter_cfg::e, GeV, reco::isLepton(), metsig::jet, AlCaHLTBitMon_ParallelJobs::p, muons2muons_cfi::photon, DiDispStaMuonMonitor_cfi::pt, and PVValHelper::pT.

◆ fatjets()

Jets Rivet::RivetAnalysis::fatjets ( ) const

inline

Definition at line 29 of file RivetAnalysis.h.

29 { return _fatjets; }

References _fatjets.

Referenced by ParticleLevelProducer::produce().

◆ finalize()

void Rivet::RivetAnalysis::finalize ( void )

inlineoverride

Definition at line 224 of file RivetAnalysis.h.

224 {};

◆ init()

void Rivet::RivetAnalysis::init ( void )

inlineoverride

Definition at line 78 of file RivetAnalysis.h.

                          {
       // Cuts
       Cut particle_cut = (Cuts::abseta < _particleMaxEta) and (Cuts::pT > _particleMinPt * GeV);
       Cut lepton_cut = (Cuts::abseta < _lepMaxEta) and (Cuts::pT > _lepMinPt * GeV);
  
       // Generic final state
       FinalState fs(particle_cut);
  
       declare(fs, "FS");
  
       // Dressed leptons
       ChargedLeptons charged_leptons(fs);
       IdentifiedFinalState photons(fs);
       photons.acceptIdPair(PID::PHOTON);
  
       PromptFinalState prompt_leptons(charged_leptons);
       prompt_leptons.acceptMuonDecays(true);
       prompt_leptons.acceptTauDecays(true);
       declare(prompt_leptons, "PromptLeptons");
  
       PromptFinalState prompt_photons(photons);
       prompt_photons.acceptMuonDecays(true);
       prompt_photons.acceptTauDecays(true);
  
       // useDecayPhotons=true allows for photons with tau ancestor,
       // photons from hadrons are vetoed by the PromptFinalState;
       // will be default DressedLeptons behaviour for Rivet >= 2.5.4
       DressedLeptons dressed_leptons(
           prompt_photons, prompt_leptons, _lepConeSize, lepton_cut, /*useDecayPhotons*/ true);
       if (not _usePromptFinalStates)
         dressed_leptons = DressedLeptons(photons, charged_leptons, _lepConeSize, lepton_cut, /*useDecayPhotons*/ true);
       declare(dressed_leptons, "DressedLeptons");
  
       declare(photons, "Photons");
  
       // Jets
       VetoedFinalState fsForJets(fs);
       if (_usePromptFinalStates and _excludePromptLeptonsFromJetClustering)
         fsForJets.addVetoOnThisFinalState(dressed_leptons);
       JetAlg::Invisibles invisiblesStrategy = JetAlg::Invisibles::DECAY;
       if (_excludeNeutrinosFromJetClustering)
         invisiblesStrategy = JetAlg::Invisibles::NONE;
       declare(FastJets(fsForJets, FastJets::ANTIKT, _jetConeSize, JetAlg::Muons::ALL, invisiblesStrategy), "Jets");
  
       // FatJets
       declare(FastJets(fsForJets, FastJets::ANTIKT, _fatJetConeSize), "FatJets");
  
       // Neutrinos
       IdentifiedFinalState neutrinos(fs);
       neutrinos.acceptNeutrinos();
       if (_usePromptFinalStates) {
         PromptFinalState prompt_neutrinos(neutrinos);
         prompt_neutrinos.acceptMuonDecays(true);
         prompt_neutrinos.acceptTauDecays(true);
         declare(prompt_neutrinos, "Neutrinos");
       } else
         declare(neutrinos, "Neutrinos");
  
       // MET
       declare(MissingMomentum(fs), "MET");
     };

References _excludeNeutrinosFromJetClustering, _excludePromptLeptonsFromJetClustering, _fatJetConeSize, _jetConeSize, _lepConeSize, _lepMaxEta, _lepMinPt, _particleMaxEta, _particleMinPt, _usePromptFinalStates, diffTwoXMLs::ALL(), GeV, neutrinos(), NONE, pat::PHOTON, photons(), and PVValHelper::pT.

◆ jets()

Jets Rivet::RivetAnalysis::jets ( ) const

inline

Definition at line 28 of file RivetAnalysis.h.

28 { return _jets; }

References _jets.

Referenced by ParticleLevelProducer::produce().

◆ leptons()

Particles Rivet::RivetAnalysis::leptons ( void ) const

inline

Definition at line 25 of file RivetAnalysis.h.

25 { return _leptons; }

References _leptons.

Referenced by ParticleLevelProducer::produce().

◆ met()

Vector3 Rivet::RivetAnalysis::met ( ) const

inline

Definition at line 30 of file RivetAnalysis.h.

30 { return _met; }

References _met.

Referenced by objects.METAnalyzer.METAnalyzer::applyDeltaMet(), and ParticleLevelProducer::produce().

◆ neutrinos()

Particles Rivet::RivetAnalysis::neutrinos ( ) const

inline

Definition at line 27 of file RivetAnalysis.h.

27 { return _neutrinos; }

References _neutrinos.

Referenced by init(), and ParticleLevelProducer::produce().

◆ photons()

Particles Rivet::RivetAnalysis::photons ( ) const

inline

Definition at line 26 of file RivetAnalysis.h.

26 { return _photons; }

References _photons.

Referenced by init(), and ParticleLevelProducer::produce().

◆ status()

std::string Rivet::RivetAnalysis::status ( void ) const

inlineoverride

Definition at line 226 of file RivetAnalysis.h.

226 { return "VALIDATED"; }

Member Data Documentation

◆ _doJetClustering

bool Rivet::RivetAnalysis::_doJetClustering

private

Definition at line 36 of file RivetAnalysis.h.

Referenced by analyze().

◆ _excludeNeutrinosFromJetClustering

bool Rivet::RivetAnalysis::_excludeNeutrinosFromJetClustering

private

Definition at line 35 of file RivetAnalysis.h.

Referenced by init().

◆ _excludePromptLeptonsFromJetClustering

bool Rivet::RivetAnalysis::_excludePromptLeptonsFromJetClustering

private

Definition at line 34 of file RivetAnalysis.h.

Referenced by init().

◆ _fatJetConeSize

double Rivet::RivetAnalysis::_fatJetConeSize

private

Definition at line 41 of file RivetAnalysis.h.

Referenced by init().

◆ _fatJetMaxEta

double Rivet::RivetAnalysis::_fatJetMaxEta

private

Definition at line 41 of file RivetAnalysis.h.

Referenced by analyze().

◆ _fatJetMinPt

double Rivet::RivetAnalysis::_fatJetMinPt

private

Definition at line 41 of file RivetAnalysis.h.

Referenced by analyze().

◆ _fatjets

Jets Rivet::RivetAnalysis::_fatjets

private

Definition at line 46 of file RivetAnalysis.h.

Referenced by analyze(), and fatjets().

◆ _jetConeSize

double Rivet::RivetAnalysis::_jetConeSize

private

Definition at line 40 of file RivetAnalysis.h.

Referenced by init().

◆ _jetMaxEta

double Rivet::RivetAnalysis::_jetMaxEta

private

Definition at line 40 of file RivetAnalysis.h.

Referenced by analyze().

◆ _jetMinPt

double Rivet::RivetAnalysis::_jetMinPt

private

Definition at line 40 of file RivetAnalysis.h.

Referenced by analyze().

◆ _jets

Jets Rivet::RivetAnalysis::_jets

private

Definition at line 46 of file RivetAnalysis.h.

Referenced by analyze(), and jets().

◆ _lepConeSize

double Rivet::RivetAnalysis::_lepConeSize

private

Definition at line 39 of file RivetAnalysis.h.

Referenced by init().

◆ _lepMaxEta

double Rivet::RivetAnalysis::_lepMaxEta

private

Definition at line 39 of file RivetAnalysis.h.

Referenced by init().

◆ _lepMinPt

double Rivet::RivetAnalysis::_lepMinPt

private

Definition at line 39 of file RivetAnalysis.h.

Referenced by init().

◆ _leptons

Particles Rivet::RivetAnalysis::_leptons

private

Definition at line 45 of file RivetAnalysis.h.

Referenced by analyze(), and leptons().

◆ _met

Vector3 Rivet::RivetAnalysis::_met

private

Definition at line 47 of file RivetAnalysis.h.

Referenced by analyze(), and met().

◆ _neutrinos

Particles Rivet::RivetAnalysis::_neutrinos

private

Definition at line 45 of file RivetAnalysis.h.

Referenced by analyze(), and neutrinos().

◆ _particleMaxEta

double Rivet::RivetAnalysis::_particleMaxEta

private

Definition at line 38 of file RivetAnalysis.h.

Referenced by init().

◆ _particleMinPt

double Rivet::RivetAnalysis::_particleMinPt

private

Definition at line 38 of file RivetAnalysis.h.

Referenced by init().

◆ _phoIsoConeSize

double Rivet::RivetAnalysis::_phoIsoConeSize

private

Definition at line 43 of file RivetAnalysis.h.

Referenced by analyze().

◆ _phoMaxEta

double Rivet::RivetAnalysis::_phoMaxEta

private

Definition at line 43 of file RivetAnalysis.h.

Referenced by analyze().

◆ _phoMaxRelIso

double Rivet::RivetAnalysis::_phoMaxRelIso

private

Definition at line 43 of file RivetAnalysis.h.

Referenced by analyze().

◆ _phoMinPt

double Rivet::RivetAnalysis::_phoMinPt

private

Definition at line 43 of file RivetAnalysis.h.

Referenced by analyze().

◆ _photons

Particles Rivet::RivetAnalysis::_photons

private

Definition at line 45 of file RivetAnalysis.h.

Referenced by analyze(), and photons().

◆ _usePromptFinalStates

bool Rivet::RivetAnalysis::_usePromptFinalStates

private

Definition at line 33 of file RivetAnalysis.h.

Referenced by init().

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ RivetAnalysis()

Member Function Documentation

◆ analyze()

◆ fatjets()

◆ finalize()

◆ init()

◆ jets()

◆ leptons()

◆ met()

◆ neutrinos()

◆ photons()

◆ status()

Member Data Documentation

◆ _doJetClustering

◆ _excludeNeutrinosFromJetClustering

◆ _excludePromptLeptonsFromJetClustering

◆ _fatJetConeSize

◆ _fatJetMaxEta

◆ _fatJetMinPt

◆ _fatjets

◆ _jetConeSize

◆ _jetMaxEta

◆ _jetMinPt

◆ _jets

◆ _lepConeSize

◆ _lepMaxEta

◆ _lepMinPt

◆ _leptons

◆ _met

◆ _neutrinos

◆ _particleMaxEta

◆ _particleMinPt

◆ _phoIsoConeSize

◆ _phoMaxEta

◆ _phoMaxRelIso

◆ _phoMinPt

◆ _photons

◆ _usePromptFinalStates