26 double eta1 = p.
eta();
27 double phi1 = p.
phi();
29 std::vector<double> Rxy;
31 for(reco::TrackCollection::const_iterator
32 recTrack = recCollection->begin(); recTrack!= recCollection->end(); recTrack++)
34 double pt = recTrack->pt();
35 double eta2 = recTrack->eta();
36 double phi2 = recTrack->phi();
41 double dR = dRDistance(eta1,phi1,eta2,phi2);
46 sort(Rxy.begin(), Rxy.end());
50 if(dR < Rxy[Rxy.size()-1])
52 Rxy[Rxy.size()-1] =
dR;
53 sort(Rxy.begin(), Rxy.end());
58 if (Rxy.size()<
y) rxy=4;
else rxy = Rxy[(int)y-1];
dRxyCalculator(const edm::Event &iEvent, const edm::EventSetup &iSetup, edm::InputTag trackLabel)
virtual double eta() const
momentum pseudorapidity
double getDRxy(const reco::Photon cluster, double x, double y)
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
virtual double phi() const
momentum azimuthal angle