Aligner.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_MFT_ALIGNER_H
#define ALICEO2_MFT_ALIGNER_H
 
#include <array>
#include <vector>
 
#include <Rtypes.h>
#include <TString.h>
#include <TFile.h>
 
#include "ITSMFTReconstruction/ChipMappingMFT.h"
#include "ForwardAlign/MillePede2.h"
 
namespace o2
{
namespace mft
{
 
class Aligner
{
 public:
  Aligner();
 
  virtual ~Aligner();
 
  virtual void init() = 0;
 
  // simple setters
 
  void setChi2CutNStdDev(const Int_t value) { mChi2CutNStdDev = value; }
  void setResidualCutInitial(const Double_t value) { mResCutInitial = value; }
  void setResidualCut(const Double_t value) { mResCut = value; }
  void setAllowedVariationDeltaX(const double value) { mAllowVar[0] = value; }
  void setAllowedVariationDeltaY(const double value) { mAllowVar[1] = value; }
  void setAllowedVariationDeltaZ(const double value) { mAllowVar[3] = value; }
  void setAllowedVariationDeltaRz(const double value) { mAllowVar[2] = value; }
  void setChi2CutFactor(const double value) { mStartFac = value; }
 
  int getNDofPerSensor() const { return mNDofPerSensor; }
 
 protected:
  static constexpr int mNumberOfTrackParam = 4;                                  
  static constexpr int mNDofPerSensor = 4;                                       
  static o2::itsmft::ChipMappingMFT mChipMapping;                                
  static constexpr int mNumberOfSensors = mChipMapping.getNChips();              
  static constexpr int mNumberOfGlobalParam = mNDofPerSensor * mNumberOfSensors; 
  std::array<double, mNDofPerSensor> mAllowVar;                                  
  double mStartFac;                                                              
  int mChi2CutNStdDev;                                                           
  double mResCutInitial;                                                         
  double mResCut;                                                                
  TString mMilleRecordsFileName;                                                 
  TString mMilleConstraintsRecFileName;                                          
  bool mIsInitDone = false;                                                      
  std::vector<int> mGlobalParameterStatus;                                       
 
  // used to fix some degrees of freedom
 
  static constexpr int mFixedParId = -1;
  static constexpr int mFreeParId = mFixedParId - 1;
 
  ClassDef(Aligner, 0);
};
 
} // namespace mft
} // namespace o2
 
#endif