188 std::set<unsigned int> whiteList;
189 std::set<reco::TrackRef> whiteListTk;
194 for (
unsigned int i = 0;
i < svWhiteList.
size();
i++) {
196 for (
unsigned int j = 0;
j < svWhiteList[
i].numberOfSourceCandidatePtrs();
j++) {
198 if (c.
id() == cands.
id())
199 whiteList.insert(c.
key());
202 for (
auto dau = svWhiteList[
i].
begin(); dau != svWhiteList[
i].
end(); dau++) {
204 if (chCand !=
nullptr) {
205 whiteListTk.insert(chCand->
track());
224 std::vector<pat::HcalDepthEnergyFractions> hcalDepthEnergyFractions;
225 hcalDepthEnergyFractions.reserve(cands->size());
226 std::vector<pat::HcalDepthEnergyFractions> hcalDepthEnergyFractions_Ordered;
227 hcalDepthEnergyFractions_Ordered.reserve(cands->size());
231 auto outPtrP = std::make_unique<std::vector<pat::PackedCandidate>>();
232 std::vector<int> mapping(cands->size());
233 std::vector<int> mappingReverse(cands->size());
234 std::vector<int> mappingTk(TKOrigs->size(), -1);
236 for (
unsigned int ic = 0, nc = cands->size(); ic < nc; ++ic) {
247 for (
size_t ii = 0;
ii < PVs->size();
ii++) {
270 float ptTrk = ctrack->
pt();
271 float etaAtVtx = ctrack->
eta();
272 float phiAtVtx = ctrack->
phi();
292 outPtrP->back().setLostInnerHits(lostHits);
294 whiteList.find(ic) != whiteList.end() ||
297 if (
abs(outPtrP->back().pdgId()) == 22) {
301 outPtrP->back().setTrackProperties(*ctrack,
312 outPtrP->back().setTrackProperties(*ctrack,
316 outPtrP->back().setTrackProperties(*ctrack,
326 outPtrP->back().setMuonID(cand.
muonRef()->isStandAloneMuon(), cand.
muonRef()->isGlobalMuon());
331 PVpos = PV->position();
336 outPtrP->back().setAssociationQuality(
342 bool isIsolatedChargedHadron =
false;
345 isIsolatedChargedHadron =
347 outPtrP->back().setIsIsolatedChargedHadron(isIsolatedChargedHadron);
356 outPtrP->back().setHcalFraction(0);
357 outPtrP->back().setCaloFraction(0);
360 if (isIsolatedChargedHadron) {
364 outPtrP->back().setRawCaloFraction(0);
365 outPtrP->back().setRawHcalFraction(0);
368 std::vector<float> dummyVector;
376 fabs(outPtrP->back().eta()) > 1.3) {
381 hcalDepthEFrac.reset(hcalDepthEnergyFractionTmp);
384 hcalDepthEnergyFractions.push_back(hcalDepthEFrac);
390 outPtrP->back().setGoodEgamma();
396 float puppiWeightVal = (*puppiWeight)[pkref];
397 float puppiWeightNoLepVal = (*puppiWeightNoLep)[pkref];
398 outPtrP->back().setPuppiWeight(puppiWeightVal, puppiWeightNoLepVal);
405 auto t0Err = (*t0ErrMap)[cand.
trackRef()];
406 outPtrP->back().setTime(
t0, t0Err);
421 auto outPtrPSorted = std::make_unique<std::vector<pat::PackedCandidate>>();
423 std::vector<size_t> reverseOrder(order.size());
424 for (
size_t i = 0, nc = cands->size();
i < nc;
i++) {
425 outPtrPSorted->push_back((*outPtrP)[order[
i]]);
426 reverseOrder[order[
i]] =
i;
427 mappingReverse[order[
i]] =
i;
428 hcalDepthEnergyFractions_Ordered.push_back(hcalDepthEnergyFractions[order[i]]);
432 for (
size_t i = 0, ntk = mappingTk.size(); i < ntk; i++) {
433 if (mappingTk[i] >= 0)
434 mappingTk[
i] = reverseOrder[mappingTk[
i]];
440 auto pf2pc = std::make_unique<edm::Association<pat::PackedCandidateCollection>>(oh);
441 auto pc2pf = std::make_unique<edm::Association<reco::PFCandidateCollection>>(
cands);
444 pf2pcFiller.insert(cands, mappingReverse.begin(), mappingReverse.end());
445 pc2pfFiller.insert(oh, order.begin(), order.end());
447 pf2pcFiller.insert(TKOrigs, mappingTk.begin(), mappingTk.end());
455 auto hcalDepthEnergyFractionsV = std::make_unique<edm::ValueMap<HcalDepthEnergyFractions>>();
457 fillerHcalDepthEnergyFractions.insert(
458 cands, hcalDepthEnergyFractions_Ordered.begin(), hcalDepthEnergyFractions_Ordered.end());
459 fillerHcalDepthEnergyFractions.fill();
462 iEvent.
put(
std::move(hcalDepthEnergyFractionsV),
"hcalDepthEnergyFractions");
double ecalEnergy() const
return corrected Ecal energy
const Point & referencePoint() const
Reference point on the track.
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
const edm::EventSetup & c
bool isNonnull() const
Checks for non-null.
double rawEcalEnergy() const
return corrected Ecal energy
double pt() const final
transverse momentum
bool hasValidHitInPixelLayer(enum PixelSubdetector::SubDetector, uint16_t layer) const
float mva_nothing_gamma() const
mva for gamma detection
bool getByToken(EDGetToken token, Handle< PROD > &result) const
double phi() const
azimuthal angle of momentum vector
static const int qualityMap[8]
conversion map from quality flags used in PV association and miniAOD one
const edm::EDGetTokenT< reco::PFCandidateCollection > Cands_
string chargedHadronIsolation
key_type key() const
Accessor for product key.
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
ProductID id() const
Accessor for product ID.
const edm::EDGetTokenT< edm::ValueMap< int > > PVAssoQuality_
const edm::EDGetTokenT< edm::ValueMap< float > > PuppiWeightNoLep_
reco::TrackRef trackRef() const
int pdgId() const final
PDG identifier.
const double minPtForTrackProperties_
bool isTimeValid() const
do we have a valid time information
double eta() const
pseudorapidity of momentum vector
const edm::EDGetTokenT< edm::ValueMap< float > > t0ErrMap_
reco::TrackRef track() const override
reference to a track
double pt() const
track transverse momentum
const PolarLorentzVector & polarP4() const final
four-momentum Lorentz vector
Abs< T >::type abs(const T &t)
const bool storeChargedHadronIsolation_
const bool timeFromValueMap_
static std::vector< size_t > sort_indexes(const std::vector< T > &v)
const edm::EDGetTokenT< reco::TrackCollection > TKOrigs_
const bool storeHcalDepthEndcapOnly_
edm::Ref< VertexCollection > VertexRef
persistent reference to a Vertex
const int covarianceVersion_
reco::MuonRef muonRef() const
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
LostInnerHits
Enumerator specifying the.
const std::array< float, 7 > & hcalDepthEnergyFractions() const
fraction of hcal energy at a given depth (index 0..6 for depth 1..7)
const edm::EDGetTokenT< edm::ValueMap< bool > > ChargedHadronIsolation_
edm::Ptr< Candidate > CandidatePtr
persistent reference to an object in a collection of Candidate objects
T const * product() const
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
const std::vector< int > covariancePackingSchemas_
XYZPointD XYZPoint
point in space with cartesian internal representation
ProductID id() const
Accessor for product ID.
const edm::EDGetTokenT< edm::ValueMap< float > > PuppiWeight_
const edm::EDGetTokenT< reco::VertexCollection > PVOrigs_
int numberOfLostHits(HitCategory category) const
const Point & vertex() const override
vertex position (overwritten by PF...)
Particle reconstructed by the particle flow algorithm.
std::vector< edm::EDGetTokenT< edm::View< reco::Candidate > > > SVWhiteLists_
reco::GsfTrackRef gsfTrackRef() const
const std::vector< int > pfCandidateTypesForHcalDepth_
double hcalEnergy() const
return corrected Hcal energy
int numberOfValidPixelHits() const
const_iterator end() const
const edm::EDGetTokenT< edm::Association< reco::VertexCollection > > PVAsso_
const edm::EDGetTokenT< reco::VertexCollection > PVs_
const double minPtForLowQualityTrackProperties_
const edm::EDGetTokenT< edm::ValueMap< float > > t0Map_
virtual ParticleType particleId() const
uint16_t getHitPattern(HitCategory category, int position) const
double phi() const final
momentum azimuthal angle
const double minPtForChargedHadronProperties_
int charge() const final
electric charge
double rawHcalEnergy() const
return raw Hcal energy
double energy() const final
energy
double eta() const final
momentum pseudorapidity
Static Public Member Functions inherited from edm::global::EDProducerBase
1.8.5
