Add optional residuals from IB/OB tracks extrapolations to OB/IB clusters

Detectors/GlobalTrackingWorkflow/study/include/GlobalTrackingStudy/CheckResidConfig.h

@@ -28,6 +28,7 @@ struct CheckResidConfig : o2::conf::ConfigurableParamHelper<CheckResidConfig> {
   bool addPVAsCluster = true;
   bool useStableRef = true;
   bool doIBOB = true;
+  bool doIBOBresiduals = true; ///< fill residuals of OB(IB) points wrt extrapolation of the IB(OB) sub-track, w/o intermediate updates
   bool doResid = true;
   bool useITSGloContributors = false;
 
@@ -48,6 +49,7 @@ struct CheckResidConfig : o2::conf::ConfigurableParamHelper<CheckResidConfig> {
   float zranges[8] = {10.f, 15.f, 15.f, 15.f, 40.f, 40.f, 74.f, 74.f};
   float maxDYZ[8] = {0.03, 0.015, 0.01, 0.01, 0.08, 0.08, 0.12, 0.1};
   float maxDPar[5] = {0.15, 0.15, 0.015, 0.015, 1.};
+  float IOBScaleFactor = 5.f; ///< multiplier applied to maxDYZ, resMMLrY, resMMLrZ when booking/drawing the (larger range) IOB residuals
   float maxHPVXY = 0.5;
   float maxHPVZ = 20.;
   float maxHPVN = 200.;

Detectors/GlobalTrackingWorkflow/study/include/GlobalTrackingStudy/CheckResidTypes.h

@@ -21,13 +21,17 @@ struct Point {
   float dz = 0.f;
   float sig2y = 0.f;
   float sig2z = 0.f;
+  float dyIOB = 0.f; ///< residual wrt extrapolation of the IB (for OB points) or OB (for IB points) sub-track, no intermediate updates
+  float dzIOB = 0.f;
+  float sig2yIOB = 0.f; ///< <= 0 means IOB residual was not computed for this point
+  float sig2zIOB = 0.f;
   float phi = 0.f;
   float z = 0.f;
   float snp = 0.f;
   int16_t clSize = 0;
   int16_t sens = -1;
   int8_t lr = -1; // -1 = vtx
-  ClassDefNV(Point, 2)
+  ClassDefNV(Point, 3)
 };
 
 struct Track {