#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

std::vector< DressedLepton >	leptons () const

Vector3	met () const

ParticleVector	neutrinos () const

ParticleVector	photons () const

	RivetAnalysis (const edm::ParameterSet &pset)

Private Attributes
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

std::vector< DressedLepton >	_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")),
 
       _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 145 of file RivetAnalysis.h.

References _phoIsoConeSize, _phoMaxRelIso, _phoMinPt, funct::abs(), GetRecoTauVFromDQM_MC_cff::cl, deltaR(), MillePedeFileConverter_cfg::e, event(), GeV, AlCaHLTBitMon_ParallelJobs::p, HadronAndPartonSelector_cfi::particles, and muons2muons_cfi::photon.

                                                 {
         _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").dressedLeptons();
         // search tau ancestors
         for ( auto & lepton : _leptons ) {
           const auto & cl = lepton.constituentLepton();
           for ( auto & p : cl.ancestors()) {
             if (p.abspid() == 15) {
               p.setMomentum(p.momentum()*10e-20);
               lepton.addPhoton(p, false);
             }
           }
         }
 
         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);
 
         }
 
         _jets      = applyProjection<FastJets>(event, "Jets").jetsByPt(jet_cut);
         _fatjets   = applyProjection<FastJets>(event, "FatJets").jetsByPt(fatjet_cut);
         _neutrinos = applyProjection<FinalState>(event, "Neutrinos").particlesByPt();
         _met       = applyProjection<MissingMomentum>(event, "MET").missingMomentum().p3();
       };

Jets Rivet::RivetAnalysis::fatjets ( ) const

inline

Definition at line 30 of file RivetAnalysis.h.

References _fatjets.

Referenced by ParticleLevelProducer::produce().

30 {return _fatjets;}

Rivet::RivetAnalysis::_fatjets

Jets _fatjets

Definition: RivetAnalysis.h:47

void Rivet::RivetAnalysis::finalize ( void )

inlineoverride

Definition at line 205 of file RivetAnalysis.h.

205 {};

void Rivet::RivetAnalysis::init ( void )

inlineoverride

Definition at line 79 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);
         
         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, /*cluster*/ true, /*useDecayPhotons*/ true);
         if (not _usePromptFinalStates)
           dressed_leptons = DressedLeptons(photons, charged_leptons, _lepConeSize, 
                             lepton_cut, /*cluster*/ true, /*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().

29 {return _jets;}

Rivet::RivetAnalysis::_jets

Jets _jets

Definition: RivetAnalysis.h:47

std::vector<DressedLepton> Rivet::RivetAnalysis::leptons ( void ) const

inline

Definition at line 26 of file RivetAnalysis.h.

References _leptons.

Referenced by ParticleLevelProducer::produce().

26 {return _leptons;}

Rivet::RivetAnalysis::_leptons

std::vector< DressedLepton > _leptons

Definition: RivetAnalysis.h:45

Vector3 Rivet::RivetAnalysis::met ( ) const

inline

Definition at line 31 of file RivetAnalysis.h.

References _met.

Referenced by ParticleLevelProducer::produce().

31 {return _met;}

Rivet::RivetAnalysis::_met

Vector3 _met

Definition: RivetAnalysis.h:48

ParticleVector Rivet::RivetAnalysis::neutrinos ( ) const

inline

Definition at line 28 of file RivetAnalysis.h.

References _neutrinos.

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

28 {return _neutrinos;}

Rivet::RivetAnalysis::_neutrinos

ParticleVector _neutrinos

Definition: RivetAnalysis.h:46

ParticleVector Rivet::RivetAnalysis::photons ( ) const

inline

Definition at line 27 of file RivetAnalysis.h.

References _photons.

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

27 {return _photons;}

Rivet::RivetAnalysis::_photons

ParticleVector _photons

Definition: RivetAnalysis.h:46

Member Data Documentation

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 41 of file RivetAnalysis.h.

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.

Jets Rivet::RivetAnalysis::_fatjets

private

Definition at line 47 of file RivetAnalysis.h.

Referenced by fatjets().

double Rivet::RivetAnalysis::_jetConeSize

private

Definition at line 40 of file RivetAnalysis.h.

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.

Jets Rivet::RivetAnalysis::_jets

private

Definition at line 47 of file RivetAnalysis.h.

Referenced by jets().

double Rivet::RivetAnalysis::_lepConeSize

private

Definition at line 39 of file RivetAnalysis.h.

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.

std::vector<DressedLepton> Rivet::RivetAnalysis::_leptons

private

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

double Rivet::RivetAnalysis::_particleMinPt

private

Definition at line 38 of file RivetAnalysis.h.

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.

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().

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.

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation