Skip to content
Snippets Groups Projects
Commit 40221904 authored by Christopher Rob Jones's avatar Christopher Rob Jones
Browse files

Tighten sanity check a little in RichCKResolutionFitter.h

parent 7db362db
No related branches found
No related tags found
No related merge requests found
Pipeline #7415372 passed
......@@ -159,11 +159,13 @@ namespace Rich::Rec {
public: // Fit Quality options
/// Maximum abs fitted mean for OK fit
std::array<double, NRads> RichMaxMean{{0.006, 0.003}};
std::array<double, NRads> RichMaxMean{{0.0040, 0.0020}};
/// Minimum abs fitted sigma for OK fit
std::array<double, NRads> RichMinSigma{{0.0003, 0.0001}};
/// Maximum abs fitted sigma for OK fit
std::array<double, NRads> RichMaxSigma{{0.010, 0.005}};
std::array<double, NRads> RichMaxSigma{{0.0030, 0.0020}};
/// Maximum shift error for OK fit
std::array<double, NRads> MaxShiftError{{0.001, 0.001}};
std::array<double, NRads> MaxShiftError{{0.0010, 0.0010}};
/// Maximum resolution error for OK fit
std::array<double, NRads> MaxSigmaError{{0.0002, 0.0001}};
};
......@@ -245,8 +247,11 @@ namespace Rich::Rec {
fitOK &= fFitF->GetParError( 2 ) < params().MaxSigmaError[rIn];
// Final sanity check on the fit results to detect complete failures
// such as when there is no Cherenkov signal at all
fitOK &= fabs( fFitF->GetParameter( 1 ) ) < params().RichMaxMean[rIn];
fitOK &= fabs( fFitF->GetParameter( 2 ) ) < params().RichMaxSigma[rIn];
const auto ckMean = fFitF->GetParameter( 1 );
const auto ckRes = fFitF->GetParameter( 2 );
fitOK &= fabs( ckMean ) < params().RichMaxMean[rIn];
fitOK &= fabs( ckRes ) > params().RichMinSigma[rIn];
fitOK &= fabs( ckRes ) < params().RichMaxSigma[rIn];
// return the final status
return fitOK;
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment