18 TrackAnalysis(
name), m_pTCut(pTCut), m_etaCut(etaCut), m_d0Cut(d0Cut), m_z0Cut(z0Cut) {
46 50, -3.0, 3.0, 500, -30000.0, 30000.0);
48 50, -3.0, 3.0, 500, -30000.0, 30000.0);
50 50, -3.0, 3.0, 500, -300.0, 300.0);
52 50, -3.0, 3.0, 500, -3.0, 3.0);
54 std::string(
m_name+
" resolution vs pixelHits (w b-layer)").c_str(), 50, -3, 3, 10, 0, 10);
56 std::string(
m_name+
" resolution vs pixelHits (w/o b-layer)").c_str(), 50, -3, 3, 10, 0, 10);
66 50, -5.0, 5.0, 500, -3.0, 3.0);
68 50, -5.0, 5.0, 500, -3.0, 3.0);
70 std::string(
m_name+
" d0VsPixelHits pull (without b-layer)").c_str(), 50, -5.0, 5.0, 10, 0.0, 10.0);
72 std::string(
m_name+
" z0VsPixelHits pull (without b-layer)").c_str(), 50, -5.0, 5.0, 10, 0.0, 10.0);
74 std::string(
m_name+
" d0VsPixelHits pull (with b-layer)").c_str(), 50, -5.0, 5.0, 10, 0.0, 10.0);
76 std::string(
m_name+
" z0VsPixelHits pull (with b-layer)").c_str(), 50, -5.0, 5.0, 10, 0.0, 10.0);
94 const std::vector<TIDA::Track*>& ,
98 std::vector<TIDA::Track*>::const_iterator
reference, referenceEnd=referenceTracks.end();
102 double referenceEta = (*reference)->eta();
103 double referencePhi =
phi((*reference)->phi());
104 double referenceZ0 = (*reference)->z0();
105 double referenceD0 = (*reference)->a0();
106 double referencePT = (*reference)->pT();
108 if (fabs(referencePT)<
m_pTCut)
continue;
109 if (fabs(referenceEta)>
m_etaCut)
continue;
110 if (fabs(referenceD0)>
m_d0Cut)
continue;
111 if (fabs(referenceZ0)>
m_z0Cut)
continue;
121 double testEta =
test->eta();
122 double testPhi =
phi(
test->phi());
123 double testZ0 =
test->z0();
124 double testD0 =
test->a0();
125 double testPT =
test->pT();
128 if(referencePT==0 || testPT==0)
continue;
131 double eeta = sqrt(
test->deta()*
test->deta() + (*reference)->deta()*(*reference)->deta() );
132 double ephi = sqrt(
test->dphi()*
test->dphi() + (*reference)->dphi()*(*reference)->dphi() );
133 double ez0 = sqrt(
test->dz0() *
test->dz0() + (*reference)->dz0() *(*reference)->dz0() );
134 double ed0 = sqrt(
test->da0() *
test->da0() + (*reference)->da0() *(*reference)->da0() );
135 double einvpT = sqrt(
test->dpT() *
test->dpT() + (*reference)->dpT() *(*reference)->dpT() );
151 double referenceBLayerHits = (*reference)->bLayerHits();
152 double referencePixelHits = (*reference)->pixelHits();
155 if (referenceBLayerHits>=1) {
160 if (referencePT <= 10000 || testPT <= 10000)
continue;
162 if (referenceBLayerHits>=1) {
167 if (referenceBLayerHits>=1)
173 if (referenceBLayerHits>=1) {
178 if (referenceBLayerHits>=1) {