Inheritance diagram for HTXSRivetProducer:

Public Member Functions
	HTXSRivetProducer (const edm::ParameterSet &cfg)

Public Member Functions inherited from edm::one::EDProducer< edm::one::WatchRuns, edm::one::SharedResources >
	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 () noexcept(false) 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	beginRun (edm::Run const &iRun, edm::EventSetup const &es) override

void	endRun (edm::Run const &iRun, edm::EventSetup const &es) override

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
Rivet::AnalysisHandler	_analysisHandler

edm::EDGetTokenT< edm::HepMCProduct >	_hepmcCollection

Rivet::HiggsTemplateCrossSections *	_HTXS

bool	_isFirstEvent

edm::EDGetTokenT< LHERunInfoProduct >	_lheRunInfo

std::string	_prodMode

HTXS::HiggsClassification	cat_

HTXS::HiggsProdMode	m_HiggsProdMode

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

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 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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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)

Detailed Description

Definition at line 30 of file HTXSRivetProducer.cc.

Constructor & Destructor Documentation

HTXSRivetProducer::HTXSRivetProducer ( const edm::ParameterSet & cfg )

inlineexplicit

Definition at line 32 of file HTXSRivetProducer.cc.

References Ecal2004TBTDCRanges_v1_cff::endRun, edm::ParameterSet::getParameter(), and HTXS::UNKNOWN.

       : _hepmcCollection(consumes<HepMCProduct>(cfg.getParameter<edm::InputTag>("HepMCCollection"))),
         _lheRunInfo(consumes<LHERunInfoProduct, edm::InRun>(cfg.getParameter<edm::InputTag>("LHERunInfo"))) {
     usesResource("Rivet");
 
     _HTXS = new Rivet::HiggsTemplateCrossSections();
 
     _isFirstEvent = true;
     _prodMode = cfg.getParameter<string>("ProductionMode");
     m_HiggsProdMode = HTXS::UNKNOWN;
 
     produces<HTXS::HiggsClassification>("HiggsClassification").setBranchAlias("HiggsClassification");
   }

Member Function Documentation

void HTXSRivetProducer::beginRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	es
	)

overrideprivate

Definition at line 171 of file HTXSRivetProducer.cc.

References DEFINE_FWK_MODULE, edm::Run::getByLabel(), HTXS::GG2ZH, HTXS::GGF, LHERunInfoProduct::headers_begin(), LHERunInfoProduct::headers_end(), mps_splice::line, groupFilesInBlocks::lines, edm::Handle< T >::product(), HTXS::QQ2ZH, writedatasetfile::run, AlCaHLTBitMon_QueryRunRegistry::string, HTXS::UNKNOWN, and HTXS::VBF.

                                                                             {
   if (_prodMode == "AUTO") {
     edm::Handle<LHERunInfoProduct> run;
     bool product_exists = iRun.getByLabel(edm::InputTag("externalLHEProducer"), run);
     if (product_exists) {
       typedef std::vector<LHERunInfoProduct::Header>::const_iterator headers_const_iterator;
       LHERunInfoProduct myLHERunInfoProduct = *(run.product());
       for (headers_const_iterator iter = myLHERunInfoProduct.headers_begin(); iter != myLHERunInfoProduct.headers_end();
            iter++) {
         std::vector<std::string> lines = iter->lines();
         for (unsigned int iLine = 0; iLine < lines.size(); iLine++) {
           std::string line = lines.at(iLine);
           // POWHEG
           if (line.find("gg_H_quark-mass-effects") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::GGF;
             break;
           }
           if (line.find("Process: HJ") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::GGF;
             break;
           }
           if (line.find("Process: HJJ") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::GGF;
             break;
           }
           if (line.find("VBF_H") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::VBF;
             break;
           }
           if (line.find("HZJ") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::QQ2ZH;
             break;
           }
           if (line.find("ggHZ") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::GG2ZH;
             break;
           }
           // MC@NLO
           if (line.find("ggh012j") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::GGF;
             break;
           }
           if (line.find("vbfh") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::VBF;
             break;
           }
           if (line.find("zh012j") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::QQ2ZH;
             break;
           }
           if (line.find("ggzh01j") != std::string::npos) {
             edm::LogInfo("HTXSRivetProducer") << iLine << " " << line << std::endl;
             m_HiggsProdMode = HTXS::GG2ZH;
             break;
           }
         }
 
         if (m_HiggsProdMode != HTXS::UNKNOWN)
           break;
       }
     }
   }
 }

void HTXSRivetProducer::endRun	(	edm::Run const &	iRun,
		edm::EventSetup const &	es
	)

overrideprivate

Definition at line 169 of file HTXSRivetProducer.cc.

169 { _HTXS->printClassificationSummary(); }

Rivet::HiggsTemplateCrossSections::printClassificationSummary

void printClassificationSummary()

Definition: HiggsTemplateCrossSections.cc:1045

HTXSRivetProducer::_HTXS

Rivet::HiggsTemplateCrossSections * _HTXS

Definition: HTXSRivetProducer.cc:56

void HTXSRivetProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &
	)

overrideprivate

Definition at line 65 of file HTXSRivetProducer.cc.

