TofQElecSvc Class Reference

#include <TofQElecSvc.h>

Inheritance diagram for TofQElecSvc:

List of all members.

Public Member Functions
	TofQElecSvc (const std::string &name, ISvcLocator *svcloc)
	~TofQElecSvc ()
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvUnknown)
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	handle (const Incident &)
const unsigned int	Board (bool barrel, int id, bool eastEnd)
const unsigned int	Crate (bool barrel, int id, bool eastEnd)
const unsigned int	Fee (bool barrel, int id, bool eastEnd)
const unsigned int	Channel (bool barrel, int id, bool eastEnd)
const double	BQTC1 (int id, double q)
const double	BQTC2 (int id, double q)
const double	EQTC (int id, double q)
const double	BQChannel1 (int id, double qtc)
const double	BQChannel2 (int id, double qtc)
const double	EQChannel (int id, double qtc)
Static Public Member Functions
static const InterfaceID &	interfaceID ()
Private Member Functions
StatusCode	FillfromDatabase ()
Private Attributes
bool	m_hasbeeninitialized
IDataProviderSvc *	m_pCalibDataSvc
IntegerProperty	m_run
CalibData::TofElecData *	m_data

---

Detailed Description

Definition at line 17 of file TofQElecSvc.h.

---

Constructor & Destructor Documentation

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

Definition at line 22 of file TofQElecSvc.cxx.

References m_run.

                                                                     : Service(name, svcloc) {
  declareProperty("Run",m_run=1);
}

TofQElecSvc::~TofQElecSvc

(

)

Definition at line 76 of file TofQElecSvc.cxx.

References m_data.

                         {
  delete m_data;
  return;
}

---

Member Function Documentation

const unsigned int TofQElecSvc::Board	(	bool	barrel,
		int	id,
		bool	eastEnd
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 151 of file TofQElecSvc.cxx.

References CalibData::TofElecData::getBTof(), CalibData::TofElecData::getETof(), m_data, msgSvc(), and Bes_Common::WARNING.

                                                                         {
  MsgStream log(msgSvc(), name());    

  if( barrel ) {
    if( id<0 || id>175 ) {
      log << MSG::WARNING <<"Board : TofId is out of Range, tofid=" << id << endreq;  
      return 1000;
    }
  }
  else {
    if( id<0 || id>95 ) {
      log << MSG::WARNING <<"Board : TofId is out of Range, tofid=" << id << endreq;  
      return 1000;
    }
  }

  unsigned int board;
  if( barrel ) {
    if( eastEnd ) {
      board = (m_data->getBTof(id)).getNumEast(0);
    }
    else {
      board = (m_data->getBTof(id)).getNumWest(0);
    }
  }
  else {
    board = (m_data->getETof(id)).getNum(0);
  }

  return board;
}

const double TofQElecSvc::BQChannel1	(	int	id,
		double	qtc
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 361 of file TofQElecSvc.cxx.

References calibUtil::ERROR, CalibData::TofElecData::getBTof(), genRecEmupikp::i, m_data, and msgSvc().

                                                       {
  double qChannel = -999.0;
  MsgStream log(msgSvc(), name());
  if( id<0 || id>175 ) {
    log << MSG::ERROR <<"BChannel1 : TofId is out of Range, tofid=" << id << endreq;  
    return qChannel;
  }

  double ratio = (m_data->getBTof(id)).getSimP1(0);
  double p[10];
  for(unsigned int i=0; i<10; i++ ) {
    p[i] = (m_data->getBTof(id)).getSimP1(i+1);
  }
  qChannel = ( p[0] + p[1]*qtc + p[2]*qtc*qtc + p[3]*qtc*qtc*qtc )*( 1.0 - TMath::Erf((qtc-p[4])/p[5]) ) + ( p[6] + p[7]*qtc )*( 1.0 + TMath::Erf((qtc-p[8])/p[9]) );
  qChannel = ratio*qChannel;

  //  cout << "Sim East tofid=" << id << " ratio=" << ratio;
  //  for( unsigned int i=0; i<11; i++ ) { cout << " P" << i << "= " << p[i]; }
  //  cout << endl;

  return qChannel;
}

const double TofQElecSvc::BQChannel2	(	int	id,
		double	qtc
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 385 of file TofQElecSvc.cxx.

References calibUtil::ERROR, CalibData::TofElecData::getBTof(), genRecEmupikp::i, m_data, and msgSvc().

                                                       {
  double qChannel = -999.0;
  MsgStream log(msgSvc(), name());    
  if( id<0 || id>175 ) {
    log << MSG::ERROR <<"BChannel2 : TofId is out of Range, tofid=" << id << endreq;  
    return qChannel;
  }

  double ratio = (m_data->getBTof(id)).getSimP2(0);
  double p[10];
  for(unsigned int i=0; i<10; i++ ) {
    p[i] = (m_data->getBTof(id)).getSimP2(i+1);
  }
  qChannel = ( p[0] + p[1]*qtc + p[2]*qtc*qtc + p[3]*qtc*qtc*qtc )*( 1.0 - TMath::Erf((qtc-p[4])/p[5]) ) + ( p[6] + p[7]*qtc )*( 1.0 + TMath::Erf((qtc-p[8])/p[9]) );
  qChannel = ratio*qChannel;

  //  cout << "Sim West tofid=" << id << " ratio=" << ratio;
  //  for( unsigned int i=0; i<11; i++ ) { cout << " P" << i << "= " << p[i]; }
  //  cout << endl;

  return qChannel;
}

const double TofQElecSvc::BQTC1	(	int	id,
		double	q
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 283 of file TofQElecSvc.cxx.

References exp(), CalibData::TofElecData::getBTof(), genRecEmupikp::i, m_data, msgSvc(), and Bes_Common::WARNING.

                                                {
  double qtcNew = -999.0;
  MsgStream log(msgSvc(), name());    
  if( id<0 || id>175 ) {
    log << MSG::WARNING <<"BQTC1 : TofId is out of Range, tofid=" << id << endreq;
    return qtcNew;
  }
  if( fabs(q-10000.0)<1.0e-6 ) { return 10000.0; }
  else if( fabs(q+999.0)<1.0e-6 ) { return -999.0; }

  double ratio = (m_data->getBTof(id)).getP1(0);
  double p[10];
  for(unsigned int i=0; i<10; i++ ) {
    p[i] = (m_data->getBTof(id)).getP1(i+1);
  }
  q = q*ratio;
  qtcNew = p[0] + p[1]*q + p[2]*q*q + p[3]*q*q*q + p[4]*q*q*q*q + p[7]*exp(p[5]*(q-p[6]));

  //  cout << "East tofid=" << id << " ratio=" << ratio;
  //  for( unsigned int i=0; i<8; i++ ) { cout << " P" << i << "= " << p[i]; }
  //  cout << endl;

  return qtcNew;
}

const double TofQElecSvc::BQTC2	(	int	id,
		double	q
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 309 of file TofQElecSvc.cxx.

References exp(), CalibData::TofElecData::getBTof(), genRecEmupikp::i, m_data, msgSvc(), and Bes_Common::WARNING.

                                                {
  double qtcNew = -999.0;
  MsgStream log(msgSvc(), name());    
  if( id<0 || id>175 ) {
    log << MSG::WARNING <<"BQTC2 : TofId is out of Range, tofid=" << id << endreq;
    return qtcNew;
  }
  if( fabs(q-10000.0)<1.0e-6 ) { return 10000.0; }
  else if( fabs(q+999.0)<1.0e-6 ) { return -999.0; }

  double ratio = (m_data->getBTof(id)).getP2(0);
  double p[10];
  for(unsigned int i=0; i<10; i++ ) {
    p[i] = (m_data->getBTof(id)).getP2(i+1);
  }
  q = q*ratio;
  qtcNew = p[0] + p[1]*q + p[2]*q*q + p[3]*q*q*q + p[4]*q*q*q*q + p[7]*exp(p[5]*(q-p[6]));

  //  cout << "West tofid=" << id << " ratio=" << ratio;
  //  for( unsigned int i=0; i<8; i++ ) { cout << " P" << i << "= " << p[i]; }
  //  cout << endl;

  return qtcNew;
}

const unsigned int TofQElecSvc::Channel	(	bool	barrel,
		int	id,
		bool	eastEnd
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 250 of file TofQElecSvc.cxx.

References CalibData::TofElecData::getBTof(), CalibData::TofElecData::getETof(), m_data, msgSvc(), and Bes_Common::WARNING.

                                                                           {
  MsgStream log(msgSvc(), name());    

  if( barrel ) {
    if( id<0 || id>175 ) {
      log << MSG::WARNING <<"Channel : TofId is out of Range, tofid=" << id << endreq;  
      return 1000;
    }
  }
  else {
    if( id<0 || id>95 ) {
      log << MSG::WARNING <<"Channel : TofId is out of Range, tofid=" << id << endreq;  
      return 1000;
    }
  }

  unsigned int channel;
  if( barrel ) {
    if( eastEnd ) {
      channel = (m_data->getBTof(id)).getNumEast(1);
    }
    else {
      channel = (m_data->getBTof(id)).getNumWest(1);
    }
  }
  else {
    channel = (m_data->getETof(id)).getNum(1);
  }

  return channel;
}

const unsigned int TofQElecSvc::Crate	(	bool	barrel,
		int	id,
		bool	eastEnd
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 184 of file TofQElecSvc.cxx.

References CalibData::TofElecData::getBTof(), CalibData::TofElecData::getETof(), m_data, msgSvc(), and Bes_Common::WARNING.

                                                                         {
  MsgStream log(msgSvc(), name());    

  if( barrel ) {
    if( id<0 || id>175 ) {
      log << MSG::WARNING <<"Crate : TofId is out of Range, tofid=" << id << endreq;  
      return 1000;
    }
  }
  else {
    if( id<0 || id>95 ) {
      log << MSG::WARNING <<"Crate : TofId is out of Range, tofid=" << id << endreq;
      return 1000;
    }
  }

  unsigned int crate;
  if( barrel ) {
    if( eastEnd ) {
      crate = (m_data->getBTof(id)).getNumEast(2);
    }
    else {
      crate = (m_data->getBTof(id)).getNumWest(2);
    }
  }
  else {
    crate = (m_data->getETof(id)).getNum(2);
  }

  return crate;
}

const double TofQElecSvc::EQChannel	(	int	id,
		double	qtc
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 409 of file TofQElecSvc.cxx.

References calibUtil::ERROR, CalibData::TofElecData::getETof(), genRecEmupikp::i, m_data, and msgSvc().

                                                      {
  double qChannel = -999.0;
  MsgStream log(msgSvc(), name());    
  if( id<0 || id>95 ) {
    log << MSG::ERROR <<"EQTC : TofId is out of Range, tofid=" << id << endreq;
    return qChannel;
  }

  double ratio = (m_data->getETof(id)).getSimP(0);
  double p[10];
  for(unsigned int i=0; i<10; i++ ) {
    p[i] = (m_data->getETof(id)).getSimP(i+1);
  }
  qChannel = ( p[0] + p[1]*qtc + p[2]*qtc*qtc + p[3]*qtc*qtc*qtc )*( 1.0 - TMath::Erf((qtc-p[4])/p[5]) ) + ( p[6] + p[7]*qtc )*( 1.0 + TMath::Erf((qtc-p[8])/p[9]) );
  qChannel = ratio*qChannel;

  //  cout << "Endcap tofid=" << id << " ratio=" << ratio;
  //  for( unsigned int i=0; i<11; i++ ) { cout << " P" << i << "= " << p[i]; }
  //  cout << endl;

  return qChannel;
}

const double TofQElecSvc::EQTC	(	int	id,
		double	q
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 335 of file TofQElecSvc.cxx.

References exp(), CalibData::TofElecData::getETof(), genRecEmupikp::i, m_data, msgSvc(), and Bes_Common::WARNING.

                                               {
  double qtcNew = -999.0;
  MsgStream log(msgSvc(), name());    
  if( id<0 || id>95 ) {
    log << MSG::WARNING <<"EQTC : TofId is out of Range, tofid=" << id << endreq;
    return qtcNew;
  }
  if( fabs(q-10000.0)<1.0e-6 ) { return 10000.0; }
  else if( fabs(q+999.0)<1.0e-6 ) { return -999.0; }

  double ratio = (m_data->getETof(id)).getP(0);
  double p[10];
  for(unsigned int i=0; i<10; i++ ) {
    p[i] = (m_data->getETof(id)).getP(i+1);
  }
  q = q*ratio;
  qtcNew = p[0] + p[1]*q + p[2]*q*q + p[3]*q*q*q + p[4]*q*q*q*q + p[7]*exp(p[5]*(q-p[6]));

  //  cout << "Endcap tofid=" << id << " ratio=" << ratio;
  //  for( unsigned int i=0; i<8; i++ ) { cout << " P" << i << "= " << p[i]; }
  //  cout << endl;

  return qtcNew;
}

const unsigned int TofQElecSvc::Fee	(	bool	barrel,
		int	id,
		bool	eastEnd
	)			[virtual]

Implements ITofQElecSvc.

Definition at line 217 of file TofQElecSvc.cxx.

References CalibData::TofElecData::getBTof(), CalibData::TofElecData::getETof(), m_data, msgSvc(), and Bes_Common::WARNING.

                                                                       {
  MsgStream log(msgSvc(), name());    

  if( barrel ) {
    if( id<0 || id>175 ) {
      log << MSG::WARNING <<"Fee : TofId is out of Range, tofid=" << id << endreq;  
      return 1000;
    }
  }
  else {
    if( id<0 || id>95 ) {
      log << MSG::WARNING <<"Fee : TofId is out of Range, tofid=" << id << endreq;  
      return 1000;
    }
  }

  unsigned int crate;
  if( barrel ) {
    if( eastEnd ) {
      crate = (m_data->getBTof(id)).getNumEast(3);
    }
    else {
      crate = (m_data->getBTof(id)).getNumWest(3);
    }
  }
  else {
    crate = (m_data->getETof(id)).getNum(3);
  }

  return crate;
}

StatusCode TofQElecSvc::FillfromDatabase

(

)

[private]

Definition at line 82 of file TofQElecSvc.cxx.

References Bes_Common::FATAL, Bes_Common::INFO, m_data, m_pCalibDataSvc, msgSvc(), CalibData::TofElecData::sim(), deljobs::string, and test.

Referenced by handle().

                                        {

  MsgStream log(msgSvc(), name());    
  std::string fullPath = "/Calib/TofQElec";
  log << MSG::INFO<<" Tof calib fullPath = "<<fullPath<< endreq;
  SmartDataPtr<CalibData::TofElecData> test(m_pCalibDataSvc, fullPath);
  if(!test){
    log << MSG::FATAL << "TofQElecSvc could not find TofQElecData in TCDS!!" << endreq;
    return StatusCode::FAILURE;
  }       

  if(m_data) delete m_data;
  m_data = new CalibData::TofElecData;
  m_data->sim(test);
   
  /*
  for( vector<CalibData::BarrelElec*>::iterator iter1 = fBarrelElec.begin(); iter1 != fBarrelElec.end(); iter1++ ) {
    delete *iter1;
  }
  fBarrelElec.clear();
  for( vector<CalibData::EndcapElec*>::iterator iter2 = fEndcapElec.begin(); iter2 != fEndcapElec.end(); iter2++ ) {
    delete *iter2;
  }
  fEndcapElec.clear();

  for( unsigned ib=0; ib<176; ib++ ) {
    BarrelElec* bElec = new BarrelElec;
    for( unsigned int j=0; j<4; j++ ) {
      bElec->setNumEast( j, ib );
      bElec->setNumWest( j, ib );
    }
    for( unsigned int k=0; k<11; k++ ) {
      bElec->setP1( k, ib );
      bElec->setP2( k, ib );
      bElec->setSimP1( k, ib );
      bElec->setSimP2( k, ib );
    }
    fBarrelElec.push_back( bElec );
  }

  for( unsigned ie=0; ie<96; ie++ ) {
    EndcapElec* eElec = new EndcapElec;
    for( unsigned int j=0; j<4; j++ ) {
      eElec->setNum( j, ie );
    }
    for( unsigned int k=0; k<11; k++ ) {
      eElec->setP( k, ie );
      eElec->setSimP( k, ie );
    }
    fEndcapElec.push_back( eElec );
  }
  */

  return StatusCode::SUCCESS;
}

StatusCode TofQElecSvc::finalize

(

)

[virtual]

Definition at line 69 of file TofQElecSvc.cxx.

References Bes_Common::INFO, and msgSvc().

                                   {
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << name() << ": End of Run" << endreq;
  return StatusCode::SUCCESS;
}

void TofQElecSvc::handle

(

const Incident &

)

Definition at line 139 of file TofQElecSvc.cxx.

References Bes_Common::DEBUG, and FillfromDatabase().

                                           {
  MsgStream log( messageService(), name() );
  log << MSG::DEBUG << "handle: " << inc.type() << endreq;

  if ( inc.type() == "NewRun" ){
    log << MSG::DEBUG << "New Run" << endreq;
    StatusCode sc= FillfromDatabase();
  }
  return;
}

StatusCode TofQElecSvc::initialize

(

)

[virtual]

Definition at line 37 of file TofQElecSvc.cxx.

References Bes_Common::DEBUG, calibUtil::ERROR, Bes_Common::INFO, m_data, m_hasbeeninitialized, m_pCalibDataSvc, and msgSvc().

                                     {
  m_hasbeeninitialized=false;

  MsgStream log(msgSvc(), name());
  log << MSG::INFO << name() << ": Start of run initialisation" << 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);
  }

  StatusCode scc;
  
  log<<MSG::INFO << "setProperties()" << endreq;
  scc = service("CalibDataSvc", m_pCalibDataSvc, true);
  if ( !scc.isSuccess() ) {
    log<<MSG::ERROR<<"Could not get IDataProviderSvc interface of QElecXmlCnvSvc"<<endreq;
    return scc;
  } else {
    log << MSG::DEBUG<< "Retrieved IDataProviderSvc interface of QElecXmlCnvSvc"<<endreq;
  }
  // Get properties from the JobOptionsSvc
  scc = setProperties();
  m_data=0;
  return sc;
}  

static const InterfaceID& ITofQElecSvc::interfaceID

(

)

[inline, static, inherited]

Definition at line 14 of file ITofQElecSvc.h.

References IID_ITofQElecSvc().

{ return IID_ITofQElecSvc; }

StatusCode TofQElecSvc::queryInterface	(	const InterfaceID &	riid,
		void **	ppvUnknown
	)			[virtual]

Definition at line 27 of file TofQElecSvc.cxx.

References IID_ITofQElecSvc().

                                                                                    {
  
  if ( IID_ITofQElecSvc.versionMatch(riid) ) { 
    *ppvInterface = static_cast<ITofQElecSvc*> (this); 
  } else { 
    return Service::queryInterface(riid, ppvInterface) ; 
  }
  return StatusCode::SUCCESS;
}

---

Member Data Documentation

CalibData::TofElecData* TofQElecSvc::m_data [private]

Definition at line 49 of file TofQElecSvc.h.

Referenced by Board(), BQChannel1(), BQChannel2(), BQTC1(), BQTC2(), Channel(), Crate(), EQChannel(), EQTC(), Fee(), FillfromDatabase(), initialize(), and ~TofQElecSvc().

bool TofQElecSvc::m_hasbeeninitialized [private]

Definition at line 43 of file TofQElecSvc.h.

Referenced by initialize().

IDataProviderSvc* TofQElecSvc::m_pCalibDataSvc [private]

Definition at line 46 of file TofQElecSvc.h.

Referenced by FillfromDatabase(), and initialize().

IntegerProperty TofQElecSvc::m_run [private]

Definition at line 47 of file TofQElecSvc.h.

Referenced by TofQElecSvc().

---