IsolatedParticlesGeneratedJets Class Reference

#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 ()
	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 &	itemsToGetFromEvent () 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
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

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

Definition at line 34 of file IsolatedParticlesGeneratedJets.cc.

References funct::abs(), ALCARECOTkAlJpsiMuMu_cff::charge, clearTreeVectors(), edm::EDConsumerBase::consumesCollector(), gather_cfg::cout, debug, reco::LeafCandidate::energy(), stringResolutionProvider_cfi::eta, reco::LeafCandidate::eta(), patPFMETCorrections_cff::genJets, GenHFHadronMatcher_cfi::genParticles, edm::Event::getByToken(), JetMatchingTools::getGenParticles(), AlCaHLTBitMon_ParallelJobs::p, cosmictrackingParticleSelector_cfi::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_, and tok_parts_.

                                                                                                {
  
  //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, and t_jetTrkPt.

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();
}

void IsolatedParticlesGeneratedJets::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 170 of file IsolatedParticlesGeneratedJets.cc.

References DEFINE_FWK_MODULE.

Referenced by o2olib.O2ORunMgr::executeJob().

                                            {
}

Member Data Documentation

bool IsolatedParticlesGeneratedJets::debug

private

Definition at line 49 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), rrapi.RRApi::dprint(), util.rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), util.rrapi.RRApi::get(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), IsolatedParticlesGeneratedJets(), and pkg.AbstractPkg::write().

edm::Service<TFileService> IsolatedParticlesGeneratedJets::fs

private

Definition at line 50 of file IsolatedParticlesGeneratedJets.h.

Referenced by BookHistograms().

std::vector<double>* IsolatedParticlesGeneratedJets::t_gjetE

private

Definition at line 57 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

std::vector<double> * IsolatedParticlesGeneratedJets::t_gjetEta

private

Definition at line 57 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

std::vector<int>* IsolatedParticlesGeneratedJets::t_gjetN

private

Definition at line 56 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

std::vector<double> * IsolatedParticlesGeneratedJets::t_gjetPhi

private

Definition at line 57 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

std::vector<double> * IsolatedParticlesGeneratedJets::t_gjetPt

private

Definition at line 57 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

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

private

Definition at line 63 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

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

private

Definition at line 60 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

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

private

Definition at line 58 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

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

private

Definition at line 62 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

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

private

Definition at line 61 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

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

private

Definition at line 59 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), BookHistograms(), and clearTreeVectors().

edm::EDGetTokenT<reco::GenJetCollection> IsolatedParticlesGeneratedJets::tok_jets_

private

Definition at line 53 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), and IsolatedParticlesGeneratedJets().

edm::EDGetTokenT<reco::GenParticleCollection> IsolatedParticlesGeneratedJets::tok_parts_

private

Definition at line 54 of file IsolatedParticlesGeneratedJets.h.

Referenced by analyze(), and IsolatedParticlesGeneratedJets().

TTree* IsolatedParticlesGeneratedJets::tree

private

Definition at line 51 of file IsolatedParticlesGeneratedJets.h.

Referenced by eventstfile.Events::__iter__(), svgfig.SVG::__str__(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::beginLoop(), MetTreeProducer.MetTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::fillTree(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::process(), MetTreeProducer.MetTreeProducer::process(), python.cmstools.EventTree::SetAlias(), and eventstfile.Events::to().