47 inline double ptFast(
const double energy,
50 const auto v = position - origin;
65 if ( reconstructionStep_ ==
"final" ) {
86 if (not pfEg.label().empty())
88 consumes<reco::PFCandidateCollection>(pfEg);
99 if (not hcTow.label().empty())
101 consumes<CaloTowerCollection>(hcTow);
123 const std::vector<std::string> flagnamesEB =
124 config.
getParameter<std::vector<std::string> >(
"RecHitFlagToBeExcludedEB");
126 const std::vector<std::string> flagnamesEE =
127 config.
getParameter<std::vector<std::string> >(
"RecHitFlagToBeExcludedEE");
130 StringToEnumValue<EcalRecHit::Flags>(flagnamesEB);
133 StringToEnumValue<EcalRecHit::Flags>(flagnamesEE);
135 const std::vector<std::string> severitynamesEB =
136 config.
getParameter<std::vector<std::string> >(
"RecHitSeverityToBeExcludedEB");
139 StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEB);
141 const std::vector<std::string> severitynamesEE =
142 config.
getParameter<std::vector<std::string> >(
"RecHitSeverityToBeExcludedEE");
145 StringToEnumValue<EcalSeverityLevel::SeverityLevel>(severitynamesEE);
151 thePhotonEnergyCorrector_->gedRegression()->setConsumes(sumes);
194 if ( reconstructionStep_ !=
"final"){
242 auto outputPhotonCollection_p = std::make_unique<reco::PhotonCollection>();
247 bool validPhotonCoreHandle=
false;
249 bool validPhotonHandle=
false;
263 validPhotonHandle=
true;
270 if (photonCoreHandle.
isValid()) {
271 validPhotonCoreHandle=
true;
279 bool validEcalRecHits=
true;
283 if (!barrelHitHandle.
isValid()) {
285 <<
"Error! Can't get the barrelEcalHits";
292 if (!endcapHitHandle.
isValid()) {
294 <<
"Error! Can't get the endcapEcalHits";
298 bool validPreshowerRecHits=
true;
302 if (!preshowerHitHandle.
isValid()) {
304 <<
"Error! Can't get the preshowerEcalHits";
306 if( validPreshowerRecHits ) preshowerRecHits = *(preshowerHitHandle.
product());
314 if (!pfEGCandidateHandle.
isValid()) {
316 <<
"Error! Can't get the pfEgammaCandidates";
325 if (!pfCandidateHandle.
isValid()) {
327 <<
"Error! Can't get the pfCandidates";
358 bool validVertex=
true;
363 <<
"Error! Can't get the product primary Vertex Collection";
381 if ( validPhotonCoreHandle)
392 outputPhotonCollection,
404 outputPhotonCollection,
406 phoChargedIsolationMap_CITK,
407 phoNeutralHadronIsolationMap_CITK,
408 phoPhotonIsolationMap_CITK);
413 edm::LogInfo(
"GEDPhotonProducer") <<
" Put in the event " << iSC <<
" Photon Candidates \n";
414 outputPhotonCollection_p->assign(outputPhotonCollection.begin(),outputPhotonCollection.end());
420 auto pfEGCandToPhotonMap_p = std::make_unique<edm::ValueMap<reco::PhotonRef>>();
422 unsigned nObj = pfEGCandidateHandle->size();
423 std::vector<reco::PhotonRef>
values(nObj);
425 for(
unsigned int lCand=0; lCand < nObj; lCand++) {
429 for(
unsigned int lSC=0; lSC < photonOrphHandle->size(); lSC++) {
432 if ( pfScRef != scRef )
continue;
433 values[lCand] = photonRef;
438 filler.insert(pfEGCandidateHandle,values.begin(),values.end());
466 std::vector<double> preselCutValues;
467 std::vector<int> flags_, severitiesexcl_;
469 for(
unsigned int lSC=0; lSC < photonCoreHandle->size(); lSC++) {
480 int thedet = scRef->seed()->hitsAndFractions()[0].first.det();
481 int subdet = scRef->seed()->hitsAndFractions()[0].first.subdetId();
484 hits = ecalBarrelHits;
489 hits = ecalEndcapHits;
498 edm::LogWarning(
"")<<
"GEDPhotonProducer: do not know if it is a barrel or endcap SuperCluster" << thedet <<
' ' << subdet;
506 ptFast(parentSCRef->energy(),parentSCRef->position(),
math::XYZPoint(0,0,0)) <= preselCutValues[0] )
continue;
512 std::vector<CaloTowerDetId> TowersBehindClus;
513 float hcalDepth1OverEcalBc,hcalDepth2OverEcalBc;
514 hcalDepth1OverEcalBc=hcalDepth2OverEcalBc=0.f;
525 TowersBehindClus = towerIsoBehindClus.
towersOf(*scRef);
526 hcalDepth1OverEcalBc = towerIsoBehindClus.
getDepth1HcalESum(TowersBehindClus)/scRef->energy();
527 hcalDepth2OverEcalBc = towerIsoBehindClus.
getDepth2HcalESum(TowersBehindClus)/scRef->energy();
537 float e2x5 = ( hits !=
nullptr ? EcalClusterTools::e2x5Max( *(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
538 float e3x3 = ( hits !=
nullptr ? EcalClusterTools::e3x3( *(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
540 std::vector<float> cov = ( hits !=
nullptr ? EcalClusterTools::covariances( *(scRef->seed()), &(*hits), &(*topology),
geometry) : std::vector<float>( {0.f,0.f,0.f} ) );
541 std::vector<float> locCov = ( hits !=
nullptr ? EcalClusterTools::localCovariances( *(scRef->seed()), &(*hits), &(*topology)) : std::vector<float>( {0.f,0.f,0.f} ) );
544 float sigmaIetaIeta =
sqrt(locCov[0]);
556 float full5x5_sigmaEtaEta =
sqrt(full5x5_cov[0]);
557 float full5x5_sigmaIetaIeta =
sqrt(full5x5_locCov[0]);
564 double photonEnergy=1.;
566 if (vertexCollection.size()>0) vtx = vertexCollection.begin()->position();
585 newCandidate.setFiducialVolumeFlags( fiducialFlags );
586 newCandidate.setIsolationVariables(isolVarR04, isolVarR03 );
593 showerShape.
e3x3= e3x3;
605 const float sep = locCov[1];
608 showerShape.
e2nd = ( hits !=
nullptr ? EcalClusterTools::e2nd(*(scRef->seed()),&(*hits)) : 0.
f );
609 showerShape.
eTop = ( hits !=
nullptr ? EcalClusterTools::eTop(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
610 showerShape.
eLeft = ( hits !=
nullptr ? EcalClusterTools::eLeft(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
611 showerShape.
eRight = ( hits !=
nullptr ? EcalClusterTools::eRight(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
612 showerShape.
eBottom = ( hits !=
nullptr ? EcalClusterTools::eBottom(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
613 showerShape.
e1x3 = ( hits !=
nullptr ? EcalClusterTools::e1x3(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
614 showerShape.
e2x2 = ( hits !=
nullptr ? EcalClusterTools::e2x2(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
615 showerShape.
e2x5Max = ( hits !=
nullptr ? EcalClusterTools::e2x5Max(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
616 showerShape.
e2x5Left = ( hits !=
nullptr ? EcalClusterTools::e2x5Left(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
617 showerShape.
e2x5Right = ( hits !=
nullptr ? EcalClusterTools::e2x5Right(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
618 showerShape.
e2x5Top = ( hits !=
nullptr ? EcalClusterTools::e2x5Top(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
619 showerShape.
e2x5Bottom = ( hits !=
nullptr ? EcalClusterTools::e2x5Bottom(*(scRef->seed()), &(*hits), &(*topology)) : 0.
f );
622 const float sigmaRR = toolsforES.eseffsirir( *scRef );
624 newCandidate.setShowerShapeVariables ( showerShape );
629 int nSaturatedXtals = 0;
630 bool isSeedSaturated =
false;
631 if (hits !=
nullptr) {
632 const auto hitsAndFractions = scRef->hitsAndFractions();
633 for (
auto&& hitFractionPair : hitsAndFractions) {
636 if (
ecalRecHit->checkFlag(EcalRecHit::Flags::kSaturated)) {
639 isSeedSaturated =
true;
645 newCandidate.setSaturationInfo(saturationInfo);
649 full5x5_showerShape.
e1x5= full5x5_e1x5;
650 full5x5_showerShape.
e2x5= full5x5_e2x5;
651 full5x5_showerShape.
e3x3= full5x5_e3x3;
652 full5x5_showerShape.
e5x5= full5x5_e5x5;
654 full5x5_showerShape.
sigmaEtaEta = full5x5_sigmaEtaEta;
657 const float full5x5_spp = (!
edm::isFinite(full5x5_locCov[2]) ? 0. :
sqrt(full5x5_locCov[2]));
658 const float full5x5_sep = full5x5_locCov[1];
675 newCandidate.full5x5_setShowerShapeVariables ( full5x5_showerShape );
713 newCandidate.setMIPVariables(mipVar);
720 if ( newCandidate.pt() <
highEt_) {
721 if ( newCandidate.hadronicOverEm() >= preselCutValues[1] ) isLooseEM=
false;
722 if ( newCandidate.ecalRecHitSumEtConeDR04() > preselCutValues[2]+ preselCutValues[3]*newCandidate.pt() ) isLooseEM=
false;
723 if ( newCandidate.hcalTowerSumEtConeDR04() > preselCutValues[4]+ preselCutValues[5]*newCandidate.pt() ) isLooseEM=
false;
724 if ( newCandidate.nTrkSolidConeDR04() >
int(preselCutValues[6]) ) isLooseEM=
false;
725 if ( newCandidate.nTrkHollowConeDR04() >
int(preselCutValues[7]) ) isLooseEM=
false;
726 if ( newCandidate.trkSumPtSolidConeDR04() > preselCutValues[8] ) isLooseEM=
false;
727 if ( newCandidate.trkSumPtHollowConeDR04() > preselCutValues[9] ) isLooseEM=
false;
728 if ( newCandidate.sigmaIetaIeta() > preselCutValues[10] ) isLooseEM=
false;
734 outputPhotonCollection.push_back(newCandidate);
754 std::vector<double> preselCutValues;
757 for(
unsigned int lSC=0; lSC < photonHandle->size(); lSC++) {
761 int thedet = scRef->seed()->hitsAndFractions()[0].first.det();
762 int subdet = scRef->seed()->hitsAndFractions()[0].first.subdetId();
770 edm::LogWarning(
"")<<
"GEDPhotonProducer: do not know if it is a barrel or endcap SuperCluster" << thedet <<
' ' << subdet;
777 ptFast(parentSCRef->energy(),parentSCRef->position(),
math::XYZPoint(0,0,0)) <= preselCutValues[0] )
continue;
791 pfIso.
photonIso = (*photons_)[photonPtr];
819 outputPhotonCollection.push_back(newCandidate);
edm::InputTag photonProducer_
void setPflowIsolationVariables(const PflowIsolationVariables &pfisol)
Set Particle Flow Isolation variables.
T getParameter(std::string const &) const
PhotonEnergyCorrector * thePhotonEnergyCorrector_
PhotonMIPHaloTagger * thePhotonMIPHaloTagger_
virtual void beginRun(edm::Run const &r, edm::EventSetup const &es) override final
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
void calculate(const reco::Photon *, const edm::Event &, const edm::EventSetup &es, reco::Photon::FiducialFlags &phofid, reco::Photon::IsolationVariables &phoisolR03, reco::Photon::IsolationVariables &phoisolR04) const
bool isNonnull() const
Checks for non-null.
edm::EDGetTokenT< reco::PFCandidateCollection > pfCandidates_
edm::ESHandle< CaloGeometry > theCaloGeom_
bool existsAs(std::string const ¶meterName, bool trackiness=true) const
checks if a parameter exists as a given type
std::vector< CaloTowerDetId > hcalTowersBehindClusters
ROOT::Math::LorentzVector< ROOT::Math::PxPyPzE4D< double > > XYZTLorentzVectorD
Lorentz vector with cylindrical internal representation using pseudorapidity.
void setCandidateP4type(const P4type type)
edm::EDGetTokenT< CaloTowerCollection > hcalTowers_
isolationSumsCalculatorSet
CaloTopology const * topology(0)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
virtual void produce(edm::Event &evt, const edm::EventSetup &es) override
float hcalDepth2OverEcalBc
edm::EDGetTokenT< edm::ValueMap< float > > phoPhotonIsolationToken_CITK
std::vector< int > flagsexclEB_
std::vector< Vertex > VertexCollection
collection of Vertex objects
edm::EDGetTokenT< EcalRecHitCollection > endcapEcalHits_
std::string reconstructionStep_
PhotonIsolationCalculator * thePhotonIsolationCalculator_
std::vector< int > severitiesexclEE_
std::unique_ptr< ModifyObjectValueBase > & gedRegression()
double ptFast(const double energy, const math::XYZPoint &position, const math::XYZPoint &origin)
std::string photonCollection_
virtual void endRun(edm::Run const &, edm::EventSetup const &) override final
edm::ESHandle< CaloTopology > theCaloTopo_
void setTowerCollection(const CaloTowerCollection *towercollection)
edm::EDGetTokenT< EcalRecHitCollection > barrelEcalHits_
PositionCalc posCalculator_
double getDepth1HcalESum(const reco::SuperCluster &sc) const
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
void setup(const edm::ParameterSet &conf, std::vector< int > const &flagsEB_, std::vector< int > const &flagsEE_, std::vector< int > const &severitiesEB_, std::vector< int > const &severitiesEE_, edm::ConsumesCollector &&iC)
void setPflowIDVariables(const PflowIDVariables &pfid)
TypeLabelItem const & produces()
declare what type of product will make and with which optional label
std::vector< int > flagsexclEE_
edm::EDGetTokenT< reco::PhotonCoreCollection > photonCoreProducerT_
edm::EDGetTokenT< edm::ValueMap< float > > phoChargedIsolationToken_CITK
void init(const edm::EventSetup &theEventSetup)
void setup(const edm::ParameterSet &conf, edm::ConsumesCollector &&iC)
double getTowerESum(const reco::Candidate *cand, const std::vector< CaloTowerDetId > *detIdToExclude=0) const
edm::EDGetTokenT< reco::PFCandidateCollection > pfEgammaCandidates_
std::vector< int > severitiesexclEB_
const_iterator end() const
void calculate(edm::Event &evt, reco::Photon &, int subdet, const reco::VertexCollection &vtxcol, const edm::EventSetup &iSetup)
DetId seed() const
return DetId of seed
GEDPhotonProducer(const edm::ParameterSet &ps)
float hcalDepth1OverEcalBc
T const * product() const
std::vector< double > preselCutValuesBarrel_
edm::EDGetTokenT< reco::PhotonCollection > photonProducerT_
XYZVectorD XYZVector
spatial vector with cartesian internal representation
XYZPointD XYZPoint
point in space with cartesian internal representation
std::vector< Photon > PhotonCollection
collectin of Photon objects
std::vector< double > preselCutValuesEndcap_
const LorentzVector & p4(P4type type) const
std::string candidateP4type_
void fillPhotonCollection(edm::Event &evt, edm::EventSetup const &es, const edm::Handle< reco::PhotonCoreCollection > &photonCoreHandle, const CaloTopology *topology, const EcalRecHitCollection *ecalBarrelHits, const EcalRecHitCollection *ecalEndcapHits, const EcalRecHitCollection *preshowerHits, const edm::Handle< CaloTowerCollection > &hcalTowersHandle, const reco::VertexCollection &pvVertices, reco::PhotonCollection &outputCollection, int &iSC)
std::vector< CaloTowerDetId > towersOf(const reco::SuperCluster &sc) const
edm::EDGetTokenT< edm::ValueMap< float > > phoNeutralHadronIsolationToken_CITK
double getDepth2HcalESum(const reco::SuperCluster &sc) const
ESHandle< TrackerGeometry > geometry
iterator find(key_type k)
static int position[264][3]
bool isUninitialized() const
edm::EDGetTokenT< reco::VertexCollection > vertexProducer_
T const * product() const
edm::EDGetTokenT< EcalRecHitCollection > preshowerHits_
void setP4(P4type type, const LorentzVector &p4, float p4Error, bool setToRecoCandidate)
std::string valueMapPFCandPhoton_