MdcPrintSvc Class Reference

#include <MdcPrintSvc.h>

Inheritance diagram for MdcPrintSvc:

Public Member Functions
	MdcPrintSvc (const std::string &name, ISvcLocator *svcloc)
	~MdcPrintSvc ()
StatusCode	initialize ()
StatusCode	finalize ()
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvUnknown)
void	printRecMdcTrackCol () const
void	printRecMdcTrack (const RecMdcTrack *tk) const
void	printDigi (uint32_t getDigiFlag=0) const
void	printMdcMcHitCol () const
virtual void	printRecMdcTrack (const RecMdcTrack *tk) const =0
virtual void	printDigi (uint32_t getDigiFlag=0) const =0
virtual void	printMdcMcHitCol () const =0

Additional Inherited Members
Static Public Member Functions inherited from IMdcPrintSvc
static const InterfaceID &	interfaceID ()

Detailed Description

Definition at line 21 of file MdcPrintSvc.h.

Constructor & Destructor Documentation

MdcPrintSvc()

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

Definition at line 30 of file MdcPrintSvc.cxx.

                                                                 :
  Service (name, svcloc) {
 
    // declare properties
    declareProperty("printHit",     m_printHit=1);
    declareProperty("printTrackDetail", m_printTrackDetail=1);
    declareProperty("printErrMat",  m_printErrMat=0);
    declareProperty("printTdcAdc",  m_printTdcAdc=0);
    declareProperty("printMcTk",    m_printMcTk=0);
  }

~MdcPrintSvc()

MdcPrintSvc::~MdcPrintSvc

(

)

Definition at line 41 of file MdcPrintSvc.cxx.

                         {
}

Member Function Documentation

finalize()

StatusCode MdcPrintSvc::finalize

(

)

Definition at line 70 of file MdcPrintSvc.cxx.

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

initialize()

StatusCode MdcPrintSvc::initialize

(

)

Definition at line 44 of file MdcPrintSvc.cxx.

                                  {
  MsgStream log(messageService(), name());
  log << MSG::INFO << "MdcPrintSvc::initialize()" << endreq;
 
  StatusCode sc = Service::initialize();
  if( sc.isFailure() ) return sc;
 
  ISvcLocator* svcLocator = Gaudi::svcLocator();
  svcLocator->service("EventDataSvc",eventSvc,true);
  if (NULL == eventSvc){
    std::cout << "ERROR Could not load service DataSvc" << std::endl;
    return StatusCode::FAILURE;
  }
 
  IRawDataProviderSvc* irawDataProviderSvc;
  sc = service ("RawDataProviderSvc", irawDataProviderSvc);
  m_rawDataProviderSvc = dynamic_cast<RawDataProviderSvc*> (irawDataProviderSvc);
  if ( sc.isFailure() ){
    log << MSG::FATAL << "Could not load RawDataProviderSvc!" << endreq;
    return StatusCode::FAILURE;
  }
 
 
  return StatusCode::SUCCESS;
}

printDigi()

void MdcPrintSvc::printDigi ( uint32_t  getDigiFlag = 0 ) const

virtual

Implements IMdcPrintSvc.

Definition at line 86 of file MdcPrintSvc.cxx.

                                                    {
  MdcDigiVec mdcDigiVec = m_rawDataProviderSvc->getMdcDigiVec(getDigiFlag);
  MdcDigiCol::iterator iter = mdcDigiVec.begin();
  std::cout<<"//==== "<<name()<<" print MdcDigiVec, nDigi="<<mdcDigiVec.size()<<" :"<<std::endl;
  for (int iDigi=0;iter!= mdcDigiVec.end(); iter++,iDigi++ ) {
    long l = MdcID::layer((*iter)->identify());
    long w = MdcID::wire((*iter)->identify());
    cout<<"("<<l<<","<<w;
    if(m_printMcTk) {
      int tkTruth = (*iter)->getTrackIndex();
      cout<<";"<<tkTruth;
    }
    if(m_printTdcAdc) {
      double tdc = RawDataUtil::MdcTime((*iter)->getTimeChannel());
      double adc = RawDataUtil::MdcCharge((*iter)->getChargeChannel());
      cout<<",t "<<tdc<<",c "<<adc;
    }
    cout<<")";
    if(iDigi%4==0) std::cout<<std::endl;
  }
  std::cout<<std::endl;
}

printMdcMcHitCol()

void MdcPrintSvc::printMdcMcHitCol ( ) const

virtual

Implements IMdcPrintSvc.

Definition at line 200 of file MdcPrintSvc.cxx.

                                       {
  MsgStream log(messageService(), name());
 
  SmartDataPtr<MdcMcHitCol> mdcMcHitCol(eventSvc,"/Event/MC/MdcMcHitCol");
  if (!mdcMcHitCol) {
    log << MSG::WARNING << "Could not find MdcMcHitCol" << endreq;
    return;
  }
 
  int ihit= 0;
  cout<<  "//==== "<<name()<<" print MdcMcHitCol nMcHit="<<mdcMcHitCol->size()<< endl;
  if (mdcMcHitCol->size()>0){
    cout<<"No. (layer,wire) mcTk  (x,y,z)   lr  driftdist(cm)  energy"<<endl;
  }
  MdcMcHitCol::iterator iter_mchit = mdcMcHitCol->begin();
  for (;iter_mchit != mdcMcHitCol->end(); iter_mchit++,ihit++ ) {
    const Identifier id= (*iter_mchit)->identify();
    int layer = MdcID::layer(id);
    int wire = MdcID::wire(id);
    int iMcTk =  (*iter_mchit)->getTrackIndex();
    double mcX = (*iter_mchit)->getPositionX()/10.;
    double mcY = (*iter_mchit)->getPositionY()/10.;
    double mcZ = (*iter_mchit)->getPositionZ()/10.;
    double mcLR = (*iter_mchit)->getPositionFlag();
    double mcDrift = (*iter_mchit)->getDriftDistance()/10.;
    double mcEnergy = (*iter_mchit)->getDepositEnergy()/10.;
    if (mcLR == 0) mcLR = -1;
    cout<<ihit <<" ("<<layer<<","<<wire<<") "<< iMcTk
      <<" ("<< mcX<<"," <<mcY<<","<< mcZ
      <<")  " <<mcLR
      <<"  " <<mcDrift
      <<"  " <<mcEnergy
      <<endl;
  }
}

printRecMdcTrack()

void MdcPrintSvc::printRecMdcTrack ( const RecMdcTrack *  tk ) const

virtual

Implements IMdcPrintSvc.

Definition at line 126 of file MdcPrintSvc.cxx.

                                                            {
  std::cout<< "//==== "<<name()<<" print RecMdcTrack No."<<tk->trackId()<<" :"<< std::endl;
  cout <<" dr "<<tk->helix(0)
    <<" phi0 "<<tk->helix(1)
    <<" cpa "<<tk->helix(2)
    <<" z0 "<<tk->helix(3)
    <<" tanl "<<tk->helix(4)
    <<endl;
  if(m_printTrackDetail){
    std::cout<<" q "<<tk->charge() 
      <<" theta "<<tk->theta()
      <<" phi "<<tk->phi()
      <<" x0 "<<tk->x()
      <<" y0 "<<tk->y()
      <<" z0 "<<tk->z()
      <<" r0 "<<tk->r()
      <<endl;
    std::cout <<" p "<<tk->p()
      <<" pt "<<tk->pxy()
      <<" pxyz("<<tk->px()
      <<","<<tk->py()
      <<","<<tk->pz()
      <<")"<<endl;
    std::cout<<" tkStat "<<tk->stat()
      <<" chi2 "<<tk->chi2()
      <<" ndof "<<tk->ndof()
      <<" nhit "<<tk->getNhits()
      <<" nst "<<tk->nster()
      <<endl;
    if(m_printErrMat){
      std::cout<< "errmat   " << std::endl;
      for (int i=0; i<15; i++){ std::cout<< " "<<tk->err(i); }
      std::cout<< "   " << std::endl;
    }
  }
 
  if(m_printHit){
    int haveDigi[43][288];
    if(m_printMcTk) {
      for(int ii=0;ii<43;ii++){
    for(int jj=0;jj<43;jj++){
      haveDigi[ii][jj]=-9999;
    }
      }
      MdcDigiVec mdcDigiVec = m_rawDataProviderSvc->getMdcDigiVec();
      MdcDigiCol::iterator iter = mdcDigiVec.begin();
      std::cout<<name()<<"//==== "<<name()<<" print MdcDigiVec, nDigi="<<mdcDigiVec.size()<<" :"<<std::endl;
      for (int iDigi=0;iter!= mdcDigiVec.end(); iter++,iDigi++ ) {
    long l = MdcID::layer((*iter)->identify());
    long w = MdcID::wire((*iter)->identify());
    haveDigi[l][w]=(*iter)->getTrackIndex();
      }
    }
 
    int nhits = tk->getVecHits().size();
    std::cout<<"nHits=" <<nhits<< std::endl;
    cout<<"No. ";
    if(m_printMcTk) cout<<"mcTkId ";
    cout<<"(layer,wire,lr) stat z"<<endl;
    for(int ii=0; ii <nhits ; ii++){
      Identifier id(tk->getVecHits()[ii]->getMdcId());
      int layer = MdcID::layer(id);
      int wire = MdcID::wire(id);
      cout<<ii<<":";
      if(m_printMcTk) { cout<<haveDigi[layer][wire]; }
      cout<<" ("<<layer<<","<<wire
    <<","<<tk->getVecHits()[ii]->getFlagLR()
    <<")  "<<tk->getVecHits()[ii]->getStat()
    <<"  "<<tk->getVecHits()[ii]->getZhit()<< " "<<std::endl;
    }//end of hit list
    std::cout << "  "<< std::endl;
  }
}

Referenced by printRecMdcTrackCol().

printRecMdcTrackCol()

void MdcPrintSvc::printRecMdcTrackCol

(

)

const

Definition at line 109 of file MdcPrintSvc.cxx.

                                          {
  MsgStream log(messageService(), name());
 
  SmartDataPtr<RecMdcTrackCol> recMdcTrackCol(eventSvc,"/Event/Recon/RecMdcTrackCol");
  if (!recMdcTrackCol) {
    log << MSG::WARNING << "Could not find RecMdcTrackCol" << endreq;
    return;
  }    
 
  std::cout<< "//======= "<<name()<<" RecMdcTrackCol size = "<<recMdcTrackCol->size() << std::endl;
  RecMdcTrackCol::iterator it = recMdcTrackCol->begin();
  for (;it!= recMdcTrackCol->end();it++){
    printRecMdcTrack(*it);
  }
 
}

Referenced by MdcHoughFinder::execute(), MdcTrkRecon::execute(), and MdcxCosmicSewer::execute().

queryInterface()

StatusCode MdcPrintSvc::queryInterface ( const InterfaceID &  riid, void **  ppvUnknown  )

virtual

Definition at line 77 of file MdcPrintSvc.cxx.

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

---

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

• MdcPrintSvc.h
• MdcPrintSvc.cxx