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
ParticleVector	leptons () const
Vector3	met () const
ParticleVector	neutrinos () const
ParticleVector	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
ParticleVector	_leptons
Vector3	_met
ParticleVector	_neutrinos
double	_particleMaxEta
double	_particleMinPt
double	_phoIsoConeSize
double	_phoMaxEta
double	_phoMaxRelIso
double	_phoMinPt
ParticleVector	_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 51 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 147 of file RivetAnalysis.h.

References _phoIsoConeSize, _phoMaxRelIso, _phoMinPt, funct::abs(), haddnano::cl, boostedElectronIsolation_cff::deltaR, MillePedeFileConverter_cfg::e, event(), GeV, reco::isLepton(), metsig::jet, AlCaHLTBitMon_ParallelJobs::p, HadronAndPartonSelector_cfi::particles, muons2muons_cfi::photon, and EnergyCorrector::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   = applyProjection<DressedLeptons>(event, "DressedLeptons").particlesByPt();
        // search tau ancestors
        Particles promptleptons = applyProjection<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;
            }
          }
        }
        
        Particles fsparticles = applyProjection<FinalState>(event,"FS").particles();
        
        for ( auto & photon : applyProjection<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 = applyProjection<FinalState>(event, "Neutrinos").particlesByPt();
        _met       = applyProjection<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 30 of file RivetAnalysis.h.

References _fatjets.

Referenced by ParticleLevelProducer::produce().

{return _fatjets;}

void Rivet::RivetAnalysis::finalize ( void  )

inlineoverride

Definition at line 231 of file RivetAnalysis.h.

{};

void Rivet::RivetAnalysis::init ( void  )

inlineoverride

Definition at line 80 of file RivetAnalysis.h.

References GeV, neutrinos(), 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);
        
        addProjection(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);
        addProjection(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);
        addProjection(dressed_leptons, "DressedLeptons");
        
        // Photons
        addProjection(photons, "Photons");
        
        // Jets
        VetoedFinalState fsForJets(fs);
        if (_usePromptFinalStates and _excludePromptLeptonsFromJetClustering)
          fsForJets.addVetoOnThisFinalState(dressed_leptons);
        JetAlg::InvisiblesStrategy invisiblesStrategy = JetAlg::DECAY_INVISIBLES;
        if (_excludeNeutrinosFromJetClustering)
          invisiblesStrategy = JetAlg::NO_INVISIBLES;
        addProjection(FastJets(fsForJets, FastJets::ANTIKT, _jetConeSize,
                               JetAlg::ALL_MUONS, invisiblesStrategy), "Jets");
        
        // FatJets
        addProjection(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);
          addProjection(prompt_neutrinos, "Neutrinos");
        }
        else
          addProjection(neutrinos, "Neutrinos");
        
        // MET
        addProjection(MissingMomentum(fs), "MET");
      };

Jets Rivet::RivetAnalysis::jets ( ) const

inline

Definition at line 29 of file RivetAnalysis.h.

References _jets.

Referenced by ParticleLevelProducer::produce().

{return _jets;}

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

inline

Definition at line 26 of file RivetAnalysis.h.

References _leptons.

Referenced by ParticleLevelProducer::produce().

{return _leptons;}

Vector3 Rivet::RivetAnalysis::met ( ) const

inline

Definition at line 31 of file RivetAnalysis.h.

References _met.

Referenced by ParticleLevelProducer::produce().

{return _met;}

ParticleVector Rivet::RivetAnalysis::neutrinos ( ) const

inline

Definition at line 28 of file RivetAnalysis.h.

References _neutrinos.

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

{return _neutrinos;}

ParticleVector Rivet::RivetAnalysis::photons ( ) const

inline

Definition at line 27 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 233 of file RivetAnalysis.h.

{ return "VALIDATED"; }

Member Data Documentation

bool Rivet::RivetAnalysis::_doJetClustering

private

Definition at line 37 of file RivetAnalysis.h.

bool Rivet::RivetAnalysis::_excludeNeutrinosFromJetClustering

private

Definition at line 36 of file RivetAnalysis.h.

bool Rivet::RivetAnalysis::_excludePromptLeptonsFromJetClustering

private

Definition at line 35 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_fatJetConeSize

private

Definition at line 42 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_fatJetMaxEta

private

Definition at line 42 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_fatJetMinPt

private

Definition at line 42 of file RivetAnalysis.h.

Jets Rivet::RivetAnalysis::_fatjets

private

Definition at line 47 of file RivetAnalysis.h.

Referenced by fatjets().

double Rivet::RivetAnalysis::_jetConeSize

private

Definition at line 41 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_jetMaxEta

private

Definition at line 41 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_jetMinPt

private

Definition at line 41 of file RivetAnalysis.h.

Jets Rivet::RivetAnalysis::_jets

private

Definition at line 47 of file RivetAnalysis.h.

Referenced by jets().

double Rivet::RivetAnalysis::_lepConeSize

private

Definition at line 40 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_lepMaxEta

private

Definition at line 40 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_lepMinPt

private

Definition at line 40 of file RivetAnalysis.h.

ParticleVector Rivet::RivetAnalysis::_leptons

private

Definition at line 46 of file RivetAnalysis.h.

Referenced by leptons().

Vector3 Rivet::RivetAnalysis::_met

private

Definition at line 48 of file RivetAnalysis.h.

Referenced by met().

ParticleVector Rivet::RivetAnalysis::_neutrinos

private

Definition at line 46 of file RivetAnalysis.h.

Referenced by neutrinos().

double Rivet::RivetAnalysis::_particleMaxEta

private

Definition at line 39 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_particleMinPt

private

Definition at line 39 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_phoIsoConeSize

private

Definition at line 44 of file RivetAnalysis.h.

Referenced by analyze().

double Rivet::RivetAnalysis::_phoMaxEta

private

Definition at line 44 of file RivetAnalysis.h.

double Rivet::RivetAnalysis::_phoMaxRelIso

private

Definition at line 44 of file RivetAnalysis.h.

Referenced by analyze().

double Rivet::RivetAnalysis::_phoMinPt

private

Definition at line 44 of file RivetAnalysis.h.

Referenced by analyze().

ParticleVector Rivet::RivetAnalysis::_photons

private

Definition at line 46 of file RivetAnalysis.h.

Referenced by photons().

bool Rivet::RivetAnalysis::_usePromptFinalStates

private

Definition at line 34 of file RivetAnalysis.h.