Inheritance diagram for PhotonProducer:

Public Member Functions
	PhotonProducer (const edm::ParameterSet &ps)

void	produce (edm::Event &evt, const edm::EventSetup &es) 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

Private Member Functions
void	fillPhotonCollection (edm::Event &evt, edm::EventSetup const &es, const edm::Handle< reco::PhotonCoreCollection > &photonCoreHandle, const CaloTopology topology, const HcalPFCuts hcalCuts, const EcalRecHitCollection ecalBarrelHits, const EcalRecHitCollection ecalEndcapHits, ElectronHcalHelper const &hcalHelperCone, ElectronHcalHelper const &hcalHelperBc, reco::VertexCollection &pvVertices, reco::PhotonCollection &outputCollection, int &iSC)

Private Attributes
edm::EDGetTokenT< EcalRecHitCollection >	barrelEcalHits_

const edm::ESGetToken< CaloGeometry, CaloGeometryRecord >	caloGeomToken_

std::string	candidateP4type_

bool	cutsFromDB_

edm::EDGetTokenT< EcalRecHitCollection >	endcapEcalHits_

std::vector< int >	flagsexclEB_

std::vector< int >	flagsexclEE_

edm::EDGetTokenT< HBHERecHitCollection >	hbheRecHits_

edm::ESGetToken< HcalPFCuts, HcalPFCutsRcd >	hcalCutsToken_

std::unique_ptr< ElectronHcalHelper >	hcalHelperBc_

std::unique_ptr< ElectronHcalHelper >	hcalHelperCone_

bool	hcalRun2EffDepth_

double	highEt_

double	hOverEConeSize_

double	minR9Barrel_

double	minR9Endcap_

std::string	PhotonCollection_

edm::EDGetTokenT< reco::PhotonCoreCollection >	photonCoreProducer_

PhotonEnergyCorrector	photonEnergyCorrector_

PhotonIsolationCalculator	photonIsolationCalculator_

PhotonMIPHaloTagger	photonMIPHaloTagger_

std::unique_ptr< PhotonMVABasedHaloTagger >	photonMVABasedHaloTagger_ = nullptr

PositionCalc	posCalculator_

std::vector< double >	preselCutValuesBarrel_

std::vector< double >	preselCutValuesEndcap_

bool	runMIPTagger_

bool	runMVABasedHaloTagger_

std::vector< int >	severitiesexclEB_

std::vector< int >	severitiesexclEE_

const edm::ESGetToken< CaloTopology, CaloTopologyRecord >	topologyToken_

bool	usePrimaryVertex_

bool	validConversions_

bool	validPixelSeeds_

edm::EDGetTokenT< reco::VertexCollection >	vertexProducer_

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

Author: Nancy Marinelli, U. of Notre Dame, US

Definition at line 48 of file PhotonProducer.cc.

Constructor & Destructor Documentation

◆ PhotonProducer()

PhotonProducer::PhotonProducer ( const edm::ParameterSet & ps )

Definition at line 125 of file PhotonProducer.cc.

