EnergyRingsTableProducer Class Reference

Inheritance diagram for EnergyRingsTableProducer:

Public Member Functions
	EnergyRingsTableProducer (const edm::ParameterSet &iConfig)
	~EnergyRingsTableProducer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Private Member Functions
void	produce (edm::Event &, edm::EventSetup const &) override

Private Attributes
std::string	name_
edm::EDGetTokenT< edm::View< pat::Jet > >	srcJet_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 22 of file EnergyRingsTableProducer.cc.

Constructor & Destructor Documentation

EnergyRingsTableProducer()

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

inlineexplicit

Definition at line 24 of file EnergyRingsTableProducer.cc.

      : name_(iConfig.getParameter<std::string>("name")),
        srcJet_(consumes<edm::View<pat::Jet>>(iConfig.getParameter<edm::InputTag>("src"))) {
    produces<nanoaod::FlatTable>();
  }

~EnergyRingsTableProducer()

EnergyRingsTableProducer::~EnergyRingsTableProducer ( )

inlineoverride

Definition at line 30 of file EnergyRingsTableProducer.cc.

{};

Member Function Documentation

fillDescriptions()

static void EnergyRingsTableProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions )

inlinestatic

Definition at line 32 of file EnergyRingsTableProducer.cc.

                                                                           {
    edm::ParameterSetDescription desc;
    desc.add<edm::InputTag>("src")->setComment("input Jet collection");
    desc.add<std::string>("name")->setComment("name of the Jet FlatTable we are extending with energy rings");
    descriptions.add("EnergyRingsTable", desc);
  }

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, and AlCaHLTBitMon_QueryRunRegistry::string.

produce()

void EnergyRingsTableProducer::produce ( edm::Event &  iEvent, edm::EventSetup const &  iSetup  )

overrideprivate

