L1Validator Class Reference

#include <L1Validator.h>

Inheritance diagram for L1Validator:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	L1Validator (const edm::ParameterSet &)
	~L1Validator () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
const reco::LeafCandidate *	FindBest (const reco::GenParticle *, const std::vector< l1extra::L1EmParticle > *, const std::vector< l1extra::L1EmParticle > *)
const reco::LeafCandidate *	FindBest (const reco::GenParticle *, const std::vector< l1extra::L1JetParticle > *, const std::vector< l1extra::L1JetParticle > *)
const reco::LeafCandidate *	FindBest (const reco::GenParticle *, const std::vector< l1extra::L1MuonParticle > *)

Private Attributes
std::string	_dirName
std::string	_fileName
edm::EDGetTokenT< reco::GenParticleCollection >	_GenSource
L1ValidatorHists	_Hists
edm::EDGetTokenT< l1t::EGammaBxCollection >	_L1EGammaBXSource
edm::EDGetTokenT< reco::GenJetCollection >	_L1GenJetSource
edm::EDGetTokenT< l1t::JetBxCollection >	_L1JetBXSource
edm::EDGetTokenT< l1t::MuonBxCollection >	_L1MuonBXSource
edm::EDGetTokenT< l1t::TauBxCollection >	_L1TauBXSource
edm::EDGetTokenT< GenEventInfoProduct >	_srcToken

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 67 of file L1Validator.h.

Constructor & Destructor Documentation

L1Validator()

L1Validator::L1Validator ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 54 of file L1Validator.cc.

                                                       {
  _dirName = iConfig.getParameter<std::string>("dirName");
  _GenSource = consumes<reco::GenParticleCollection>(iConfig.getParameter<edm::InputTag>("GenSource"));
 
  _L1MuonBXSource = consumes<l1t::MuonBxCollection>(iConfig.getParameter<edm::InputTag>("L1MuonBXSource"));
  _L1EGammaBXSource = consumes<l1t::EGammaBxCollection>(iConfig.getParameter<edm::InputTag>("L1EGammaBXSource"));
  _L1TauBXSource = consumes<l1t::TauBxCollection>(iConfig.getParameter<edm::InputTag>("L1TauBXSource"));
  _L1JetBXSource = consumes<l1t::JetBxCollection>(iConfig.getParameter<edm::InputTag>("L1JetBXSource"));
  _srcToken = mayConsume<GenEventInfoProduct>(iConfig.getParameter<edm::InputTag>("srcToken"));
  _L1GenJetSource = consumes<reco::GenJetCollection>(iConfig.getParameter<edm::InputTag>("L1GenJetSource"));
 
  //_fileName = iConfig.getParameter<std::string>("fileName");
}

References _dirName, _GenSource, _L1EGammaBXSource, _L1GenJetSource, _L1JetBXSource, _L1MuonBXSource, _L1TauBXSource, _srcToken, edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

~L1Validator()

L1Validator::~L1Validator ( )

override

Definition at line 68 of file L1Validator.cc.

{}

Member Function Documentation

analyze()

void L1Validator::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

overridevirtual

select the final state (i.e status==1) muons (pdg==+/-13)

select the final state (i.e status==1) electrons (pdg==+/-11) and photons (pdg==22)

select the matrix element (i.e status==2) taus (pdg==+/-15) before decay

Reimplemented from DQMEDAnalyzer.

Definition at line 75 of file L1Validator.cc.

                                                                             {
  using namespace edm;
  using namespace std;
  using namespace l1extra;
  using namespace reco;
 
  Handle<GenParticleCollection> GenParticles;
  Handle<l1t::MuonBxCollection> MuonsBX;
  Handle<l1t::EGammaBxCollection> EGammasBX;
  Handle<l1t::TauBxCollection> TausBX;
  Handle<l1t::JetBxCollection> JetsBX;
  Handle<GenEventInfoProduct> genEvtInfoProduct;
  Handle<reco::GenJetCollection> GenJets;
 
  bool GotEverything = true;
 
  if (!iEvent.getByToken(_GenSource, GenParticles))
    GotEverything = false;
  if (!iEvent.getByToken(_L1MuonBXSource, MuonsBX))
    GotEverything = false;
  if (!iEvent.getByToken(_L1EGammaBXSource, EGammasBX))
    GotEverything = false;
  if (!iEvent.getByToken(_L1TauBXSource, TausBX))
    GotEverything = false;
  if (!iEvent.getByToken(_L1JetBXSource, JetsBX))
    GotEverything = false;
  if (!iEvent.getByToken(_srcToken, genEvtInfoProduct))
    GotEverything = false;
  if (!iEvent.getByToken(_L1GenJetSource, GenJets))
    GotEverything = false;
 
  if (!GotEverything)
    return;
 
  /*
  std::string moduleName = "";
  if( genEvtInfoProduct.isValid() ) {
          const edm::Provenance& prov =
  iEvent.getProvenance(genEvtInfoProduct.id()); moduleName =
  edm::moduleName(prov);
          //cout<<" generator name: "<<moduleName<<endl;
  }
  */
 
  _Hists.NEvents++;
 
  int nL1Muons = 0, nL1EGammas = 0, nL1Taus = 0, nL1Jets = 0;
  if (MuonsBX->getFirstBX() >= 0)
    nL1Muons = MuonsBX->size(0);
  if (EGammasBX->getFirstBX() >= 0)
    nL1EGammas = EGammasBX->size(0);
  if (TausBX->getFirstBX() >= 0)
    nL1Taus = TausBX->size(0);
  if (JetsBX->getFirstBX() >= 0)
    nL1Jets = JetsBX->size(0);
 
  _Hists.FillNumber(L1ValidatorHists::Type::Muon, nL1Muons);
  _Hists.FillNumber(L1ValidatorHists::Type::Egamma, nL1EGammas);
  _Hists.FillNumber(L1ValidatorHists::Type::Tau, nL1Taus);
  _Hists.FillNumber(L1ValidatorHists::Type::Jet, nL1Jets);
 
  // For gen jet
 
  for (auto &Genjet : *GenJets) {
    // eta within calorimeter acceptance 4.7
    if (fabs((&Genjet)->eta()) > 4.7)
      continue;
 
    // only consider the gen jet with pt greater than 10 GeV
    if ((&Genjet)->pt() < 10.0)
      continue;
 
    double minDR = 999.0;
 
    // match L1T object
    const l1t::Jet *L1Part = nullptr;
    for (int iBx = JetsBX->getFirstBX(); iBx <= JetsBX->getLastBX(); ++iBx) {
      if (iBx > 0)
        continue;
      for (std::vector<l1t::Jet>::const_iterator jet = JetsBX->begin(iBx); jet != JetsBX->end(iBx); ++jet) {
        double idR = reco::deltaR((&Genjet)->eta(), (&Genjet)->phi(), jet->eta(), jet->phi());
        if (idR < minDR) {
          minDR = idR;
          L1Part = &(*jet);
        }
      }
    }
    _Hists.Fill(L1ValidatorHists::Type::Jet, &Genjet, L1Part);
  }
 
  for (uint i = 0; i < GenParticles->size(); i++) {
    const GenParticle *GenPart = &GenParticles->at(i);
 
    int pdg = GenPart->pdgId(), status = GenPart->status();
 
    double minDR = 999.0;
 
    // only consider the gen particle with pt greater than 10 GeV
    if (GenPart->pt() < 10.0)
      continue;
 
    if (status == 1 && abs(pdg) == 13) {  // Muon
 
      // eta within tracker acceptance 2.4
      if (fabs(GenPart->eta()) > 2.4)
        continue;
 
      // match L1T object
      const l1t::Muon *L1Part = nullptr;
      for (int iBx = MuonsBX->getFirstBX(); iBx <= MuonsBX->getLastBX(); ++iBx) {
        if (iBx > 0)
          continue;
        for (std::vector<l1t::Muon>::const_iterator mu = MuonsBX->begin(iBx); mu != MuonsBX->end(iBx); ++mu) {
          double idR = reco::deltaR(GenPart->eta(), GenPart->phi(), mu->eta(), mu->phi());
          if (idR < minDR) {
            minDR = idR;
            L1Part = &(*mu);
          }
        }
        _Hists.Fill(L1ValidatorHists::Type::Muon, GenPart, L1Part);
      }
 
    } else if (status == 1 && (abs(pdg) == 11 || pdg == 22)) {  // Egamma
 
      // eta within EM calorimeter acceptance 2.5
      if (fabs(GenPart->eta()) > 2.5)
        continue;
 
      // exclude the calorimeter barrel and endcap overlap region
      if (fabs(GenPart->eta()) > 1.4442 && fabs(GenPart->eta()) < 1.5660)
        continue;
 
      // match L1T object
      const l1t::EGamma *L1Part = nullptr;
      for (int iBx = EGammasBX->getFirstBX(); iBx <= EGammasBX->getLastBX(); ++iBx) {
        if (iBx > 0)
          continue;
        for (std::vector<l1t::EGamma>::const_iterator eg = EGammasBX->begin(iBx); eg != EGammasBX->end(iBx); ++eg) {
          double idR = reco::deltaR(GenPart->eta(), GenPart->phi(), eg->eta(), eg->phi());
          if (idR < minDR) {
            minDR = idR;
            L1Part = &(*eg);
          }
        }
      }
      _Hists.Fill(L1ValidatorHists::Type::Egamma, GenPart, L1Part);
 
    } else if (status == 2 && abs(pdg) == 15) {  // Tau
 
      // eta within tracker acceptance 2.4
      if (fabs(GenPart->eta()) > 2.4)
        continue;
 
      // match L1T object
      const l1t::Tau *L1Part = nullptr;
      for (int iBx = TausBX->getFirstBX(); iBx <= TausBX->getLastBX(); ++iBx) {
        if (iBx > 0)
          continue;
        for (std::vector<l1t::Tau>::const_iterator tau = TausBX->begin(iBx); tau != TausBX->end(iBx); ++tau) {
          double idR = reco::deltaR(GenPart->eta(), GenPart->phi(), tau->eta(), tau->phi());
          if (idR < minDR) {
            minDR = idR;
            L1Part = &(*tau);
          }
        }
      }
      _Hists.Fill(L1ValidatorHists::Type::Tau, GenPart, L1Part);
    }
  }
}

