21 m_h_res_eta =
new TH1D(std::string(
m_name+
"-Res-eta").c_str(), std::string(
m_name+
" eta resolution").c_str(), 100, -0.02, 0.02);
22 m_h_res_phi =
new TH1D(std::string(
m_name+
"-Res-phi").c_str(), std::string(
m_name+
" phi resolution").c_str(), 100, -0.01, 0.01);
23 m_h_res_z0 =
new TH1D(std::string(
m_name+
"-Res-z0").c_str(), std::string(
m_name+
" z0 resolution").c_str(), 100, -3.0, 3.0);
24 m_h_res_d0 =
new TH1D(std::string(
m_name+
"-Res-d0").c_str(), std::string(
m_name+
" d0 resolution").c_str(), 100, -3.0, 3.0);
25 m_h_res_invpT =
new TH1D(std::string(
m_name+
"-Res-pT").c_str(), std::string(
m_name+
" inv-pT resolution").c_str(), 100, -0.0004, 0.0004);
33 m_h_pull_eta =
new TH1D(std::string(
m_name+
"-Pull-eta").c_str(), std::string(
m_name+
" eta pull").c_str(), 100, -20.0, 20.0);
34 m_h_pull_phi =
new TH1D(std::string(
m_name+
"-Pull-phi").c_str(), std::string(
m_name+
" phi pull").c_str(), 100, -20.0, 20.0);
35 m_h_pull_z0 =
new TH1D(std::string(
m_name+
"-Pull-z0").c_str(), std::string(
m_name+
" z0 pull").c_str(), 100, -20.0, 20.0);
36 m_h_pull_d0 =
new TH1D(std::string(
m_name+
"-Pull-d0").c_str(), std::string(
m_name+
" d0 pull").c_str(), 100, -20.0, 20.0);
37 m_h_pull_invpT =
new TH1D(std::string(
m_name+
"-Pull-pT").c_str(), std::string(
m_name+
" inv-pT pull").c_str(), 100, -20.0, 20.0);
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() );
138 m_h_res_eta->Fill(fabs(referenceEta)-fabs(testEta));
140 m_h_res_z0->Fill(fabs(referenceZ0)-fabs(testZ0));
141 m_h_res_d0->Fill(fabs(referenceD0)-fabs(testD0));
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) {