IsolatedParticlesGeneratedJets Class Reference

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

Inheritance diagram for IsolatedParticlesGeneratedJets:

Public Member Functions
	IsolatedParticlesGeneratedJets (const edm::ParameterSet &)
	~IsolatedParticlesGeneratedJets () override
Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () 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
ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESResolverIndex > 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::array< std::vector< ModuleDescription const *> *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)

Private Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	bookHistograms ()
void	clearTreeVectors ()
void	endJob () override

Private Attributes
const bool	debug_
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
const edm::EDGetTokenT< reco::GenJetCollection >	tok_jets_
const edm::EDGetTokenT< reco::GenParticleCollection >	tok_parts_
TTree *	tree_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
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)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: Studies properties of jets at generator level in context of isolated particles

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

Definition at line 49 of file IsolatedParticlesGeneratedJets.cc.

Constructor & Destructor Documentation

IsolatedParticlesGeneratedJets()

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

explicit

Definition at line 80 of file IsolatedParticlesGeneratedJets.cc.

References TFileService::kSharedResource.

    : 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"))) {
  usesResource(TFileService::kSharedResource);
}

~IsolatedParticlesGeneratedJets()

IsolatedParticlesGeneratedJets::~IsolatedParticlesGeneratedJets ( )

inlineoverride

Definition at line 52 of file IsolatedParticlesGeneratedJets.cc.

{}

Member Function Documentation

analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 95 of file IsolatedParticlesGeneratedJets.cc.

References funct::abs(), ALCARECOTkAlJpsiMuMu_cff::charge, clearTreeVectors(), edm::EDConsumerBase::consumesCollector(), debug_, reco::LeafCandidate::energy(), PVValHelper::eta, reco::LeafCandidate::eta(), l1tCaloJetHTTProducer_cfi::genJets, AJJGenJetFilter_cfi::genParticles, iEvent, generator_cfi::jetMatching, BTaggingMonitoring_cff::njets, AlCaHLTBitMon_ParallelJobs::p, EgammaValidation_cff::pdgId, phi, reco::LeafCandidate::phi(), DiDispStaMuonMonitor_cfi::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_)
        edm::LogVerbatim("IsoTrack") << "Jet(pt,Eta,Phi) " << genJetPt << " " << genJetEta << " " << genJetPhi;
      for (unsigned int ic = 0; ic < genJetConstituents.size(); ic++) {
        if (debug_)
          edm::LogVerbatim("IsoTrack") << "p,pt,eta,phi " << genJetConstituents[ic]->p() << " "
                                       << genJetConstituents[ic]->pt() << " " << genJetConstituents[ic]->eta() << " "
                                       << genJetConstituents[ic]->phi();

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

  if (debug_) {
    unsigned int indx = 0;
    reco::GenParticleCollection::const_iterator ig = genParticles->begin();
    for (; ig != genParticles->end(); ++ig, ++indx) {
      edm::LogVerbatim("IsoTrack") << "Track " << indx << " Status " << ig->status() << " charge " << ig->charge()
                                   << " pdgId " << ig->pdgId() << " mass " << ig->mass() << " P " << ig->momentum()
                                   << " E " << ig->energy() << " Origin " << ig->vertex();
    }
  }

  tree_->Fill();
}

beginJob()

void IsolatedParticlesGeneratedJets::beginJob ( )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 173 of file IsolatedParticlesGeneratedJets.cc.

References bookHistograms().

{ bookHistograms(); }

bookHistograms()

void IsolatedParticlesGeneratedJets::bookHistograms ( )

private

Definition at line 190 of file IsolatedParticlesGeneratedJets.cc.

References compareTotals::fs, 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().

                                                    {
  edm::Service<TFileService> fs;
  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", "std::vector<int>", &t_gjetN);
  tree_->Branch("t_gjetE", "std::vector<double>", &t_gjetE);
  tree_->Branch("t_gjetPt", "std::vector<double>", &t_gjetPt);
  tree_->Branch("t_gjetEta", "std::vector<double>", &t_gjetEta);
  tree_->Branch("t_gjetPhi", "std::vector<double>", &t_gjetPhi);

  tree_->Branch("t_jetTrkP", "std::vector<vector<double> >", &t_jetTrkP);
  tree_->Branch("t_jetTrkPt", "std::vector<vector<double> >", &t_jetTrkPt);
  tree_->Branch("t_jetTrkEta", "std::vector<vector<double> >", &t_jetTrkEta);
  tree_->Branch("t_jetTrkPhi", "std::vector<vector<double> >", &t_jetTrkPhi);
  tree_->Branch("t_jetTrkPdg", "std::vector<vector<double> >", &t_jetTrkPdg);
  tree_->Branch("t_jetTrkCharge", "std::vector<vector<double> >", &t_jetTrkCharge);
}

clearTreeVectors()

void IsolatedParticlesGeneratedJets::clearTreeVectors ( )

private

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

endJob()

void IsolatedParticlesGeneratedJets::endJob ( void  )

inlineoverrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 59 of file IsolatedParticlesGeneratedJets.cc.

Referenced by o2olib.O2ORunMgr::executeJob().

{}

fillDescriptions()

void IsolatedParticlesGeneratedJets::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

static

Definition at line 87 of file IsolatedParticlesGeneratedJets.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and ProducerED_cfi::InputTag.

                                                                                                {
  edm::ParameterSetDescription desc;
  desc.addUntracked<bool>("Debug", true);
  desc.add<edm::InputTag>("JetSource", edm::InputTag("ak5GenJets"));
  desc.add<edm::InputTag>("ParticleSource", edm::InputTag("genParticles"));
  descriptions.add("isolatedParticlesGeneratedJets", desc);
}

Member Data Documentation

debug_

const bool IsolatedParticlesGeneratedJets::debug_

private

Definition at line 64 of file IsolatedParticlesGeneratedJets.cc.

Referenced by analyze().

t_gjetE

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

private

Definition at line 71 of file IsolatedParticlesGeneratedJets.cc.

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

t_gjetEta

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

private

Definition at line 71 of file IsolatedParticlesGeneratedJets.cc.

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

t_gjetN

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

private

Definition at line 70 of file IsolatedParticlesGeneratedJets.cc.

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

t_gjetPhi

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

private

Definition at line 71 of file IsolatedParticlesGeneratedJets.cc.

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

t_gjetPt

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

private

Definition at line 71 of file IsolatedParticlesGeneratedJets.cc.

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

t_jetTrkCharge

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

private

Definition at line 77 of file IsolatedParticlesGeneratedJets.cc.

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

t_jetTrkEta

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

private

Definition at line 74 of file IsolatedParticlesGeneratedJets.cc.

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

t_jetTrkP

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

private

Definition at line 72 of file IsolatedParticlesGeneratedJets.cc.

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

t_jetTrkPdg

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

private

Definition at line 76 of file IsolatedParticlesGeneratedJets.cc.

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

t_jetTrkPhi

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

private

Definition at line 75 of file IsolatedParticlesGeneratedJets.cc.

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

t_jetTrkPt

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

private

Definition at line 73 of file IsolatedParticlesGeneratedJets.cc.

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

tok_jets_

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

private

Definition at line 67 of file IsolatedParticlesGeneratedJets.cc.

Referenced by analyze().

tok_parts_

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

private

Definition at line 68 of file IsolatedParticlesGeneratedJets.cc.

Referenced by analyze().

tree_

TTree* IsolatedParticlesGeneratedJets::tree_

private

Definition at line 65 of file IsolatedParticlesGeneratedJets.cc.

Referenced by analyze(), and bookHistograms().