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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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 >
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::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 ( const edm::ParameterSet &  iConfig )

explicit

Definition at line 54 of file L1Validator.cc.

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

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

L1Validator::~L1Validator ( )

override

Definition at line 68 of file L1Validator.cc.

{}

Member Function Documentation

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.

References _GenSource, _Hists, _L1EGammaBXSource, _L1GenJetSource, _L1JetBXSource, _L1MuonBXSource, _L1TauBXSource, _srcToken, funct::abs(), reco::deltaR(), PVValHelper::eta, reco::LeafCandidate::eta(), L1ValidatorHists::Fill(), L1ValidatorHists::FillNumber(), jetmet_cfg::GenJets, edm::Event::getByToken(), mps_fire::i, configurableAnalysis::Jet, metsig::jet, HLT_FULL_cff::minDR, RPCpg::mu, dumpRecoGeometry_cfg::Muon, L1ValidatorHists::NEvents, reco::LeafCandidate::pdgId(), phi, reco::LeafCandidate::phi(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), dt_dqm_sourceclient_common_cff::reco, mps_update::status, reco::LeafCandidate::status(), metsig::tau, HLT_FULL_cff::Tau, and parallelization::uint().

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

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 70 of file L1Validator.cc.

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

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

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

static

Definition at line 332 of file L1Validator.cc.

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

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

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.

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

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

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.

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

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

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

private

Definition at line 313 of file L1Validator.cc.

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

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

Member Data Documentation

std::string L1Validator::_dirName

private

Definition at line 80 of file L1Validator.h.

Referenced by bookHistograms(), and L1Validator().

std::string L1Validator::_fileName

private

Definition at line 81 of file L1Validator.h.

Referenced by trackingPlots.TrackingSeedingLayerTable::draw().

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

private

Definition at line 83 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

L1ValidatorHists L1Validator::_Hists

private

Definition at line 91 of file L1Validator.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 85 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

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

private

Definition at line 89 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

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

private

Definition at line 87 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

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

private

Definition at line 84 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

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

private

Definition at line 86 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

edm::EDGetTokenT<GenEventInfoProduct> L1Validator::_srcToken

private

Definition at line 88 of file L1Validator.h.

Referenced by analyze(), and L1Validator().