#include <AnalysisRootpleProducerOnlyMC.h>

Inheritance diagram for AnalysisRootpleProducerOnlyMC:

Public Member Functions
	AnalysisRootpleProducerOnlyMC (const edm::ParameterSet &)

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

void	beginJob () override

void	endJob () override

void	fillChargedJet (float, float, float, float)

void	fillEventInfo (int)

void	fillInclusiveJet (float, float, float, float)

void	fillMCParticles (float, float, float, float)

void	store ()

	~AnalysisRootpleProducerOnlyMC () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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 &	operator= (EDConsumerBase &&)=default

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

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 Attributes
TTree *	AnalysisTree

edm::Handle< std::vector< reco::GenParticle > >	CandHandleMC

TClonesArray *	ChargedJet

edm::Handle< reco::GenJetCollection >	ChgGenJetsHandle

edm::EDGetTokenT< std::vector< reco::GenParticle > >	chgGenPartCollToken

edm::EDGetTokenT< reco::GenJetCollection >	chgJetCollToken

float	EtaCJ [NCJMAX]

float	EtaEHJ [NEHJMAX]

float	EtaIJ [NIJMAX]

float	EtaMC [NMCPMAX]

float	EtaTJ [NTJMAX]

float	EtaTK [NTKMAX]

int	EventKind

edm::Handle< edm::HepMCProduct >	EvtHandle

edm::Service< TFileService >	fs

edm::EDGetTokenT< reco::GenJetCollection >	genJetCollToken

edm::Handle< reco::GenJetCollection >	GenJetsHandle

TClonesArray *	InclusiveJet

edm::EDGetTokenT< edm::HepMCProduct >	mcEventToken

float	MomentumCJ [NCJMAX]

float	MomentumEHJ [NEHJMAX]

float	MomentumIJ [NIJMAX]

float	MomentumMC [NMCPMAX]

float	MomentumTJ [NTJMAX]

float	MomentumTK [NTKMAX]

TClonesArray *	MonteCarlo

int	NumberCaloJet

int	NumberChargedJet

int	NumberInclusiveJet

int	NumberMCParticles

int	NumberTracks

int	NumberTracksJet

float	PhiCJ [NCJMAX]

float	PhiEHJ [NEHJMAX]

float	PhiIJ [NIJMAX]

float	PhiMC [NMCPMAX]

float	PhiTJ [NTJMAX]

float	PhiTK [NTKMAX]

float	piG

float	TransverseMomentumCJ [NCJMAX]

float	TransverseMomentumEHJ [NEHJMAX]

float	TransverseMomentumIJ [NIJMAX]

float	TransverseMomentumMC [NMCPMAX]

float	TransverseMomentumTJ [NTJMAX]

float	TransverseMomentumTK [NTKMAX]

Static Private Attributes
static const int	NCJMAX = 10000

static const int	NEHJMAX = 10000

static const int	NIJMAX = 10000

static const int	NMCPMAX = 10000

static const int	NTJMAX = 10000

static const int	NTKMAX = 10000

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

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

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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

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

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

Detailed Description

Definition at line 31 of file AnalysisRootpleProducerOnlyMC.h.

Constructor & Destructor Documentation

◆ AnalysisRootpleProducerOnlyMC()

AnalysisRootpleProducerOnlyMC::AnalysisRootpleProducerOnlyMC ( const edm::ParameterSet & pset )

explicit

Definition at line 52 of file AnalysisRootpleProducerOnlyMC.cc.

                                                                                      {
   mcEventToken = consumes<edm::HepMCProduct>(pset.getUntrackedParameter<InputTag>("MCEvent", std::string("")));
   genJetCollToken =
       consumes<reco::GenJetCollection>(pset.getUntrackedParameter<InputTag>("GenJetCollectionName", std::string("")));
   chgJetCollToken = consumes<reco::GenJetCollection>(
       pset.getUntrackedParameter<InputTag>("ChgGenJetCollectionName", std::string("")));
   chgGenPartCollToken = consumes<std::vector<reco::GenParticle> >(
       pset.getUntrackedParameter<InputTag>("ChgGenPartCollectionName", std::string("")));
  
   piG = acos(-1.);
   NumberMCParticles = 0;
   NumberInclusiveJet = 0;
   NumberChargedJet = 0;
 }

References muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

◆ ~AnalysisRootpleProducerOnlyMC()

AnalysisRootpleProducerOnlyMC::~AnalysisRootpleProducerOnlyMC ( )

inlineoverride

Definition at line 34 of file AnalysisRootpleProducerOnlyMC.h.

34 {}

Member Function Documentation

◆ analyze()

void AnalysisRootpleProducerOnlyMC::analyze	(	const edm::Event &	e,
		const edm::EventSetup &
	)

overridevirtual

Implements edm::EDAnalyzer.

Definition at line 111 of file AnalysisRootpleProducerOnlyMC.cc.

                                                                              {
   e.getByToken(mcEventToken, EvtHandle);
   e.getByToken(chgGenPartCollToken, CandHandleMC);
   e.getByToken(chgJetCollToken, ChgGenJetsHandle);
   e.getByToken(genJetCollToken, GenJetsHandle);
  
   const HepMC::GenEvent* Evt = EvtHandle->GetEvent();
  
   EventKind = Evt->signal_process_id();
  
   std::vector<math::XYZTLorentzVector> GenPart;
   std::vector<GenJet> ChgGenJetContainer;
   std::vector<GenJet> GenJetContainer;
  
   GenPart.clear();
   ChgGenJetContainer.clear();
   GenJetContainer.clear();
  
   ChargedJet->Clear();
   InclusiveJet->Clear();
   MonteCarlo->Clear();
  
   if (!ChgGenJetsHandle->empty()) {
     for (GenJetCollection::const_iterator it(ChgGenJetsHandle->begin()), itEnd(ChgGenJetsHandle->end()); it != itEnd;
          ++it) {
       ChgGenJetContainer.push_back(*it);
     }
  
     std::stable_sort(ChgGenJetContainer.begin(), ChgGenJetContainer.end(), GenJetSort());
  
     std::vector<GenJet>::const_iterator it(ChgGenJetContainer.begin()), itEnd(ChgGenJetContainer.end());
     for (int iChargedJet(0); it != itEnd; ++it, ++iChargedJet) {
       fillChargedJet(it->p(), it->pt(), it->eta(), it->phi());
       new ((*ChargedJet)[iChargedJet]) TLorentzVector(it->px(), it->py(), it->pz(), it->energy());
     }
   }
  
   if (!GenJetsHandle->empty()) {
     for (GenJetCollection::const_iterator it(GenJetsHandle->begin()), itEnd(GenJetsHandle->end()); it != itEnd; ++it) {
       GenJetContainer.push_back(*it);
     }
  
     std::stable_sort(GenJetContainer.begin(), GenJetContainer.end(), GenJetSort());
  
     std::vector<GenJet>::const_iterator it(GenJetContainer.begin()), itEnd(GenJetContainer.end());
     for (int iInclusiveJet(0); it != itEnd; ++it, ++iInclusiveJet) {
       fillInclusiveJet(it->p(), it->pt(), it->eta(), it->phi());
       new ((*InclusiveJet)[iInclusiveJet]) TLorentzVector(it->px(), it->py(), it->pz(), it->energy());
     }
   }
  
   if (!CandHandleMC->empty()) {
     for (vector<GenParticle>::const_iterator it(CandHandleMC->begin()), itEnd(CandHandleMC->end()); it != itEnd; it++) {
       GenPart.push_back(it->p4());
     }
  
     std::stable_sort(GenPart.begin(), GenPart.end(), GreaterPt());
  
     std::vector<math::XYZTLorentzVector>::const_iterator it(GenPart.begin()), itEnd(GenPart.end());
     for (int iMonteCarlo(0); it != itEnd; ++it, ++iMonteCarlo) {
       fillMCParticles(it->P(), it->Pt(), it->Eta(), it->Phi());
       new ((*MonteCarlo)[iMonteCarlo]) TLorentzVector(it->Px(), it->Py(), it->Pz(), it->E());
     }
   }
  
   store();
 }

References MillePedeFileConverter_cfg::e, nanoDQM_cfi::GenPart, and HLT_FULL_cff::MonteCarlo.

◆ beginJob()

void AnalysisRootpleProducerOnlyMC::beginJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 67 of file AnalysisRootpleProducerOnlyMC.cc.

                                              {
   // use TFileService for output to root file
   AnalysisTree = fs->make<TTree>("AnalysisTree", "MBUE Analysis Tree ");
  
   // process type
   AnalysisTree->Branch("EventKind", &EventKind, "EventKind/I");
  
   // store p, pt, eta, phi for particles and jets
  
   // GenParticles at hadron level
   AnalysisTree->Branch("NumberMCParticles", &NumberMCParticles, "NumberMCParticles/I");
   AnalysisTree->Branch("MomentumMC", MomentumMC, "MomentumMC[NumberMCParticles]/F");
   AnalysisTree->Branch("TransverseMomentumMC", TransverseMomentumMC, "TransverseMomentumMC[NumberMCParticles]/F");
   AnalysisTree->Branch("EtaMC", EtaMC, "EtaMC[NumberMCParticles]/F");
   AnalysisTree->Branch("PhiMC", PhiMC, "PhiMC[NumberMCParticles]/F");
  
   // GenJets
   AnalysisTree->Branch("NumberInclusiveJet", &NumberInclusiveJet, "NumberInclusiveJet/I");
   AnalysisTree->Branch("MomentumIJ", MomentumIJ, "MomentumIJ[NumberInclusiveJet]/F");
   AnalysisTree->Branch("TrasverseMomentumIJ", TransverseMomentumIJ, "TransverseMomentumIJ[NumberInclusiveJet]/F");
   AnalysisTree->Branch("EtaIJ", EtaIJ, "EtaIJ[NumberInclusiveJet]/F");
   AnalysisTree->Branch("PhiIJ", PhiIJ, "PhiIJ[NumberInclusiveJet]/F");
  
   // jets from charged GenParticles
   AnalysisTree->Branch("NumberChargedJet", &NumberChargedJet, "NumberChargedJet/I");
   AnalysisTree->Branch("MomentumCJ", MomentumCJ, "MomentumCJ[NumberChargedJet]/F");
   AnalysisTree->Branch("TrasverseMomentumCJ", TransverseMomentumCJ, "TransverseMomentumCJ[NumberChargedJet]/F");
   AnalysisTree->Branch("EtaCJ", EtaCJ, "EtaCJ[NumberChargedJet]/F");
   AnalysisTree->Branch("PhiCJ", PhiCJ, "PhiCJ[NumberChargedJet]/F");
  
   // alternative storage method:
   // save TClonesArrays of TLorentzVectors
   // i.e. store 4-vectors of particles and jets
  
   MonteCarlo = new TClonesArray("TLorentzVector", 10000);
   AnalysisTree->Branch("MonteCarlo", "TClonesArray", &MonteCarlo, 128000, 0);
  
   InclusiveJet = new TClonesArray("TLorentzVector", 10000);
   AnalysisTree->Branch("InclusiveJet", "TClonesArray", &InclusiveJet, 128000, 0);
  
   ChargedJet = new TClonesArray("TLorentzVector", 10000);
   AnalysisTree->Branch("ChargedJet", "TClonesArray", &ChargedJet, 128000, 0);
 }

References HLT_FULL_cff::MonteCarlo.

◆ endJob()

void AnalysisRootpleProducerOnlyMC::endJob ( void )

overridevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 179 of file AnalysisRootpleProducerOnlyMC.cc.

179 {}

◆ fillChargedJet()

void AnalysisRootpleProducerOnlyMC::fillChargedJet	(	float	p,
		float	pt,
		float	eta,
		float	phi
	)

Definition at line 44 of file AnalysisRootpleProducerOnlyMC.cc.

                                                                                           {
   MomentumCJ[NumberChargedJet] = p;
   TransverseMomentumCJ[NumberChargedJet] = pt;
   EtaCJ[NumberChargedJet] = eta;
   PhiCJ[NumberChargedJet] = phi;
   NumberChargedJet++;
 }

References PVValHelper::eta, AlCaHLTBitMon_ParallelJobs::p, and DiDispStaMuonMonitor_cfi::pt.

◆ fillEventInfo()

void AnalysisRootpleProducerOnlyMC::fillEventInfo ( int e )

Definition at line 26 of file AnalysisRootpleProducerOnlyMC.cc.

26 { EventKind = e; }

References MillePedeFileConverter_cfg::e.

◆ fillInclusiveJet()

void AnalysisRootpleProducerOnlyMC::fillInclusiveJet	(	float	p,
		float	pt,
		float	eta,
		float	phi
	)

Definition at line 36 of file AnalysisRootpleProducerOnlyMC.cc.

                                                                                             {
   MomentumIJ[NumberInclusiveJet] = p;
   TransverseMomentumIJ[NumberInclusiveJet] = pt;
   EtaIJ[NumberInclusiveJet] = eta;
   PhiIJ[NumberInclusiveJet] = phi;
   NumberInclusiveJet++;
 }

References PVValHelper::eta, AlCaHLTBitMon_ParallelJobs::p, and DiDispStaMuonMonitor_cfi::pt.

◆ fillMCParticles()

void AnalysisRootpleProducerOnlyMC::fillMCParticles	(	float	p,
		float	pt,
		float	eta,
		float	phi
	)

Definition at line 28 of file AnalysisRootpleProducerOnlyMC.cc.

                                                                                            {
   MomentumMC[NumberMCParticles] = p;
   TransverseMomentumMC[NumberMCParticles] = pt;
   EtaMC[NumberMCParticles] = eta;
   PhiMC[NumberMCParticles] = phi;
   NumberMCParticles++;
 }

References PVValHelper::eta, AlCaHLTBitMon_ParallelJobs::p, and DiDispStaMuonMonitor_cfi::pt.

◆ store()

void AnalysisRootpleProducerOnlyMC::store ( )

Definition at line 18 of file AnalysisRootpleProducerOnlyMC.cc.

                                           {
   AnalysisTree->Fill();
  
   NumberMCParticles = 0;
   NumberInclusiveJet = 0;
   NumberChargedJet = 0;
 }