Rivet::RivetAnalysis Class Reference

#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

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

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

inlineoverride

Definition at line 141 of file RivetAnalysis.h.

References _doJetClustering, _fatJetMinPt, _fatjets, _jetMinPt, _jets, _leptons, _met, _neutrinos, _phoIsoConeSize, _phoMaxEta, _phoMaxRelIso, _phoMinPt, _photons, funct::abs(), haddnano::cl, HLT_FULL_cff::deltaR, alignCSCRings::e, edmPickEvents::event, GeV, reco::isLepton(), metsig::jet, AlCaHLTBitMon_ParallelJobs::p, and DiDispStaMuonMonitor_cfi::pt.

                                              {
      _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);
          }
        }
      }
    };

Jets Rivet::RivetAnalysis::fatjets ( ) const

inline

Definition at line 29 of file RivetAnalysis.h.

References _fatjets.

Referenced by ParticleLevelProducer::produce().

{ return _fatjets; }

void Rivet::RivetAnalysis::finalize ( void  )

inlineoverride

Definition at line 224 of file RivetAnalysis.h.

{};

void Rivet::RivetAnalysis::init ( void  )

inlineoverride

Definition at line 78 of file RivetAnalysis.h.

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

                         {
      // 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");
    };

Jets Rivet::RivetAnalysis::jets ( ) const

inline

Definition at line 28 of file RivetAnalysis.h.

References _jets.

Referenced by ParticleLevelProducer::produce().

{ return _jets; }

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

inline

Definition at line 25 of file RivetAnalysis.h.

References _leptons.

Referenced by ParticleLevelProducer::produce().

{ return _leptons; }

Vector3 Rivet::RivetAnalysis::met ( ) const

inline

Definition at line 30 of file RivetAnalysis.h.

References _met.

Referenced by ParticleLevelProducer::produce().

{ return _met; }

Particles Rivet::RivetAnalysis::neutrinos ( ) const

inline

Definition at line 27 of file RivetAnalysis.h.

References _neutrinos.

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

{ return _neutrinos; }

Particles Rivet::RivetAnalysis::photons ( ) const

inline

Definition at line 26 of file RivetAnalysis.h.

References _photons.

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

{ return _photons; }

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

inlineoverride

Definition at line 226 of file RivetAnalysis.h.

{ return "VALIDATED"; }

Member Data Documentation

bool Rivet::RivetAnalysis::_doJetClustering

private

Definition at line 36 of file RivetAnalysis.h.

Referenced by analyze().

bool Rivet::RivetAnalysis::_excludeNeutrinosFromJetClustering

private

Definition at line 35 of file RivetAnalysis.h.

Referenced by init().

bool Rivet::RivetAnalysis::_excludePromptLeptonsFromJetClustering

private

Definition at line 34 of file RivetAnalysis.h.

Referenced by init().

double Rivet::RivetAnalysis::_fatJetConeSize

private

Definition at line 41 of file RivetAnalysis.h.

Referenced by init().

double Rivet::RivetAnalysis::_fatJetMaxEta

private

Definition at line 41 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_fatJetMinPt

private

Definition at line 41 of file RivetAnalysis.h.

Referenced by analyze().

Jets Rivet::RivetAnalysis::_fatjets

private

Definition at line 46 of file RivetAnalysis.h.

Referenced by analyze(), and fatjets().

double Rivet::RivetAnalysis::_jetConeSize

private

Definition at line 40 of file RivetAnalysis.h.

Referenced by init().

double Rivet::RivetAnalysis::_jetMaxEta

private

Definition at line 40 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_jetMinPt

private

Definition at line 40 of file RivetAnalysis.h.

Referenced by analyze().

Jets Rivet::RivetAnalysis::_jets

private

Definition at line 46 of file RivetAnalysis.h.

Referenced by analyze(), and jets().

double Rivet::RivetAnalysis::_lepConeSize

private

Definition at line 39 of file RivetAnalysis.h.

Referenced by init().

double Rivet::RivetAnalysis::_lepMaxEta

private

Definition at line 39 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_lepMinPt

private

Definition at line 39 of file RivetAnalysis.h.

Referenced by init().

Particles Rivet::RivetAnalysis::_leptons

private

Definition at line 45 of file RivetAnalysis.h.

Referenced by analyze(), and leptons().

Vector3 Rivet::RivetAnalysis::_met

private

Definition at line 47 of file RivetAnalysis.h.

Referenced by analyze(), and met().

Particles Rivet::RivetAnalysis::_neutrinos

private

Definition at line 45 of file RivetAnalysis.h.

Referenced by analyze(), and neutrinos().

double Rivet::RivetAnalysis::_particleMaxEta

private

Definition at line 38 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_particleMinPt

private

Definition at line 38 of file RivetAnalysis.h.

Referenced by init().

double Rivet::RivetAnalysis::_phoIsoConeSize

private

Definition at line 43 of file RivetAnalysis.h.

Referenced by analyze().

double Rivet::RivetAnalysis::_phoMaxEta

private

Definition at line 43 of file RivetAnalysis.h.

Referenced by analyze().

double Rivet::RivetAnalysis::_phoMaxRelIso

private

Definition at line 43 of file RivetAnalysis.h.

Referenced by analyze().

double Rivet::RivetAnalysis::_phoMinPt

private

Definition at line 43 of file RivetAnalysis.h.

Referenced by analyze().

Particles Rivet::RivetAnalysis::_photons

private

Definition at line 45 of file RivetAnalysis.h.

Referenced by analyze(), and photons().

bool Rivet::RivetAnalysis::_usePromptFinalStates

private

Definition at line 33 of file RivetAnalysis.h.

Referenced by init().