T0Fit.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 O2_TRD_T0FIT_H
#define O2_TRD_T0FIT_H
 
#include "DataFormatsTRD/T0FitHistos.h"
#include "DetectorsCalibration/TimeSlotCalibration.h"
#include "DetectorsCalibration/TimeSlot.h"
#include "DataFormatsTRD/Constants.h"
#include "CCDB/CcdbObjectInfo.h"
#include "DataFormatsTRD/CalT0.h"
#include "TRDCalibration/CalibrationParams.h"
 
#include "Rtypes.h"
#include "TProfile.h"
#include "TF1.h"
#include "Fit/Fitter.h"
#include "TFile.h"
#include "TTree.h"
 
#include <array>
#include <cstdlib>
#include <memory>
 
namespace o2
{
namespace trd
{
//______________________________________________________________________________________________
struct ErfLandauChi2Functor {
  double operator()(const double* par) const;
  std::vector<float> x; 
  std::vector<float> y; 
  float lowerBoundFit;  
  float upperBoundFit;  
};
 
//______________________________________________________________________________________________
class T0Fit final : public o2::calibration::TimeSlotCalibration<o2::trd::T0FitHistos>
{
  using Slot = o2::calibration::TimeSlot<o2::trd::T0FitHistos>;
 
 public:
  T0Fit() = default;
  ~T0Fit() final = default;
 
  bool hasEnoughData(const Slot& slot) const final { return slot.getContainer()->getNEntries() >= mParams.minEntriesTotalT0Fit; }
  void initOutput() final;
  void finalizeSlot(Slot& slot) final;
  Slot& emplaceNewSlot(bool front, TFType tStart, TFType tEnd) final;
 
  void createOutputFile();
 
  void closeOutputFile();
 
  const std::vector<o2::trd::CalT0>& getCcdbObjectVector() const { return mObjectVector; }
  std::vector<o2::ccdb::CcdbObjectInfo>& getCcdbObjectInfoVector() { return mInfoVector; }
 
  void initProcessing();
 
 private:
  bool mInitDone{false};                                         
  const TRDCalibParams& mParams{TRDCalibParams::Instance()};     
  bool mEnableOutput{false};                                     
  std::unique_ptr<TFile> mOutFile{nullptr};                      
  std::unique_ptr<TTree> mOutTree{nullptr};                      
  ErfLandauChi2Functor mFitFunctor;                              
  ROOT::Fit::Fitter mFitter;                                     
  std::array<double, 4> mParamsStart;                            
  std::unique_ptr<TF1> mFuncErfLandau;                           
  float mDummyT0{-5};                                            
  std::array<float, o2::trd::constants::MAXCHAMBER> t0_chambers; 
  float t0_average{-5};                                          
 
  std::vector<o2::ccdb::CcdbObjectInfo> mInfoVector; 
  std::vector<o2::trd::CalT0> mObjectVector;         
 
  std::unique_ptr<TProfile> adcProfIncl;                                            
  std::array<std::unique_ptr<TProfile>, o2::trd::constants::MAXCHAMBER> adcProfDet; 
 
  ClassDefNV(T0Fit, 1);
};
 
} // namespace trd
} // namespace o2
 
#endif // O2_TRD_T0FIT_H