TofGeomSvc.cxx

Go to the documentation of this file.

#include "TofGeomSvc.h"
#include "GaudiKernel/IInterface.h"
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/StatusCode.h"
 
DECLARE_COMPONENT( TofGeomSvc )
 
TofGeomSvc::TofGeomSvc( const std::string& name, ISvcLocator* svcloc )
    : base_class( name, svcloc ) {}
 
/*TofGeomSvc::TofGeomSvc( const std::string& name, ISvcLocator* svcloc ) : Service(name,
svcloc) {}
 
StatusCode TofGeomSvc::queryInterface (const InterfaceID& riid, void** ppvInterface ){
 
  if ( IID_ITofGeomSvc.versionMatch(riid) ) {
    *ppvInterface = static_cast<ITofGeomSvc*> (this);
  } else {
    return Service::queryInterface(riid, ppvInterface) ;
  }
  return StatusCode::SUCCESS;
}*/
 
StatusCode TofGeomSvc::initialize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << name() << ": Start of run initialisation" << endmsg;
 
  StatusCode sc = Service::initialize();
  if ( sc.isFailure() ) return sc;
 
  //  get geometry data
  // Fill();
  return StatusCode::SUCCESS;
}
 
StatusCode TofGeomSvc::finalize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << name() << ": End of Run" << endmsg;
  return StatusCode::SUCCESS;
}
 
TofGeomSvc::~TofGeomSvc() {
  for ( vector<BTofGeo*>::iterator it1 = fBTofGeo.begin(); it1 != fBTofGeo.end(); it1++ )
    delete *it1;
  for ( vector<ETofGeo*>::iterator it2 = fETofGeo.begin(); it2 != fETofGeo.end(); it2++ )
    delete *it2;
  fBTofGeo.clear();
  fETofGeo.clear();
}
/*
void TofGeomSvc::Fill(){
    DB2BesGeoBarTOF obj1;
    DB2BesGeoEndTOF obj2;
    vector<DBBesGeoBarTOF> VecA;
    vector<DBBesGeoEndTOF> VecB;
    obj1.get_DBBesGeoBarTOF(VecA);
    obj2.get_DBBesGeoEndTOF(VecB);
    double inradm = VecA[0].get_inrad();
    double length = VecA[0].get_length();
    double thick  = VecA[0].get_thick();
   for(int ib=0;ib<176;ib++){
        BTofGeo* btof= new BTofGeo;
        double phi=0,phiMin=0,phiMax=0,inrad=0;
        double dphi = 0.0714;
        if(ib<88){
           inrad=inradm;
           phi  = dphi/2.+ib*dphi;      //phi at center of the module
           phiMin = phi-dphi/2.;
           phiMax = phi+dphi/2.;
        }
        else if(ib>87) {
           int ibb=ib-88;
           inrad=inradm+5.1;
           phi    = ibb*dphi;      //phi at center of the module
           phiMin = phi-dphi/2.;
           phiMax = phi+dphi/2.;
        }
        btof->setPhiMax(phiMax);
        btof->setPhiMin(phiMin);
        btof->setInrad(inrad);
        btof->setLength(length);
        btof->setThick(thick);
        fBTofGeo.push_back(btof);
 
    }
    for(int ie=0;ie<96;ie++){
     double inrad_e  = VecB[0].get_inrad();
    double outrad_e = VecB[0].get_outrad();
    double thick_e  = VecB[0].get_thick();
        double phi_e=0,phiMin_e=0,phiMax_e=0,thetaMax_e=0,thetaMin_e=0;
        ETofGeo* etof=new ETofGeo;
        double dphi = 0.1309;
        if(ie<48){
           phi_e  = dphi/2.+ie*dphi;      //phi at center of the counter
           phiMin_e = phi_e-dphi/2.;
           phiMax_e = phi_e+dphi/2.;
           thetaMax_e=0.93;
           thetaMin_e=0.83;
        }else{
           int ie2=ie-48;
           phi_e  = dphi/2.+ie2*dphi;      //phi at center of the counter
           phiMin_e = phi_e-dphi/2.;
           phiMax_e = phi_e+dphi/2.;
           thetaMax_e=-0.83;
           thetaMin_e=-0.93;
        }
        etof->setPhiMin(phiMin_e);
        etof->setPhiMax(phiMax_e);
        etof->setThetaMax(thetaMax_e);
        etof->setThetaMin(thetaMin_e);
        etof->setInrad(inrad_e);
        etof->setOutrad(outrad_e);
        etof->setThick(thick_e);
        fETofGeo.push_back(etof);
    }
}
*/
 
void TofGeomSvc::Dump() {
  std::cout << "Now can get the TOF Geometry Service!!" << std::endl;
  /*
    std::cout<<"Barrel TOF Counter Number = "<<fBTofGeo.size()<<std::endl;
    std::cout<<"Endcap TOF Counter Number = "<<fETofGeo.size()<<std::endl;
    for(int ib=0; ib<176; ib++){
       cout<<"The"<<" "<<ib<<" "<<"Barrel TOF phiMax is"<<" "<<BTof(ib)->getPhiMax()<<endl;
       cout<<"The"<<" "<<ib<<" "<<"Barrel TOF phiMin is"<<" "<<BTof(ib)->getPhiMin()<<endl;
    }
    for(int ie=0; ie<96; ie++){
       cout<<"The"<<" "<<ie<<" "<<"Endcap TOF phiMax is"<<" "<<ETof(ie)->getPhiMax()<<endl;
       cout<<"The"<<" "<<ie<<" "<<"Endcap TOF phiMin is"<<" "<<ETof(ie)->getPhiMin()<<endl;
    }
    cout<<"The 33rd  Barrel TOF Counter  inrad is"<<" "<<BTof(33)->getInrad()<<endl;
    cout<<"The 133rd Barrel TOF Counter  inrad is"<<" "<<BTof(133)->getInrad()<<endl;
    cout<<"The 133rd Barrel TOF Counter  length is"<<" "<<BTof(133)->getLength()<<endl;
    cout<<"The 133rd Barrel TOF Counter  thick  is"<<" "<<BTof(133)->getThick()<<endl;
    cout<<"The 22rd  Endcap TOF Counter  inrad is"<<" "<<ETof(22)->getInrad()<<endl;
    cout<<"The 22rd  Endcap TOF Counter  outrad is"<<" "<<ETof(22)->getOutrad()<<endl;
    cout<<"The 22rd  Endcap TOF Counter  thick  is"<<" "<<ETof(22)->getThick()<<endl;
    cout<<"The 22rd  Endcap TOF Counter  thetaMax is"<<" "<<ETof(22)->getThetaMax()<<endl;
    cout<<"The 22rd  Endcap TOF Counter  thetaMin is"<<" "<<ETof(22)->getThetaMin()<<endl;
    cout<<"The 90rd  Endcap TOF Counter  thetaMax is"<<" "<<ETof(90)->getThetaMax()<<endl;
    cout<<"The 90rd  Endcap TOF Counter  thetaMin is"<<" "<<ETof(90)->getThetaMin()<<endl;
  */
}
 
const double TofGeomSvc::getBPhiMax( unsigned id ) {
  if ( id < fBTofGeo.size() ) return fBTofGeo[id]->getPhiMax();
  return 0;
}
 
const double TofGeomSvc::getBPhiMin( unsigned id ) {
  if ( id < fBTofGeo.size() ) return fBTofGeo[id]->getPhiMin();
  return 0;
}
 
const double TofGeomSvc::getEPhiMax( unsigned id ) {
  if ( id < fETofGeo.size() ) return fETofGeo[id]->getPhiMax();
  return 0;
}
 
const double TofGeomSvc::getEPhiMin( unsigned id ) {
  if ( id < fETofGeo.size() ) return fETofGeo[id]->getPhiMin();
  return 0;
}
 
BTofGeo* TofGeomSvc::BTof( unsigned id ) const {
  if ( id < fBTofGeo.size() ) return fBTofGeo[id];
  return 0;
}
 
ETofGeo* TofGeomSvc::ETof( unsigned id ) const {
  if ( id < fETofGeo.size() ) return fETofGeo[id];
  return 0;
}