References barrelEcalHits_, candidateP4type_, cutsFromDB_, endcapEcalHits_, flagsexclEB_, flagsexclEE_, hbheRecHits_, hcalCutsToken_, hcalHelperBc_, hcalHelperCone_, hcalRun2EffDepth_, highEt_, ElectronHcalHelper::Configuration::hOverEConeSize, hOverEConeSize_, gedPhotons_cfi::isolationSumsCalculatorSet, minR9Barrel_, minR9Endcap_, gedPhotons_cfi::mipVariableSet, PhotonCollection_, photonCoreProducer_, photonIsolationCalculator_, photonMIPHaloTagger_, ecalPiZeroTask_cfi::posCalcParameters, posCalculator_, preselCutValuesBarrel_, preselCutValuesEndcap_, runMIPTagger_, runMVABasedHaloTagger_, PhotonMIPHaloTagger::setup(), PhotonIsolationCalculator::setup(), severitiesexclEB_, severitiesexclEE_, AlCaHLTBitMon_QueryRunRegistry::string, usePrimaryVertex_, and vertexProducer_.

     : caloGeomToken_(esConsumes()), topologyToken_(esConsumes()), photonEnergyCorrector_(config, consumesCollector()) {
   // use configuration file to setup input/output collection names
 
   photonCoreProducer_ = consumes<reco::PhotonCoreCollection>(config.getParameter<edm::InputTag>("photonCoreProducer"));
   barrelEcalHits_ = consumes<EcalRecHitCollection>(config.getParameter<edm::InputTag>("barrelEcalHits"));
   endcapEcalHits_ = consumes<EcalRecHitCollection>(config.getParameter<edm::InputTag>("endcapEcalHits"));
   vertexProducer_ = consumes<reco::VertexCollection>(config.getParameter<edm::InputTag>("primaryVertexProducer"));
   hbheRecHits_ = consumes<HBHERecHitCollection>(config.getParameter<edm::InputTag>("hbheRecHits"));
   hOverEConeSize_ = config.getParameter<double>("hOverEConeSize");
   highEt_ = config.getParameter<double>("highEt");
   // R9 value to decide converted/unconverted
   minR9Barrel_ = config.getParameter<double>("minR9Barrel");
   minR9Endcap_ = config.getParameter<double>("minR9Endcap");
   usePrimaryVertex_ = config.getParameter<bool>("usePrimaryVertex");
   runMIPTagger_ = config.getParameter<bool>("runMIPTagger");
   runMVABasedHaloTagger_ = config.getParameter<bool>("runMVABasedHaloTagger");
 
   candidateP4type_ = config.getParameter<std::string>("candidateP4type");
 
   edm::ParameterSet posCalcParameters = config.getParameter<edm::ParameterSet>("posCalcParameters");
   posCalculator_ = PositionCalc(posCalcParameters);
 
   //Retrieve HCAL PF thresholds - from config or from DB
   cutsFromDB_ = config.getParameter<bool>("usePFThresholdsFromDB");
   if (cutsFromDB_) {
     hcalCutsToken_ = esConsumes<HcalPFCuts, HcalPFCutsRcd>(edm::ESInputTag("", "withTopo"));
   }
 
   //AA
   //Flags and Severities to be excluded from photon calculations
   const std::vector<std::string> flagnamesEB =
       config.getParameter<std::vector<std::string>>("RecHitFlagToBeExcludedEB");
 
   const std::vector<std::string> flagnamesEE =
       config.getParameter<std::vector<std::string>>("RecHitFlagToBeExcludedEE");
 
   flagsexclEB_ = StringToEnumValue<EcalRecHit::Flags>(flagnamesEB);
 
   flagsexclEE_ = StringToEnumValue<EcalRecHit::Flags>(flagnamesEE);
 
   const std::vector<std::string> severitynamesEB =
       config.getParameter<std::vector<std::string>>("RecHitSeverityToBeExcludedEB");
 
   severitiesexclEB_ = StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEB);
 
   const std::vector<std::string> severitynamesEE =
       config.getParameter<std::vector<std::string>>("RecHitSeverityToBeExcludedEE");
 
   severitiesexclEE_ = StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEE);
 
   ElectronHcalHelper::Configuration cfgCone, cfgBc;
   cfgCone.hOverEConeSize = hOverEConeSize_;
   if (cfgCone.hOverEConeSize > 0) {
     cfgCone.onlyBehindCluster = false;
     cfgCone.checkHcalStatus = false;
 
     cfgCone.hbheRecHits = hbheRecHits_;
 
     cfgCone.eThresHB = config.getParameter<EgammaHcalIsolation::arrayHB>("recHitEThresholdHB");
     cfgCone.maxSeverityHB = config.getParameter<int>("maxHcalRecHitSeverity");
     cfgCone.eThresHE = config.getParameter<EgammaHcalIsolation::arrayHE>("recHitEThresholdHE");
     cfgCone.maxSeverityHE = cfgCone.maxSeverityHB;
   }
 
   cfgBc.hOverEConeSize = 0.;
   cfgBc.onlyBehindCluster = true;
   cfgBc.checkHcalStatus = false;
 
   cfgBc.hbheRecHits = hbheRecHits_;
 
   cfgBc.eThresHB = config.getParameter<EgammaHcalIsolation::arrayHB>("recHitEThresholdHB");
   cfgBc.maxSeverityHB = config.getParameter<int>("maxHcalRecHitSeverity");
   cfgBc.eThresHE = config.getParameter<EgammaHcalIsolation::arrayHE>("recHitEThresholdHE");
   cfgBc.maxSeverityHE = cfgBc.maxSeverityHB;
 
   hcalHelperCone_ = std::make_unique<ElectronHcalHelper>(cfgCone, consumesCollector());
   hcalHelperBc_ = std::make_unique<ElectronHcalHelper>(cfgBc, consumesCollector());
 
   hcalRun2EffDepth_ = config.getParameter<bool>("hcalRun2EffDepth");
 
   //AA
 
   //
 
   // Parameters for the position calculation:
   //  std::map<std::string,double> providedParameters;
   // providedParameters.insert(std::make_pair("LogWeighted",config.getParameter<bool>("posCalc_logweight")));
   //providedParameters.insert(std::make_pair("T0_barl",config.getParameter<double>("posCalc_t0_barl")));
   //providedParameters.insert(std::make_pair("T0_endc",config.getParameter<double>("posCalc_t0_endc")));
   //providedParameters.insert(std::make_pair("T0_endcPresh",config.getParameter<double>("posCalc_t0_endcPresh")));
   //providedParameters.insert(std::make_pair("W0",config.getParameter<double>("posCalc_w0")));
   //providedParameters.insert(std::make_pair("X0",config.getParameter<double>("posCalc_x0")));
   //posCalculator_ = PositionCalc(providedParameters);
   // cut values for pre-selection
   preselCutValuesBarrel_.push_back(config.getParameter<double>("minSCEtBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("maxHoverEBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("ecalRecHitSumEtOffsetBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("ecalRecHitSumEtSlopeBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("hcalRecHitSumEtOffsetBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("hcalRecHitSumEtSlopeBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("nTrackSolidConeBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("nTrackHollowConeBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("trackPtSumSolidConeBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("trackPtSumHollowConeBarrel"));
   preselCutValuesBarrel_.push_back(config.getParameter<double>("sigmaIetaIetaCutBarrel"));
   //
   preselCutValuesEndcap_.push_back(config.getParameter<double>("minSCEtEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("maxHoverEEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("ecalRecHitSumEtOffsetEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("ecalRecHitSumEtSlopeEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("hcalRecHitSumEtOffsetEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("hcalRecHitSumEtSlopeEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("nTrackSolidConeEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("nTrackHollowConeEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("trackPtSumSolidConeEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("trackPtSumHollowConeEndcap"));
   preselCutValuesEndcap_.push_back(config.getParameter<double>("sigmaIetaIetaCutEndcap"));
   //
 
   edm::ParameterSet isolationSumsCalculatorSet = config.getParameter<edm::ParameterSet>("isolationSumsCalculatorSet");
   photonIsolationCalculator_.setup(
       isolationSumsCalculatorSet, flagsexclEB_, flagsexclEE_, severitiesexclEB_, severitiesexclEE_, consumesCollector());
 
   edm::ParameterSet mipVariableSet = config.getParameter<edm::ParameterSet>("mipVariableSet");
   photonMIPHaloTagger_.setup(mipVariableSet, consumesCollector());
 
   // Register the product
   produces<reco::PhotonCollection>(PhotonCollection_);
 }

Member Function Documentation

◆ fillPhotonCollection()

void PhotonProducer::fillPhotonCollection	(	edm::Event &	evt,
		edm::EventSetup const &	es,
		const edm::Handle< reco::PhotonCoreCollection > &	photonCoreHandle,
		const CaloTopology *	topology,
		const HcalPFCuts *	hcalCuts,
		const EcalRecHitCollection *	ecalBarrelHits,
		const EcalRecHitCollection *	ecalEndcapHits,
		ElectronHcalHelper const &	hcalHelperCone,
		ElectronHcalHelper const &	hcalHelperBc,
		reco::VertexCollection &	pvVertices,
		reco::PhotonCollection &	outputCollection,
		int &	iSC
	)

private

fill shower shape block

fill full5x5 shower shape block

get ecal photon specific corrected energy plus values from regressions and store them in the Photon
Pre-selection loose isolation cuts

Definition at line 348 of file PhotonProducer.cc.

References PhotonEnergyCorrector::calculate(), PhotonIsolationCalculator::calculate(), PositionCalc::Calculate_Location(), caloGeomToken_, candidateP4type_, EcalClusterToolsT< noZS >::covariances(), EgHLTOffHistBins_cfi::e1x5, EcalClusterToolsT< noZS >::e1x5(), reco::Photon::ShowerShape::e1x5, EgHLTOffHistBins_cfi::e2x5, reco::Photon::ShowerShape::e2x5, EcalClusterToolsT< noZS >::e2x5Max(), EcalClusterToolsT< noZS >::e3x3(), reco::Photon::ShowerShape::e3x3, EcalClusterToolsT< noZS >::e5x5(), reco::Photon::ShowerShape::e5x5, electrons_cff::e5x5, DetId::Ecal, reco::Photon::ecal_photons, EcalBarrel, EcalEndcap, EcalPreshower, cosmicPhotonAnalyzer_cfi::eMax, EcalClusterToolsT< noZS >::eMax(), flagsexclEB_, flagsexclEE_, HLT_2024v12_cff::full5x5_sigmaIetaIeta, relativeConstraints::geometry, edm::EventSetup::getData(), ElectronHcalHelper::hcalESum(), reco::Photon::ShowerShape::hcalOverEcal, reco::Photon::ShowerShape::hcalOverEcalBc, ElectronHcalHelper::hcalTowersBehindClusters(), reco::Photon::ShowerShape::hcalTowersBehindClusters, highEt_, hfClusterShapes_cfi::hits, l1ctLayer2EG_cff::id, PhotonEnergyCorrector::init(), createfilelist::int, EcalClusterToolsT< noZS >::localCovariances(), reco::Photon::ShowerShape::maxEnergyXtal, EgHLTOffEleSelection_cfi::minR9, minR9Barrel_, minR9Endcap_, PhotonMIPHaloTagger::MIPcalculate(), gedPhotons_cfi::outputPhotonCollection, MaterialEffects_cfi::photonEnergy, photonEnergyCorrector_, photonIsolationCalculator_, photonMIPHaloTagger_, posCalculator_, reco::Photon::ShowerShape::pre7DepthHcal, preselCutValuesBarrel_, preselCutValuesEndcap_, electrons_cff::r9, reco::Photon::regression1, runMIPTagger_, severitiesexclEB_, severitiesexclEE_, reco::Photon::ShowerShape::sigmaEtaEta, HLT_2024v12_cff::sigmaEtaEta, run3scouting_cff::sigmaIetaIeta, reco::Photon::ShowerShape::sigmaIetaIeta, mathSSE::sqrt(), parallelization::uint, spclusmultinvestigator_cfi::vertexCollection, and L1BJetProducer_cff::vtx.

Referenced by produce().

                                                     {
   // get the geometry from the event setup:
   const CaloGeometry* geometry = &es.getData(caloGeomToken_);
   const CaloSubdetectorGeometry* subDetGeometry = nullptr;
   const CaloSubdetectorGeometry* geometryES = geometry->getSubdetectorGeometry(DetId::Ecal, EcalPreshower);
   const EcalRecHitCollection* hits = nullptr;
   std::vector<double> preselCutValues;
   float minR9 = 0;
 
   photonEnergyCorrector_.init(es);
 
   std::vector<int> flags_, severitiesexcl_;
 
   for (unsigned int lSC = 0; lSC < photonCoreHandle->size(); lSC++) {
     reco::PhotonCoreRef coreRef(reco::PhotonCoreRef(photonCoreHandle, lSC));
     reco::SuperClusterRef scRef = coreRef->superCluster();
     iSC++;
 
     int subdet = scRef->seed()->hitsAndFractions()[0].first.subdetId();
     subDetGeometry = geometry->getSubdetectorGeometry(DetId::Ecal, subdet);
 
     if (subdet == EcalBarrel) {
       preselCutValues = preselCutValuesBarrel_;
       minR9 = minR9Barrel_;
       hits = ecalBarrelHits;
       flags_ = flagsexclEB_;
       severitiesexcl_ = severitiesexclEB_;
     } else if (subdet == EcalEndcap) {
       preselCutValues = preselCutValuesEndcap_;
       minR9 = minR9Endcap_;
       hits = ecalEndcapHits;
       flags_ = flagsexclEE_;
       severitiesexcl_ = severitiesexclEE_;
     } else {
       edm::LogWarning("") << "PhotonProducer: do not know if it is a barrel or endcap SuperCluster";
     }
     if (hits == nullptr)
       continue;
 
     // SC energy preselection
     if (scRef->energy() / cosh(scRef->eta()) <= preselCutValues[0])
       continue;
 
     // recalculate position of seed BasicCluster taking shower depth for unconverted photon
     math::XYZPoint unconvPos =
         posCalculator_.Calculate_Location(scRef->seed()->hitsAndFractions(), hits, subDetGeometry, geometryES);
 
     float maxXtal = EcalClusterTools::eMax(*(scRef->seed()), &(*hits));
     //AA
     //Change these to consider severity level of hits
     float e1x5 = EcalClusterTools::e1x5(*(scRef->seed()), &(*hits), &(*topology));
     float e2x5 = EcalClusterTools::e2x5Max(*(scRef->seed()), &(*hits), &(*topology));
     float e3x3 = EcalClusterTools::e3x3(*(scRef->seed()), &(*hits), &(*topology));
     float e5x5 = EcalClusterTools::e5x5(*(scRef->seed()), &(*hits), &(*topology));
     const auto& cov = EcalClusterTools::covariances(*(scRef->seed()), &(*hits), &(*topology), geometry);
     const auto& locCov = EcalClusterTools::localCovariances(*(scRef->seed()), &(*hits), &(*topology));
 
     float sigmaEtaEta = sqrt(cov[0]);
     float sigmaIetaIeta = sqrt(locCov[0]);
     float r9 = e3x3 / (scRef->rawEnergy());
 
     float full5x5_maxXtal = noZS::EcalClusterTools::eMax(*(scRef->seed()), &(*hits));
     //AA
     //Change these to consider severity level of hits
     float full5x5_e1x5 = noZS::EcalClusterTools::e1x5(*(scRef->seed()), &(*hits), &(*topology));
     float full5x5_e2x5 = noZS::EcalClusterTools::e2x5Max(*(scRef->seed()), &(*hits), &(*topology));
     float full5x5_e3x3 = noZS::EcalClusterTools::e3x3(*(scRef->seed()), &(*hits), &(*topology));
     float full5x5_e5x5 = noZS::EcalClusterTools::e5x5(*(scRef->seed()), &(*hits), &(*topology));
     const auto& full5x5_cov = noZS::EcalClusterTools::covariances(*(scRef->seed()), &(*hits), &(*topology), geometry);
     const auto& full5x5_locCov = noZS::EcalClusterTools::localCovariances(*(scRef->seed()), &(*hits), &(*topology));
 
     float full5x5_sigmaEtaEta = sqrt(full5x5_cov[0]);
     float full5x5_sigmaIetaIeta = sqrt(full5x5_locCov[0]);
 
     // compute position of ECAL shower
     math::XYZPoint caloPosition;
     if (r9 > minR9) {
       caloPosition = unconvPos;
     } else {
       caloPosition = scRef->position();
     }
 
     double photonEnergy = 1.;
     math::XYZPoint vtx(0., 0., 0.);
     if (!vertexCollection.empty())
       vtx = vertexCollection.begin()->position();
     // compute momentum vector of photon from primary vertex and cluster position
     math::XYZVector direction = caloPosition - vtx;
     math::XYZVector momentum = direction.unit();
 
     // Create dummy candidate with unit momentum and zero energy to allow setting of all variables. The energy is set for last.
     math::XYZTLorentzVectorD p4(momentum.x(), momentum.y(), momentum.z(), photonEnergy);
     reco::Photon newCandidate(p4, caloPosition, coreRef, vtx);
 
     // Calculate fiducial flags and isolation variable. Blocked are filled from the isolationCalculator
     reco::Photon::FiducialFlags fiducialFlags;
     reco::Photon::IsolationVariables isolVarR03, isolVarR04;
     photonIsolationCalculator_.calculate(&newCandidate, evt, es, fiducialFlags, isolVarR04, isolVarR03, hcalCuts);
     newCandidate.setFiducialVolumeFlags(fiducialFlags);
     newCandidate.setIsolationVariables(isolVarR04, isolVarR03);
 
     reco::Photon::ShowerShape showerShape;
     showerShape.e1x5 = e1x5;
     showerShape.e2x5 = e2x5;
     showerShape.e3x3 = e3x3;
     showerShape.e5x5 = e5x5;
     showerShape.maxEnergyXtal = maxXtal;
     showerShape.sigmaEtaEta = sigmaEtaEta;
     showerShape.sigmaIetaIeta = sigmaIetaIeta;
     for (uint id = 0; id < showerShape.hcalOverEcal.size(); ++id) {
       showerShape.hcalOverEcal[id] = hcalHelperCone.hcalESum(*scRef, id + 1, hcalCuts) / scRef->energy();
       showerShape.hcalOverEcalBc[id] = hcalHelperBc.hcalESum(*scRef, id + 1, hcalCuts) / scRef->energy();
     }
     showerShape.hcalTowersBehindClusters = hcalHelperBc.hcalTowersBehindClusters(*scRef);
     showerShape.pre7DepthHcal = false;
     newCandidate.setShowerShapeVariables(showerShape);
 
     reco::Photon::ShowerShape full5x5_showerShape;
     full5x5_showerShape.e1x5 = full5x5_e1x5;
     full5x5_showerShape.e2x5 = full5x5_e2x5;
     full5x5_showerShape.e3x3 = full5x5_e3x3;
     full5x5_showerShape.e5x5 = full5x5_e5x5;
     full5x5_showerShape.maxEnergyXtal = full5x5_maxXtal;
     full5x5_showerShape.sigmaEtaEta = full5x5_sigmaEtaEta;
     full5x5_showerShape.sigmaIetaIeta = full5x5_sigmaIetaIeta;
     for (uint id = 0; id < full5x5_showerShape.hcalOverEcal.size(); ++id) {
       full5x5_showerShape.hcalOverEcal[id] = hcalHelperCone.hcalESum(*scRef, id + 1, hcalCuts) / full5x5_e5x5;
       full5x5_showerShape.hcalOverEcalBc[id] = hcalHelperBc.hcalESum(*scRef, id + 1, hcalCuts) / full5x5_e5x5;
     }
     full5x5_showerShape.hcalTowersBehindClusters = hcalHelperBc.hcalTowersBehindClusters(*scRef);
     full5x5_showerShape.pre7DepthHcal = false;
     newCandidate.full5x5_setShowerShapeVariables(full5x5_showerShape);
 
     // Photon candidate takes by default (set in photons_cfi.py)  a 4-momentum derived from the ecal photon-specific corrections.
     photonEnergyCorrector_.calculate(evt, newCandidate, subdet, vertexCollection, es);
     if (candidateP4type_ == "fromEcalEnergy") {
       newCandidate.setP4(newCandidate.p4(reco::Photon::ecal_photons));
       newCandidate.setCandidateP4type(reco::Photon::ecal_photons);
     } else if (candidateP4type_ == "fromRegression") {
       newCandidate.setP4(newCandidate.p4(reco::Photon::regression1));
       newCandidate.setCandidateP4type(reco::Photon::regression1);
     }
 
     // fill MIP Vairables for Halo: Block for MIP are filled from PhotonMIPHaloTagger
     reco::Photon::MIPVariables mipVar;
     if (subdet == EcalBarrel && runMIPTagger_) {
       photonMIPHaloTagger_.MIPcalculate(&newCandidate, evt, es, mipVar);
       newCandidate.setMIPVariables(mipVar);
     }
 
     bool isLooseEM = true;
     if (newCandidate.pt() < highEt_) {
       if (newCandidate.hadronicOverEm() >= preselCutValues[1])
         isLooseEM = false;
       if (newCandidate.ecalRecHitSumEtConeDR04() > preselCutValues[2] + preselCutValues[3] * newCandidate.pt())
         isLooseEM = false;
       if (newCandidate.hcalTowerSumEtConeDR04() > preselCutValues[4] + preselCutValues[5] * newCandidate.pt())
         isLooseEM = false;
       if (newCandidate.nTrkSolidConeDR04() > int(preselCutValues[6]))
         isLooseEM = false;
       if (newCandidate.nTrkHollowConeDR04() > int(preselCutValues[7]))
         isLooseEM = false;
       if (newCandidate.trkSumPtSolidConeDR04() > preselCutValues[8])
         isLooseEM = false;
       if (newCandidate.trkSumPtHollowConeDR04() > preselCutValues[9])
         isLooseEM = false;
       if (newCandidate.sigmaIetaIeta() > preselCutValues[10])
         isLooseEM = false;
     }
 
     if (isLooseEM)
       outputPhotonCollection.push_back(newCandidate);
   }
 }

◆ produce()

void PhotonProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 256 of file PhotonProducer.cc.

References barrelEcalHits_, HLT_2024v12_cff::barrelRecHits, cutsFromDB_, endcapEcalHits_, HLT_2024v12_cff::endcapRecHits, fillPhotonCollection(), edm::Event::getByToken(), edm::EventSetup::getData(), hcalCutsToken_, hcalHelperBc_, hcalHelperCone_, hcalRun2EffDepth_, edm::HandleBase::isValid(), eostools::move(), gedPhotons_cfi::outputPhotonCollection, PhotonCollection_, photonCoreProducer_, edm::Handle< T >::product(), edm::Event::put(), HLT_2024v12_cff::topology, topologyToken_, usePrimaryVertex_, spclusmultinvestigator_cfi::vertexCollection, and vertexProducer_.

                                                                                    {
   HcalPFCuts const* hcalCuts = nullptr;
   if (cutsFromDB_) {
     hcalCuts = &theEventSetup.getData(hcalCutsToken_);
   }
   using namespace edm;
   //  nEvt_++;
 
   reco::PhotonCollection outputPhotonCollection;
   auto outputPhotonCollection_p = std::make_unique<reco::PhotonCollection>();
 
   // Get the PhotonCore collection
   bool validPhotonCoreHandle = true;
   Handle<reco::PhotonCoreCollection> photonCoreHandle;
   theEvent.getByToken(photonCoreProducer_, photonCoreHandle);
   if (!photonCoreHandle.isValid()) {
     edm::LogError("PhotonProducer") << "Error! Can't get the photonCoreProducer";
     validPhotonCoreHandle = false;
   }
 
   // Get EcalRecHits
   bool validEcalRecHits = true;
   Handle<EcalRecHitCollection> barrelHitHandle;
   EcalRecHitCollection barrelRecHits;
   theEvent.getByToken(barrelEcalHits_, barrelHitHandle);
   if (!barrelHitHandle.isValid()) {
     edm::LogError("PhotonProducer") << "Error! Can't get the barrelEcalHits";
     validEcalRecHits = false;
   }
   if (validEcalRecHits)
     barrelRecHits = *(barrelHitHandle.product());
 
   Handle<EcalRecHitCollection> endcapHitHandle;
   theEvent.getByToken(endcapEcalHits_, endcapHitHandle);
   EcalRecHitCollection endcapRecHits;
   if (!endcapHitHandle.isValid()) {
     edm::LogError("PhotonProducer") << "Error! Can't get the endcapEcalHits";
     validEcalRecHits = false;
   }
   if (validEcalRecHits)
     endcapRecHits = *(endcapHitHandle.product());
 
   const CaloTopology* topology = &theEventSetup.getData(topologyToken_);
 
   // prepare access to hcal data
   hcalHelperCone_->beginEvent(theEvent, theEventSetup);
   hcalHelperBc_->beginEvent(theEvent, theEventSetup);
 
   // Get the primary event vertex
   Handle<reco::VertexCollection> vertexHandle;
   reco::VertexCollection vertexCollection;
   bool validVertex = true;
   if (usePrimaryVertex_) {
     theEvent.getByToken(vertexProducer_, vertexHandle);
     if (!vertexHandle.isValid()) {
       edm::LogWarning("PhotonProducer") << "Error! Can't get the product primary Vertex Collection "
                                         << "\n";
       validVertex = false;
     }
     if (validVertex)
       vertexCollection = *(vertexHandle.product());
   }
 
   int iSC = 0;  // index in photon collection
   // Loop over barrel and endcap SC collections and fill the  photon collection
   if (validPhotonCoreHandle)
     fillPhotonCollection(theEvent,
                          theEventSetup,
                          photonCoreHandle,
                          topology,
                          hcalCuts,
                          &barrelRecHits,
                          &endcapRecHits,
                          *hcalHelperCone_,
                          *hcalHelperBc_,
                          vertexCollection,
                          outputPhotonCollection,
                          iSC);
 
   // put the product in the event
   edm::LogInfo("PhotonProducer") << " Put in the event " << iSC << " Photon Candidates \n";
   outputPhotonCollection_p->assign(outputPhotonCollection.begin(), outputPhotonCollection.end());
 
   // go back to run2-like 2 effective depths if desired - depth 1 is the normal depth 1, depth 2 is the sum over the rest
   if (hcalRun2EffDepth_) {
     for (auto& pho : *outputPhotonCollection_p)
       pho.hcalToRun2EffDepth();
   }
 
   theEvent.put(std::move(outputPhotonCollection_p), PhotonCollection_);
 }