TGainMatch.h

Go to the documentation of this file.

#ifndef TGAINMATCH_H
#define TGAINMATCH_H
 
/** \addtogroup Calibration
 *  @{
 */
 
#include "TCal.h"
#include "TCalManager.h"
#include "TPeak.h"
#include "TSpectrum.h"
#include "TH2.h"
#include "TF1.h"
 
class TGainMatch : public TCal {
public:
   TGainMatch() : fCoarseRange(fDefaultCoarseRange) {}
   TGainMatch(const char* name, const char* title) : TCal(name, title), fCoarseRange(fDefaultCoarseRange) { Clear(); }
   TGainMatch(const TGainMatch&);
   TGainMatch(TGainMatch&&) noexcept            = default;
   TGainMatch& operator=(const TGainMatch&)     = default;
   TGainMatch& operator=(TGainMatch&&) noexcept = default;
   ~TGainMatch()                                = default;
 
   void Copy(TObject& obj) const override;
 
   void CalculateGain(Double_t cent1, Double_t cent2, Double_t eng1, Double_t eng2);
 
   static Bool_t CoarseMatchAll(TCalManager* cm, TH2* mat, Double_t energy1 = 1173.228, Double_t energy2 = 1332.492);
   static Bool_t FineMatchFastAll(TCalManager* cm, TH2* mat, Double_t energy1, Double_t energy2);
   static Bool_t FineMatchFastAll(TCalManager* cm, TH2* mat1, Double_t energy1, TH2* mat2, Double_t energy2);
   static Bool_t FineMatchFastAll(TCalManager* cm, TH2* mat1, TPeak* peak1, TH2* mat2, TPeak* peak2);
   static Bool_t FineMatchFastAll(TCalManager* cm, TH2* mat, TPeak* peak1, TPeak* peak2);
 
   static Bool_t FineMatchAll(TCalManager* cm, TH2* charge_mat, TH2* eng_mat, Int_t testchan, Double_t energy1,
                              Double_t energy2, Int_t low_range = 100, Int_t high_range = 600);
 
   static Bool_t AlignAll(TCalManager* cm, TH1* hist, TH2* mat, Int_t low_range = 100, Int_t high_range = 600);
 
   Bool_t CoarseMatch(TH1* hist, Int_t chanNum = 9999, Double_t energy1 = 1173.228, Double_t energy2 = 1332.492);
   Bool_t FineMatchFast(TH1* hist1, TPeak* peak1, TH1* hist2, TPeak* peak2, Int_t channelNum = 9999);
   Bool_t FineMatchFast(TH1* hist, TPeak* peak1, TPeak* peak2, Int_t channelNum = 9999);
   Bool_t FineMatchFast(TH1* hist, Double_t energy1, Double_t energy2, Int_t channelNum = 9999);
   Bool_t FineMatchFast(TH1* hist1, Double_t energy1, TH1* hist2, Double_t energy2, Int_t channelNum = 9999);
 
   Bool_t FineMatch(TH1* energyHist, TH1* testhist, TH1* chargeHist, Double_t energy1, Double_t energy2,
                    Int_t low_range = 100, Int_t high_range = 600, Int_t channelNum = 9999);
 
   Bool_t Align(TH1* test, TH1* hist, Int_t low_range = 100, Int_t high_range = 600);
 
   void Clear(Option_t* opt = "") override;
   void Print(Option_t* opt = "") const override;
 
   Bool_t IsGroupable() const override { return false; }
   void   WriteToChannel() const override;
 
   void      SetNucleus(TNucleus*, Option_t* = "") override { Warning("SetNucleus", "Is not used in TGainMatching"); }
   TNucleus* GetNucleus() const override
   {
      Warning("GetNucleus", "Is not used in TGainMatching");
      return nullptr;
   }
 
   void SetHist(TH1*) override { Warning("SetHist", "Is not used in TGainMatching"); }
   TH1* GetHist() const
   {
      Warning("GetHist", "Is not used in TGainMatching");
      return nullptr;
   }
 
   void        SetCoarseRange(Double_t coarseRange) { fCoarseRange = coarseRange; }
   Double_t    GetCoarseRange() const { return fCoarseRange; }
   static void SetDefaultCoarseRange(Double_t coarseRange)
   {
      std::cout << "All new TGainMatch objects will have their range set to " << coarseRange << std::endl;
      fDefaultCoarseRange = coarseRange;
   }
   static Double_t GetDefaultCoarseRange() { return fDefaultCoarseRange; }
 
private:
   Bool_t                  fCoarseMatch{false};
   Bool_t                  fAligned{false};
   TH1*                    fHist{nullptr};
   std::array<Double_t, 2> fAlignCoeffs{0.};
   std::array<Double_t, 2> fGainCoeffs{0.};
   Double_t                fCoarseRange{0.};
   Double_t                HistCompare(Double_t* x, Double_t* par);
 
   static Double_t fDefaultCoarseRange;
 
   /// \cond CLASSIMP
   ClassDefOverride(TGainMatch, 1)   // NOLINT(readability-else-after-return)
   /// \endcond
};
/*! @} */
#endif