144 std::vector<reco::Candidate>::const_iterator
cand;
151 std::vector<reco::CandidatePtr> puppiCandsNoLepPtrs;
158 for (
auto pup : *puppiCandsNoLep){
159 puppiCandsNoLepPtrs.push_back(pup.sourceCandidatePtr(0));
162 std::vector<int> mappingPuppi(
usePuppi_ ? puppiCands->size() : 0);
178 std::set<unsigned int> whiteList;
179 std::set<reco::TrackRef> whiteListTk;
184 for(
unsigned int i=0;
i<svWhiteList.
size();
i++) {
186 for(
unsigned int j=0; j< svWhiteList[
i].numberOfSourceCandidatePtrs(); j++) {
188 if(c.
id() == cands.
id()) whiteList.insert(c.
key());
192 for(
auto dau = svWhiteList[
i].
begin(); dau != svWhiteList[
i].
end() ; dau++){
194 if(chCand!=
nullptr) {
195 whiteListTk.insert(chCand->
track());
211 auto outPtrP = std::make_unique<std::vector<pat::PackedCandidate>>();
212 std::vector<int>
mapping(cands->size());
213 std::vector<int> mappingReverse(cands->size());
214 std::vector<int> mappingTk(TKOrigs->size(), -1);
216 for(
unsigned int ic=0, nc = cands->size(); ic < nc; ++ic) {
227 for(
size_t ii=0;
ii<PVs->size();
ii++){
229 if(dz<dist) {pvi=
ii;dist=
dz; }
242 float ptTrk = ctrack->
pt();
243 float etaAtVtx = ctrack->
eta();
244 float phiAtVtx = ctrack->
phi();
263 outPtrP->back().setLostInnerHits( lostHits );
266 whiteList.find(ic)!=whiteList.end() ||
270 if(
abs(outPtrP->back().pdgId())==22) {
281 if(outPtrP->back().pt() > 0.5 ){
290 outPtrP->back().setMuonID(cand.
muonRef()->isStandAloneMuon(), cand.
muonRef()->isGlobalMuon());
296 PVpos = PV->position();
305 bool isIsolatedChargedHadron =
false;
309 outPtrP->back().setIsIsolatedChargedHadron(isIsolatedChargedHadron);
314 }
else if(isIsolatedChargedHadron) {
318 outPtrP->back().setHcalFraction(0);
326 float puppiWeightVal = (*puppiWeight)[pkref];
327 float puppiWeightNoLepVal = 0.0;
331 if ( puppiWeightNoLep.
isValid() ) {
333 auto pkrefPtr = pkref->sourceCandidatePtr(0);
335 bool foundNoLep =
false;
336 for (
size_t ipcnl = 0; ipcnl < puppiCandsNoLepPtrs.size(); ipcnl++){
337 if (puppiCandsNoLepPtrs[ipcnl] == pkrefPtr){
339 puppiWeightNoLepVal = puppiCandsNoLep->at(ipcnl).pt()/cand.
pt();
343 if ( !foundNoLep || puppiWeightNoLepVal > 1 ) {
344 puppiWeightNoLepVal = 1.0;
347 outPtrP->back().setPuppiWeight( puppiWeightVal, puppiWeightNoLepVal );
349 mappingPuppi[((*puppiCandsMap)[pkref]).
key()]=ic;
359 auto outPtrPSorted = std::make_unique<std::vector<pat::PackedCandidate>>();
361 std::vector<size_t> reverseOrder(order.size());
362 for(
size_t i=0,nc=cands->size();
i<nc;
i++) {
363 outPtrPSorted->push_back((*outPtrP)[order[
i]]);
364 reverseOrder[order[
i]] =
i;
365 mappingReverse[order[
i]]=
i;
369 for(
size_t i=0,ntk=mappingTk.size();
i<ntk;
i++){
370 if(mappingTk[
i] >= 0)
371 mappingTk[
i]=reverseOrder[mappingTk[
i]];
374 for(
size_t i=0,ntk=mappingPuppi.size();
i<ntk;
i++){
375 mappingPuppi[
i]=reverseOrder[mappingPuppi[
i]];
381 auto pf2pc = std::make_unique<edm::Association<pat::PackedCandidateCollection>>(oh);
382 auto pc2pf = std::make_unique<edm::Association<reco::PFCandidateCollection>>(cands);
385 pf2pcFiller.insert(cands, mappingReverse.begin(), mappingReverse.end());
386 pc2pfFiller.insert(oh , order.begin(), order.end());
388 pf2pcFiller.insert(TKOrigs, mappingTk.begin(), mappingTk.end());
389 if(
usePuppi_) pf2pcFiller.insert(puppiCands, mappingPuppi.begin(), mappingPuppi.end());
virtual double pt() const final
transverse momentum
double ecalEnergy() const
return corrected Ecal energy
const Point & referencePoint() const
Reference point on the track.
const edm::EDGetTokenT< edm::ValueMap< float > > PuppiWeight_
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
bool isNonnull() const
Checks for non-null.
double rawEcalEnergy() const
return corrected Ecal energy
bool hasValidHitInPixelLayer(enum PixelSubdetector::SubDetector, uint16_t layer) const
virtual double eta() const final
momentum pseudorapidity
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_
key_type key() const
Accessor for product key.
ProductID id() const
Accessor for product ID.
const edm::EDGetTokenT< edm::ValueMap< int > > PVAssoQuality_
reco::TrackRef trackRef() const
virtual double phi() const final
momentum azimuthal angle
const double minPtForTrackProperties_
double eta() const
pseudorapidity of momentum vector
double pt() const
track transverse momentum
virtual int pdgId() const final
PDG identifier.
virtual double energy() const final
energy
Abs< T >::type abs(const T &t)
const bool storeChargedHadronIsolation_
virtual const Point & vertex() const
vertex position (overwritten by PF...)
const edm::EDGetTokenT< reco::TrackCollection > TKOrigs_
const edm::EDGetTokenT< std::vector< reco::PFCandidate > > PuppiCands_
std::vector< size_t > sort_indexes(const std::vector< T > &v) const
const edm::EDGetTokenT< edm::ValueMap< reco::CandidatePtr > > PuppiCandsMap_
const edm::EDGetTokenT< edm::ValueMap< float > > PuppiWeightNoLep_
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 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_
const edm::EDGetTokenT< std::vector< reco::PFCandidate > > PuppiCandsNoLep_
XYZPointD XYZPoint
point in space with cartesian internal representation
virtual reco::TrackRef track() const
reference to a track
ProductID id() const
Accessor for product ID.
const edm::EDGetTokenT< reco::VertexCollection > PVOrigs_
int numberOfLostHits(HitCategory category) const
Particle reconstructed by the particle flow algorithm.
reco::GsfTrackRef gsfTrackRef() const
std::vector< edm::EDGetTokenT< edm::View< reco::Candidate > > > SVWhiteLists_
virtual const PolarLorentzVector & polarP4() const final
four-momentum Lorentz vector
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_
uint16_t getHitPattern(HitCategory category, int position) const
const double minPtForChargedHadronProperties_
double rawHcalEnergy() const
return raw Hcal energy
Public Member Functions inherited from edm::global::EDProducer<>
1.8.11
