#include <AnalysisRootpleProducerOnlyMC.h>

Inheritance diagram for AnalysisRootpleProducerOnlyMC:

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

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

virtual void	beginJob ()

virtual void	endJob ()

void	fillChargedJet (float, float, float, float)

void	fillEventInfo (int)

void	fillInclusiveJet (float, float, float, float)

void	fillMCParticles (float, float, float, float)

void	store ()

virtual	~AnalysisRootpleProducerOnlyMC ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
	EDConsumerBase ()

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

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

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

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Attributes
TTree *	AnalysisTree

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

TClonesArray *	ChargedJet

edm::Handle < reco::GenJetCollection >	ChgGenJetsHandle

edm::InputTag	chgGenPartCollName

edm::InputTag	chgJetCollName

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::InputTag	genJetCollName

edm::Handle < reco::GenJetCollection >	GenJetsHandle

TClonesArray *	InclusiveJet

edm::InputTag	mcEvent

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

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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)

Detailed Description

Definition at line 31 of file AnalysisRootpleProducerOnlyMC.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 68 of file AnalysisRootpleProducerOnlyMC.cc.

References edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

 {
   mcEvent = pset.getUntrackedParameter<InputTag>("MCEvent",std::string(""));
   genJetCollName = pset.getUntrackedParameter<InputTag>("GenJetCollectionName",std::string(""));
   chgJetCollName = pset.getUntrackedParameter<InputTag>("ChgGenJetCollectionName",std::string(""));
   chgGenPartCollName = pset.getUntrackedParameter<InputTag>("ChgGenPartCollectionName",std::string(""));
 
   piG = acos(-1.);
   NumberMCParticles=0;
   NumberInclusiveJet=0;
   NumberChargedJet=0;
 }

virtual AnalysisRootpleProducerOnlyMC::~AnalysisRootpleProducerOnlyMC ( )

inlinevirtual

Definition at line 37 of file AnalysisRootpleProducerOnlyMC.h.

37 {}

Member Function Documentation

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

virtual

Implements edm::EDAnalyzer.

Definition at line 129 of file AnalysisRootpleProducerOnlyMC.cc.

References edm::Event::getByLabel().

 {
 
   e.getByLabel( mcEvent           , EvtHandle        ) ;
   e.getByLabel( chgGenPartCollName, CandHandleMC     );
   e.getByLabel( chgJetCollName    , ChgGenJetsHandle );
   e.getByLabel( genJetCollName    , 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->size()){
 
     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->size()){
     
     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->size()){
 
     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();
 }

void AnalysisRootpleProducerOnlyMC::beginJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 82 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);
 }

void AnalysisRootpleProducerOnlyMC::endJob ( void )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 210 of file AnalysisRootpleProducerOnlyMC.cc.

211 {

212 }

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

Definition at line 59 of file AnalysisRootpleProducerOnlyMC.cc.

References eta(), AlCaHLTBitMon_ParallelJobs::p, phi, and RecoTauCleanerPlugins::pt.

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

void AnalysisRootpleProducerOnlyMC::fillEventInfo ( int e )

Definition at line 36 of file AnalysisRootpleProducerOnlyMC.cc.

References alignCSCRings::e.

 {
   EventKind = e;
 }

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

Definition at line 50 of file AnalysisRootpleProducerOnlyMC.cc.

References eta(), AlCaHLTBitMon_ParallelJobs::p, phi, and RecoTauCleanerPlugins::pt.

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

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

Definition at line 41 of file AnalysisRootpleProducerOnlyMC.cc.

References eta(), AlCaHLTBitMon_ParallelJobs::p, phi, and RecoTauCleanerPlugins::pt.

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

void AnalysisRootpleProducerOnlyMC::store ( )

Definition at line 27 of file AnalysisRootpleProducerOnlyMC.cc.

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

Member Data Documentation

TTree* AnalysisRootpleProducerOnlyMC::AnalysisTree

private

Definition at line 66 of file AnalysisRootpleProducerOnlyMC.h.

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

private

Definition at line 57 of file AnalysisRootpleProducerOnlyMC.h.

TClonesArray* AnalysisRootpleProducerOnlyMC::ChargedJet

private

Definition at line 86 of file AnalysisRootpleProducerOnlyMC.h.

edm::Handle< reco::GenJetCollection > AnalysisRootpleProducerOnlyMC::ChgGenJetsHandle

private

Definition at line 59 of file AnalysisRootpleProducerOnlyMC.h.

edm::InputTag AnalysisRootpleProducerOnlyMC::chgGenPartCollName

private

Definition at line 54 of file AnalysisRootpleProducerOnlyMC.h.

edm::InputTag AnalysisRootpleProducerOnlyMC::chgJetCollName

private

Definition at line 53 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::EtaCJ[NCJMAX]

private

Definition at line 80 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::EtaEHJ[NEHJMAX]

private

Definition at line 82 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::EtaIJ[NIJMAX]

private

Definition at line 79 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::EtaMC[NMCPMAX]

private

Definition at line 77 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::EtaTJ[NTJMAX]

private

Definition at line 81 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::EtaTK[NTKMAX]

private

Definition at line 78 of file AnalysisRootpleProducerOnlyMC.h.

int AnalysisRootpleProducerOnlyMC::EventKind

private

Definition at line 75 of file AnalysisRootpleProducerOnlyMC.h.

edm::Handle< edm::HepMCProduct > AnalysisRootpleProducerOnlyMC::EvtHandle

private

Definition at line 56 of file AnalysisRootpleProducerOnlyMC.h.

edm::Service<TFileService> AnalysisRootpleProducerOnlyMC::fs

private

Definition at line 64 of file AnalysisRootpleProducerOnlyMC.h.

edm::InputTag AnalysisRootpleProducerOnlyMC::genJetCollName

private

Definition at line 52 of file AnalysisRootpleProducerOnlyMC.h.

edm::Handle< reco::GenJetCollection > AnalysisRootpleProducerOnlyMC::GenJetsHandle

private

Definition at line 58 of file AnalysisRootpleProducerOnlyMC.h.

TClonesArray* AnalysisRootpleProducerOnlyMC::InclusiveJet

private

Definition at line 85 of file AnalysisRootpleProducerOnlyMC.h.

edm::InputTag AnalysisRootpleProducerOnlyMC::mcEvent

private

Definition at line 51 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::MomentumCJ[NCJMAX]

private

Definition at line 80 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::MomentumEHJ[NEHJMAX]

private

Definition at line 82 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::MomentumIJ[NIJMAX]

private

Definition at line 79 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::MomentumMC[NMCPMAX]

private

Definition at line 77 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::MomentumTJ[NTJMAX]

private

Definition at line 81 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::MomentumTK[NTKMAX]

private

Definition at line 78 of file AnalysisRootpleProducerOnlyMC.h.

TClonesArray* AnalysisRootpleProducerOnlyMC::MonteCarlo

private

Definition at line 84 of file AnalysisRootpleProducerOnlyMC.h.

const int AnalysisRootpleProducerOnlyMC::NCJMAX = 10000

staticprivate

Definition at line 71 of file AnalysisRootpleProducerOnlyMC.h.

const int AnalysisRootpleProducerOnlyMC::NEHJMAX = 10000

staticprivate

Definition at line 73 of file AnalysisRootpleProducerOnlyMC.h.

const int AnalysisRootpleProducerOnlyMC::NIJMAX = 10000

staticprivate

Definition at line 70 of file AnalysisRootpleProducerOnlyMC.h.

const int AnalysisRootpleProducerOnlyMC::NMCPMAX = 10000

staticprivate

Definition at line 68 of file AnalysisRootpleProducerOnlyMC.h.

const int AnalysisRootpleProducerOnlyMC::NTJMAX = 10000

staticprivate

Definition at line 72 of file AnalysisRootpleProducerOnlyMC.h.

const int AnalysisRootpleProducerOnlyMC::NTKMAX = 10000

staticprivate

Definition at line 69 of file AnalysisRootpleProducerOnlyMC.h.

int AnalysisRootpleProducerOnlyMC::NumberCaloJet

private

Definition at line 75 of file AnalysisRootpleProducerOnlyMC.h.

int AnalysisRootpleProducerOnlyMC::NumberChargedJet

private

Definition at line 75 of file AnalysisRootpleProducerOnlyMC.h.

int AnalysisRootpleProducerOnlyMC::NumberInclusiveJet

private

Definition at line 75 of file AnalysisRootpleProducerOnlyMC.h.

int AnalysisRootpleProducerOnlyMC::NumberMCParticles

private

Definition at line 75 of file AnalysisRootpleProducerOnlyMC.h.

int AnalysisRootpleProducerOnlyMC::NumberTracks

private

Definition at line 75 of file AnalysisRootpleProducerOnlyMC.h.

int AnalysisRootpleProducerOnlyMC::NumberTracksJet

private

Definition at line 75 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::PhiCJ[NCJMAX]

private

Definition at line 80 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::PhiEHJ[NEHJMAX]

private

Definition at line 82 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::PhiIJ[NIJMAX]

private

Definition at line 79 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::PhiMC[NMCPMAX]

private

Definition at line 77 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::PhiTJ[NTJMAX]

private

Definition at line 81 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::PhiTK[NTKMAX]

private

Definition at line 78 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::piG

private

Definition at line 62 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::TransverseMomentumCJ[NCJMAX]

private

Definition at line 80 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::TransverseMomentumEHJ[NEHJMAX]

private

Definition at line 82 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::TransverseMomentumIJ[NIJMAX]

private

Definition at line 79 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::TransverseMomentumMC[NMCPMAX]

private

Definition at line 77 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::TransverseMomentumTJ[NTJMAX]

private

Definition at line 81 of file AnalysisRootpleProducerOnlyMC.h.

float AnalysisRootpleProducerOnlyMC::TransverseMomentumTK[NTKMAX]

private

Definition at line 78 of file AnalysisRootpleProducerOnlyMC.h.

Public Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation