EmcBhaCalib.h

Go to the documentation of this file.

//--------------------------------------------------------------------------
// Environment:
//      This software was developed for the BESIII collaboration.  If you
//      use all or part of it, please give an appropriate acknowledgement.
//
// Copyright Information:
//  Copyright (C) 2005              IHEP
//
//------------------------------------------------------------------------
 
#ifndef EMCBHACALIB_H
#define EMCBHACALIB_H
 
//-------------
// C Headers --
//-------------
#include <ctime>
 
#include <cstring>
 
#include "GaudiKernel/Algorithm.h"
 
#include "EmcBhaCalibData.h"
#include "EmcLSSMatrix.h"
 
#include "EmcCalibConstSvc/IEmcCalibConstSvc.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/ISvcLocator.h"
 
// Ntuple
#include "GaudiKernel/NTuple.h"
 
using namespace std;
 
//------------------------------------
// Collaborating Class Declarations --
//------------------------------------
 
//      ---------------------
//      -- Class Interface --
//      ---------------------
/**
 *  class  EmcBhaCalib Algorithm - performs calibration of EMC Xtals with Bhabha
 *  events and Chi^2 fit, the resulting system of linear equations
 *  of the fit is solved with the SLAP Algebra package
 *
 *  The matrix and the vector can either be fetched from file(s).
 *  If they are read from files the file names are supposed to have the name
 *  CalibVector.dat and calibMatrix.dat with an extension at the beginning.
 *  Only these extensions (usually the runnumbers) need to be listed in the read-in-file
 *  that is specified with runNumberFile.
 *
 *  @author Chunxiu Liu        (originator/contributor etc.);
 */
 
class EmcBhaCalib : public Algorithm {
 
public:
  EmcBhaCalib( const std::string& name, ISvcLocator* pSvcLocator );
 
  //--------------
  // Destructor
  //--------------
  virtual ~EmcBhaCalib();
 
  StatusCode initialize();
  StatusCode execute();
  StatusCode finalize();
 
  virtual void help();
 
  //------------------
  // Helper functions
  //------------------
 
private:
  // read in calibration constants from file or initialize them to 1
  void initCalibConst();
 
  // solve system of linear equations with SLAP package
  bool solveEqua();
 
  // output constants to file
  void writeOutConst();
 
  // output matrix and vector to file
  void writeOut();
 
  // read in matrix and vector from file
  bool readInFromFile();
 
  // read in matrix and vector from the database
  bool readInFromDB();
 
  // make the absolute constants and fit constants vs theta to get values
  // where we did not find some
  bool prepareConstants();
 
  // fill nutuple
  void ntupleOut();
 
  // print data information and only used in testing
  void printInfo( std::string xtalHitsDirFile );
  //
 
  void digiConstCor();
  //--------------
  //  Data
  //--------------
 
private:
  double m_peakCor[6240];
 
  // cut on minimum number of direct hits needed per xtal
  int m_dirHitsCut;
 
  // convergence criterion
  double m_convCrit;
 
  // ask before doing the matrix inversion
  bool m_askForMatrixInversion;
 
  // fit a polynom second order to the constants theta distribution
  // to get constants where we did not find one
  bool m_fitPolynom; // now not been used
 
  // write the whole matrix, the vector and the constants to file ?
  bool m_writeToFile;
 
  // read calibration data (matrix) from the database (or from files)
  bool m_readDataFromDB;
 
  // which method for solving the matrix equation
  std::string m_equationSolver;
 
  // file name extensions for different output files
  std::string m_fileExt;
 
  // file directory
  std::string m_fileDir;
 
  // collected calibration data
  EmcBhaCalibData* m_myCalibData;
 
  // was calibration successful ?
  bool m_calibrationSuccessful;
 
  // array of relative calibration constants, will be reduced to only
  // living xtals for the matrix inversion (no zero hit entries)
  double* m_calibConst;
 
  // relative calibration constants, full list including zero entries
  double* m_calibConstUnred;
 
  // absolute calibration constants (new relative * old constants)
  double* m_absoluteConstants;
 
  // the old absolute constants the new once are relative to
  double* m_oldConstants;
 
  // number of non zeros in matrix M
  long int m_nrNonZeros;
 
  // number of xtals with minimum required direct hits
  int m_nrXtalsEnoughHits;
 
  // file with a list of runnumbers to read in for which runs matrix file
  //  to load (Default: ./runnumber.dat)
  std::string m_runNumberFile;
 
  // message output in the tool classes
  int m_MsgFlag;
 
  // ntuple for analysis
  NTuple::Tuple*      m_tuple1;
  NTuple::Item<long>  ixtal;
  NTuple::Item<float> gi0;
  NTuple::Item<float> calibConst;
  NTuple::Item<float> err;
  NTuple::Item<long>  nhitxtal;
 
  double DigiCalibConst;
 
  IEmcCalibConstSvc* m_emcCalibConstSvc;
};
 
#endif // EMCBHACALIB_H