XTCalibration.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
#include "string"
#include "TH2D.h"
#include "TH1D.h"
#include "TF1.h"
namespace Belle2 {
  namespace CDC {
    class XTCalibration {
    public:
      XTCalibration();
      virtual ~XTCalibration() {}
      virtual void BField(bool bfield) {m_BField = bfield;}
      virtual void setDebug(bool debug = false) {m_debug = debug; }
      virtual void setUseDB(bool useDB = false) {m_useDB = useDB; }
      virtual void setMinimumNDF(double minndf) {m_ndfmin = minndf;}
      virtual void setMinimumPval(double minPval) {m_Pvalmin = minPval;}
      virtual void inputFileNames(std::string inputname) {m_inputRootFileNames.assign(inputname);}
      virtual void profileFileNames(std::string profileFileName) {m_profileFileName.assign(profileFileName);}
      virtual void useProfileFromInputXT(bool useProfileXTFromInputXT) {m_useProfileXTFromInputXT = useProfileXTFromInputXT;}
      virtual void setXTFileName(std::string name) {m_xtfile.assign(name);}
      virtual void setMode(unsigned short mode = 1) {m_xtmode = mode;}
      virtual void setStoreHisto(bool storeHist = false) {m_storeHisto = storeHist;}
      void setLRSeparate(bool lr = true) {m_LRseparate = lr;}
 
 
      void execute()
      {
        calibrate();
      }
 
    protected:
      virtual bool calibrate();
      virtual void readXTFromDB();
      virtual void readXTFromText();
      virtual void CreateHisto();
      virtual void readProfile();
      virtual void readXT();
      virtual void Write();
      virtual void storeHisto();
 
    private:
      double m_ndfmin = 5;    
      double m_Pvalmin = 0.;  
      bool m_debug = false;   
      bool m_storeHisto = false;  
      bool m_useDB = false;       
      bool m_useProfileXTFromInputXT = true; 
      bool m_LRseparate = true; 
      bool m_useSliceFit = false; 
      bool m_BField = true;  
      //double m_XT_fit[56][2][18][7][8];  /**< Fitted parameter*/
      double xtold[56][2][18][7][8];     
      int fitflag[56][2][20][10];         
      TF1* xtf5r[56][2][20][10];          
      TProfile* hprof[56][2][20][10];     
      TH2D* hist2d[56][2][20][10];        
      TH2D* hist2d_draw[56][20][10];       
      TH1D* hist2d_1[56][2][20][10];       
      /*********************************
      Fit Flag
      =-1: low statitic
       =1: good
      =0: Fit failure
      =2: Error Outer
      =3: Error Inner part;
      **********************************/
 
      std::string m_OutputXTFileName = "xt_new.dat"; 
      std::string m_inputRootFileNames = "rootfile/output*"; 
      std::string m_profileFileName = "xt_profile";      
      std::string m_xtfile = "cdc/data/xt.dat";          
      //int m_firstExperiment; /**< First experiment. */
      //int m_firstRun; /**< First run. */
      //int m_lastExperiment; /**< Last experiment */
      //int m_lastRun; /**< Last run. */
 
      //      double m_alpha[18][3];//alpha bin, 18bin, 0=low, 1 = up, 2 = alpha
      //double m_theta[10][3];//alpha bin, 18bin, 0=low, 1 = up, 2 = alpha
      //      unsigned short m_xtParamMode;    /*!< Mode for xt parameterization */
      //unsigned short m_nXtParams;      /*!< no. of xt parameters per bin */
      int m_nalpha; 
      int m_ntheta;
      double l_alpha[18];
      double u_alpha[18];
      double ialpha[18]; 
      double l_theta[7]; 
      double u_theta[7];
      double itheta[7]; 
      int nalpha_old; 
      int ntheta_old;
      double l_alpha_old[18];
      double u_alpha_old[18];
      double ialpha_old[18]; 
      double l_theta_old[7]; 
      double u_theta_old[7];
      double itheta_old[7]; 
      unsigned short xtmode_old; 
      int m_MAXalpha = 18;        
      int m_MAXtheta = 7;         
      unsigned short m_xtmode = 1; 
      int m_smallestEntryRequire = 1000; 
      double m_par6[56] = {89,   91,  94,  99,  104, 107, 110, 117,
                           126, 144, 150, 157, 170, 180,
                           160, 167, 183, 205, 200, 194,
                           177, 189, 192, 206, 224, 234,
                           193, 206, 209, 215, 222, 239,
                           204, 212, 217, 227, 235, 240,
                           215, 222, 230, 239, 246, 253,
                           227, 232, 239, 243, 253, 258,
                           231, 243, 246, 256, 263, 300
                          };
 
    };
  }
}