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
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

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<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 22 of file EnergyRingsTableProducer.cc.

Constructor & Destructor Documentation

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

inlineoverride

Definition at line 30 of file EnergyRingsTableProducer.cc.

{};

Member Function Documentation

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

inlinestatic

Definition at line 32 of file EnergyRingsTableProducer.cc.

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

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

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

overrideprivate

Definition at line 47 of file EnergyRingsTableProducer.cc.

References funct::abs(), ztail::d, HLT_FULL_cff::deltaR, edm::Event::getByToken(), metsig::jet, cuda_std::lower_bound(), eostools::move(), name_, edm::Event::put(), and srcJet_.

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

Member Data Documentation

std::string EnergyRingsTableProducer::name_

private

Definition at line 42 of file EnergyRingsTableProducer.cc.

Referenced by produce().

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

private

Definition at line 43 of file EnergyRingsTableProducer.cc.

Referenced by produce().