27 for(reco::TrackCollection::const_iterator
28 recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
30 double pt = recTrack->pt();
45 double eta1 = p.
eta();
46 double phi1 = p.
phi();
48 for(reco::TrackCollection::const_iterator
49 recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
51 double pt = recTrack->pt();
52 if ( (pt > ptCut) && ( fabs(eta1 - recTrack->eta()) < detaCut) && ( fabs(calcDphi(recTrack->phi(),phi1)) < 3.141592/2. ) )
70 double eta1 = p.
eta();
71 double phi1 = p.
phi();
75 for(reco::TrackCollection::const_iterator
76 recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
78 double pt = recTrack->pt();
79 double eta2 = recTrack->eta();
80 double phi2 = recTrack->phi();
82 if(dRDistance(eta1,phi1,eta2,phi2) >= 0.1 * x)
97 double eta1 = p.
eta();
98 double phi1 = p.
phi();
102 for(reco::TrackCollection::const_iterator
103 recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
105 double pt = recTrack->pt();
106 double eta2 = recTrack->eta();
107 double phi2 = recTrack->phi();
108 if (dRDistance(eta1,phi1,eta2,phi2) >= outerR)
110 if (dRDistance(eta1,phi1,eta2,phi2) <= innerR)
double getTxy(const reco::Photon &p, double x, double y)
virtual float phi() const
momentum azimuthal angle
double getHollSxy(const reco::Photon &p, double thePtCut, double outerR, double innerR)
virtual float eta() const
momentum pseudorapidity
int getNumAllTracks(double ptCut)
TxyCalculator(const edm::Event &iEvent, const edm::EventSetup &iSetup, const edm::InputTag &trackLabel)
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
int getNumLocalTracks(const reco::Photon &p, double detaCut, double ptCut)