References funct::abs(), HTXS::BBH, eostools::cat(), class-composition::children, Exception, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), HTXS::GG2ZH, HTXS::GGF, eostools::move(), HadronAndPartonSelector_cfi::particles, edm::Event::put(), HTXS::QQ2ZH, HTXS::Rivet2Root(), HTXS::TH, HTXS::TTH, HTXS::UNKNOWN, HTXS::VBF, and HTXS::WH.

                                                                       {
   //get the hepmc product from the event
   edm::Handle<HepMCProduct> evt;
 
   bool product_exists = iEvent.getByToken(_hepmcCollection, evt);
   if (product_exists) {
     // get HepMC GenEvent
     const HepMC::GenEvent* myGenEvent = evt->GetEvent();
 
     if (_prodMode == "AUTO") {
       // for these prod modes, don't change what is set in BeginRun
       if (m_HiggsProdMode != HTXS::GGF && m_HiggsProdMode != HTXS::VBF && m_HiggsProdMode != HTXS::GG2ZH) {
         unsigned nWs = 0;
         unsigned nZs = 0;
         unsigned nTs = 0;
         unsigned nBs = 0;
         unsigned nHs = 0;
 
         HepMC::GenVertex* HSvtx = myGenEvent->signal_process_vertex();
 
         if (HSvtx) {
           for (auto ptcl : particles(HSvtx, HepMC::children)) {
             if (std::abs(ptcl->pdg_id()) == 24)
               ++nWs;
             if (ptcl->pdg_id() == 23)
               ++nZs;
             if (abs(ptcl->pdg_id()) == 6)
               ++nTs;
             if (abs(ptcl->pdg_id()) == 5)
               ++nBs;
             if (ptcl->pdg_id() == 25)
               ++nHs;
           }
         }
 
         if (nZs == 1 && nHs == 1 && (nWs + nTs) == 0) {
           m_HiggsProdMode = HTXS::QQ2ZH;
         } else if (nWs == 1 && nHs == 1 && (nZs + nTs) == 0) {
           m_HiggsProdMode = HTXS::WH;
         } else if (nTs == 2 && nHs == 1 && nZs == 0) {
           m_HiggsProdMode = HTXS::TTH;
         } else if (nTs == 1 && nHs == 1 && nZs == 0) {
           m_HiggsProdMode = HTXS::TH;
         } else if (nBs == 2 && nHs == 1 && nZs == 0) {
           m_HiggsProdMode = HTXS::BBH;
         }
 
         _HTXS->setHiggsProdMode(m_HiggsProdMode);
       }
     }
 
     if (_isFirstEvent) {
       _analysisHandler.addAnalysis(_HTXS);
 
       // set the production mode if not done already
       if (_prodMode == "GGF")
         m_HiggsProdMode = HTXS::GGF;
       else if (_prodMode == "VBF")
         m_HiggsProdMode = HTXS::VBF;
       else if (_prodMode == "WH")
         m_HiggsProdMode = HTXS::WH;
       else if (_prodMode == "ZH")
         m_HiggsProdMode = HTXS::QQ2ZH;
       else if (_prodMode == "QQ2ZH")
         m_HiggsProdMode = HTXS::QQ2ZH;
       else if (_prodMode == "GG2ZH")
         m_HiggsProdMode = HTXS::GG2ZH;
       else if (_prodMode == "TTH")
         m_HiggsProdMode = HTXS::TTH;
       else if (_prodMode == "BBH")
         m_HiggsProdMode = HTXS::BBH;
       else if (_prodMode == "TH")
         m_HiggsProdMode = HTXS::TH;
       else if (_prodMode == "AUTO") {
         edm::LogInfo("HTXSRivetProducer")
             << "Using AUTO for HiggsProdMode, found it to be: " << m_HiggsProdMode << "\n";
         edm::LogInfo("HTXSRivetProducer")
             << "(UNKNOWN=0, GGF=1, VBF=2, WH=3, QQ2ZH=4, GG2ZH=5, TTH=6, BBH=7, TH=8)" << endl;
       } else {
         throw cms::Exception("HTXSRivetProducer")
             << "ProductionMode must be one of: GGF,VBF,WH,ZH,QQ2ZH,GG2ZH,TTH,BBH,TH,AUTO ";
       }
       _HTXS->setHiggsProdMode(m_HiggsProdMode);
 
       // at this point the production mode must be known
       if (m_HiggsProdMode == HTXS::UNKNOWN) {
         edm::LogInfo("HTXSRivetProducer") << "HTXSRivetProducer WARNING: HiggsProduction mode is UNKNOWN" << endl;
       }
 
       // initialize rivet analysis
       _analysisHandler.init(*myGenEvent);
       _isFirstEvent = false;
     }
 
     // classify the event
     Rivet::HiggsClassification rivet_cat = _HTXS->classifyEvent(*myGenEvent, m_HiggsProdMode);
     cat_ = HTXS::Rivet2Root(rivet_cat);
 
     unique_ptr<HTXS::HiggsClassification> cat(new HTXS::HiggsClassification(cat_));
 
     iEvent.put(std::move(cat), "HiggsClassification");
   }
 }

Member Data Documentation

Rivet::AnalysisHandler HTXSRivetProducer::_analysisHandler

private

Definition at line 55 of file HTXSRivetProducer.cc.

edm::EDGetTokenT<edm::HepMCProduct> HTXSRivetProducer::_hepmcCollection

private

Definition at line 52 of file HTXSRivetProducer.cc.

Rivet::HiggsTemplateCrossSections* HTXSRivetProducer::_HTXS

private

Definition at line 56 of file HTXSRivetProducer.cc.

bool HTXSRivetProducer::_isFirstEvent

private

Definition at line 58 of file HTXSRivetProducer.cc.

edm::EDGetTokenT<LHERunInfoProduct> HTXSRivetProducer::_lheRunInfo

private

Definition at line 53 of file HTXSRivetProducer.cc.

std::string HTXSRivetProducer::_prodMode

private

Definition at line 59 of file HTXSRivetProducer.cc.

HTXS::HiggsClassification HTXSRivetProducer::cat_

private

Definition at line 62 of file HTXSRivetProducer.cc.

HTXS::HiggsProdMode HTXSRivetProducer::m_HiggsProdMode

private

Definition at line 60 of file HTXSRivetProducer.cc.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation