EmcCalibTagDeadOrHot.cxx

Go to the documentation of this file.

//--------------------------------------------------------------------------
// File and Version Information:
//
// tag hot channel and dead channel in the Emc Calibration constant file
//
// Environment:
//  Software developed for the BESIII Detector at the BEPCII.
//
// Author List:
//      Chunxiu Liu
//
// Copyright Information:
//  Copyright (C) 2012-04-19               IHEP
//
//------------------------------------------------------------------------
 
//-----------------------
// This Class's Header --
//-----------------------
#include "EmcCalibTagDeadOrHot.h"
 
#include <cassert>
#include <cmath>
#include <cstdlib>
#include <fstream>
#include <iostream>
//-------------------------------
// Collaborating Class Headers --
//-------------------------------
 
#include "GaudiKernel/DataObject.h"
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/ISvcLocator.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/PropertyMgr.h"
#include "GaudiKernel/SmartDataPtr.h"
 
#include "CLHEP/Vector/ThreeVector.h"
using namespace std;
 
using CLHEP::Hep3Vector;
 
//--------------------
#include "GaudiKernel/MsgStream.h"
 
#include "TFile.h"
#include "TH1F.h"
#include "TObjArray.h"
#include "TROOT.h"
#include "TTree.h"
 
DECLARE_COMPONENT( EmcCalibTagDeadOrHot )
//----------------
// Constructors --
//----------------
EmcCalibTagDeadOrHot::EmcCalibTagDeadOrHot( const std::string& name, ISvcLocator* pSvcLocator )
    : Algorithm( name, pSvcLocator )
    , m_fileExt( "" )
    , m_fileDir( "/home/besdata/public/liucx/Calib/" )
    , m_deadIxtal( -999 )
    , m_hotIxtal( -999 ) {
 
  // Declare the properties
  declareProperty( "fileExt", m_fileExt );
  declareProperty( "fileDir", m_fileDir );
  declareProperty( "deadIxtal", m_deadIxtal );
  declareProperty( "hotIxtal", m_hotIxtal );
 
  m_oldConstants = new double[6240];
  m_IxtalNumber  = new int[6240];
}
 
//--------------
// Destructor --
//--------------
EmcCalibTagDeadOrHot::~EmcCalibTagDeadOrHot() {
 
  if ( 0 != m_oldConstants )
  {
    delete[] m_oldConstants;
    m_oldConstants = 0;
  }
  if ( 0 != m_IxtalNumber )
  {
    delete[] m_IxtalNumber;
    m_IxtalNumber = 0;
  }
}
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode EmcCalibTagDeadOrHot::initialize() {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "in initialize()" << endmsg;
 
  // use EmcCalibConstSvc
  StatusCode scCalib;
  scCalib = Gaudi::svcLocator()->service( "EmcCalibConstSvc", m_emcCalibConstSvc );
  if ( scCalib != StatusCode::SUCCESS )
  { log << MSG::ERROR << "can not use EmcCalibConstSvc" << endmsg; }
  else
  {
    std::cout << "Test EmcCalibConstSvc   DigiCalibConst(0)=  "
              << m_emcCalibConstSvc->getDigiCalibConst( 0 ) << std::endl;
  }
 
  // init starting values for calibration constants from file or set to 1
 
  initCalibConst();
 
  SetDeadOrHotChannel();
 
  return StatusCode::SUCCESS;
}
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode EmcCalibTagDeadOrHot::execute() {
 
  MsgStream log( msgSvc(), name() );
  log << MSG::DEBUG << "in execute()" << endmsg;
 
  return StatusCode::SUCCESS;
}
 
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
StatusCode EmcCalibTagDeadOrHot::finalize() {
 
  MsgStream log( msgSvc(), name() );
 
  log << MSG::INFO << "in endRun()" << endmsg;
 
  return StatusCode::SUCCESS;
}
 
//* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
void EmcCalibTagDeadOrHot::initCalibConst() {
 
  MsgStream log( msgSvc(), name() );
 
  int nConstEmc;
 
  nConstEmc = m_emcCalibConstSvc->getDigiCalibConstNo();
 
  if ( nConstEmc != 6240 ) cout << "number of calibconst=" << nConstEmc << endl;
 
  for ( int i = 0; i < nConstEmc; i++ )
  {
 
    m_oldConstants[i] = m_emcCalibConstSvc->getDigiCalibConst( i );
    m_IxtalNumber[i]  = m_emcCalibConstSvc->getIxtalNumber( i );
  }
  /*
  ifstream In;
  In.open("/home/bes/liucx/GrelConstIxtal2013.txt");
  ofstream out;
  out.open("thephi.dat");
    int ixtal,ithe,iphi;
  double rel;
  while( !In.eof() ) {
 
    In >> ixtal>>rel;
    ithe= m_emcCalibConstSvc -> getThetaIndex(ixtal);
    iphi= m_emcCalibConstSvc -> getPhiIndex(ixtal);
    out<<ixtal<<"\t"<<ithe<<"\t"<<iphi<<"\t"<<rel<<endl;
  }
  In.close();
 
  out.close();
 
  cout<<"1,10,101"<<"\t"<<m_emcCalibConstSvc -> getIndex(1,10,101)<<endl;
  cout<<"1,30,105"<<"\t"<<m_emcCalibConstSvc -> getIndex(1,30,105)<<endl;
  cout<<"1,35,40"<<"\t"<<m_emcCalibConstSvc -> getIndex(1,35,40)<<endl;
  cout<<"0,2,40"<<"\t"<<m_emcCalibConstSvc -> getIndex(0,2,40)<<endl;
  cout<<"2,2,40"<<"\t"<<m_emcCalibConstSvc -> getIndex(2,2,40)<<endl;
  */
}
 
//* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * *
void EmcCalibTagDeadOrHot::SetDeadOrHotChannel() {
 
  int Nxtals = 6240;
 
  double DigiConst[6240];
  int    IxtalNumber[6240];
 
  for ( int ind = 0; ind < Nxtals; ind++ )
  {
 
    DigiConst[ind]   = m_oldConstants[ind];
    IxtalNumber[ind] = m_IxtalNumber[ind];
    if ( ind == m_deadIxtal ) IxtalNumber[ind] = -9;
    if ( ind == m_hotIxtal ) IxtalNumber[ind] = -99;
    // cout<<"ind="<<ind<<"\t"<< DigiConst[ind]<<endl;
  }
 
  TFile fconst( "EmcCalibConst_DeadOrHot.root", "recreate" );
 
  // define tree fill the absolute digicalibconsts into the root-file
  TTree* constr = new TTree( "DigiCalibConst", "DigiCalibConst" );
  constr->Branch( "DigiCalibConst", DigiConst, "DigiConst[6240]/D" );
  constr->Branch( "IxtalNumber", IxtalNumber, "IxtalNumber[6240]/I" );
 
  constr->Fill();
 
  constr->Write();
 
  delete constr;
 
  fconst.Close();
 
  // end tree
}