eclGammaGammaEAlgorithm.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2018 - Belle II Collaboration                             *
 *                                                                        *
 * Analyze histograms of amplitudes for each ECL crystal from gamma pair  *
 * events. Code can either find most-likely energy deposit for each       *
 * crystal (MC) or calibration constant for each crystal (data)           *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Christopher Hearty                                       *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
#include <ecl/calibration/eclGammaGammaEAlgorithm.h>
#include <calibration/CalibrationAlgorithm.h>
 
namespace Belle2 {
  namespace ECL {
 
    class eclGammaGammaEAlgorithm : public CalibrationAlgorithm {
    public:
 
      eclGammaGammaEAlgorithm();
 
      virtual ~eclGammaGammaEAlgorithm() {}
 
      void setOutputName(const std::string& outputName) {m_outputName = outputName;}
 
      std::string getOutputName() {return m_outputName;}
 
      void setCellIDLo(int cellIDLo) {m_cellIDLo = cellIDLo;}
 
      int getCellIDLo() {return m_cellIDLo;}
 
      void setCellIDHi(int cellIDHi) {m_cellIDHi = cellIDHi;}
 
      int getCellIDHi() {return m_cellIDHi;}
 
      void setMinEntries(int minEntries) {m_minEntries = minEntries;}
 
      int getMinEntries() {return m_minEntries;}
 
      void setMaxIterations(int maxIterations) {m_maxIterations = maxIterations;}
 
      int getMaxIterations() {return m_maxIterations;}
 
      void setTRatioMin(double tRatioMin) {m_tRatioMin = tRatioMin;}
 
      double getTRatioMin() {return m_tRatioMin;}
 
      void setTRatioMax(double tRatioMax) {m_tRatioMax = tRatioMax;}
 
      double getTRatioMax() {return m_tRatioMax;}
 
      void setUpperEdgeThresh(double upperEdgeThresh) {m_upperEdgeThresh = upperEdgeThresh;}
 
      double getUpperEdgeThresh() {return m_upperEdgeThresh;}
 
      void setPerformFits(bool performFits) {m_performFits = performFits;}
 
      bool getPerformFits() {return m_performFits;}
 
      void setFindExpValues(bool findExpValues) {m_findExpValues = findExpValues;}
 
      bool getFindExpValues() {return m_findExpValues;}
 
      void setStoreConst(int storeConst) {m_storeConst = storeConst;}
 
      int getStoreConst() {return m_storeConst;}
 
 
    protected:
 
      virtual EResult calibrate() override;
 
    private:
 
      std::string m_outputName = "eclGammaGammaEAlgorithm.root"; 
      int m_cellIDLo = 1;  
      int m_cellIDHi = 8736;  
      int m_minEntries = 150;  
      int m_maxIterations = 10; 
      double m_tRatioMin = 0.45;  
      double m_tRatioMax = 0.60; 
      double m_upperEdgeThresh = 0.02; 
      bool m_performFits = true;  
      bool m_findExpValues =
        false;  
      int m_storeConst = 0; 
      int fitOK = 16; 
      int iterations = 8; 
      int atLimit = 4; 
      int poorFit = 3; 
      int noPeak = 2; 
      int notFit = -1; 
    };
  }
} // namespace Belle2