16 |
const float TofDedxHCut::_S31_SA_params[] = { 1.6, 1.2, -0.15, -2., 0.8, 1.6, 1.2, -0.2, 0., 0.8 }; |
const float TofDedxHCut::_S31_SA_params[] = { 1.6, 1.2, -0.15, -2., 0.8, 1.6, 1.2, -0.2, 0., 0.8 }; |
17 |
const float TofDedxHCut::_S32_SA_params[] = { 1.6, 1.2, -0.15, -2., 0.8, 1.6, 1.2, -0.2, 0., 0.8 }; |
const float TofDedxHCut::_S32_SA_params[] = { 1.6, 1.2, -0.15, -2., 0.8, 1.6, 1.2, -0.2, 0., 0.8 }; |
18 |
|
|
19 |
TofDedxHCut::TofDedxHCut(const char *cutName, unsigned int layer, bool standAlone) : |
const float TofDedxHCut::_S11_SA_params_strict[] = { 1.5, 1.2, -0.15, -1.7, 0.8, 1.6, 1.2, -0.2, -0.8, 0.8 }; |
20 |
PamCut(cutName), _layer(layer), _standAlone(standAlone), _p(NULL) { |
const float TofDedxHCut::_S12_SA_params_strict[] = { 1.5, 1.2, -0.15, -1.7, 0.8, 1.6, 1.2, -0.2, -0.8, 0.8 }; |
21 |
|
const float TofDedxHCut::_S21_SA_params_strict[] = { 1.6, 1.2, -0.15, -1.7, 0.8, 1.6, 1.2, -0.2, -0.8, 0.8 }; |
22 |
switch (_layer) { |
const float TofDedxHCut::_S22_SA_params_strict[] = { 1.6, 1.2, -0.15, -1.7, 0.8, 1.6, 1.2, -0.2, -0.8, 0.8 }; |
23 |
case S11: |
const float TofDedxHCut::_S31_SA_params_strict[] = { 1.6, 1.2, -0.15, -1.7, 0.8, 1.6, 1.2, -0.2, -0.8, 0.8 }; |
24 |
_p = _S11_SA_params; |
const float TofDedxHCut::_S32_SA_params_strict[] = { 1.6, 1.2, -0.15, -1.7, 0.8, 1.6, 1.2, -0.2, -0.8, 0.8 }; |
25 |
break; |
|
26 |
case S12: |
TofDedxHCut::TofDedxHCut(const char *cutName, unsigned int layer, bool standAlone, bool strict) : |
27 |
_p = _S12_SA_params; |
PamCut(cutName), _layer(layer), _standAlone(standAlone), _p(NULL) { |
28 |
break; |
|
29 |
case S21: |
if (strict) { |
30 |
_p = _S21_SA_params; |
switch (_layer) { |
31 |
break; |
case S11: |
32 |
case S22: |
_p = _S11_SA_params_strict; |
33 |
_p = _S22_SA_params; |
break; |
34 |
break; |
case S12: |
35 |
case S31: |
_p = _S12_SA_params_strict; |
36 |
_p = _S31_SA_params; |
break; |
37 |
break; |
case S21: |
38 |
case S32: |
_p = _S21_SA_params_strict; |
39 |
_p = _S32_SA_params; |
break; |
40 |
break; |
case S22: |
41 |
default: |
_p = _S22_SA_params_strict; |
42 |
_p = _S21_SA_params; |
break; |
43 |
|
case S31: |
44 |
|
_p = _S31_SA_params_strict; |
45 |
|
break; |
46 |
|
case S32: |
47 |
|
_p = _S32_SA_params_strict; |
48 |
|
break; |
49 |
|
default: |
50 |
|
_p = _S21_SA_params_strict; |
51 |
|
} |
52 |
|
} |
53 |
|
else { |
54 |
|
switch (_layer) { |
55 |
|
case S11: |
56 |
|
_p = _S11_SA_params; |
57 |
|
break; |
58 |
|
case S12: |
59 |
|
_p = _S12_SA_params; |
60 |
|
break; |
61 |
|
case S21: |
62 |
|
_p = _S21_SA_params; |
63 |
|
break; |
64 |
|
case S22: |
65 |
|
_p = _S22_SA_params; |
66 |
|
break; |
67 |
|
case S31: |
68 |
|
_p = _S31_SA_params; |
69 |
|
break; |
70 |
|
case S32: |
71 |
|
_p = _S32_SA_params; |
72 |
|
break; |
73 |
|
default: |
74 |
|
_p = _S21_SA_params; |
75 |
|
} |
76 |
} |
} |
|
|
|
77 |
} |
} |
78 |
|
|
79 |
int TofDedxHCut::Check(PamLevel2 *event) { |
int TofDedxHCut::Check(PamLevel2 *event) { |
80 |
|
|
81 |
static int trkSeqNo; |
static int trkSeqNo; |
82 |
|
static int idx; |
83 |
if (_standAlone) |
if (_standAlone) |
84 |
trkSeqNo = 0; |
trkSeqNo = -1; |
85 |
else { |
else |
86 |
for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) { |
trkSeqNo = event->GetTrack(0)->GetTrkTrack()->GetSeqNo(); |
87 |
if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == event->GetTrack(0)->GetTrkTrack()->GetSeqNo()) |
|
88 |
trkSeqNo = i; |
for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) { |
89 |
} |
if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == trkSeqNo) |
90 |
|
idx = i; |
91 |
} |
} |
92 |
|
|
93 |
float beta = event->GetToFLevel2()->GetToFTrkVar(trkSeqNo)->beta[12]; |
float beta = event->GetToFLevel2()->GetToFTrkVar(idx)->beta[12]; |
94 |
|
|
95 |
static unsigned int auxLayer; |
static unsigned int auxLayer; |
96 |
static float dEdx; |
static float dEdx; |
103 |
nLayers = 0; |
nLayers = 0; |
104 |
while (auxLayer > 0) { |
while (auxLayer > 0) { |
105 |
if ((auxLayer & 0x1) == 1) { |
if ((auxLayer & 0x1) == 1) { |
106 |
dEdx += event->GetToFLevel2()->GetdEdx(trkSeqNo, currLayer, 100); |
dEdx += event->GetToFLevel2()->GetdEdx(idx, currLayer, 100); |
107 |
nLayers++; |
nLayers++; |
108 |
|
// Patch to exclude S115b |
109 |
|
if (currLayer == 0) { |
110 |
|
ToFTrkVar *tofTrk = event->GetToFLevel2()->GetToFTrkVar(idx); |
111 |
|
bool has8 = false, has9 = false; |
112 |
|
float dedx9 = 0.; |
113 |
|
for (Int_t ib = 0; ib < tofTrk->npmtadc; ib++) { |
114 |
|
if (tofTrk->pmtadc[ib] == 8) { |
115 |
|
has8 = true; |
116 |
|
} |
117 |
|
if (tofTrk->pmtadc[ib] == 9) { |
118 |
|
has9 = true; |
119 |
|
dedx9 = tofTrk->dedx[ib]; |
120 |
|
} |
121 |
|
} |
122 |
|
|
123 |
|
if (has9) { |
124 |
|
dEdx -= (dedx9 / ((int) has8 + 1)); |
125 |
|
dEdx *= 2.; |
126 |
|
} |
127 |
|
} |
128 |
|
//End of patch |
129 |
} |
} |
130 |
currLayer++; |
currLayer++; |
131 |
auxLayer >>= 1; |
auxLayer >>= 1; |
133 |
|
|
134 |
dEdx /= nLayers; |
dEdx /= nLayers; |
135 |
|
|
136 |
if (dEdx < 1 / (_p[0] * beta + _p[1]) + _p[2] + _p[3]*beta) |
if (dEdx < _p[0] / (_p[1] * beta + _p[2]) + _p[3] + _p[4] * beta) |
137 |
return 0; |
return 0; |
138 |
|
|
139 |
if (dEdx > 1 / (_p[4] * beta + _p[5]) + _p[6] + _p[7]*beta) |
if (dEdx > _p[5] / (_p[6] * beta + _p[7]) + _p[8] + _p[9] * beta) |
140 |
return 0; |
return 0; |
141 |
|
|
142 |
return CUTOK; |
return CUTOK; |