#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 &ev, edm::EventSetup const &es) final

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) override

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer ()

	DQMEDAnalyzer (DQMEDAnalyzer const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer &&)=delete

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) override

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) override

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from edm::one::EDProducer< edm::Accumulator, edm::EndLuminosityBlockProducer, edm::EndRunProducer, edm::one::WatchLuminosityBlocks, edm::one::WatchRuns >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

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

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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 edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 68 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 69 of file L1Validator.cc.

69 {

70 }

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 77 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_cfi::Jet, electronTrackIsolations_cfi::minDR, RPCpg::mu, nanoDQM_cfi::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 &
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 72 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 308 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 236 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 263 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 290 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 81 of file L1Validator.h.

Referenced by bookHistograms(), and L1Validator().

std::string L1Validator::_fileName

private

Definition at line 82 of file L1Validator.h.

Referenced by trackingPlots.TrackingSeedingLayerTable::draw().

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

private

Definition at line 84 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

L1ValidatorHists L1Validator::_Hists

private

Definition at line 92 of file L1Validator.h.

Referenced by analyze(), and bookHistograms().

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

private

Definition at line 86 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

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

private

Definition at line 90 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

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

private

Definition at line 88 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

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

private

Definition at line 85 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

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

private

Definition at line 87 of file L1Validator.h.

Referenced by analyze(), and L1Validator().

edm::EDGetTokenT<GenEventInfoProduct> L1Validator::_srcToken

private

Definition at line 89 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

Member Function Documentation

Member Data Documentation