References _GenSource, _Hists, _L1EGammaBXSource, _L1GenJetSource, _L1JetBXSource, _L1MuonBXSource, _L1TauBXSource, _srcToken, funct::abs(), BXVector< T >::begin(), reco::deltaR(), BXVector< T >::end(), PVValHelper::eta, L1ValidatorHists::Fill(), L1ValidatorHists::FillNumber(), jetmet_cfg::GenJets, nanoDQM_cfi::GenPart, heavyFlavorValidation_cfi::GenParticles, BXVector< T >::getFirstBX(), BXVector< T >::getLastBX(), mps_fire::i, iEvent, L1TRate_cfi::Jet, metsig::jet, PDWG_EXOHSCP_cff::minDR, amptDefaultParameters_cff::mu, dumpRecoGeometry_cfg::Muon, L1ValidatorHists::NEvents, phi, DiDispStaMuonMonitor_cfi::pt, BXVector< T >::size(), mps_update::status, L1TRate_Offline_cfi::Tau, metsig::tau, and parallelization::uint.

bookHistograms()

void L1Validator::bookHistograms ( DQMStore::IBooker &  iBooker, edm::Run const &  , edm::EventSetup const &    )

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 70 of file L1Validator.cc.

                                                                                                {
  // iBooker.setCurrentFolder(_dirName);
  _Hists.Book(iBooker, _dirName);
};

References _dirName, _Hists, and L1ValidatorHists::Book().

fillDescriptions()

void L1Validator::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 332 of file L1Validator.cc.

                                                                             {
  // The following says we do not know what parameters are allowed so do no
  // validation
  // Please change this to state exactly what you do use, even if it is no
  // parameters
  edm::ParameterSetDescription desc;
  desc.setUnknown();
  descriptions.addDefault(desc);
}

References edm::ConfigurationDescriptions::addDefault(), and edm::ParameterSetDescription::setUnknown().

FindBest() [1/3]

const reco::LeafCandidate * L1Validator::FindBest ( const reco::GenParticle *  GenPart, const std::vector< l1extra::L1EmParticle > *  Collection1, const std::vector< l1extra::L1EmParticle > *  Collection2 = nullptr  )

private

Definition at line 253 of file L1Validator.cc.

                                                                                                              {
  const reco::LeafCandidate *BestPart = nullptr;
  double BestDR = 999.;
 
  for (uint i = 0; i < Collection1->size(); i++) {
    const reco::LeafCandidate *ThisPart = &Collection1->at(i);
    double ThisDR = reco::deltaR(GenPart->eta(), GenPart->phi(), ThisPart->eta(), ThisPart->phi());
    if (ThisDR < BestDR) {
      BestDR = ThisDR;
      BestPart = ThisPart;
    }
  }
 
  if (Collection2 == nullptr)
    return BestPart;
 
  for (uint i = 0; i < Collection2->size(); i++) {
    const reco::LeafCandidate *ThisPart = &Collection2->at(i);
    double ThisDR = reco::deltaR(GenPart->eta(), GenPart->phi(), ThisPart->eta(), ThisPart->phi());
    if (ThisDR < BestDR) {
      BestDR = ThisDR;
      BestPart = ThisPart;
    }
  }
 
  return BestPart;
}

References reco::deltaR(), reco::LeafCandidate::eta(), nanoDQM_cfi::GenPart, mps_fire::i, reco::LeafCandidate::phi(), and parallelization::uint.

FindBest() [2/3]

const reco::LeafCandidate * L1Validator::FindBest ( const reco::GenParticle *  GenPart, const std::vector< l1extra::L1JetParticle > *  Collection1, const std::vector< l1extra::L1JetParticle > *  Collection2 = nullptr  )

private

Definition at line 283 of file L1Validator.cc.

                                                                                                               {
  const reco::LeafCandidate *BestPart = nullptr;
  double BestDR = 999.;
 
  for (uint i = 0; i < Collection1->size(); i++) {
    const reco::LeafCandidate *ThisPart = &Collection1->at(i);
    double ThisDR = reco::deltaR(GenPart->eta(), GenPart->phi(), ThisPart->eta(), ThisPart->phi());
    if (ThisDR < BestDR) {
      BestDR = ThisDR;
      BestPart = ThisPart;
    }
  }
 
  if (Collection2 == nullptr)
    return BestPart;
 
  for (uint i = 0; i < Collection2->size(); i++) {
    const reco::LeafCandidate *ThisPart = &Collection2->at(i);
    double ThisDR = reco::deltaR(GenPart->eta(), GenPart->phi(), ThisPart->eta(), ThisPart->phi());
    if (ThisDR < BestDR) {
      BestDR = ThisDR;
      BestPart = ThisPart;
    }
  }
 
  return BestPart;
}

References reco::deltaR(), reco::LeafCandidate::eta(), nanoDQM_cfi::GenPart, mps_fire::i, reco::LeafCandidate::phi(), and parallelization::uint.

FindBest() [3/3]

const reco::LeafCandidate * L1Validator::FindBest ( const reco::GenParticle *  GenPart, const std::vector< l1extra::L1MuonParticle > *  Collection1  )

private

Definition at line 313 of file L1Validator.cc.

                                                                                                      {
  const reco::LeafCandidate *BestPart = nullptr;
  double BestDR = 999.;
 
  for (uint i = 0; i < Collection1->size(); i++) {
    const reco::LeafCandidate *ThisPart = &Collection1->at(i);
    double ThisDR = reco::deltaR(GenPart->eta(), GenPart->phi(), ThisPart->eta(), ThisPart->phi());
    if (ThisDR < BestDR) {
      BestDR = ThisDR;
      BestPart = ThisPart;
    }
  }
 
  return BestPart;
}

References reco::deltaR(), reco::LeafCandidate::eta(), nanoDQM_cfi::GenPart, mps_fire::i, reco::LeafCandidate::phi(), and parallelization::uint.

Member Data Documentation

_dirName

std::string L1Validator::_dirName

private

Definition at line 80 of file L1Validator.h.

Referenced by bookHistograms(), and L1Validator().

_fileName

std::string L1Validator::_fileName

private

Definition at line 81 of file L1Validator.h.

Referenced by trackingPlots.TrackingSeedingLayerTable::draw().

_GenSource

edm::EDGetTokenT<reco::GenParticleCollection> L1Validator::_GenSource

private

Definition at line 83 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

_Hists

L1ValidatorHists L1Validator::_Hists

private

Definition at line 91 of file L1Validator.h.

Referenced by analyze(), and bookHistograms().

_L1EGammaBXSource

edm::EDGetTokenT<l1t::EGammaBxCollection> L1Validator::_L1EGammaBXSource

private

Definition at line 85 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

_L1GenJetSource

edm::EDGetTokenT<reco::GenJetCollection> L1Validator::_L1GenJetSource

private

Definition at line 89 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

_L1JetBXSource

edm::EDGetTokenT<l1t::JetBxCollection> L1Validator::_L1JetBXSource

private

Definition at line 87 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

_L1MuonBXSource

edm::EDGetTokenT<l1t::MuonBxCollection> L1Validator::_L1MuonBXSource

private

Definition at line 84 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

_L1TauBXSource

edm::EDGetTokenT<l1t::TauBxCollection> L1Validator::_L1TauBXSource

private

Definition at line 86 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

_srcToken

edm::EDGetTokenT<GenEventInfoProduct> L1Validator::_srcToken

private

Definition at line 88 of file L1Validator.h.

Referenced by analyze(), and L1Validator().