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 noexcept final
bool	wantsGlobalRuns () const noexcept final
bool	wantsInputProcessBlocks () const noexcept final
bool	wantsProcessBlocks () const noexcept 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 noexcept
bool	wantsStreamRuns () const noexcept
	~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().