#include <RecoEgamma/EgammaTools/plugins/HGCalElectronIDValueMapProducer.cc>

Inheritance diagram for HGCalElectronIDValueMapProducer:

Public Member Functions
	HGCalElectronIDValueMapProducer (const edm::ParameterSet &)

	~HGCalElectronIDValueMapProducer () 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	beginStream (edm::StreamID) override

void	endStream () override

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

Private Attributes
std::unique_ptr< HGCalEgammaIDHelper >	eIDHelper_

edm::EDGetTokenT< edm::View< reco::GsfElectron > >	electronsToken_

std::map< const std::string, std::vector< float > >	maps_

float	radius_

Static Private Attributes
static const std::vector< std::string >	valuesProduced_

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

Description: [one line class summary]

Implementation: [Notes on implementation]

Definition at line 41 of file HGCalElectronIDValueMapProducer.cc.

Constructor & Destructor Documentation

◆ HGCalElectronIDValueMapProducer()

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

explicit

Definition at line 82 of file HGCalElectronIDValueMapProducer.cc.

References eIDHelper_, crabWrapper::key, maps_, and valuesProduced_.

     : electronsToken_(consumes(iConfig.getParameter<edm::InputTag>("electrons"))),
       radius_(iConfig.getParameter<double>("pcaRadius")) {
   for (const auto& key : valuesProduced_) {
     maps_[key] = {};
     produces<edm::ValueMap<float>>(key);
   }
 
   eIDHelper_ = std::make_unique<HGCalEgammaIDHelper>(iConfig, consumesCollector());
 }

◆ ~HGCalElectronIDValueMapProducer()

HGCalElectronIDValueMapProducer::~HGCalElectronIDValueMapProducer ( )

override

Definition at line 93 of file HGCalElectronIDValueMapProducer.cc.

93 {}

Member Function Documentation

◆ beginStream()

void HGCalElectronIDValueMapProducer::beginStream ( edm::StreamID )

overrideprivate

Definition at line 214 of file HGCalElectronIDValueMapProducer.cc.

214 {}

◆ endStream()

void HGCalElectronIDValueMapProducer::endStream ( )

overrideprivate

Definition at line 217 of file HGCalElectronIDValueMapProducer.cc.

217 {}

◆ fillDescriptions()

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

static

Definition at line 220 of file HGCalElectronIDValueMapProducer.cc.

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

                                                                                                  {
   // Auto-generate hgcalElectronIDValueMap_cfi
   edm::ParameterSetDescription desc;
   desc.add<edm::InputTag>("electrons", edm::InputTag("ecalDrivenGsfElectronsHGC"));
   desc.add<double>("pcaRadius", 3.0);
   desc.add<std::vector<std::string>>("variables", valuesProduced_);
   desc.add<std::vector<double>>("dEdXWeights")
       ->setComment("This must be copied from dEdX_weights in RecoLocalCalo.HGCalRecProducers.HGCalRecHit_cfi");
   desc.add<unsigned int>("isoNRings", 5);
   desc.add<double>("isoDeltaR", 0.15);
   desc.add<double>("isoDeltaRmin", 0.0);
   desc.add<edm::InputTag>("EERecHits", edm::InputTag("HGCalRecHit", "HGCEERecHits"));
   desc.add<edm::InputTag>("FHRecHits", edm::InputTag("HGCalRecHit", "HGCHEFRecHits"));
   desc.add<edm::InputTag>("BHRecHits", edm::InputTag("HGCalRecHit", "HGCHEBRecHits"));
   desc.add<edm::InputTag>("hitMapTag", edm::InputTag("hgcalRecHitMapProducer"));
   descriptions.add("hgcalElectronIDValueMap", desc);
 }

◆ produce()

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

overrideprivate

Definition at line 96 of file HGCalElectronIDValueMapProducer.cc.

References eIDHelper_, HPSPFTauProducerPuppi_cfi::electron, electronsToken_, Exception, trigObjTnPSource_cfi::filler, iEvent, maps_, eostools::move(), MillePedeFileConverter_cfg::out, radius_, and valuesProduced_.

                                                                                            {
   using namespace edm;
 
   auto electronsH = iEvent.getHandle(electronsToken_);
 
   // Clear previous map
   for (auto&& kv : maps_) {
     kv.second.clear();
     kv.second.reserve(electronsH->size());
   }
 
   // Set up helper tool
   eIDHelper_->eventInit(iEvent, iSetup);
 
   for (const auto& electron : *electronsH) {
     if (electron.isEB()) {
       // Fill some dummy value
       for (auto&& kv : maps_) {
         kv.second.push_back(0.);
       }
     } else {
       eIDHelper_->computeHGCAL(electron, radius_);
 
       // check the PCA has worked out
       if (eIDHelper_->sigmaUU() == -1) {
         for (auto&& kv : maps_) {
           kv.second.push_back(0.);
         }
         continue;
       }
 
       hgcal::LongDeps ld(eIDHelper_->energyPerLayer(radius_, true));
       float measuredDepth, expectedDepth, expectedSigma;
       float depthCompatibility = eIDHelper_->clusterDepthCompatibility(ld, measuredDepth, expectedDepth, expectedSigma);
 
       // Fill here all the ValueMaps from their appropriate functions
 
       // Energies / PT
       const auto* eleCluster = electron.electronCluster().get();
       const double sinTheta = eleCluster->position().rho() / eleCluster->position().r();
       maps_["ecOrigEt"].push_back(eleCluster->energy() * sinTheta);
       maps_["ecOrigEnergy"].push_back(eleCluster->energy());
 
       // energies calculated in an cylinder around the axis of the electron cluster
       float ec_tot_energy = ld.energyEE() + ld.energyFH() + ld.energyBH();
       maps_["ecEt"].push_back(ec_tot_energy * sinTheta);
       maps_["ecEnergy"].push_back(ec_tot_energy);
       maps_["ecEnergyEE"].push_back(ld.energyEE());
       maps_["ecEnergyFH"].push_back(ld.energyFH());
       maps_["ecEnergyBH"].push_back(ld.energyBH());
 
       // Cluster shapes
       // PCA related
       maps_["pcaEig1"].push_back(eIDHelper_->eigenValues()(0));
       maps_["pcaEig2"].push_back(eIDHelper_->eigenValues()(1));
       maps_["pcaEig3"].push_back(eIDHelper_->eigenValues()(2));
       maps_["pcaSig1"].push_back(eIDHelper_->sigmas()(0));
       maps_["pcaSig2"].push_back(eIDHelper_->sigmas()(1));
       maps_["pcaSig3"].push_back(eIDHelper_->sigmas()(2));
       maps_["pcaAxisX"].push_back(eIDHelper_->axis().x());
       maps_["pcaAxisY"].push_back(eIDHelper_->axis().y());
       maps_["pcaAxisZ"].push_back(eIDHelper_->axis().z());
       maps_["pcaPositionX"].push_back(eIDHelper_->barycenter().x());
       maps_["pcaPositionY"].push_back(eIDHelper_->barycenter().y());
       maps_["pcaPositionZ"].push_back(eIDHelper_->barycenter().z());
 
       // transverse shapes
       maps_["sigmaUU"].push_back(eIDHelper_->sigmaUU());
       maps_["sigmaVV"].push_back(eIDHelper_->sigmaVV());
       maps_["sigmaEE"].push_back(eIDHelper_->sigmaEE());
       maps_["sigmaPP"].push_back(eIDHelper_->sigmaPP());
 
       // long profile
       maps_["nLayers"].push_back(ld.nLayers());
       maps_["firstLayer"].push_back(ld.firstLayer());
       maps_["lastLayer"].push_back(ld.lastLayer());
       maps_["e4oEtot"].push_back(ld.e4oEtot());
       maps_["layerEfrac10"].push_back(ld.layerEfrac10());
       maps_["layerEfrac90"].push_back(ld.layerEfrac90());
 
       // depth
       maps_["measuredDepth"].push_back(measuredDepth);
       maps_["expectedDepth"].push_back(expectedDepth);
       maps_["expectedSigma"].push_back(expectedSigma);
       maps_["depthCompatibility"].push_back(depthCompatibility);
 
       // Isolation
       maps_["caloIsoRing0"].push_back(eIDHelper_->getIsolationRing(0));
       maps_["caloIsoRing1"].push_back(eIDHelper_->getIsolationRing(1));
       maps_["caloIsoRing2"].push_back(eIDHelper_->getIsolationRing(2));
       maps_["caloIsoRing3"].push_back(eIDHelper_->getIsolationRing(3));
       maps_["caloIsoRing4"].push_back(eIDHelper_->getIsolationRing(4));
     }
   }
 
   // Check we didn't make up a new variable and forget it in valuesProduced_
   if (maps_.size() != valuesProduced_.size()) {
     throw cms::Exception("HGCalElectronIDValueMapProducer")
         << "We have a miscoded value map producer, since map size changed";
   }
 
   for (auto&& kv : maps_) {
     // Check we didn't forget any values
     if (kv.second.size() != electronsH->size()) {
       throw cms::Exception("HGCalElectronIDValueMapProducer")
           << "We have a miscoded value map producer, since the variable " << kv.first << " wasn't filled.";
     }
     // Do the filling
     auto out = std::make_unique<edm::ValueMap<float>>();
     edm::ValueMap<float>::Filler filler(*out);
     filler.insert(electronsH, kv.second.begin(), kv.second.end());
     filler.fill();
     // and put it into the event
     iEvent.put(std::move(out), kv.first);
   }
 }

