#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 one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

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

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.

68 {}

Member Function Documentation

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

override

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

Definition at line 75 of file L1Validator.cc.

References _GenSource, _Hists, _L1EGammaBXSource, _L1GenJetSource, _L1JetBXSource, _L1MuonBXSource, _L1TauBXSource, _srcToken, funct::abs(), BXVector< T >::begin(), reco::deltaR(), BXVector< T >::end(), PVValHelper::eta, reco::LeafCandidate::eta(), L1ValidatorHists::Fill(), L1ValidatorHists::FillNumber(), nanoDQM_cfi::GenPart, TauGenJets_cfi::GenParticles, edm::Event::getByToken(), BXVector< T >::getFirstBX(), BXVector< T >::getLastBX(), mps_fire::i, metsig::jet, nanoDQM_cff::Jet, electronTrackIsolations_cfi::minDR, RPCpg::mu, nanoDQM_cff::Muon, L1ValidatorHists::NEvents, reco::LeafCandidate::pdgId(), phi, reco::LeafCandidate::phi(), EnergyCorrector::pt, reco::LeafCandidate::pt(), BXVector< T >::size(), mps_update::status, reco::LeafCandidate::status(), 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);
     }
   }
 }

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

overrideprotected

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(), DEFINE_FWK_MODULE, 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 = `0`
	)

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 = `0`
	)

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

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

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation