EmcCalibConstSvc Class Reference

#include <EmcCalibConstSvc.h>

Inheritance diagram for EmcCalibConstSvc:

Public Member Functions
	EmcCalibConstSvc (const std::string &name, ISvcLocator *svcloc)
	~EmcCalibConstSvc ()
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	handle (const Incident &)
double	getDigiCalibConst (int No) const
int	getIxtalNumber (int No) const
int	getDigiCalibConstNo () const
int	getIxtalNumberNo () const
int	getIndex (unsigned int PartId, unsigned int ThetaIndex, unsigned int PhiIndex) const
unsigned int	getPartID (int Index) const
unsigned int	getThetaIndex (int Index) const
unsigned int	getPhiIndex (int Index) const
double	getCrystalEmaxData (int Index) const
void	setCrystalEmaxData (double CrystalEmaxData[6240])
double	getCrystalDeadEcut (int Index) const
void	setCrystalDeadEcut (double CrystalDeadEcut[6240])
void	Dump ()

Public Attributes
IDataProviderSvc *	m_calDataSvc
EmcStructure *	m_theEmcStruc
double	m_CrystalEmaxData [6240]
double	m_CrystalEmaxData_1 [6240]
double	m_CrystalEmaxData_2 [6240]
double	m_CrystalEmaxData_4740 [6240]
double	m_CrystalEmaxData_4750 [6240]
double	m_CrystalEmaxData_4780 [6240]
double	m_CrystalEmaxData_4840 [6240]
double	m_CrystalEmaxData_4914 [6240]
double	m_CrystalEmaxData_4946 [6240]
bool	m_EmaxVersion1
bool	m_EmaxVersion2
bool	m_EmaxVersion4740
bool	m_EmaxVersion4750
bool	m_EmaxVersion4780
bool	m_EmaxVersion4840
bool	m_EmaxVersion4914
bool	m_EmaxVersion4946
double	m_CrystalDeadEcut [6240]
double	m_CrystalDeadEcut_2020 [6240]
bool	m_DeadEcutVersion2020

Detailed Description

Definition at line 16 of file EmcCalibConstSvc.h.

Constructor & Destructor Documentation

EmcCalibConstSvc()

EmcCalibConstSvc::EmcCalibConstSvc	(	const std::string &	name,
		ISvcLocator *	svcloc )

for version 1/////////

Definition at line 35 of file EmcCalibConstSvc.cxx.

                                                                              :
   base_class (name, svcloc) {
  declareProperty("ReadSatuDeadEnDb", m_readSatuDeadEnDb = true);
 
  for(int i=0;i<6240;i++){
    m_CrystalEmaxData[i]=0;
    m_CrystalDeadEcut[i]=-1; 
  }
  if (m_readSatuDeadEnDb==true){
    m_runFrom = 0;
    m_runTo = 0;
    m_inFlag=false;
  }
  /////////////////////////
  if (m_readSatuDeadEnDb==false){
    m_EmaxVersion1=false;
    m_EmaxVersion2=false;
    m_EmaxVersion4740=false;
    m_EmaxVersion4750=false;
    m_EmaxVersion4780=false;
    m_EmaxVersion4840=false;
    m_EmaxVersion4914=false;
    m_EmaxVersion4946=false;
    
    m_DeadEcutVersion2020=false;
    
    for(int i=0;i<6240;i++){
      m_CrystalEmaxData_1[i]=0;
      m_CrystalEmaxData_2[i]=0;
      m_CrystalEmaxData_4740[i]=0;
      m_CrystalEmaxData_4750[i]=0;
      m_CrystalEmaxData_4780[i]=0;
      m_CrystalEmaxData_4840[i]=0;
      m_CrystalEmaxData_4914[i]=0;
      m_CrystalEmaxData_4946[i]=0;
 
      m_CrystalDeadEcut_2020[i]=-1;
    }
 
  
    int ixtal;
    double emaxData;
    int Nixt=0;
    
    ///for version 1/////////
    string paraPath = getenv("EMCCALIBCONSTSVCROOT");
    paraPath += "/share/emax_data.dat";  // version=1
    ifstream in;
    in.open(paraPath.c_str());
    assert(in);
    
    while (in.peek() != EOF) 
      { 
    
    in >> ixtal>>emaxData;
    
    m_CrystalEmaxData_1[ixtal]=emaxData;
    Nixt++; 
    // cout<<"readFile==============="<<emaxData<<endl;
      }
    
    in.close();
    //at the end of version 1///// 
    
    ////version 2-4680data@2020///////////
    string paraPath2 = getenv("EMCCALIBCONSTSVCROOT");
    paraPath2 += "/share/emax_data_2.dat";  // version=2
    
    ifstream in2;
    in2.open(paraPath2.c_str());
    assert(in2);
    
    Nixt=0;
    while (in2.peek() != EOF) 
      { 
    
    in2 >> ixtal>>emaxData;
    
    m_CrystalEmaxData_2[ixtal]=emaxData;
    
    Nixt++; 
    // cout<<"readFile2==============="<<emaxData<<endl;
      }
    
    in2.close();
    //at the end of version 2/////
    
    ////version 4740data@2021///////////
    string paraPath4740 = getenv("EMCCALIBCONSTSVCROOT");
    paraPath4740 += "/share/emax_data_4740.dat";  // version=4740
    
    ifstream in4740;
    in4740.open(paraPath4740.c_str());
    assert(in4740);
    
    Nixt=0;
    while (in4740.peek() != EOF) 
      { 
    
    in4740 >> ixtal>>emaxData;
    
    m_CrystalEmaxData_4740[ixtal]=emaxData;
     
    Nixt++; 
    // cout<<"readFile4740==============="<<emaxData<<endl;
      }
    
    in4740.close();
    //at the end of version 4740/////
    
    ////version 4750data@2021///////////
    string paraPath4750 = getenv("EMCCALIBCONSTSVCROOT");
    paraPath4750 += "/share/emax_data_4750.dat";  // version=4750
    
    ifstream in4750;
    in4750.open(paraPath4750.c_str());
    assert(in4750);
    
    Nixt=0;
    while (in4750.peek() != EOF) 
      { 
    
    in4750 >> ixtal>>emaxData;
    
    m_CrystalEmaxData_4750[ixtal]=emaxData;
    
    Nixt++; 
    // cout<<"readFile4750==============="<<emaxData<<endl;
      }
    
    in4750.close();
    //at the end of version 4750/////
    
    ////version 4780data@2021///////////
    string paraPath4780 = getenv("EMCCALIBCONSTSVCROOT");
    paraPath4780 += "/share/emax_data_4780.dat";  // version=4780
    
    ifstream in4780;
    in4780.open(paraPath4780.c_str());
    assert(in4780);
    
    Nixt=0;
    while (in4780.peek() != EOF) 
      { 
    
    in4780 >> ixtal>>emaxData;
    
    m_CrystalEmaxData_4780[ixtal]=emaxData;
     
    Nixt++; 
    // cout<<"readFile4780==============="<<emaxData<<endl;
      }
    
    in4780.close();
    //at the end of version 4780/////
    
    ////version 4840data@2021///////////
    string paraPath4840 = getenv("EMCCALIBCONSTSVCROOT");
    paraPath4840 += "/share/emax_data_4840.dat";  // version=4840
    
    ifstream in4840;
    in4840.open(paraPath4840.c_str());
    assert(in4840);
    
    Nixt=0;
    while (in4840.peek() != EOF) 
      { 
    
    in4840 >> ixtal>>emaxData;
    
    m_CrystalEmaxData_4840[ixtal]=emaxData;
    
    Nixt++; 
    // cout<<"readFile4840==============="<<emaxData<<endl;
      }
    
    in4840.close();
    //at the end of version 4840/////
    
 
    ////version 4914data@2021///////////
    string paraPath4914 = getenv("EMCCALIBCONSTSVCROOT");
    paraPath4914 += "/share/emax_data_4914.dat";  // version=4914
    
    ifstream in4914;
    in4914.open(paraPath4914.c_str());
    assert(in4914);
    
    Nixt=0;
    while (in4914.peek() != EOF) 
      { 
    
    in4914 >> ixtal>>emaxData;
    
    m_CrystalEmaxData_4914[ixtal]=emaxData;
    
    Nixt++; 
    // cout<<"readFile4914==============="<<emaxData<<endl;
      }
    
    in4914.close();
    //at the end of version 4914/////
    
    ////version 4946data@2021///////////
    string paraPath4946 = getenv("EMCCALIBCONSTSVCROOT");
    paraPath4946 += "/share/emax_data_4946.dat";  // version=4946
    
    ifstream in4946;
    in4946.open(paraPath4946.c_str());
    assert(in4946);
    
    Nixt=0;
    while (in4946.peek() != EOF) 
      { 
    
    in4946 >> ixtal>>emaxData;
    
    m_CrystalEmaxData_4946[ixtal]=emaxData;
    
    Nixt++; 
    // cout<<"readFile4946==============="<<emaxData<<endl;
      }
    
    in4946.close();
    //at the end of version 4946/////
    
    
    ////read dead channel energy threshould Ecut///////////
    string paraPath3 = getenv("EMCCALIBCONSTSVCROOT");
    paraPath3 += "/share/DeadThresholdEcut_2020.conf";  
    
    ifstream in3;
    in3.open(paraPath3.c_str());
    assert(in3);
    
    Nixt=0;
    double satuEn,Ecut;
    double DeltaEcut=0.0; 
    
    while (in3.peek() != EOF) 
      { 
    
    in3 >> ixtal>>satuEn>>Ecut;
    if ((Ecut-satuEn)>DeltaEcut) {
      m_CrystalDeadEcut_2020[ixtal]=Ecut-DeltaEcut;
    } else {
      m_CrystalDeadEcut_2020[ixtal]=Ecut;
    }
    Nixt++; 
    
      }
    
    in3.close();
    
  }  //end of m_readSatuDeadEnDb==false
 
}

Referenced by EmcCalibConstSvc().

~EmcCalibConstSvc()

EmcCalibConstSvc::~EmcCalibConstSvc

(

)

Definition at line 293 of file EmcCalibConstSvc.cxx.

                                   {
 
}

Member Function Documentation

Dump()

void EmcCalibConstSvc::Dump

(

)

Definition at line 726 of file EmcCalibConstSvc.cxx.

{
  
  for(int iNo=0;iNo<6;iNo++){
    cout<<"getDigiCalibConst "<<getDigiCalibConst(iNo)<<endl;
  }
  /*
  cout<<"getDigiCalibConstNo "<<getDigiCalibConstNo()<<endl;
 
  
   cout<<"ind"<<"  "<<"getThetaIndex(ind)"<<"  "
       <<"getPhiIndex(ind)"<<"getPartID"<<endl;
 
  for(int ind=0; ind<6240;ind++){
 
    cout<<ind<<"  "<<getThetaIndex(ind)<<"  "
    <<getPhiIndex(ind)<<"  "<<getPartID(ind)<<endl;  
  
  }
 
  cout<<"getIndex(0,5,95)="<<getIndex(0,5,95)<<endl;
  cout<<"getIndex(1,43,119)="<<getIndex(1,43,119)<<endl;
  cout<<"getIndex(2,5,95)="<<getIndex(2,5,95)<<endl;
  cout<<"getIndex(0,5,96)="<<getIndex(0,5,96)<<endl;
  cout<<"getIndex(1,43,120)="<<getIndex(1,43,120)<<endl;
  cout<<"getIndex(2,5,96)="<<getIndex(2,5,96)<<endl;
  */
  for(int i=480;i<600;i++){
    cout<<i<<"\t"<<getCrystalEmaxData(i) <<"\t"<<getCrystalDeadEcut(i)<<endl;
  }
  
 
}

finalize()

StatusCode EmcCalibConstSvc::finalize ( )

virtual

Definition at line 365 of file EmcCalibConstSvc.cxx.

                                     {
 
     delete m_theEmcStruc;
 
     MsgStream log(messageService(), name());
     log << MSG::INFO << "EmcCalibConstSvc::finalize()" << endreq;
     return StatusCode::SUCCESS;
 
}

getCrystalDeadEcut()

double EmcCalibConstSvc::getCrystalDeadEcut

(

int

Index

)

const

Definition at line 717 of file EmcCalibConstSvc.cxx.

{
  //using the percent m_CrystalDeadEcut[ixtal]=deadEn/satuEn, the deadEcut of the corresponding data are calulated by  m_CrystalDeadEcut[ixtal]*m_CrystalEmaxData[ixtal]. 
  double CrystalDeadEcut;
  CrystalDeadEcut=m_CrystalDeadEcut[Index]*m_CrystalEmaxData[Index];
  return CrystalDeadEcut;
}

Referenced by Dump().

getCrystalEmaxData()

double EmcCalibConstSvc::getCrystalEmaxData

(

int

Index

)

const

Definition at line 712 of file EmcCalibConstSvc.cxx.

{
  return m_CrystalEmaxData[Index];
}

Referenced by Dump().

getDigiCalibConst()

double EmcCalibConstSvc::getDigiCalibConst

(

int

No

)

const

Definition at line 617 of file EmcCalibConstSvc.cxx.

 {
     double digiCalibConst = 0.0;
     MsgStream log(messageService(), name());
     
     std::string fullPath = "/Calib/EmcCal";
     SmartDataPtr<CalibData::EmcCalibData> calConst(m_calDataSvc, fullPath);
     if( ! calConst ){
      log << MSG::ERROR << "can not access to EmcDigi CalibData via SmartPtr" 
          << endreq;
     }else {
         digiCalibConst = calConst->getDigiCalibConst(No);
     }
     
     return digiCalibConst;
}

Referenced by Dump().

getDigiCalibConstNo()

int EmcCalibConstSvc::getDigiCalibConstNo

(

)

const

Definition at line 651 of file EmcCalibConstSvc.cxx.

 {
     int digiCalibConstNo = 0;
     MsgStream log(messageService(), name());
     
     std::string fullPath = "/Calib/EmcCal";
     SmartDataPtr<CalibData::EmcCalibData> calConst(m_calDataSvc, fullPath);
     if( ! calConst ){
      log << MSG::ERROR << "can not access to EmcDigiNo CalibData via SmartPtr" 
          << endreq;
     }else {
        digiCalibConstNo = calConst->getDigiCalibConstNo();
     }
     
     return digiCalibConstNo;
}

getIndex()

int EmcCalibConstSvc::getIndex	(	unsigned int	PartId,
		unsigned int	ThetaIndex,
		unsigned int	PhiIndex ) const

Definition at line 687 of file EmcCalibConstSvc.cxx.

{
  return m_theEmcStruc->getGeomIndex(PartId, ThetaIndex, PhiIndex);
}

getIxtalNumber()

int EmcCalibConstSvc::getIxtalNumber

(

int

No

)

const

Definition at line 634 of file EmcCalibConstSvc.cxx.

 {
     int IxtalNumber = 999999;
     MsgStream log(messageService(), name());
     
     std::string fullPath = "/Calib/EmcCal";
     SmartDataPtr<CalibData::EmcCalibData> calConst(m_calDataSvc, fullPath);
     if( ! calConst ){
      log << MSG::ERROR << "can not access to EmcDigi CalibData via SmartPtr" 
          << endreq;
     }else {
         IxtalNumber = calConst->getIxtalNumber(No);
     }
     
     return IxtalNumber;
}

getIxtalNumberNo()

int EmcCalibConstSvc::getIxtalNumberNo

(

)

const

Definition at line 668 of file EmcCalibConstSvc.cxx.

 {
     int IxtalNumberNo = 0;
     MsgStream log(messageService(), name());
     
     std::string fullPath = "/Calib/EmcCal";
     SmartDataPtr<CalibData::EmcCalibData> calConst(m_calDataSvc, fullPath);
     if( ! calConst ){
      log << MSG::ERROR << "can not access to EmcDigiNo CalibData via SmartPtr" 
          << endreq;
     }else {
        IxtalNumberNo = calConst->getIxtalNumberNo();
     }
     
     return IxtalNumberNo;
}

getPartID()

unsigned int EmcCalibConstSvc::getPartID

(

int

Index

)

const

Definition at line 695 of file EmcCalibConstSvc.cxx.

{
  return m_theEmcStruc->getPartId(Index);
}

getPhiIndex()

unsigned int EmcCalibConstSvc::getPhiIndex

(

int

Index

)

const

Definition at line 705 of file EmcCalibConstSvc.cxx.

{
  return m_theEmcStruc->getPhi(Index);
}

getThetaIndex()

unsigned int EmcCalibConstSvc::getThetaIndex

(

int

Index

)

const

Definition at line 700 of file EmcCalibConstSvc.cxx.

{
  return m_theEmcStruc->getTheta(Index);
}

handle()

void EmcCalibConstSvc::handle

(

const Incident &

inc

)

Definition at line 375 of file EmcCalibConstSvc.cxx.

                                                {
  MsgStream log( messageService(), name() );
  log << MSG::DEBUG << "handle: " << inc.type() << endreq;
  
  if ( inc.type() == "NewRun" ){
    log << MSG::DEBUG << "NewRun" << endreq;
    
    SmartDataPtr<Event::EventHeader> eventHeader(m_eventSvc,"/Event/EventHeader");
    int  run = eventHeader->runNumber();
    
    log << MSG::DEBUG << "handle: " <<"Run in handle is: "<<run<< endreq;
      
    if (m_readSatuDeadEnDb==true&&run<0){  //only for MC  
      if (run<0) run = -run;
 
      if (run>=m_runFrom&&run<=m_runTo) {
    m_inFlag=true;
      } else {
    m_inFlag=false;
      }
      if (m_inFlag==false){
    bool getSatuDead;
    getSatuDead=m_EmcSatuDeadEnSvc ->getEmcSatuDeadEnSvcInfo();
    if(!getSatuDead){
      log << MSG::ERROR << "can not get EmcSatuDeadEnSvc" << endreq;
    } else {
      m_runFrom=m_EmcSatuDeadEnSvc -> getRunFrom();
      m_runTo=m_EmcSatuDeadEnSvc -> getRunTo();   
      ////to read EMC saturation energy and the energy threshold of 'special' dead channel from DataBase///////
      cout << "current run=" << run<<endl;
      cout <<"RunFrom="<< m_runFrom<<",RunTo="<< m_runTo<<endl;
      cout << "in EmcCalibConstSvc  open getSatuEnFile()= " << m_EmcSatuDeadEnSvc -> getSatuEnFile()<<endl;
      cout << "open getDeadEnFile()= " << m_EmcSatuDeadEnSvc -> getDeadEnFile()<<endl;        
      string aEmcSatuEnFile;
      string aEmcDeadEnFile;      
      aEmcSatuEnFile= m_EmcSatuDeadEnSvc -> getSatuEnFile();
      aEmcDeadEnFile= m_EmcSatuDeadEnSvc -> getDeadEnFile();
      
      ifstream inSatuEn;
      inSatuEn.open(aEmcSatuEnFile.c_str());
 
      int ixtal;
      if (inSatuEn.is_open()){
      
        assert(inSatuEn);     
        double emaxData;      
        while (inSatuEn.peek() != EOF) 
          {     
        inSatuEn >> ixtal>>emaxData;    
        m_CrystalEmaxData[ixtal]=emaxData;          
          }   
        inSatuEn.close();   
      
      } else {
        std::cout << "EmcCalibConstSvc ERROR::Could  not open the file of " <<aEmcSatuEnFile<< endl;
        exit(1);
        
      }
 
      ifstream inDeadEn;
      inDeadEn.open(aEmcDeadEnFile.c_str());
 
      if (inDeadEn.is_open()){
 
        assert(inDeadEn);     
        double satuEn,deadEn;
        while (inDeadEn.peek() != EOF) 
          {         
        inDeadEn >> ixtal>>satuEn>>deadEn;
        //using the percent deadEn/satuEn, the deadEcut of the corresponding data are calulated by the emaxData*(deadEn/satuEn),i.e. m_CrystalDeadEcut[ixtal]*m_CrystalEmaxData[ixtal].
        m_CrystalDeadEcut[ixtal]=deadEn/satuEn;             
          }   
        inDeadEn.close();
      } else {
 
        std::cout << "EmcCalibConstSvc ERROR::Could  not open the file of " <<aEmcDeadEnFile<< endl;
        exit(1);       
      }
 
 
    }
    //Dump();
      } // end of m_inFlag==false   
    } //end of m_readSatuDeadEnDb==ture
 
 
    /////////////////////////////
 
    if (m_readSatuDeadEnDb==false&&run<0){  //only for MC  
      
      //cout<<"&&&&&&&&&&&&&&:&&&&&&&&&&&&&&&&&&&& EmcCalibConstSvc handle,run="<<run<<endl;
      if (run<0) run = -run;
      
      
      log << MSG::DEBUG << "handle: " <<"Run in handle is: "<<run<< endreq;
      //run <63075: EmaxVersion=1; run>=63075: EmaxVersion=2;
      cout<< "handle: " <<"Run in handle is: "<<run<<endl;
      
      
      if (run<63075&& !m_EmaxVersion1) {
    setCrystalEmaxData(m_CrystalEmaxData_1); 
    
    m_EmaxVersion1=true;
    m_EmaxVersion2=false;
    m_EmaxVersion4740=false;
    m_EmaxVersion4750=false;
    m_EmaxVersion4780=false;
    m_EmaxVersion4840=false;
    m_EmaxVersion4914=false;
    m_EmaxVersion4946=false;
    
    cout<< "handle: " <<"EmaxVersion="<<"111111111111111"<<endl;
      }
      
      
      if (run>=63075&&run<=65207&& !m_EmaxVersion2) {  // 2020 topup data~4680
    setCrystalEmaxData(m_CrystalEmaxData_2); 
    
    m_EmaxVersion1=false;
    m_EmaxVersion2=true;
    m_EmaxVersion4740=false;
    m_EmaxVersion4750=false;
    m_EmaxVersion4780=false;
    m_EmaxVersion4840=false;
    m_EmaxVersion4914=false;
    m_EmaxVersion4946=false;
    cout<< "handle: " <<"EmaxVersion="<<"emax data from topup data 4680@2020"<<endl;
      }
      
      
      if (run>=65208&&run<=65321&& !m_EmaxVersion4740) {  // 2021 data~4740
    setCrystalEmaxData(m_CrystalEmaxData_4740);         
    
    m_EmaxVersion1=false;
    m_EmaxVersion2=false;
    m_EmaxVersion4740=true;
    m_EmaxVersion4750=false;
    m_EmaxVersion4780=false;
    m_EmaxVersion4840=false;
    m_EmaxVersion4914=false;
    m_EmaxVersion4946=false;
    
    cout<< "handle: " <<"EmaxVersion="<<"emax data from 4740data@2021"<<endl;
      }
      
      if (run>=65322&&run<=65494&& !m_EmaxVersion4750) {  // 2021 data~4750
    setCrystalEmaxData(m_CrystalEmaxData_4750);         
    
    m_EmaxVersion1=false;
    m_EmaxVersion2=false;
    m_EmaxVersion4740=false;
    m_EmaxVersion4750=true;
    m_EmaxVersion4780=false;
    m_EmaxVersion4840=false;
    m_EmaxVersion4914=false;
    m_EmaxVersion4946=false;
    
    cout<< "handle: " <<"EmaxVersion="<<"emax data from 4750data@2021"<<endl;
      }
      
      if (run>=65495&&run<=65646&& !m_EmaxVersion4780) {  // 2021 data~4780
    setCrystalEmaxData(m_CrystalEmaxData_4780);         
    
    m_EmaxVersion1=false;
    m_EmaxVersion2=false;
    m_EmaxVersion4740=false;
    m_EmaxVersion4750=false;
    m_EmaxVersion4780=true;
    m_EmaxVersion4840=false;
    m_EmaxVersion4914=false;
    m_EmaxVersion4946=false;
    
    cout<< "handle: " <<"EmaxVersion="<<"emax data from 4780data@2021"<<endl;
      }
      
      
      if (run>=65647&&run<=65866&& !m_EmaxVersion4840) {  // 2021 data~4840
    setCrystalEmaxData(m_CrystalEmaxData_4840);         
    
    m_EmaxVersion1=false;
    m_EmaxVersion2=false;
    m_EmaxVersion4740=false;
    m_EmaxVersion4750=false;
    m_EmaxVersion4780=false;
    m_EmaxVersion4840=true;
    m_EmaxVersion4914=false;
    m_EmaxVersion4946=false;
    
    cout<< "handle: " <<"EmaxVersion="<<"emax data from 4840data@2021"<<endl;
      }
      
      if (run>=65867&&run<=65937&& !m_EmaxVersion4914) {  // 2021 data~4914
    setCrystalEmaxData(m_CrystalEmaxData_4914);         
    
    m_EmaxVersion1=false;
    m_EmaxVersion2=false;
    m_EmaxVersion4740=false;
    m_EmaxVersion4750=false;
    m_EmaxVersion4780=false;
    m_EmaxVersion4840=false;
    m_EmaxVersion4914=true;
    m_EmaxVersion4946=false;
    
    cout<< "handle: " <<"EmaxVersion="<<"emax data from 4914data@2021"<<endl;
      }
      
      if (run>=65938&& !m_EmaxVersion4946) {  // 2021 data~4946
    setCrystalEmaxData(m_CrystalEmaxData_4946);         
    
    m_EmaxVersion1=false;
    m_EmaxVersion2=false;
    m_EmaxVersion4740=false;
    m_EmaxVersion4750=false;
    m_EmaxVersion4780=false;
    m_EmaxVersion4840=false;
    m_EmaxVersion4914=false;
    m_EmaxVersion4946=true;
    
    cout<< "handle: " <<"EmaxVersion="<<"emax data from 4946data@2021"<<endl;
      }
     
      if (!m_DeadEcutVersion2020){
    setCrystalDeadEcut(m_CrystalDeadEcut_2020);
    m_DeadEcutVersion2020=true;
    cout<< "handle: " <<"DeadEcutVersion="<<"2020"<<endl;
      }
      
      //if(!getCrystalEmaxDataFileInfo()){
      //    log << MSG::ERROR << "can not find the crystal EmaxData file" << endreq;
      //}
      
      //Dump();
    } //end of m_readSatuDeadEnDb==false
 
 
    
  } //end of NewRun 
 
}

initialize()

StatusCode EmcCalibConstSvc::initialize ( )

virtual

Definition at line 306 of file EmcCalibConstSvc.cxx.

                                       {
  MsgStream log(messageService(), name());
  log << MSG::INFO << "EmcCalibConstSvc::initialize()" << endreq;
  
  StatusCode sc = Service::initialize();
  if( sc.isFailure() ) return sc;
 
  
  IIncidentSvc* incsvc;
  sc = service("IncidentSvc", incsvc);
  int priority = 100;
  if( sc.isSuccess() ){
    incsvc -> addListener(this, "NewRun", priority);
  }
  
  sc = serviceLocator()->service("EventDataSvc", m_eventSvc, true);
  if (sc .isFailure() ) {
    log << MSG::ERROR << "Unable to find EventDataSvc " << endreq;
    return sc;
  }
 
  sc = service("CalibDataSvc", m_calDataSvc, true);
  if( sc == StatusCode::SUCCESS ){
    log << MSG::INFO << "Retrieve IDataProviderSvc" << endreq;
  }else{
    log << MSG::FATAL << "can not get IDataProviderSvc" << endreq;
  }
 
  //     sc = service("EmcRecGeoSvc", m_emcGeomSvc);
  //     if( sc != StatusCode::SUCCESS ){
  //      log << MSG::ERROR << "can not use EmcRecGeoSvc" << endreq;
  // }
  
  m_theEmcStruc= new EmcStructure() ; 
  m_theEmcStruc->setEmcStruc();
 
  /////////////////
  ISvcLocator* svcLocator = Gaudi::svcLocator();
  sc = svcLocator->service("EmcSatuDeadEnSvc", m_EmcSatuDeadEnSvc);
 
  //sc = serviceLocator()->service("EmcSatuDeadEnSvc", m_EmcSatuDeadEnSvc,true);
  if( sc == StatusCode::SUCCESS){
    //log << MSG::INFO 
      std::cout<<"Retrieve EmcSatuDeadEnSvc" << endl;
 
 
  } 
  else {
    log << MSG::FATAL << "can not get EmcSatuDeadEnSvc" << endreq; 
  }
 
 
 
 
  //Dump();
  //cout<<"&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& EmcCalibConstSvc initialize"<<endl;
  return StatusCode::SUCCESS;
}

setCrystalDeadEcut()

void EmcCalibConstSvc::setCrystalDeadEcut ( double CrystalDeadEcut[6240] )

inline

Definition at line 61 of file EmcCalibConstSvc.h.

  {
    for(int i=0;i<6240;i++){
      m_CrystalDeadEcut[i]=CrystalDeadEcut[i];
      
    }  
  };

Referenced by handle().

setCrystalEmaxData()

void EmcCalibConstSvc::setCrystalEmaxData ( double CrystalEmaxData[6240] )

inline

Definition at line 51 of file EmcCalibConstSvc.h.

  {
    for(int i=0;i<6240;i++){
      m_CrystalEmaxData[i]=CrystalEmaxData[i];
      
    }  
  };

Referenced by handle().

Member Data Documentation

m_calDataSvc

IDataProviderSvc* EmcCalibConstSvc::m_calDataSvc

Definition at line 31 of file EmcCalibConstSvc.h.

Referenced by getDigiCalibConst(), getDigiCalibConstNo(), getIxtalNumber(), getIxtalNumberNo(), and initialize().

m_CrystalDeadEcut

double EmcCalibConstSvc::m_CrystalDeadEcut[6240]

Definition at line 94 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), getCrystalDeadEcut(), handle(), and setCrystalDeadEcut().

m_CrystalDeadEcut_2020

double EmcCalibConstSvc::m_CrystalDeadEcut_2020[6240]

Definition at line 95 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData

double EmcCalibConstSvc::m_CrystalEmaxData[6240]

Definition at line 75 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), getCrystalDeadEcut(), getCrystalEmaxData(), handle(), and setCrystalEmaxData().

m_CrystalEmaxData_1

double EmcCalibConstSvc::m_CrystalEmaxData_1[6240]

Definition at line 76 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_2

double EmcCalibConstSvc::m_CrystalEmaxData_2[6240]

Definition at line 77 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4740

double EmcCalibConstSvc::m_CrystalEmaxData_4740[6240]

Definition at line 78 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4750

double EmcCalibConstSvc::m_CrystalEmaxData_4750[6240]

Definition at line 79 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4780

double EmcCalibConstSvc::m_CrystalEmaxData_4780[6240]

Definition at line 80 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4840

double EmcCalibConstSvc::m_CrystalEmaxData_4840[6240]

Definition at line 81 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4914

double EmcCalibConstSvc::m_CrystalEmaxData_4914[6240]

Definition at line 82 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4946

double EmcCalibConstSvc::m_CrystalEmaxData_4946[6240]

Definition at line 83 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_DeadEcutVersion2020

bool EmcCalibConstSvc::m_DeadEcutVersion2020

Definition at line 97 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion1

bool EmcCalibConstSvc::m_EmaxVersion1

Definition at line 85 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion2

bool EmcCalibConstSvc::m_EmaxVersion2

Definition at line 86 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4740

bool EmcCalibConstSvc::m_EmaxVersion4740

Definition at line 87 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4750

bool EmcCalibConstSvc::m_EmaxVersion4750

Definition at line 88 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4780

bool EmcCalibConstSvc::m_EmaxVersion4780

Definition at line 89 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4840

bool EmcCalibConstSvc::m_EmaxVersion4840

Definition at line 90 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4914

bool EmcCalibConstSvc::m_EmaxVersion4914

Definition at line 91 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4946

bool EmcCalibConstSvc::m_EmaxVersion4946

Definition at line 92 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_theEmcStruc

EmcStructure* EmcCalibConstSvc::m_theEmcStruc

Definition at line 73 of file EmcCalibConstSvc.h.

Referenced by finalize(), getIndex(), getPartID(), getPhiIndex(), getThetaIndex(), and initialize().

---

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

• EmcCalibConstSvc.h
• EmcCalibConstSvc.cxx