Member Data Documentation

◆ eIDHelper_

std::unique_ptr<HGCalEgammaIDHelper> HGCalElectronIDValueMapProducer::eIDHelper_

private

Definition at line 60 of file HGCalElectronIDValueMapProducer.cc.

Referenced by HGCalElectronIDValueMapProducer(), and produce().

◆ electronsToken_

edm::EDGetTokenT<edm::View<reco::GsfElectron> > HGCalElectronIDValueMapProducer::electronsToken_

private

Definition at line 54 of file HGCalElectronIDValueMapProducer.cc.

Referenced by produce().

◆ maps_

std::map<const std::string, std::vector<float> > HGCalElectronIDValueMapProducer::maps_

private

Definition at line 58 of file HGCalElectronIDValueMapProducer.cc.

Referenced by HGCalElectronIDValueMapProducer(), and produce().

◆ radius_

float HGCalElectronIDValueMapProducer::radius_

private

Definition at line 55 of file HGCalElectronIDValueMapProducer.cc.

Referenced by produce().

◆ valuesProduced_

const std::vector< std::string > HGCalElectronIDValueMapProducer::valuesProduced_

staticprivate

Initial value:

= {
    "ecOrigEt",      "ecOrigEnergy",  "ecEt",
    "ecEnergy",      "ecEnergyEE",    "ecEnergyFH",
    "ecEnergyBH",    "pcaEig1",       "pcaEig2",
    "pcaEig3",       "pcaSig1",       "pcaSig2",
    "pcaSig3",       "pcaAxisX",      "pcaAxisY",
    "pcaAxisZ",      "pcaPositionX",  "pcaPositionY",
    "pcaPositionZ",  "sigmaUU",       "sigmaVV",
    "sigmaEE",       "sigmaPP",       "nLayers",
    "firstLayer",    "lastLayer",     "e4oEtot",
    "layerEfrac10",  "layerEfrac90",  "measuredDepth",
    "expectedDepth", "expectedSigma", "depthCompatibility",
    "caloIsoRing0",  "caloIsoRing1",  "caloIsoRing2",
    "caloIsoRing3",  "caloIsoRing4",
}

Definition at line 57 of file HGCalElectronIDValueMapProducer.cc.

Referenced by fillDescriptions(), HGCalElectronIDValueMapProducer(), and produce().

Public Member Functions

Static Public Member Functions

Private Member Functions

Private Attributes

Static Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HGCalElectronIDValueMapProducer()

◆ ~HGCalElectronIDValueMapProducer()

Member Function Documentation

◆ beginStream()

◆ endStream()

◆ fillDescriptions()

◆ produce()

Member Data Documentation

◆ eIDHelper_

◆ electronsToken_

◆ maps_

◆ radius_

◆ valuesProduced_