#include <Calibration/IsolatedParticles/plugins/IsolatedParticlesGeneratedJets.cc>

Inheritance diagram for IsolatedParticlesGeneratedJets:

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

	~IsolatedParticlesGeneratedJets ()

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
std::vector< ConsumesInfo >	consumesInfo () const

	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

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

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

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

virtual void	beginJob ()

void	BookHistograms ()

void	clearTreeVectors ()

virtual void	endJob ()

Private Attributes
bool	debug

edm::Service< TFileService >	fs

std::vector< double > *	t_gjetE

std::vector< double > *	t_gjetEta

std::vector< int > *	t_gjetN

std::vector< double > *	t_gjetPhi

std::vector< double > *	t_gjetPt

std::vector< std::vector < double > > *	t_jetTrkCharge

std::vector< std::vector < double > > *	t_jetTrkEta

std::vector< std::vector < double > > *	t_jetTrkP

std::vector< std::vector < double > > *	t_jetTrkPdg

std::vector< std::vector < double > > *	t_jetTrkPhi

std::vector< std::vector < double > > *	t_jetTrkPt

edm::EDGetTokenT < reco::GenJetCollection >	tok_jets_

edm::EDGetTokenT < reco::GenParticleCollection >	tok_parts_

TTree *	tree

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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 35 of file IsolatedParticlesGeneratedJets.h.

Constructor & Destructor Documentation

IsolatedParticlesGeneratedJets::IsolatedParticlesGeneratedJets ( const edm::ParameterSet & iConfig )

explicit

Definition at line 22 of file IsolatedParticlesGeneratedJets.cc.

References debug, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), tok_jets_, and tok_parts_.

                                                                                              {
 
   debug      = iConfig.getUntrackedParameter<bool>  ("Debug", false);
   tok_jets_  = consumes<reco::GenJetCollection>(iConfig.getParameter<edm::InputTag>("JetSource"));
   tok_parts_ = consumes<reco::GenParticleCollection>(iConfig.getParameter<edm::InputTag>("ParticleSource"));
 }

IsolatedParticlesGeneratedJets::~IsolatedParticlesGeneratedJets ( )

Definition at line 30 of file IsolatedParticlesGeneratedJets.cc.

                                                                 {
 
 }

Member Function Documentation

void IsolatedParticlesGeneratedJets::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 34 of file IsolatedParticlesGeneratedJets.cc.

References funct::abs(), RecoTauCleanerPlugins::charge, clearTreeVectors(), edm::EDConsumerBase::consumesCollector(), gather_cfg::cout, debug, reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), eta(), genParticleCandidates2GenParticles_cfi::genParticles, edm::Event::getByToken(), JetMatchingTools::getGenParticles(), HLT_25ns14e33_v1_cff::njets, AlCaHLTBitMon_ParallelJobs::p, benchmark_cfg::pdgId, phi, reco::LeafCandidate::phi(), EnergyCorrector::pt, reco::LeafCandidate::pt(), t_gjetE, t_gjetEta, t_gjetN, t_gjetPhi, t_gjetPt, t_jetTrkCharge, t_jetTrkEta, t_jetTrkP, t_jetTrkPdg, t_jetTrkPhi, t_jetTrkPt, tok_jets_, tok_parts_, and tree.

                                                                                                 {
   
   //using namespace edm;
   clearTreeVectors();
 
   //=== genJet information
   edm::Handle<reco::GenJetCollection> genJets;
   iEvent.getByToken(tok_jets_, genJets);
 
   //=== genJet information
   edm::Handle<reco::GenParticleCollection> genParticles;
   iEvent.getByToken(tok_parts_, genParticles);
 
   JetMatchingTools jetMatching (iEvent, consumesCollector());
   std::vector <std::vector <const reco::GenParticle*> > genJetConstituents (genJets->size());
 
   int njets = 0;
   for (unsigned iGenJet = 0; iGenJet < genJets->size(); ++iGenJet) {
     const reco::GenJet& genJet = (*genJets) [iGenJet];
 
     double genJetE   = genJet.energy();
     double genJetPt  = genJet.pt();
     double genJetEta = genJet.eta();
     double genJetPhi = genJet.phi();
 
     if( genJetPt> 30.0 && std::abs(genJetEta)<3.0 ) {
 
       njets++;
 
       std::vector <const reco::GenParticle*> genJetConstituents = jetMatching.getGenParticles ((*genJets) [iGenJet]);
 
       std::vector<double> v_trkP, v_trkPt, v_trkEta, v_trkPhi, v_trkPdg, v_trkCharge;
     
       if(debug) std::cout<<"Jet(pt,Eta,Phi) "<<genJetPt<<" "<<genJetEta<<" "<<genJetPhi <<std::endl;
       for(unsigned int ic=0; ic<genJetConstituents.size(); ic++) {
 
     if(debug) {
       std::cout << "p,pt,eta,phi "<<genJetConstituents[ic]->p()<<" "<<genJetConstituents[ic]->pt()
             <<" "<<genJetConstituents[ic]->eta()<<" "<<genJetConstituents[ic]->phi()
             <<std::endl;
     }
 
     v_trkP.push_back(genJetConstituents[ic]->p());
     v_trkPt.push_back(genJetConstituents[ic]->pt());
     v_trkEta.push_back(genJetConstituents[ic]->eta());
     v_trkPhi.push_back(genJetConstituents[ic]->phi());
     v_trkPdg.push_back(genJetConstituents[ic]->pdgId());
     v_trkCharge.push_back(genJetConstituents[ic]->charge());
 
       } //loop over genjet constituents
 
       t_gjetE   ->push_back(genJetE  );
       t_gjetPt  ->push_back(genJetPt );
       t_gjetEta ->push_back(genJetEta);
       t_gjetPhi ->push_back(genJetPhi);
       
       t_jetTrkP   ->push_back(v_trkP  );
       t_jetTrkPt  ->push_back(v_trkPt );
       t_jetTrkEta ->push_back(v_trkEta);
       t_jetTrkPhi ->push_back(v_trkPhi);
       t_jetTrkPdg ->push_back(v_trkPdg);
       t_jetTrkCharge ->push_back(v_trkCharge);
 
     } // if jetPt>30
 
   } //loop over genjets
 
   t_gjetN->push_back(njets);
 
   unsigned int indx = 0;
   for(reco::GenParticleCollection::const_iterator ig = genParticles->begin(); ig!= genParticles->end(); ++ig,++indx) {
  
     if (debug)
       std::cout << "Track " << indx << " Status " << ig->status() << " charge "
         << ig->charge() << " pdgId " << ig->pdgId() << " mass "
         << ig->mass() << " P " << ig->momentum() << " E "
         << ig->energy() << " Origin " << ig->vertex() << std::endl;
   }
 
 
   tree->Fill();
 }

