#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 > )

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 66 of file L1Validator.h.

Constructor & Destructor Documentation

◆ L1Validator()

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::~L1Validator ( )

override

Definition at line 68 of file L1Validator.cc.

68 {}

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.

References _GenSource, _Hists, _L1EGammaBXSource, _L1GenJetSource, _L1JetBXSource, _L1MuonBXSource, _L1TauBXSource, _srcToken, funct::abs(), BXVector< T >::begin(), reco::deltaR(), ntuple_cfi::Egamma, 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, HLT_2022v11_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.

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

◆ bookHistograms()

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

◆ fillDescriptions()

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

static

Definition at line 332 of file L1Validator.cc.

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

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

◆ 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.

References reco::deltaR(), reco::LeafCandidate::eta(), nanoDQM_cfi::GenPart, 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;
 }

◆ 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.

References reco::deltaR(), reco::LeafCandidate::eta(), nanoDQM_cfi::GenPart, 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;
 }

◆ 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.

References reco::deltaR(), reco::LeafCandidate::eta(), nanoDQM_cfi::GenPart, 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

◆ _dirName

std::string L1Validator::_dirName

private

Definition at line 79 of file L1Validator.h.

Referenced by bookHistograms(), and L1Validator().

◆ _fileName

std::string L1Validator::_fileName

private

Definition at line 80 of file L1Validator.h.

Referenced by trackingPlots.TrackingSeedingLayerTable::draw().

◆ _GenSource

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

private

Definition at line 82 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

◆ _Hists

L1ValidatorHists L1Validator::_Hists

private

Definition at line 90 of file L1Validator.h.

Referenced by analyze(), and bookHistograms().

◆ _L1EGammaBXSource

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

private

Definition at line 84 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

◆ _L1GenJetSource

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

private

Definition at line 88 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

◆ _L1JetBXSource

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

private

Definition at line 86 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

◆ _L1MuonBXSource

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

private

Definition at line 83 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

◆ _L1TauBXSource

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

private

Definition at line 85 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

◆ _srcToken

edm::EDGetTokenT<GenEventInfoProduct> L1Validator::_srcToken

private

Definition at line 87 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ L1Validator()

◆ ~L1Validator()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ fillDescriptions()

◆ FindBest() [1/3]

◆ FindBest() [2/3]

◆ FindBest() [3/3]

Member Data Documentation

◆ _dirName

◆ _fileName

◆ _GenSource

◆ _Hists

◆ _L1EGammaBXSource

◆ _L1GenJetSource

◆ _L1JetBXSource

◆ _L1MuonBXSource

◆ _L1TauBXSource

◆ _srcToken