56 bool isMiddleFailure =
true;
59 const double R_WIDTH_DEFAULT = 100;
60 const double R_WIDTH_INNER_NO_HIT = 200;
83 if (tgcHits.size()>0) {
87 if( !
prepareTgcPoints(tgcHits, tgcStripInnPoints, tgcWireInnPoints, tgcStripMidPoints, tgcWireMidPoints) ) {
89 return StatusCode::FAILURE;
94 isMiddleFailure =
false;
98 isMiddleFailure =
true;
101 ATH_MSG_DEBUG(
"Fit to TGC middle station returns only a point");
105 status =
m_tgcFit->runTgcInner(tgcStripInnPoints, tgcWireInnPoints, tgcFitResult);
110 ATH_MSG_DEBUG(
"Fit to TGC inner station returns only a point");
114 <<
"/" << tgcFitResult.
tgcInn[2] <<
"/" << tgcFitResult.
tgcInn[3]);
116 <<
"/" << tgcFitResult.
tgcMid1[2] <<
"/" << tgcFitResult.
tgcMid1[3]);
118 <<
"/" << tgcFitResult.
tgcMid2[2] <<
"/" << tgcFitResult.
tgcMid2[3]);
123 if (tgcHits.size()>0 && !isMiddleFailure){
128 const double PHI_RANGE = 12./(
M_PI/8.);
134 double PhiInOctant = std::abs(tgcFitResult.
tgcMid1[1] - Octant * (
M_PI/4.));
135 if (PhiInOctant > (
M_PI/8.)) PhiInOctant = (
M_PI/4.) - PhiInOctant;
137 int phiBin =
static_cast<int>(PhiInOctant * PHI_RANGE);
138 int etaBin =
static_cast<int>((std::abs(tgcFitResult.
tgcMid1[0]) - 1.)/0.05);
149 if (std::abs(tgcFitResult.
tgcMid1[3])>=0.) tgcFitResult.
phi = tgcFitResult.
tgcMid1[1];
150 if (std::abs(tgcFitResult.
tgcMid2[3])>=0.) tgcFitResult.
phi = tgcFitResult.
tgcMid2[1];
159 double tmp1 = (tgcFitResult.
tgcMid1[1]>0)?
186 muonRoad.
aw[endcap_middle][0] = tgcFitResult.
slope;
187 muonRoad.
bw[endcap_middle][0] = tgcFitResult.
intercept;
188 muonRoad.
aw[endcap_outer][0] = tgcFitResult.
slope;
189 muonRoad.
bw[endcap_outer][0] = tgcFitResult.
intercept;
190 muonRoad.
aw[endcap_extra][0] = tgcFitResult.
slope;
191 muonRoad.
bw[endcap_extra][0] = tgcFitResult.
intercept;
192 muonRoad.
aw[bee][0] = tgcFitResult.
slope;
194 for (
int i_layer=0; i_layer<
N_LAYER; i_layer++) {
195 muonRoad.
rWidth[endcap_middle][i_layer] = R_WIDTH_DEFAULT;
196 muonRoad.
rWidth[endcap_outer][i_layer] = R_WIDTH_DEFAULT;
197 muonRoad.
rWidth[endcap_extra][i_layer] = R_WIDTH_DEFAULT;
198 muonRoad.
rWidth[bee][i_layer] = R_WIDTH_DEFAULT;
202 muonRoad.
aw[endcap_inner][0] = tgcFitResult.
tgcInn[2]/tgcFitResult.
tgcInn[3];
203 muonRoad.
aw[barrel_inner][0] = tgcFitResult.
tgcInn[2]/tgcFitResult.
tgcInn[3];
204 muonRoad.
aw[csc][0] = tgcFitResult.
tgcInn[2]/tgcFitResult.
tgcInn[3];
205 for (
int i_layer=0; i_layer<
N_LAYER; i_layer++) muonRoad.
rWidth[endcap_inner][i_layer] = R_WIDTH_DEFAULT;
206 for (
int i_layer=0; i_layer<
N_LAYER; i_layer++) muonRoad.
rWidth[barrel_inner][i_layer] = R_WIDTH_DEFAULT;
207 for (
int i_layer=0; i_layer<
N_LAYER; i_layer++) muonRoad.
rWidth[csc][i_layer] = R_WIDTH_DEFAULT;
215 double extrInnerEta = 0;
231 = (*m_backExtrapolatorTool)->give_eta_phi_at_vertex(muonSA,
eta,sigma_eta,
phi,sigma_phi,0.);
232 if (
sc.isSuccess() ){
235 extrInnerEta = etaMiddle;
239 return StatusCode::FAILURE;
242 if (muonSA)
delete muonSA;
245 if (extrInnerEta != 0.) {
247 muonRoad.
aw[endcap_inner][0] =
std::tan(
theta)*(std::abs(extrInnerEta)/extrInnerEta);
248 muonRoad.
aw[barrel_inner][0] =
std::tan(
theta)*(std::abs(extrInnerEta)/extrInnerEta);
249 muonRoad.
aw[csc][0] =
std::tan(
theta)*(std::abs(extrInnerEta)/extrInnerEta);
251 muonRoad.
aw[endcap_inner][0] = 0;
252 muonRoad.
aw[barrel_inner][0] = 0;
253 muonRoad.
aw[csc][0] = 0;
256 for (
int i_layer=0; i_layer<
N_LAYER; i_layer++) muonRoad.
rWidth[endcap_inner][i_layer] = R_WIDTH_INNER_NO_HIT;
257 for (
int i_layer=0; i_layer<
N_LAYER; i_layer++) muonRoad.
rWidth[barrel_inner][i_layer] = R_WIDTH_INNER_NO_HIT;
258 for (
int i_layer=0; i_layer<
N_LAYER; i_layer++) muonRoad.
rWidth[csc][i_layer] = R_WIDTH_INNER_NO_HIT;
262 muonRoad.
bw[endcap_inner][0] = 0.;
263 muonRoad.
bw[barrel_inner][0] = 0.;
264 muonRoad.
bw[csc][0] = 0.;
269 ATH_MSG_DEBUG(
"Because no TGC hits are available, estimate the road from RoI");
275 muonRoad.
aw[endcap_inner][0] = aw;
276 muonRoad.
bw[endcap_inner][0] = 0;
277 muonRoad.
aw[endcap_middle][0] = aw;
278 muonRoad.
bw[endcap_middle][0] = 0;
279 muonRoad.
aw[endcap_outer][0] = aw;
280 muonRoad.
bw[endcap_outer][0] = 0;
281 muonRoad.
aw[endcap_extra][0] = aw;
282 muonRoad.
bw[endcap_extra][0] = 0;
283 muonRoad.
aw[barrel_inner][0] = aw;
284 muonRoad.
bw[barrel_inner][0] = 0;
285 muonRoad.
aw[bee][0] = aw;
286 muonRoad.
bw[bee][0] = 0;
287 muonRoad.
aw[csc][0] = aw;
288 muonRoad.
bw[csc][0] = 0;
289 for (
int i_layer=0; i_layer<
N_LAYER; i_layer++) {
291 muonRoad.
rWidth[endcap_inner][i_layer] = 300;
292 muonRoad.
rWidth[endcap_middle][i_layer] = 400;
293 muonRoad.
rWidth[endcap_outer][i_layer] = 600;
294 muonRoad.
rWidth[endcap_extra][i_layer] = 400;
295 muonRoad.
rWidth[barrel_inner][i_layer] = 250;
296 muonRoad.
rWidth[bee][i_layer] = 500;
297 muonRoad.
rWidth[csc][i_layer] = 200;
311 muonRoad.
aw[
i][1] = muonRoad.
aw[
i][0];
312 muonRoad.
bw[
i][1] = muonRoad.
bw[
i][0];
318 if( ! isMiddleFailure ) {
320 phiMiddle = tgcFitResult.
tgcMid1[1];
322 side = (p_roids->
eta()<0.)? 0 : 1;
323 phiMiddle = p_roids->
phi();
329 int sector_trigger = 99;
330 int sector_overlap = 99;
333 float tempDeltaPhi=99;
334 std::vector<Identifier> stationList;
335 std::vector<IdentifierHash> mdtHashList;
344 if(phiMax >
M_PI) phiMax -=
M_PI*2.;
345 if(phiMin <
M_PI*-1) phiMin +=
M_PI*2.;
363 if (
name.substr(0, 1) ==
"B" )
continue;
364 if (
name.substr(1, 1) !=
"M" )
continue;
367 if (
name[2]==
'S' ||
name[2]==
'E') floatPhi = floatPhi +
M_PI/8;
368 tempDeltaPhi = std::abs(floatPhi-muonRoad.
phiMiddle);
369 if (phiMiddle<0) tempDeltaPhi = std::abs(floatPhi-muonRoad.
phiMiddle-2*
M_PI);
370 if(tempDeltaPhi >
M_PI) tempDeltaPhi = std::abs(tempDeltaPhi - 2*
M_PI);
373 if(
name[2]==
'S' ||
name[2]==
'E') LargeSmall = 1;
375 if(sector_trigger == 99)
376 sector_trigger = sector;
377 else if(sector_trigger != sector)
378 sector_overlap = sector;
382 temp_sector = sector;
387 if(temp_sector != sector_trigger){
388 sector_overlap = sector_trigger;
389 sector_trigger = temp_sector;
399 for (
int i_sector=0; i_sector<
N_SECTOR; i_sector++) {
401 muonRoad.
aw[endcap_inner][i_sector] = muonRoad.
aw_ftf[3][0];
402 muonRoad.
bw[endcap_inner][i_sector] = muonRoad.
bw_ftf[3][0];
403 muonRoad.
aw[endcap_middle][i_sector] = muonRoad.
aw_ftf[4][0];
404 muonRoad.
bw[endcap_middle][i_sector] = muonRoad.
bw_ftf[4][0];
405 muonRoad.
aw[endcap_outer][i_sector] = muonRoad.
aw_ftf[5][0];
406 muonRoad.
bw[endcap_outer][i_sector] = muonRoad.
bw_ftf[5][0];
407 muonRoad.
aw[endcap_extra][i_sector] = muonRoad.
aw_ftf[6][0];
408 muonRoad.
bw[endcap_extra][i_sector] = muonRoad.
bw_ftf[6][0];
409 muonRoad.
aw[csc][i_sector] = muonRoad.
aw_ftf[7][0];
410 muonRoad.
bw[csc][i_sector] = muonRoad.
bw_ftf[7][0];
411 muonRoad.
aw[bee][i_sector] = muonRoad.
aw_ftf[8][0];
412 muonRoad.
bw[bee][i_sector] = muonRoad.
bw_ftf[8][0];
413 muonRoad.
aw[barrel_inner][i_sector] = muonRoad.
aw_ftf[0][0];
414 muonRoad.
bw[barrel_inner][i_sector] = muonRoad.
bw_ftf[0][0];
419 return StatusCode::SUCCESS;