Definition at line 47 of file EnergyRingsTableProducer.cc.

                                                                                    {
  edm::Handle<edm::View<pat::Jet>> srcJet;
  iEvent.getByToken(srcJet_, srcJet);
 
  float cone_boundaries[] = {0.05, 0.1, 0.2, 0.3, 0.4};
  size_t ncone_boundaries = sizeof(cone_boundaries) / sizeof(float);
  unsigned int nJet = srcJet->size();
  unsigned int ncand = 0;
  std::vector<float> numdaughterspt03;
  std::vector<std::vector<float>> EmFractionEnergyRings(ncone_boundaries + 1, std::vector<float>(nJet, 0.));
  std::vector<std::vector<float>> ChFractionEnergyRings(ncone_boundaries + 1, std::vector<float>(nJet, 0.));
  std::vector<std::vector<float>> NeFractionEnergyRings(ncone_boundaries + 1, std::vector<float>(nJet, 0.));
  std::vector<std::vector<float>> MuFractionEnergyRings(ncone_boundaries + 1, std::vector<float>(nJet, 0.));
 
  for (unsigned int ij = 0; ij < nJet; ij++) {
    ++ncand;
    auto jet = srcJet->ptrAt(ij);
    int numDaughtersPt03 = 0;
    for (unsigned int ijcone = 0; ijcone < ncone_boundaries; ijcone++) {
      EmFractionEnergyRings[ijcone][ij] = 0;
      MuFractionEnergyRings[ijcone][ij] = 0;
      ChFractionEnergyRings[ijcone][ij] = 0;
      NeFractionEnergyRings[ijcone][ij] = 0;
    }
    for (const auto& d : jet->daughterPtrVector()) {
      float candDr = Geom::deltaR(d->p4(), jet->p4());
      size_t icone =
          std::lower_bound(&cone_boundaries[0], &cone_boundaries[ncone_boundaries], candDr) - &cone_boundaries[0];
      float candEnergy = d->energy();
      int pdgid = abs(d->pdgId());
      if (pdgid == 22 || pdgid == 11) {
        EmFractionEnergyRings[icone][ij] += candEnergy;
      } else if (pdgid == 13) {
        MuFractionEnergyRings[icone][ij] += candEnergy;
      } else if (d->charge() != 0) {
        ChFractionEnergyRings[icone][ij] += candEnergy;
      } else {
        NeFractionEnergyRings[icone][ij] += candEnergy;
      }
      if (d->pt() > 0.3)
        numDaughtersPt03 += 1;
    }  // end of jet daughters loop
    numdaughterspt03.push_back(numDaughtersPt03);
  }                                                                            //end of jet loop
  auto tab = std::make_unique<nanoaod::FlatTable>(ncand, name_, false, true);  //extension to Jet collection set to true
  tab->addColumn<int>("numDaughtersPt03", numdaughterspt03, "number of jet daughters with pT>0.3 GeV");
 
  tab->addColumn<float>("EmFractionEnergyRing0", EmFractionEnergyRings[0], "Em energy fraction in ring in dR 0-0.05");
  tab->addColumn<float>("EmFractionEnergyRing1", EmFractionEnergyRings[1], "Em energy fraction in ring in dR 0.05-0.1");
  tab->addColumn<float>("EmFractionEnergyRing2", EmFractionEnergyRings[2], "Em energy fraction in ring in dR 0.1-0.2");
  tab->addColumn<float>("EmFractionEnergyRing3", EmFractionEnergyRings[3], "Em energy fraction in ring in dR 0.2-0.3");
  tab->addColumn<float>("EmFractionEnergyRing4", EmFractionEnergyRings[4], "Em energy fraction in ring in dR 0.3-0.4");
  tab->addColumn<float>(
      "EmFractionEnergyRing5", EmFractionEnergyRings[5], "Em energy fraction in ring in dR 0.4 overflow");
 
  tab->addColumn<float>("ChFractionEnergyRing0", ChFractionEnergyRings[0], "Ch energy fraction in ring in dR 0-0.05");
  tab->addColumn<float>("ChFractionEnergyRing1", ChFractionEnergyRings[1], "Ch energy fraction in ring in dR 0.05-0.1");
  tab->addColumn<float>("ChFractionEnergyRing2", ChFractionEnergyRings[2], "Ch energy fraction in ring in dR 0.1-0.2");
  tab->addColumn<float>("ChFractionEnergyRing3", ChFractionEnergyRings[3], "Ch energy fraction in ring in dR 0.2-0.3");
  tab->addColumn<float>("ChFractionEnergyRing4", ChFractionEnergyRings[4], "Ch energy fraction in ring in dR 0.3-0.4");
  tab->addColumn<float>(
      "ChFractionEnergyRing5", ChFractionEnergyRings[5], "Ch energy fraction in ring in dR 0.4 overflow");
 
  tab->addColumn<float>("MuFractionEnergyRing0", MuFractionEnergyRings[0], "Mu energy fraction in ring in dR 0-0.05");
  tab->addColumn<float>("MuFractionEnergyRing1", MuFractionEnergyRings[1], "Mu energy fraction in ring in dR 0.05-0.1");
  tab->addColumn<float>("MuFractionEnergyRing2", MuFractionEnergyRings[2], "Mu energy fraction in ring in dR 0.1-0.2");
  tab->addColumn<float>("MuFractionEnergyRing3", MuFractionEnergyRings[3], "Mu energy fraction in ring in dR 0.2-0.3");
  tab->addColumn<float>("MuFractionEnergyRing4", MuFractionEnergyRings[4], "Mu energy fraction in ring in dR 0.3-0.4");
  tab->addColumn<float>(
      "MuFractionEnergyRing5", MuFractionEnergyRings[5], "Mu energy fraction in ring in dR 0.4 overflow");
 
  tab->addColumn<float>("NeFractionEnergyRing0", NeFractionEnergyRings[0], "Ne energy fraction in ring in dR 0-0.05");
  tab->addColumn<float>("NeFractionEnergyRing1", NeFractionEnergyRings[1], "Ne energy fraction in ring in dR 0.05-0.1");
  tab->addColumn<float>("NeFractionEnergyRing2", NeFractionEnergyRings[2], "Ne energy fraction in ring in dR 0.1-0.2");
  tab->addColumn<float>("NeFractionEnergyRing3", NeFractionEnergyRings[3], "Ne energy fraction in ring in dR 0.2-0.3");
  tab->addColumn<float>("NeFractionEnergyRing4", NeFractionEnergyRings[4], "Ne energy fraction in ring in dR 0.3-0.4");
  tab->addColumn<float>(
      "NeFractionEnergyRing5", NeFractionEnergyRings[5], "Ne energy fraction in ring in dR 0.4 overflow");
 
  iEvent.put(std::move(tab));
}

References funct::abs(), ztail::d, PbPb_ZMuSkimMuonDPG_cff::deltaR, iEvent, metsig::jet, pfDeepBoostedJetPreprocessParams_cfi::lower_bound, eostools::move(), name_, EgammaValidation_cff::pdgid, electrons_cff::srcJet, and srcJet_.

Member Data Documentation

name_

std::string EnergyRingsTableProducer::name_

private

Definition at line 42 of file EnergyRingsTableProducer.cc.

Referenced by produce().

srcJet_

edm::EDGetTokenT<edm::View<pat::Jet> > EnergyRingsTableProducer::srcJet_

private

Definition at line 43 of file EnergyRingsTableProducer.cc.

Referenced by produce().