PHOSRunbyrunCalibrator.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_CALIBRATION_PHOSRUNBYRUN_CALIBRATOR_H
#define O2_CALIBRATION_PHOSRUNBYRUN_CALIBRATOR_H
 
 
#include "Framework/Task.h"
#include "Framework/ProcessingContext.h"
#include "DataFormatsPHOS/BadChannelsMap.h"
#include "DataFormatsPHOS/TriggerRecord.h"
#include "DetectorsCalibration/TimeSlotCalibration.h"
#include <boost/histogram.hpp>
#include "TH1.h"
#include "PHOSCalibWorkflow/RingBuffer.h"
#include "DataFormatsPHOS/Cluster.h"
#include "PHOSBase/Geometry.h"
 
using namespace o2::framework;
 
namespace o2
{
namespace phos
{
 
class PHOSRunbyrunSlot
{
 public:
  using boostHisto = boost::histogram::histogram<std::tuple<boost::histogram::axis::regular<double, boost::use_default, boost::use_default, boost::use_default>>, boost::histogram::unlimited_storage<std::allocator<char>>>;
 
  PHOSRunbyrunSlot(bool useCCDB, std::string path);
  PHOSRunbyrunSlot(const PHOSRunbyrunSlot& other);
 
  ~PHOSRunbyrunSlot() = default;
 
  void print() const;
  void fill(const gsl::span<const Cluster>& clusters, const gsl::span<const TriggerRecord>& trs);
  void fill(const gsl::span<const Cluster>& /*clusters*/){}; // not used
  void merge(const PHOSRunbyrunSlot* prev);
  void clear();
 
  boostHisto& getCollectedHistos(int m) { return mReMi[m]; }
 
  void setRunStartTime(long tf) { mRunStartTime = tf; }
 
 private:
  bool checkCluster(const Cluster& clu);
 
 private:
  bool mUseCCDB = false;
  long mRunStartTime = 0;                             
  float mPtCut = 1.5;                                 
  std::string mCCDBPath{"http://alice-ccdb.cern.ch"}; 
  std::array<boostHisto, 8> mReMi;                    
  std::unique_ptr<RingBuffer> mBuffer;                
  std::unique_ptr<BadChannelsMap> mBadMap;            
 
  ClassDefNV(PHOSRunbyrunSlot, 1);
};
 
//==========================================================================================
class PHOSRunbyrunCalibrator final : public o2::calibration::TimeSlotCalibration<o2::phos::PHOSRunbyrunSlot>
{
  using Slot = o2::calibration::TimeSlot<o2::phos::PHOSRunbyrunSlot>;
 
 public:
  PHOSRunbyrunCalibrator();
  ~PHOSRunbyrunCalibrator() final;
 
  bool hasEnoughData(const Slot& slot) const final;
  void initOutput() final;
  void finalizeSlot(Slot& slot) final;
  Slot& emplaceNewSlot(bool front, TFType tstart, TFType tend) final;
  bool process(TFType tf, const gsl::span<const Cluster>& clu, const gsl::span<const TriggerRecord>& trs);
 
  std::array<float, 8> getCalibration() { return mRunByRun; }
  void writeHistos();
 
  // Functions used in histo fittings
  double CBRatio(double* x, double* p);
  double CBSignal(double* x, double* p);
  double bg(double* x, double* p);
 
 private:
  void scanClusters(o2::framework::ProcessingContext& pc);
  bool checkCluster(const Cluster& clu);
 
 private:
  bool mUseCCDB = false;
  long mRunStartTime = 0;                             
  std::string mCCDBPath{"http://alice-ccdb.cern.ch"}; 
  std::array<float, 8> mRunByRun;                     
  std::array<TH1F*, 8> mReMi;                         
 
  ClassDefOverride(PHOSRunbyrunCalibrator, 1);
};
 
} // namespace phos
} // namespace o2
 
#endif