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 33 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 288 of file EmcCalibConstSvc.cxx.

{}

Member Function Documentation

Dump()

void EmcCalibConstSvc::Dump

(

)

Definition at line 709 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 348 of file EmcCalibConstSvc.cxx.

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

getCrystalDeadEcut()

double EmcCalibConstSvc::getCrystalDeadEcut

(

int

Index

)

const

Definition at line 701 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 697 of file EmcCalibConstSvc.cxx.

                                                             {
  return m_CrystalEmaxData[Index];
}

Referenced by Dump().

getDigiCalibConst()

double EmcCalibConstSvc::getDigiCalibConst

(

int

No

)

const

Definition at line 628 of file EmcCalibConstSvc.cxx.

                                                         {
  double    digiCalibConst = 0.0;
  MsgStream log( msgSvc(), 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" << endmsg; }
  else { digiCalibConst = calConst->getDigiCalibConst( No ); }
 
  return digiCalibConst;
}

Referenced by Dump().

getDigiCalibConstNo()

int EmcCalibConstSvc::getDigiCalibConstNo

(

)

const

Definition at line 654 of file EmcCalibConstSvc.cxx.

                                                {
  int       digiCalibConstNo = 0;
  MsgStream log( msgSvc(), 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" << endmsg; }
  else { digiCalibConstNo = calConst->getDigiCalibConstNo(); }
 
  return digiCalibConstNo;
}

getIndex()

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

Definition at line 680 of file EmcCalibConstSvc.cxx.

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

getIxtalNumber()

int EmcCalibConstSvc::getIxtalNumber

(

int

No

)

const

Definition at line 641 of file EmcCalibConstSvc.cxx.

                                                   {
  int       IxtalNumber = 999999;
  MsgStream log( msgSvc(), 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" << endmsg; }
  else { IxtalNumber = calConst->getIxtalNumber( No ); }
 
  return IxtalNumber;
}

getIxtalNumberNo()

int EmcCalibConstSvc::getIxtalNumberNo

(

)

const

Definition at line 667 of file EmcCalibConstSvc.cxx.

                                             {
  int       IxtalNumberNo = 0;
  MsgStream log( msgSvc(), 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" << endmsg; }
  else { IxtalNumberNo = calConst->getIxtalNumberNo(); }
 
  return IxtalNumberNo;
}

getPartID()

unsigned int EmcCalibConstSvc::getPartID

(

int

Index

)

const

Definition at line 685 of file EmcCalibConstSvc.cxx.

                                                          {
  return m_theEmcStruc->getPartId( Index );
}

getPhiIndex()

unsigned int EmcCalibConstSvc::getPhiIndex

(

int

Index

)

const

Definition at line 693 of file EmcCalibConstSvc.cxx.

                                                            {
  return m_theEmcStruc->getPhi( Index );
}

getThetaIndex()

unsigned int EmcCalibConstSvc::getThetaIndex

(

int

Index

)

const

Definition at line 689 of file EmcCalibConstSvc.cxx.

                                                              {
  return m_theEmcStruc->getTheta( Index );
}

handle()

void EmcCalibConstSvc::handle

(

const Incident &

inc

)

to read EMC saturation energy and the energy threshold of 'special' dead channel from DataBase///////

Definition at line 357 of file EmcCalibConstSvc.cxx.

                                                   {
  MsgStream log( msgSvc(), name() );
  log << MSG::DEBUG << "handle: " << inc.type() << endmsg;
 
  if ( inc.type() == "NewRun" )
  {
    log << MSG::DEBUG << "NewRun" << endmsg;
 
    SmartDataPtr<Event::EventHeader> eventHeader( m_eventSvc, "/Event/EventHeader" );
    int                              run = eventHeader->runNumber();
 
    log << MSG::DEBUG << "handle: "
        << "Run in handle is: " << run << endmsg;
 
    if ( m_readSatuDeadEnDb == true )
      { //for MC and data
      
      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" << endmsg; }
        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 )
    { // for MC and data
 
      // cout<<"&&&&&&&&&&&&&&:&&&&&&&&&&&&&&&&&&&& EmcCalibConstSvc handle,run="<<run<<endl;
      if ( run < 0 ) run = -run;
 
      log << MSG::DEBUG << "handle: "
          << "Run in handle is: " << run << endmsg;
      // 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" << endmsg;
      // }
 
      // Dump();
    } // end of m_readSatuDeadEnDb==false
 
  } // end of NewRun
}

initialize()

StatusCode EmcCalibConstSvc::initialize ( )

virtual

Definition at line 299 of file EmcCalibConstSvc.cxx.

                                        {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << "EmcCalibConstSvc::initialize()" << endmsg;
 
  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 " << endmsg;
    return sc;
  }
 
  sc = service( "CalibDataSvc", m_calDataSvc, true );
  if ( sc == StatusCode::SUCCESS )
  { log << MSG::INFO << "Retrieve IDataProviderSvc" << endmsg; }
  else { log << MSG::FATAL << "can not get IDataProviderSvc" << endmsg; }
 
  //     sc = service("EmcRecGeoSvc", m_emcGeomSvc);
  //     if( sc != StatusCode::SUCCESS ){
  //      log << MSG::ERROR << "can not use EmcRecGeoSvc" << endmsg;
  // }
 
  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" << endmsg; }
 
  // Dump();
  // cout<<"&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& EmcCalibConstSvc initialize"<<endl;
  return StatusCode::SUCCESS;
}

setCrystalDeadEcut()

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

inline

Definition at line 54 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 48 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 30 of file EmcCalibConstSvc.h.

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

m_CrystalDeadEcut

double EmcCalibConstSvc::m_CrystalDeadEcut[6240]

Definition at line 82 of file EmcCalibConstSvc.h.

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

m_CrystalDeadEcut_2020

double EmcCalibConstSvc::m_CrystalDeadEcut_2020[6240]

Definition at line 83 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData

double EmcCalibConstSvc::m_CrystalEmaxData[6240]

Definition at line 63 of file EmcCalibConstSvc.h.

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

m_CrystalEmaxData_1

double EmcCalibConstSvc::m_CrystalEmaxData_1[6240]

Definition at line 64 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_2

double EmcCalibConstSvc::m_CrystalEmaxData_2[6240]

Definition at line 65 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4740

double EmcCalibConstSvc::m_CrystalEmaxData_4740[6240]

Definition at line 66 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4750

double EmcCalibConstSvc::m_CrystalEmaxData_4750[6240]

Definition at line 67 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4780

double EmcCalibConstSvc::m_CrystalEmaxData_4780[6240]

Definition at line 68 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4840

double EmcCalibConstSvc::m_CrystalEmaxData_4840[6240]

Definition at line 69 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4914

double EmcCalibConstSvc::m_CrystalEmaxData_4914[6240]

Definition at line 70 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_CrystalEmaxData_4946

double EmcCalibConstSvc::m_CrystalEmaxData_4946[6240]

Definition at line 71 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_DeadEcutVersion2020

bool EmcCalibConstSvc::m_DeadEcutVersion2020

Definition at line 85 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion1

bool EmcCalibConstSvc::m_EmaxVersion1

Definition at line 73 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion2

bool EmcCalibConstSvc::m_EmaxVersion2

Definition at line 74 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4740

bool EmcCalibConstSvc::m_EmaxVersion4740

Definition at line 75 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4750

bool EmcCalibConstSvc::m_EmaxVersion4750

Definition at line 76 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4780

bool EmcCalibConstSvc::m_EmaxVersion4780

Definition at line 77 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4840

bool EmcCalibConstSvc::m_EmaxVersion4840

Definition at line 78 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4914

bool EmcCalibConstSvc::m_EmaxVersion4914

Definition at line 79 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_EmaxVersion4946

bool EmcCalibConstSvc::m_EmaxVersion4946

Definition at line 80 of file EmcCalibConstSvc.h.

Referenced by EmcCalibConstSvc(), and handle().

m_theEmcStruc

EmcStructure* EmcCalibConstSvc::m_theEmcStruc

Definition at line 61 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