113 for (
int ibe = 0; ibe < 2; ibe++) {
115 IdentifierHash idHash =
index;
127 idSide = idBarrelEndcap ? 1 : -1;
129 if (
isBarrel && (idBarrelEndcap == -1)) {
133 }
else if (ibe == 1) {
134 idSide = idBarrelEndcap ? 1 : 0;
136 if (!
isBarrel && ((idBarrelEndcap == -2) || (idBarrelEndcap == 2))) {
142 if (sectorflag == 1) {
143 if (!element)
continue;
145 for (
unsigned int istraw = 0; istraw < element->
nStraws(); istraw++) {
146 std::vector<Identifier> neighbourIDs;
153 if (idLayerWheel == 1) i_chip += 21;
155 if (idLayerWheel == 2) i_chip += 54;
157 int tempStrawNumber =
strawNumber(istraw, idStrawLayer, idLayerWheel);
159 if (tempStrawNumber < 0 || tempStrawNumber > (
s_Straw_max[ibe] - 1)) {
160 ATH_MSG_WARNING(
"Found tempStrawNumber = " << tempStrawNumber <<
" out of range.");
162 m_mat_chip_B.at(idPhiModule).at(tempStrawNumber) = i_chip;
163 m_mat_chip_B.at(idPhiModule + 32).at(tempStrawNumber) = i_chip;
165 }
else if (ibe == 1) {
167 int tempStrawNumber =
strawNumberEndCap(istraw, idStrawLayer, idLayerWheel, idPhiModule, idSide);
169 if (tempStrawNumber < 0 || tempStrawNumber > (
s_Straw_max[ibe] - 1)) {
170 ATH_MSG_WARNING(
"Found tempStrawNumber = " << tempStrawNumber <<
" out of range.");
172 m_mat_chip_E.at(idPhiModule).at(tempStrawNumber) = i_chip;
173 m_mat_chip_E.at(idPhiModule + 32).at(tempStrawNumber) = i_chip;
189 return StatusCode::SUCCESS;
196 switch (LayerNumber) {
214 int addToStrawNumber = 0;
215 int addToStrawNumberNext = 0;
217 const int numberOfStraws[75] = {
244 addToStrawNumber += numberOfStraws[
i - 1];
245 addToStrawNumberNext = addToStrawNumber + numberOfStraws[
i];
250 if (strawNumber < 0 || strawNumber >
s_Straw_max[0] - 1) {
265 static const int TripletOrientation[2][32] = {
267 1, 1, 1, 1, 1, 1, 1, 1,
268 0, 0, 0, 0, 0, 0, 0, 0,
269 1, 1, 1, 1, 1, 1, 1, 1,
270 0, 0, 0, 0, 0, 0, 0, 0
273 1, 1, 1, 1, 1, 1, 1, 1,
274 0, 0, 0, 0, 0, 0, 0, 0,
275 1, 1, 1, 1, 1, 1, 1, 1,
276 0, 0, 0, 0, 0, 0, 0, 0
281 if (
side == 2) phi1 = phi_stack,
side = 1;
282 else if (
side == -2) phi1 = 31 - phi_stack,
side = 0;
285 if (TripletOrientation[
side][phi1]) {
311 int strawNumberNew = 0;
313 if (LayerNumber < 6 && strawLayerNumber > 7) {
314 strawNumberNew = strawNumberNew + (384 * LayerNumber);
317 strawNumberNew = strawNumberNew + (384 * LayerNumber);
320 strawNumberNew = strawNumberNew + 2304 + 192 * (LayerNumber - 6);
323 strawNumberNew = strawNumberNew + 2304 + 192 * (LayerNumber - 6);
329 if (strawNumber < 0 || strawNumber >
s_Straw_max[1] - 1) {
340 float degreeValue = radValue /
M_PI * 180;
342 if (degreeValue < 0) degreeValue += 360;
359 for (; RDO_CollectionBegin != RDO_CollectionEnd; ++RDO_CollectionBegin) {
362 if (!TRT_Collection)
continue;
366 for (; p_rdo != TRT_Collection->
end(); ++p_rdo) {
369 if (!p_lolum)
continue;
480 const float timeCor = comTimeObject ? comTimeObject->
getTime() : 0;
482 auto p_trk = trackCollection.
begin();
488 int ntrackstack[2][64];
489 int nTotalTracks = 0;
490 int nTracksB[2] = {0, 0};
491 int nTracksEC[2] = {0, 0};
492 int nTracksEC_B[2] = {0, 0};
493 int nTrksperLB_B = 0;
494 int nTrksperLB_E[2] = {0, 0};
496 for (
int ibe = 0; ibe < 2; ibe++) {
497 std::fill(ntrackstack[ibe], ntrackstack[ibe] + 64, 0);
500 for (; p_trk != trackCollection.
end(); ++p_trk) {
509 AllTrkPar = ((*p_trk)->track())->trackParameters();
515 for (p_trkpariter = AllTrkPar->
begin(); p_trkpariter != AllTrkPar->
end(); ++p_trkpariter) {
517 if ((mPer =
dynamic_cast<const Trk::Perigee *
>(*p_trkpariter)))
break;
529 const Trk::TrackStates *trackStates = ((*p_trk)->track())->trackStateOnSurfaces();
531 if (trackStates ==
nullptr)
continue;
539 int n_trt_hits =
unsigned(tempHitsVariable);
541 bool is_pT_over_20GeV =
false;
544 is_pT_over_20GeV =
true;
546 is_pT_over_20GeV =
false;
549 const bool cnst_is_pT_over_20GeV = is_pT_over_20GeV;
554 if (!passed_track_preselection)
continue;
557 int checkB[2] = {0, 0};
558 int checkEC[2] = {0, 0};
559 int checkEC_B[2] = {0, 0};
561 int nTRTHits_side[2][2];
562 int nTRTHitsW_perwheel[2][18];
563 int hitontrack[2] = {0, 0};
564 int hitontrack_E_side[2] = {0, 0};
566 for (
int ibe = 0; ibe < 2; ibe++) {
567 for (
int iside = 0; iside < 2; iside++) {
568 nTRTHits_side[ibe][iside] = -1;
569 nTRTHitsW[ibe][iside] = 0;
571 std::fill(nTRTHitsW_perwheel[ibe], nTRTHitsW_perwheel[ibe] + 18, 0);
575 int layer_or_wheel = 0;
579 int nearest_straw_layer[2] = {100, 100};
580 int nearest_straw[2] = {0, 0};
581 int testLayer[2] = {100, 100};
582 float phi2D[2] = {-100, -100};
584 for (TSOSItBeginTemp = TSOSItBegin0; TSOSItBeginTemp != TSOSItEnd; ++TSOSItBeginTemp) {
585 if ((*TSOSItBeginTemp) ==
nullptr)
continue;
590 if (!trtCircle)
continue;
593 if (!aTrackParam)
continue;
596 int ibe = std::abs(barrel_ec) - 1;
602 if (layer_or_wheel >= testLayer[ibe])
continue;
603 testLayer[ibe] = layer_or_wheel;
605 if (straw_layer < nearest_straw_layer[ibe]) {
606 nearest_straw_layer[ibe] = straw_layer;
607 nearest_straw[ibe] =
straw;
611 circleElement =
nullptr;
615 if (phi2D[0] == -999) {
616 ATH_MSG_DEBUG(
"Track did not go through inner layer of Barrel.");
619 testLayer[0] <<
" m_straw_layer: " <<
620 nearest_straw_layer[0] <<
" (in the Barrel).");
623 if (phi2D[1] == -999) {
624 ATH_MSG_DEBUG(
"Track did not go through any inner layer of EndCap A or C.");
627 testLayer[1] <<
" m_straw_layer: " <<
628 nearest_straw_layer[1] <<
" (in the EndCaps).");
631 bool trackfound[2][64];
633 for (
int i = 0;
i < 2;
i++) {
634 std::fill(trackfound[
i], trackfound[
i] + 64,
false);
637 for (TSOSItBegin = TSOSItBegin0; TSOSItBegin != TSOSItEnd; ++TSOSItBegin) {
639 if ((*TSOSItBegin) ==
nullptr)
continue;
645 if (!trtCircle)
continue;
649 if (!aTrackParam)
continue;
657 int ibe = std::abs(barrel_ec) - 1;
658 int iside = barrel_ec > 0 ? 0 : 1;
659 int thisStrawNumber[2] = {-1, -1};
660 int chip[2] = {0, 0};
665 if (thisStrawNumber[ibe] >= 0 && thisStrawNumber[ibe] <
s_Straw_max[ibe]) {
666 chip[ibe] =
m_mat_chip_B.at(phi_module).at(thisStrawNumber[ibe]);
668 }
else if (ibe == 1) {
671 if (thisStrawNumber[ibe] >= 0 && thisStrawNumber[ibe] <
s_Straw_max[ibe]) {
672 chip[ibe] =
m_mat_chip_E.at(phi_module).at(thisStrawNumber[ibe]);
675 thisStrawNumber[ibe] = -1;
678 if (thisStrawNumber[ibe] < 0 || thisStrawNumber[ibe] >=
s_Straw_max[ibe])
continue;
680 if (checkB[iside] == 0 && ibe == 0) {
685 if (checkEC[iside] == 0 && ibe == 1) {
690 if (checkEC_B[iside] == 0 && checkB[iside] == 1 && ibe == 1 ) {
691 nTracksEC_B[iside]++;
692 checkEC_B[iside] = 1;
702 int iphi_module = -9999;
704 if (iside == 0) iphi_module = phi_module;
705 else if (iside == 1) iphi_module = phi_module + 32;
707 if (iphi_module >= 0 && iphi_module < 64) trackfound[ibe][iphi_module] =
true;
718 }
else if (ibe == 1) {
720 hitontrack_E_side[iside]++;
726 if (RawDriftCircle) {
727 nTRTHits_side[ibe][iside]++;
732 if (driftTimeValid) {
733 const float validRawDriftTime = RawDriftCircle->
rawDriftTime();
736 ValidRawDriftTimeonTrkS_x = thisStrawNumber[ibe];
737 ValidRawDriftTimeonTrkS_y = validRawDriftTime;
742 ValidRawDriftTimeonTrkC_x = chip[ibe] - 1;
743 ValidRawDriftTimeonTrkC_y = validRawDriftTime;
751 DriftTimeonTrkDist_B_Ar = RawDriftCircle->
rawDriftTime();
758 }
else if (ibe == 1) {
760 DriftTimeonTrkDist_E_Ar = RawDriftCircle->
rawDriftTime();
770 float locR_err = 0.0;
771 const AmgSymMatrix(5)* b_err = aTrackParam->covariance();
775 ATH_MSG_WARNING(
"Some diagonal element(s) of the covariance matrix is (are) infinite or smaller than / too close to zero or above the covariance cutoff");
781 ATH_MSG_ERROR(
"Track parameters have no covariance attached.");
792 if ((loc *
locR) < 0) loc = -loc;
798 double pull_b = -999.;
799 const double diff_loc_err = std::abs(loc_err-locR_err);
800 if ( diff_loc_err > 0 ) {
801 pull_b = (loc -
locR) /diff_loc_err ;
804 else pull_b_fill =
false;
806 const double trkdrifttime = (!rtr) ? 0 : rtr->
drifttime(std::abs(
locR));
807 const double timeresidual = RawDriftCircle->
rawDriftTime() - thist0 - trkdrifttime;
812 Pull_Biased_Barrel = pull_b;
818 Residual_B_Ar = loc -
locR;
820 Residual_noTubeHits_B_Ar = loc -
locR;
821 if (!isTubeHit)
fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), Residual_noTubeHits_B_Ar);
823 if (cnst_is_pT_over_20GeV) {
824 Residual_B_Ar_20GeV = loc -
locR;
826 Residual_noTubeHits_B_Ar_20GeV = loc -
locR;
827 if (!isTubeHit)
fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), Residual_noTubeHits_B_Ar_20GeV);
829 TimeResidual_B_Ar = timeresidual;
831 TimeResidual_noTubeHits_B_Ar = timeresidual;
832 if (!isTubeHit)
fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), TimeResidual_noTubeHits_B_Ar);
834 Residual_B = loc -
locR;
836 Residual_noTubeHits_B = loc -
locR;
837 if (!isTubeHit)
fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), Residual_noTubeHits_B);
838 TimeResidual_B = timeresidual;
840 TimeResidual_noTubeHits_B = timeresidual;
841 if (!isTubeHit)
fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), TimeResidual_noTubeHits_B);
843 if (cnst_is_pT_over_20GeV) {
844 Residual_B_20GeV = loc -
locR;
846 Residual_noTubeHits_B_20GeV = loc -
locR;
847 if (!isTubeHit)
fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), Residual_noTubeHits_B_20GeV);
850 }
else if (ibe == 1) {
853 Pull_Biased_EndCap = pull_b;
859 Residual_E_Ar = loc -
locR;
861 Residual_noTubeHits_E_Ar = loc -
locR;
863 TimeResidual_E_Ar = timeresidual;
865 TimeResidual_noTubeHits_E_Ar = timeresidual;
868 if (cnst_is_pT_over_20GeV) {
869 Residual_E_Ar_20GeV = loc -
locR;
871 Residual_noTubeHits_E_Ar_20GeV = loc -
locR;
875 Residual_E = loc -
locR;
877 Residual_noTubeHits_E = loc -
locR;
879 TimeResidual_E = timeresidual;
881 TimeResidual_noTubeHits_E = timeresidual;
884 if (cnst_is_pT_over_20GeV) {
885 Residual_E_20GeV = loc -
locR;
887 Residual_noTubeHits_E_20GeV = loc -
locR;
897 WireToTrkPosition_B_Ar =
locR;
900 WireToTrkPosition_B =
locR;
903 }
else if (ibe == 1) {
905 WireToTrkPosition_E_Ar =
locR;
908 WireToTrkPosition_E =
locR;
915 const float EP = timeCor;
921 RtRelation_B_Ar_x =
LE - EP - t0;
922 RtRelation_B_Ar_y = std::abs(
locR);
923 fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), RtRelation_B_Ar_x, RtRelation_B_Ar_y);
925 RtRelation_B_Ar_x =
LE - t0;
926 RtRelation_B_Ar_y = std::abs(
locR);
927 fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), RtRelation_B_Ar_x, RtRelation_B_Ar_y);
931 RtRelation_B_x =
LE - EP - t0;
932 RtRelation_B_y = std::abs(
locR);
935 RtRelation_B_x =
LE - t0;
936 RtRelation_B_y = std::abs(
locR);
940 }
else if (ibe == 1) {
943 RtRelation_E_Ar_x =
LE - EP - t0;
944 RtRelation_E_Ar_y = std::abs(
locR);
947 RtRelation_E_Ar_x =
LE - t0;
948 RtRelation_E_Ar_y = std::abs(
locR);
953 RtRelation_E_x =
LE - EP - t0;
954 RtRelation_E_y = std::abs(
locR);
957 RtRelation_E_x =
LE - t0;
958 RtRelation_E_y = std::abs(
locR);
978 HitTronTMapC_x = chip[ibe] - 1;
979 HitTronTMapC_y = trailingEdgeScaled;
984 HitTronTwEPCMapS_x = thisStrawNumber[ibe];
985 HitTronTwEPCMapS_y = trailingEdgeScaled - timeCor;
990 HitTronTwEPCMapC_x = chip[ibe] - 1;
991 HitTronTwEPCMapC_y = trailingEdgeScaled - timeCor;
999 TronTDist_B_Ar = trailingEdgeScaled;
1001 AvgTroTDetPhi_B_Ar_x = phi2D[ibe];
1002 AvgTroTDetPhi_B_Ar_y = trailingEdgeScaled;
1003 fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), AvgTroTDetPhi_B_Ar_x, AvgTroTDetPhi_B_Ar_y);
1005 TronTDist_B = trailingEdgeScaled;
1007 AvgTroTDetPhi_B_x = phi2D[ibe];
1008 AvgTroTDetPhi_B_y = trailingEdgeScaled;
1009 fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), AvgTroTDetPhi_B_x, AvgTroTDetPhi_B_y);
1011 }
else if (ibe == 1) {
1013 TronTDist_E_Ar = trailingEdgeScaled;
1015 AvgTroTDetPhi_E_Ar_x = phi2D[ibe];
1016 AvgTroTDetPhi_E_Ar_y = trailingEdgeScaled;
1019 TronTDist_E = trailingEdgeScaled;
1021 AvgTroTDetPhi_E_x = phi2D[ibe];
1022 AvgTroTDetPhi_E_y = trailingEdgeScaled;
1033 for (
int ibe = 0; ibe < 2; ibe++) {
1034 for (
int i = 0;
i < 64;
i++)
1035 if (trackfound[ibe][
i])
1036 ntrackstack[ibe][
i]++;
1041 NumHoTDetPhi_B_x = phi2D[ibe];
1042 NumHoTDetPhi_B_y = hitontrack[ibe];
1043 fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), NumHoTDetPhi_B_x, NumHoTDetPhi_B_y);
1049 NumHoTDetPhi_E_x = phi2D[ibe];
1050 NumHoTDetPhi_E_y = hitontrack_E_side[0];
1051 fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe)+
"0", NumHoTDetPhi_E_x, NumHoTDetPhi_E_y);
1055 NumHoTDetPhi_E_x = phi2D[ibe];
1056 NumHoTDetPhi_E_y = hitontrack_E_side[1];
1057 fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe)+
"1", NumHoTDetPhi_E_x, NumHoTDetPhi_E_y);
1062 if (phi2D[ibe] < 0)
continue;
1066 if (nTRTHitsW[ibe][0] + nTRTHitsW[ibe][1] > 0) {
1067 NumTrksDetPhi_B = phi2D[ibe];
1070 }
else if (ibe == 1) {
1071 if (nTRTHitsW[ibe][0] > 0) {
1072 NumTrksDetPhi_E = phi2D[ibe];
1076 if (nTRTHitsW[ibe][1] > 0) {
1077 NumTrksDetPhi_E = phi2D[ibe];
1084 if ((nTRTHitsW[ibe][0] + nTRTHitsW[ibe][1]) > 0) {
1087 if (comTimeObject) {
1088 if (
m_doShift && (phi2D[ibe] > 0) && (std::abs(timeCor) > 1
e-8)) {
1089 EvtPhaseDetPhi_B_x = phi2D[ibe];
1090 EvtPhaseDetPhi_B_y = timeCor;
1091 fill(
"ShiftTRTTrackHistograms"+
std::to_string(ibe), EvtPhaseDetPhi_B_x, EvtPhaseDetPhi_B_y);
1094 }
else if (ibe == 1) {
1095 for (
int iside = 0; iside < 2; iside++) {
1096 if (nTRTHitsW[ibe][iside] > 0) nTrksperLB_E[iside]++;
1097 if (comTimeObject) {
1098 if (nTRTHits_side[ibe][iside] > 5 && (std::abs(timeCor)
1101 EvtPhaseDetPhi_E_x = phi2D[ibe];
1102 EvtPhaseDetPhi_E_y = timeCor;
1112 if (comTimeObject) {
1113 if (std::abs(timeCor) > 1
e-8) {
1116 fill(
"ShiftTRTTrackHistograms0", EvtPhase);
1120 std::vector<int> trigid;
1123 std::vector<unsigned int> level1TAV = trigDecision->
tav();
1125 for (
unsigned int j = 0; j < 8 && j < level1TAV.size(); ++j) {
1126 for (
unsigned int i = 0;
i < 32; ++
i) {
1127 if ((level1TAV[j] >>
i) & 0
x1) {
1128 trigid.push_back(
i + (j % 8) * 32);
1133 for (
unsigned int j = 0; j < trigid.size(); ++j) {
1134 EvtPhaseVsTrig_x = timeCor;
1135 EvtPhaseVsTrig_y = trigid[j];
1136 fill(
"ShiftTRTTrackHistograms0", EvtPhaseVsTrig_x, EvtPhaseVsTrig_y);
1149 SummaryWeight = nTotalTracks;
1152 SummaryWeight = nTracksB[0];
1155 SummaryWeight = nTracksB[1];
1158 SummaryWeight = nTracksEC[0];
1161 SummaryWeight = nTracksEC[1];
1164 SummaryWeight = nTracksEC_B[0];
1167 SummaryWeight = nTracksEC_B[1];
1173 int lastLumiBlock = -99;
1177 float evtLumiBlock = 1.;
1178 float lumiBlockScale = (evtLumiBlock > 0) ? (1. / evtLumiBlock) : 0;
1181 NTrksperLB_x = lastLumiBlock;
1182 NTrksperLB_y = (
float)nTrksperLB_B * lumiBlockScale;
1183 fill(
"ShiftTRTTrackHistograms0", NTrksperLB_x, NTrksperLB_y);
1185 for (
int iside = 0; iside < 2; iside++) {
1186 NTrksperLB_x = lastLumiBlock;
1187 NTrksperLB_y = (
float)nTrksperLB_E[iside] * lumiBlockScale;
1193 for (
int iside = 0; iside < 2; iside++) {
1194 nTrksperLB_E[iside] = 0;
1202 return StatusCode::SUCCESS;
1217 int lumiBlockNumber;
1219 lumiBlockNumber = eventInfo.
lumiBlock();
1227 IntLumWeight = intLum;
1228 fill(
"SmryHistograms", IntLumWeight, IntLum);
1229 LBvsLum = lumiBlockNumber;
1230 LBvsLumWeight = intLum;
1231 fill(
"SmryHistograms", LBvsLumWeight, LBvsLum);
1232 LBvsTime_x = lumiBlockNumber;
1234 fill(
"SmryHistograms", LBvsTime_x, LBvsTime_y);
1243 auto p_trk = trackCollection.
begin();
1248 for (; p_trk != trackCollection.
end(); ++p_trk) {
1249 AllTrkPar = ((*p_trk)->track())->trackParameters();
1251 for (p_trkpariter = AllTrkPar->
begin(); p_trkpariter != AllTrkPar->
end(); ++p_trkpariter) {
1252 if ((perigee =
dynamic_cast<const Trk::Perigee *
>(*p_trkpariter)))
break;
1259 ATH_MSG_DEBUG(
"No perigee mearsurement found for the track. This entry will not be propogated to aging histograms.");
1263 float track_eta = perigee->eta();
1264 float track_p = (perigee->parameters()[
Trk::qOverP] != 0.) ? std::abs(1. / (perigee->parameters()[
Trk::qOverP])) : 10e7;
1265 const Trk::TrackStates *trackStates = ((*p_trk)->track())->trackStateOnSurfaces();
1267 if (trackStates ==
nullptr)
continue;
1274 int trt_hits =
unsigned(tempHitsVariable);
1278 std::abs(track_p) >= 5000.;
1280 if (!passed_track_preselection)
continue;
1283 const Trk::TrackStates *track_states = ((*p_trk)->track())->trackStateOnSurfaces();
1286 ATH_MSG_DEBUG(
"This track has " << track_states->
size() <<
" track states on surface.");
1288 ATH_MSG_DEBUG(
"This track has null track states on surface.");
1292 int barrel_ec_side = 0;
1293 int layer_or_wheel = 0;
1295 int straw_layer = 0;
1297 for (; TSOSItBegin != TSOSItEnd; ++TSOSItBegin) {
1298 if ((*TSOSItBegin) ==
nullptr)
continue;
1304 if (!trtCircle)
continue;
1305 if (!aTrackParam)
continue;
1313 int Ba_Ec = abs(barrel_ec_side) - 1;
1314 int Side = barrel_ec_side > 0 ? 0 : 1;
1318 double RPos = sqrt(xPos * xPos + yPos * yPos);
1324 if (!RawDriftCircle) {
1330 int middleHTbit = RawDriftCircle->
getWord() & 0x00020000;
1332 bool is_middleHTbit_high = (middleHTbit != 0);
1334 bool isHighLevel = is_middleHTbit_high;
1335 bool shortStraw =
false;
1338 if (std::abs(track_eta) < 2. && Ba_Ec == 0.) {
1339 if ((layer_or_wheel == 0) && (phi_module < 4 || (phi_module > 7 && phi_module < 12) || (phi_module > 15 && phi_module < 20) || (phi_module > 23 && phi_module < 28))) InputBar = 1;
1340 else if ((
runNumber >= 296939) && (layer_or_wheel == 0) && (phi_module > 27)) InputBar = 1;
1341 else if (layer_or_wheel == 0)
1343 else if ((layer_or_wheel == 1) && ((phi_module > 1 && phi_module < 6) || (phi_module > 9 && phi_module < 14) || (phi_module > 17 && phi_module < 22) || (phi_module > 25 && phi_module < 30)))
1345 else if (layer_or_wheel == 1)
1347 else if (layer_or_wheel == 2 && phi_module % 2 != 0)
1349 else if (layer_or_wheel == 2)
1356 if ((layer_or_wheel == 0) && straw_layer < 9.)
1361 if ((!shortStraw) && (Ba_Ec == 0)) {
1392 if ((Ba_Ec == 1) && (layer_or_wheel < 6) &&
1393 ((straw_layer > 3 && straw_layer < 8) ||
1394 (straw_layer > 11))) {
1397 if ((Ba_Ec == 1) && (layer_or_wheel >= 6) &&
1398 (straw_layer > 3)) {
1401 if ((Ba_Ec == 1) && (layer_or_wheel < 6) &&
1402 ((straw_layer > -1 && straw_layer < 4) ||
1403 (straw_layer > 7 && straw_layer < 12))) {
1406 if ((Ba_Ec == 1) && (layer_or_wheel >= 6) &&
1407 ((straw_layer > -1 && straw_layer < 4))) {
1411 if (WType < 0 && Ba_Ec == 1) {
1427 return StatusCode::SUCCESS;
1433 bool passEventBurst =
true;
1445 if (!xAODEventInfo.
isValid()) {
1448 return StatusCode::FAILURE;
1452 if (!trackCollection.
isValid()) {
1455 return StatusCode::FAILURE;
1460 if (!trigDecision) {
1465 const ComTime *comTimeObject=
nullptr;
1468 if (!tmp_comTimeObject.
isValid()) {
1473 comTimeObject = tmp_comTimeObject.
cptr();
1480 if (!trackCollection.
isValid()) {
1483 return StatusCode::FAILURE;
1487 if (passEventBurst) {
1493 return StatusCode::SUCCESS;