Skip to content
Snippets Groups Projects
Commit 506bcaf6 authored by Nikita Belyaev's avatar Nikita Belyaev
Browse files

Tube hits excluded from time residual

parent 8dc40852
No related merge requests found
......@@ -819,6 +819,10 @@ for (; p_trk != trackCollection.end(); ++p_trk) {
float loc_err = sqrt(Amg::error(trtCircle->localCovariance(), Trk::driftRadius));
float locR = aTrackParam->parameters()[Trk::driftRadius];
float loc = trtCircle->localParameters()[Trk::driftRadius];
const Trk::MeasurementBase* measurement = (*TSOSItBegin)->measurementOnTrack();
const Trk::RIO_OnTrack* rotp = dynamic_cast<const Trk::RIO_OnTrack*>(measurement);
float errDC = sqrt(rotp->localCovariance()(Trk::driftRadius, Trk::driftRadius));
//std::cout << "errDC = " << errDC << ";" << std::endl;
if (isTubeHit) {
bool isOK = false;
......@@ -858,12 +862,12 @@ for (; p_trk != trackCollection.end(); ++p_trk) {
fill("ShiftTRTTrackHistograms"+std::to_string(ibe), Residual_B_Ar_20GeV);
}
TimeResidual_B_Ar = timeresidual;
fill("ShiftTRTTrackHistograms"+std::to_string(ibe), TimeResidual_B_Ar);
if (errDC < 1) fill("ShiftTRTTrackHistograms"+std::to_string(ibe), TimeResidual_B_Ar);
} else {
Residual_B = loc - locR;
fill("ShiftTRTTrackHistograms"+std::to_string(ibe), Residual_B);
TimeResidual_B = timeresidual;
fill("ShiftTRTTrackHistograms"+std::to_string(ibe), TimeResidual_B);
if (errDC < 1) fill("ShiftTRTTrackHistograms"+std::to_string(ibe), TimeResidual_B);
if (cnst_is_pT_over_20GeV) {
Residual_B_20GeV = loc - locR;
......@@ -882,7 +886,7 @@ for (; p_trk != trackCollection.end(); ++p_trk) {
Residual_E_Ar = loc - locR;
fill("ShiftTRTTrackHistograms"+std::to_string(ibe)+std::to_string(iside), Residual_E_Ar);
TimeResidual_E_Ar = timeresidual;
fill("ShiftTRTTrackHistograms"+std::to_string(ibe)+std::to_string(iside), TimeResidual_E_Ar);
if (errDC < 1) fill("ShiftTRTTrackHistograms"+std::to_string(ibe)+std::to_string(iside), TimeResidual_E_Ar);
if (cnst_is_pT_over_20GeV) {
Residual_E_Ar_20GeV = loc - locR;
......@@ -892,7 +896,7 @@ for (; p_trk != trackCollection.end(); ++p_trk) {
Residual_E = loc - locR;
fill("ShiftTRTTrackHistograms"+std::to_string(ibe)+std::to_string(iside), Residual_E);
TimeResidual_E = timeresidual;
fill("ShiftTRTTrackHistograms"+std::to_string(ibe)+std::to_string(iside), TimeResidual_E);
if (errDC < 1) fill("ShiftTRTTrackHistograms"+std::to_string(ibe)+std::to_string(iside), TimeResidual_E);
if (cnst_is_pT_over_20GeV) {
Residual_E_20GeV = loc - locR;
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment