EmcBhaCalib Class Reference

#include <EmcBhaCalib.h>

Inheritance diagram for EmcBhaCalib:

Public Member Functions
	EmcBhaCalib (const std::string &name, ISvcLocator *pSvcLocator)
virtual	~EmcBhaCalib ()
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
virtual void	help ()

Detailed Description

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.);

Definition at line 58 of file EmcBhaCalib.h.

Constructor & Destructor Documentation

EmcBhaCalib()

EmcBhaCalib::EmcBhaCalib	(	const std::string &	name,
		ISvcLocator *	pSvcLocator )

Definition at line 69 of file EmcBhaCalib.cxx.

  :Algorithm(name, pSvcLocator),  
   m_dirHitsCut(200),
   m_convCrit(0.000001),
   m_askForMatrixInversion(true),
   m_fitPolynom(true),  //no used now
   m_writeToFile(false),
   m_readDataFromDB(false),
   m_equationSolver("GMR"),
   m_fileExt(""),
   m_fileDir("/home/besdata/public/liucx/Calib/"),
   m_nrNonZeros(0), 
   m_nrXtalsEnoughHits(0),
   m_runNumberFile("runnumbers.dat"),
   m_MsgFlag(0)
 
{
 
  // Declare the properties  
  declareProperty("equationSolver", m_equationSolver);
  declareProperty("dirHitsCut", m_dirHitsCut);
  declareProperty("writeToFile", m_writeToFile);
  declareProperty("fileExt", m_fileExt);
  declareProperty("fileDir", m_fileDir);
  declareProperty("readDataFromDB", m_readDataFromDB);
  declareProperty("askForMatrixInversion", m_askForMatrixInversion);
  declareProperty("fitPolynom", m_fitPolynom);
  declareProperty("convCrit",  m_convCrit);
  declareProperty("runNumberFile", m_runNumberFile);
  declareProperty("MsgFlag", m_MsgFlag);
 
  m_calibConst = new double[6240];
  m_calibConstUnred = new double[6240];
  m_absoluteConstants = new double[6240];
  m_oldConstants = new double[6240];
 
 
  //ntuple
  m_tuple1=0;
 
 
 
}

~EmcBhaCalib()

EmcBhaCalib::~EmcBhaCalib ( )

virtual

Definition at line 116 of file EmcBhaCalib.cxx.

                          {
  if ( 0 != m_calibConst) {
    delete [] m_calibConst;
    m_calibConst = 0;
  }
  if ( 0 != m_calibConstUnred) {
    delete [] m_calibConstUnred;
    m_calibConstUnred = 0;
  }
  if ( 0 != m_absoluteConstants) {
    delete [] m_absoluteConstants;
    m_absoluteConstants = 0;
  }
  if ( 0 != m_oldConstants) {
    delete [] m_oldConstants;
    m_oldConstants = 0;
  }
  if ( 0 != m_myCalibData) {
    delete m_myCalibData;
    m_myCalibData = 0;
  }
 
}

Member Function Documentation

execute()

StatusCode EmcBhaCalib::execute

(

)

Definition at line 195 of file EmcBhaCalib.cxx.

                                {
  
  MsgStream log(msgSvc(), name());
  log << MSG::DEBUG << "in execute()" << endreq;   
 //read in accumulated matrix/matrices
  if ( m_readDataFromDB ) {
    m_calibrationSuccessful = readInFromDB();
 
    log << MSG::INFO << "read in accumulated matrix from DB"<<endreq;
 
  } else {
    m_calibrationSuccessful = readInFromFile();
 
    log << MSG::INFO <<"read in accumulated matrix from file"<<endreq; 
 
  }
  
  if ( m_calibrationSuccessful == true ) {
    
    //ask first if to do a matrix inversion
    if (m_askForMatrixInversion==true) {
 
      m_calibrationSuccessful = false;
      log << MSG::INFO << " Well, we have " 
      << m_nrXtalsEnoughHits << " crystals whith more "
      << "than " <<  m_dirHitsCut 
      << " direct hits. Do you want do "
      << "a matrix inversion ? [N]" << endreq; 
 
 
      m_calibrationSuccessful = true;
 
    }
   
    if ( m_calibrationSuccessful == true ) {
      
      //write added matrix and vector to files
      if ( m_writeToFile == true) { 
        writeOut();
      }
      //cout <<"writeout"<<endl;
    m_myCalibData->printVec(2);
 
      //reduce to continious array of only non zeros elements for SLAP
      if ( m_calibrationSuccessful = m_myCalibData->reduce() ) {
 
    cout<<"m_calibrationSuccessful="<<m_calibrationSuccessful<<endl;
 
    if ( m_myCalibData->nXtalsHit() != m_myCalibData->nXtals() ){
      log << MSG::INFO << " Reduced number of Xtals for calibration: " 
          << m_myCalibData->nXtalsHit() 
          << endreq;
    }
    cout<<"m_myCalibData->nXtalsHit()="<<m_myCalibData->nXtalsHit()
        <<"m_myCalibData->nXtals()="<<m_myCalibData->nXtals()<<endl;
    m_myCalibData->printVec(10);
 
    //solve matrix equation
    if ( m_calibrationSuccessful = solveEqua() ) {
      
      for (int i=0; i<m_myCalibData->nXtalsHit(); i++){
        //fill an array of constants for storage in database
        m_calibConstUnred[m_myCalibData->xtalInd(i)] 
          = m_calibConst[i];
        
        // if (m_myCalibData->xtalHitsDir(i)>10)
          // cout<<"Index,calibconst="<<"  "<<i <<"  "<<m_myCalibData->xtalInd(i)
          //  <<"  "<<m_calibConstUnred[m_myCalibData->xtalInd(i)]
          //  <<"  "<<m_calibConstUnred[m_myCalibData->xtalInd(i)]*
          //  m_oldConstants[m_myCalibData->xtalInd(i)]<<endl;
        
      }
      //do validation, fit polynom if necessary, fill CalList     
      prepareConstants();
 
      //if wanted write constants to plain ASCII file 
      if ( m_writeToFile==true){ 
        writeOutConst();
      }
 
      ntupleOut();
 
    }
    
      } else {
    log << MSG::WARNING << " Reduction of the Emc Bhabha calibration matrix FAILED !"
        << endl
        << " Will NOT perform Emc Bhabha calibration !"
        << endreq;
    return( StatusCode::FAILURE);
      }
    }
  } else {
    log << MSG::WARNING << " ERROR in reading in Emc Bhabha calibration data !"
    << endl
    << " Will NOT perform Emc Bhabha calibration !"
    << endreq;
    return( StatusCode::FAILURE);
  }
 
  return StatusCode::SUCCESS;
}

finalize()

StatusCode EmcBhaCalib::finalize

(

)

Definition at line 299 of file EmcBhaCalib.cxx.

                                 {
 
  MsgStream log(msgSvc(), name());
 
 
  log << MSG::INFO << "in endRun()" << endreq;
 
 
  return StatusCode::SUCCESS;
}

help()

void EmcBhaCalib::help ( )

virtual

Definition at line 312 of file EmcBhaCalib.cxx.

                  {
  
  MsgStream log(msgSvc(), name());
  
  log << MSG::INFO<< "Performs the Chi square fit of the accumulated " 
      << endreq;
}

initialize()

StatusCode EmcBhaCalib::initialize

(

)

Definition at line 141 of file EmcBhaCalib.cxx.

                                  {
  
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in initialize()" << endreq;
 
  m_myCalibData = new EmcBhaCalibData(6240, m_MsgFlag);
 
  m_calibrationSuccessful = false;  
  
  StatusCode status1;
  
  NTuplePtr nt1(ntupleSvc(),"FILE302/n1");
  if ( nt1 ) m_tuple1 = nt1;
  else {
    m_tuple1=ntupleSvc()->book("FILE302/n1",CLID_ColumnWiseTuple,"Calib");
    if( m_tuple1 ) {
      
     status1 = m_tuple1->addItem ("ixtal",ixtal); 
     status1 = m_tuple1->addItem ("gi0",gi0);
     status1 = m_tuple1->addItem ("calibConst",calibConst);
     status1 = m_tuple1->addItem ("err",err);
     status1 = m_tuple1->addItem ("nhitxtal",nhitxtal);
      
    }
    else    {   // did not manage to book the N tuple....
      log << MSG::ERROR <<"Cannot book N-tuple:" << long(m_tuple1) << endmsg;
      return StatusCode::FAILURE;
    }
  }
 
  // use EmcCalibConstSvc
  StatusCode scCalib; 
  scCalib = Gaudi::svcLocator() -> service("EmcCalibConstSvc", m_emcCalibConstSvc);
  if( scCalib != StatusCode::SUCCESS){
    log << MSG::ERROR << "can not use EmcCalibConstSvc" << endreq;
  } 
  else {
    std::cout << "Test EmcCalibConstSvc   DigiCalibConst(0)=  " 
          << m_emcCalibConstSvc -> getDigiCalibConst(0) << std::endl;
  }
 
   
  //std::cout<<6239<<"\t"<< m_emcCalibConstSvc->getCrystalEmaxData(6239) <<endl;
    
 
   //init starting values for calibration constants from file or set to 1
  initCalibConst();  
 
  //digiConstCor();
 
  return StatusCode::SUCCESS;
}

---

The documentation for this class was generated from the following files:

• EmcBhaCalib.h
• EmcBhaCalib.cxx