#include <HiggsValidation.h>

Inheritance diagram for HiggsValidation:

Classes
class	MonitoredDecays

Public Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	bookHistograms (DQMStore::IBooker &i, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (const edm::Run &r, const edm::EventSetup &c) override

	HiggsValidation (const edm::ParameterSet &)

	~HiggsValidation () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &ev, edm::EventSetup const &es) final

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)

Private Member Functions
std::string	convert (int)

int	findHiggsDecayChannel (const HepMC::GenParticle , std::vector< HepMC::GenParticle > &decayprod)

Private Attributes
edm::ESHandle< HepPDT::ParticleDataTable >	fPDGTable
	PDT table. More...

edm::InputTag	hepmcCollection_

edm::EDGetTokenT< edm::HepMCProduct >	hepmcCollectionToken_

MonitorElement *	Higgs_eta

MonitorElement *	Higgs_mass

MonitorElement *	Higgs_pt

MonitorElement *	HiggsDecayChannels

std::vector< MonitorElement * >	HiggsDecayProd_eta

std::vector< MonitorElement * >	HiggsDecayProd_pt

MonitoredDecays *	monitoredDecays

MonitorElement *	nEvt

int	particle_id

std::string	particle_name

WeightManager	wmanager_

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

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 37 of file HiggsValidation.h.

Constructor & Destructor Documentation

HiggsValidation::HiggsValidation ( const edm::ParameterSet & iPSet )

explicit

Definition at line 21 of file HiggsValidation.cc.

References hepmcCollection_, hepmcCollectionToken_, and monitoredDecays.

                                                             : 
   wmanager_(iPSet,consumesCollector()),
   hepmcCollection_(iPSet.getParameter<edm::InputTag>("hepmcCollection")),
   particle_id(iPSet.getParameter<int>("pdg_id")),
   particle_name(iPSet.getParameter<std::string>("particleName"))
 {    
   monitoredDecays = new MonitoredDecays(iPSet);
   hepmcCollectionToken_=consumes<HepMCProduct>(hepmcCollection_);
 }

HiggsValidation::~HiggsValidation ( )

override

Definition at line 31 of file HiggsValidation.cc.

31 {}

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 72 of file HiggsValidation.cc.

References funct::abs(), MonitorElement::Fill(), findHiggsDecayChannel(), edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), hepmcCollectionToken_, Higgs_eta, Higgs_mass, Higgs_pt, HiggsDecayChannels, HiggsDecayProd_eta, HiggsDecayProd_pt, mps_fire::i, training_settings::idx, HiggsValidation::MonitoredDecays::isDecayParticle(), monitoredDecays, nEvt, particle_id, WeightManager::weight(), mps_merge::weight, and wmanager_.

 { 
   double weight =   wmanager_.weight(iEvent);  
   nEvt->Fill(0.5,weight);
   
   //Gathering the HepMCProduct information
   edm::Handle<HepMCProduct> evt;
   iEvent.getByToken(hepmcCollectionToken_, evt);
   
   //Get EVENT
   HepMC::GenEvent *myGenEvent = new HepMC::GenEvent(*(evt->GetEvent()));
 
   // loop over all particles
   bool filled = false;
   for(HepMC::GenEvent::particle_const_iterator iter = myGenEvent->particles_begin(); 
       iter!= myGenEvent->particles_end() && !filled; ++iter) {
     if(particle_id == std::abs((*iter)->pdg_id())){
       std::vector<HepMC::GenParticle*> decayprod;
       int channel = findHiggsDecayChannel(*iter,decayprod);
       HiggsDecayChannels->Fill(channel,weight);
       Higgs_pt->Fill((*iter)->momentum().perp(),weight); 
       Higgs_eta->Fill((*iter)->momentum().eta(),weight);
       Higgs_mass->Fill((*iter)->momentum().m(),weight);
       for(unsigned int i=0;i<decayprod.size();i++){
     int idx=monitoredDecays->isDecayParticle(decayprod.at(i)->pdg_id());
     if(0<=idx && idx<=(int)HiggsDecayProd_pt.size()){
       HiggsDecayProd_pt.at(idx)->Fill(decayprod.at(i)->momentum().perp(),weight);
       HiggsDecayProd_eta.at(idx)->Fill(decayprod.at(i)->momentum().eta(),weight);
     }
       }
       filled = true;
     }
   }
   
   delete myGenEvent;
   
 }//analyze

void HiggsValidation::bookHistograms	(	DQMStore::IBooker &	i,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 37 of file HiggsValidation.cc.

References DQMHelper::book1dHisto(), HiggsValidation::MonitoredDecays::channel(), HiggsValidation::MonitoredDecays::ConvertIndex(), dir, Higgs_eta, Higgs_mass, Higgs_pt, HiggsDecayChannels, HiggsDecayProd_eta, HiggsDecayProd_pt, training_settings::idx, monitoredDecays, HiggsValidation::MonitoredDecays::NDecayParticles(), nEvt, particle_name, MonitorElement::setBinLabel(), DQMStore::IBooker::setCurrentFolder(), HiggsValidation::MonitoredDecays::size(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                              {
 
   std::string dir="Generator/";
     dir+=particle_name;
     DQMHelper dqm(&i); i.setCurrentFolder(dir);
     
     // Number of analyzed events
     nEvt = dqm.book1dHisto("nEvt", "n analyzed Events", 1, 0., 1.,"bin","Number of Events");
     
     //decay type
     
     std::string channel = particle_name+"_DecayChannels";
     HiggsDecayChannels = dqm.book1dHisto(channel,particle_name+" decay channels",monitoredDecays->size(),0,monitoredDecays->size(),"Decay Channels","Number of Events");
     
     for(size_t j = 0; j < monitoredDecays->size(); ++j){
       HiggsDecayChannels->setBinLabel(1+j,monitoredDecays->channel(j));
     }
   
 
   //Kinematics 
     Higgs_pt = dqm.book1dHisto((particle_name+"_pt"),(particle_name+" p_{t}"),50,0,250,"P_{t}^{"+particle_name+"} (GeV)","Number of Events");
     Higgs_eta  = dqm.book1dHisto((particle_name+"_eta"),(particle_name+" #eta"),50,-5,5,"#eta^{"+particle_name+"}","Number of Events"); 
     Higgs_mass = dqm.book1dHisto((particle_name+"_m"),(particle_name+" M"),500,0,500,"M^{"+particle_name+"}","Number of Events");
 
   int idx=0;
   for(unsigned int j=0;j<monitoredDecays->NDecayParticles();j++){
     HiggsDecayProd_pt.push_back(dqm.book1dHisto((monitoredDecays->ConvertIndex(idx)+"_pt"),(monitoredDecays->ConvertIndex(idx)+" p_{t}"),50,0,250,"P_{t}^{"+monitoredDecays->ConvertIndex(idx)+"} (GeV)","Number of Events"));
     HiggsDecayProd_eta.push_back(dqm.book1dHisto((monitoredDecays->ConvertIndex(idx)+"_eta"),(monitoredDecays->ConvertIndex(idx)+" #eta"),50,-5,5,"#eta^{"+monitoredDecays->ConvertIndex(idx)+"}","Number of Events"));
     idx++;
   }
 
   return;
 }

std::string HiggsValidation::convert ( int )

private

void HiggsValidation::dqmBeginRun	(	const edm::Run &	r,
		const edm::EventSetup &	c
	)

override

Definition at line 33 of file HiggsValidation.cc.

References fPDGTable, and edm::EventSetup::getData().

                                                                          {
   c.getData( fPDGTable );
 }

int HiggsValidation::findHiggsDecayChannel	(	const HepMC::GenParticle *	genParticle,
		std::vector< HepMC::GenParticle * > &	decayprod
	)

private

Definition at line 110 of file HiggsValidation.cc.

References class-composition::children, monitoredDecays, parents, HiggsValidation::MonitoredDecays::position(), HiggsValidation::MonitoredDecays::stable(), and HiggsValidation::MonitoredDecays::undetermined().

Referenced by analyze().

                                                                                                                      {
   if(genParticle->status() == 1) return monitoredDecays->stable();
   std::vector<int> children;
   if ( genParticle->end_vertex() ) {
     HepMC::GenVertex::particle_iterator des;
     for(des = genParticle->end_vertex()->particles_begin(HepMC::descendants);
     des!= genParticle->end_vertex()->particles_end(HepMC::descendants);++des ) {
       
       if((*des)->pdg_id() == genParticle->pdg_id()) continue;
       
       HepMC::GenVertex::particle_iterator mother = (*des)->production_vertex()->particles_begin(HepMC::parents);
       if((*mother)->pdg_id() == genParticle->pdg_id()){
     children.push_back((*des)->pdg_id());
     decayprod.push_back((*des));
       }
     }
   }
   
   if(children.size() == 2 && children.at(0) != 0 && children.at(1) != 0){
     return monitoredDecays->position(children.at(0),children.at(1));
   }
   return monitoredDecays->undetermined();
 }