VDriftCorrFact.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#ifndef AliceO2_TPC_VDRIFT_CORRFACT_H
#define AliceO2_TPC_VDRIFT_CORRFACT_H
 
#include "GPUCommonRtypes.h"
#include "CommonConstants/LHCConstants.h"
 
namespace o2::tpc
{
 
struct VDriftCorrFact {
  long firstTime{};         
  long lastTime{};          
  long creationTime{};      
  float corrFact{1.0};      
  float corrFactErr{0.0};   
  float refVDrift{0.};      
  float refTimeOffset{0.};  
  float timeOffsetCorr{0.}; 
  float refTP{0.};          
 
  float getVDrift() const { return refVDrift * corrFact; }
  float getVDriftError() const { return refVDrift * corrFactErr; }
 
  float getTimeOffset() const { return refTimeOffset + timeOffsetCorr; }
 
  // renormalize VDrift reference and correction either to provided new reference (if >0) or to correction 1 wrt current reference
  void normalize(float newVRef = 0.f, float tp = 0.f)
  {
    float normVDrift = newVRef;
    if (newVRef == 0.f) {
      normVDrift = refVDrift * corrFact;
      newVRef = normVDrift;
      if ((tp > 0) && (refTP > 0)) {
        // linear scaling based on relative change of T/P
        normVDrift *= refTP / tp;
        refTP = tp; // update reference T/P
      }
    }
    float fact = refVDrift / normVDrift;
    refVDrift = newVRef;
    corrFactErr *= fact;
    corrFact *= fact;
  }
 
  // similarly, the time offset reference is set to provided newRefTimeOffset (if > -998) or modified to have timeOffsetCorr to
  // be 0 otherwise
 
  void normalizeOffset(float newRefTimeOffset = -999.)
  {
    if (newRefTimeOffset > -999.) {
      timeOffsetCorr = getTimeOffset() - newRefTimeOffset;
      refTimeOffset = newRefTimeOffset;
    } else {
      refTimeOffset = getTimeOffset();
      timeOffsetCorr = 0.;
    }
  }
 
  ClassDefNV(VDriftCorrFact, 3);
};
 
} // namespace o2::tpc
#endif