void IsolatedParticlesGeneratedJets::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 117 of file IsolatedParticlesGeneratedJets.cc.

References BookHistograms().

                                               {
 
   BookHistograms();
 
 }

void IsolatedParticlesGeneratedJets::BookHistograms ( )

private

Definition at line 138 of file IsolatedParticlesGeneratedJets.cc.

References fs, TFileService::make(), t_gjetE, t_gjetEta, t_gjetN, t_gjetPhi, t_gjetPt, t_jetTrkCharge, t_jetTrkEta, t_jetTrkP, t_jetTrkPdg, t_jetTrkPhi, t_jetTrkPt, and tree.

Referenced by beginJob().

                                                    {
 
   tree = fs->make<TTree>("tree", "tree");
 
   t_gjetN     = new std::vector<int>   ();
   t_gjetE     = new std::vector<double>();
   t_gjetPt    = new std::vector<double>();
   t_gjetEta   = new std::vector<double>();
   t_gjetPhi   = new std::vector<double>();
 
   t_jetTrkP   = new std::vector<std::vector<double> >();
   t_jetTrkPt  = new std::vector<std::vector<double> >();
   t_jetTrkEta = new std::vector<std::vector<double> >();
   t_jetTrkPhi = new std::vector<std::vector<double> >();
   t_jetTrkPdg = new std::vector<std::vector<double> >();
   t_jetTrkCharge = new std::vector<std::vector<double> >();
 
   tree->Branch("t_gjetN",     "vector<int>",             &t_gjetN);
   tree->Branch("t_gjetE",     "vector<double>",          &t_gjetE);
   tree->Branch("t_gjetPt",    "vector<double>",          &t_gjetPt);
   tree->Branch("t_gjetEta",   "vector<double>",          &t_gjetEta);
   tree->Branch("t_gjetPhi",   "vector<double>",          &t_gjetPhi);
 
   tree->Branch("t_jetTrkP",   "vector<vector<double> >", &t_jetTrkP);
   tree->Branch("t_jetTrkPt",  "vector<vector<double> >", &t_jetTrkPt);
   tree->Branch("t_jetTrkEta", "vector<vector<double> >", &t_jetTrkEta);
   tree->Branch("t_jetTrkPhi", "vector<vector<double> >", &t_jetTrkPhi);
   tree->Branch("t_jetTrkPdg", "vector<vector<double> >", &t_jetTrkPdg);
   tree->Branch("t_jetTrkCharge", "vector<vector<double> >", &t_jetTrkCharge);
 
 }

void IsolatedParticlesGeneratedJets::clearTreeVectors ( )

private

Definition at line 123 of file IsolatedParticlesGeneratedJets.cc.

References t_gjetE, t_gjetEta, t_gjetN, t_gjetPhi, t_gjetPt, t_jetTrkCharge, t_jetTrkEta, t_jetTrkP, t_jetTrkPdg, t_jetTrkPhi, and t_jetTrkPt.

Referenced by analyze().

                                                       {
   t_gjetN     ->clear();
   t_gjetE     ->clear();
   t_gjetPt    ->clear();
   t_gjetEta   ->clear();
   t_gjetPhi   ->clear();
 
   t_jetTrkP   ->clear();
   t_jetTrkPt  ->clear();
   t_jetTrkEta ->clear();
   t_jetTrkPhi ->clear();
   t_jetTrkPdg ->clear();
   t_jetTrkCharge ->clear();
 }