56 std::unique_ptr<reco::Track>
ret;
58 int nhits = hits.size();
59 if(nhits <2)
return ret;
66 for (
auto const & ih : hits)
70 points[
i] = recHit->globalPosition();
71 errors[
i] = recHit->globalPositionError();
72 isBarrel[++
i] = recHit->detUnit()->type().isBarrel();
84 float valPt = (invPt > 1.e-4) ? 1./invPt : 1.e4;
85 float errPt = 0.055*valPt + 0.017*valPt*valPt;
92 float valZip =
getZip(valTip, curvature, points[0],points[1]);
98 valPhi,valTip,valZip);
101 float errTip, errZip;
103 float errPhi = 0.002;
104 float errCotTheta = 0.002;
110 chi2 = rzLine.chi2();
tuple ret
prodAgent to be discontinued
bool isBarrel(GeomDetEnumerators::SubDetector m)
Global3DPoint GlobalPoint
void getErrTipAndErrZip(float pt, float eta, float &errZip, float &errTip) const
T mag() const
The vector magnitude. Equivalent to sqrt(vec.mag2())
T inversePt(T curvature, const edm::EventSetup &iSetup)
float getZip(float d0, float curv, const GlobalPoint &inner, const GlobalPoint &outer) const
T curvature(T InversePt, const edm::EventSetup &iSetup)
OutputIterator zip(InputIterator1 first1, InputIterator1 last1, InputIterator2 first2, InputIterator2 last2, OutputIterator result, Compare comp)
virtual RecHitPointer build(const TrackingRecHit *p) const =0
build a tracking rechit from an existing rechit
reco::Track * build(const Measurement1D &pt, const Measurement1D &phi, const Measurement1D &cotTheta, const Measurement1D &tip, const Measurement1D &zip, float chi2, int charge, const std::vector< const TrackingRecHit * > &hits, const MagneticField *mf, const GlobalPoint &reference=GlobalPoint(0, 0, 0)) const
float getCotThetaAndUpdateZip(const GlobalPoint &inner, const GlobalPoint &outer, float radius, float phi, float d0, float &zip) const
int getCharge(int phi1, int phi2, int phi3, int phi4, int mode)
T y() const
Cartesian y coordinate.
const edm::EventSetup * theES
#define declareDynArray(T, n, x)
float getPhi(float xC, float yC, int charge) const
const TransientTrackingRecHitBuilder * theTTRecHitBuilder
T x() const
Cartesian x coordinate.
const MagneticField * theField