MdcFastTrkAlg Class Reference

#include <MdcFastTrkAlg.h>

List of all members.

Public Member Functions
	MdcFastTrkAlg (const std::string &name, ISvcLocator *pSvcLocator)
StatusCode	initialize ()
StatusCode	execute ()
StatusCode	finalize ()
StatusCode	beginRun ()
Private Attributes
FTFinder *	m_ftFinder
double	m_bunchtime_MC
double	m_mdc_tcal
NTuple::Tuple *	m_tuple
IBesTimerSvc *	m_timersvc
BesTimer *	m_timer [NTIMERS]

---

Detailed Description

Definition at line 18 of file MdcFastTrkAlg.h.

---

Constructor & Destructor Documentation

MdcFastTrkAlg::MdcFastTrkAlg	(	const std::string &	name,
		ISvcLocator *	pSvcLocator
	)

Definition at line 64 of file MdcFastTrkAlg.cxx.

References m_bunchtime_MC, m_ftFinder, and m_mdc_tcal.

                                                                            :
  Algorithm(name, pSvcLocator)
{
  // Declare the properties  
        declareProperty("bunchtime",      m_bunchtime_MC=8.0);
        declareProperty("T0cal",          m_mdc_tcal=200.0);
  m_ftFinder = 0;
  
}

---

Member Function Documentation

StatusCode MdcFastTrkAlg::beginRun

(

)

Definition at line 236 of file MdcFastTrkAlg.cxx.

References FTFinder::begin_run(), Bes_Common::INFO, m_ftFinder, and msgSvc().

                                   {
 
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in beginRun()" << endreq;

  m_ftFinder->begin_run();
 
  return StatusCode::SUCCESS;
}

StatusCode MdcFastTrkAlg::execute

(

)

Definition at line 188 of file MdcFastTrkAlg.cxx.

References Bes_Common::DEBUG, BesTimer::elapsed(), calibUtil::ERROR, FTFinder::event(), g_eventtime, m_ftFinder, m_timer, m_tuple, msgSvc(), BesTimer::start(), and BesTimer::stop().

                                  {

  MsgStream log(msgSvc(), name());
  log << MSG::DEBUG << "in execute()" << endreq;

#ifndef OnlineMode
  m_timer[1]->start();
#endif

  m_ftFinder->event();

#ifndef OnlineMode
  m_timer[1]->stop();

  //cout << "m_timer[1]->elapsed()::" << m_timer[1]->elapsed() << endl;
  //cout << "m_timer[1]->mean()::" << m_timer[1]->mean() << endl;
  g_eventtime = m_timer[1]->elapsed();

  if(m_tuple){     
    StatusCode status = m_tuple->write();
    if (!status.isSuccess()) {
      log << MSG::ERROR << "Can't fill ntuple!" << endreq;
    }
  }
#endif

  return StatusCode::SUCCESS;
}

StatusCode MdcFastTrkAlg::finalize

(

)

Definition at line 218 of file MdcFastTrkAlg.cxx.

References Bes_Common::DEBUG, Bes_Common::INFO, m_ftFinder, m_timersvc, msgSvc(), num_2Dtrk, num_3Dtrk, num_finaltrk, IBesTimerSvc::print(), and FTFinder::term().

                                   {

  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in finalize()" << endreq;

  m_ftFinder->term();
#ifndef OnlineMode
  log << MSG::DEBUG <<"num_2Dtrk: " << num_2Dtrk 
      <<" num_3Dtrk: " << num_3Dtrk 
      <<" num_finaltrk: " << num_finaltrk 
      << endmsg;
  m_timersvc->print();
#endif

  return StatusCode::SUCCESS;
}

StatusCode MdcFastTrkAlg::initialize

(

)

Definition at line 75 of file MdcFastTrkAlg.cxx.

References IBesTimerSvc::addItem(), calibUtil::ERROR, g_chi2sz, g_chi2xy, g_dr, g_dz, g_estime, g_eventNo, g_eventtime, g_hitmap, g_kappa, g_naxialhit, g_ncell, g_ncellMC, g_nhit, g_nstereohit, g_ntrk, g_ntrkMC, g_p, g_phi, g_phi0, g_phi0MC, g_pt, g_ptMC, g_px, g_pxMC, g_py, g_pyMC, g_pz, g_pzMC, g_sigmaxy, g_sigmaz, g_tanl, g_theta, g_theta0MC, g_vx, g_vy, g_vz, histoSvc(), Bes_Common::INFO, FTFinder::init(), m_bunchtime_MC, m_ftFinder, m_mdc_tcal, m_timer, m_timersvc, m_tuple, msgSvc(), ntupleSvc(), num_2Dtrk, num_3Dtrk, num_finaltrk, BesTimer::propName(), FTFinder::setAlgorithmPointer(), FTFinder::setBunchtime(), FTFinder::setT0cal(), and Bes_Common::WARNING.

                                    {

  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "in initialize()" << endreq;

#ifndef OnlineMode
  NTuplePtr nt(ntupleSvc(),"FILE102/n1");
  if ( nt ) m_tuple = nt;
  else {
      m_tuple=ntupleSvc()->book("FILE102/n1",CLID_ColumnWiseTuple,"MdcRecEvent");
      if( m_tuple ) {
         m_tuple->addItem ("eventNo",  g_eventNo);
         m_tuple->addItem ("NtrackMC", g_ntrkMC,0,50);
         m_tuple->addItem ("MCtheta0", g_ntrkMC,  g_theta0MC);
         m_tuple->addItem ("MCphi0",   g_ntrkMC, g_phi0MC);
         m_tuple->addItem ("pxMC",  g_ntrkMC, g_pxMC);
         m_tuple->addItem ("pyMC",  g_ntrkMC, g_pyMC);
         m_tuple->addItem ("pzMC",  g_ntrkMC, g_pzMC);
         m_tuple->addItem ("ptMC",  g_ntrkMC, g_ptMC);
         m_tuple->addItem ("Ntrack",g_ntrk, 0, 50);
         m_tuple->addItem ("px",    g_ntrk, g_px);
         m_tuple->addItem ("py",    g_ntrk, g_py);
         m_tuple->addItem ("pz",    g_ntrk, g_pz);
         m_tuple->addItem ("pt",    g_ntrk, g_pt);
         m_tuple->addItem ("p",     g_ntrk, g_p);
         m_tuple->addItem ("phi",   g_ntrk, g_phi);
         m_tuple->addItem ("theta", g_ntrk, g_theta);
         m_tuple->addItem ("vx",    g_ntrk, g_vx);
         m_tuple->addItem ("vy",    g_ntrk, g_vy);
         m_tuple->addItem ("vz",    g_ntrk, g_vz); 
         m_tuple->addItem ("dr",    g_ntrk, g_dr);
         m_tuple->addItem ("phi0",  g_ntrk, g_phi0);
         m_tuple->addItem ("kappa", g_ntrk, g_kappa);
         m_tuple->addItem ("dz",    g_ntrk, g_dz);
         m_tuple->addItem ("tanl",  g_ntrk, g_tanl);
         m_tuple->addItem ("eventtime", g_eventtime);
         m_tuple->addItem ("Testime", g_estime);
      }
      else    {   // did not manage to book the N tuple....
         log << MSG::ERROR <<"Cannot book N-tuple:" << long(m_tuple) << endmsg;
         //return StatusCode::FAILURE;
      }
   }

   g_sigmaxy   = histoSvc()->book( "sigmaxy", "1D sigmaxy", 100, -0.2, 0.2 );
   g_sigmaz = histoSvc()->book( "sigmaz", "1D sigmaz", 100, -3.5, 3.5);
   g_chi2xy = histoSvc()->book( "chi2xy", "1D chi2xy", 100, -0.01, 0.01 );
   g_chi2sz = histoSvc()->book( "chi2sz", "1D chi2sz", 100, -8.0, 8. );
   g_ncellMC = histoSvc()->book( "ncellmc", "1D ncellMC", 100, 0.0, 7000. );
   g_ncell = histoSvc()->book( "ncell", "1D ncell", 100, 0.0, 7000. );
   g_naxialhit = histoSvc()->book( "naxialhit", "2D axial hit", 30, 0.0, 30. );
   g_nstereohit = histoSvc()->book( "nstereohit", "2D stereo hit", 30, 0.0, 30. );
   g_nhit = histoSvc()->book( "nhit", "2D hit", 50, 0.0, 50. );
   g_hitmap[0] = histoSvc()->book( "hm0", "2d hitmap1", 80, -80., 80., 80, -80., 80. );
   g_hitmap[1] = histoSvc()->book( "hm1", "2d hitmap2", 80, -80., 80., 80, -80., 80. );
   g_hitmap[2] = histoSvc()->book( "hm2", "2d hitmap3", 80, -80., 80., 80, -80., 80. );
   g_hitmap[3] = histoSvc()->book( "hm3", "2d hitmap4", 80, -80., 80., 80, -80., 80. );
   g_hitmap[4] = histoSvc()->book( "hm4", "2d hitmap5", 80, -80., 80., 80, -80., 80. );
   g_hitmap[5] = histoSvc()->book( "hm5", "2d hitmap6", 80, -80., 80., 80, -80., 80. );
   g_hitmap[6] = histoSvc()->book( "hm6", "2d hitmap7", 80, -80., 80., 80, -80., 80. );
   g_hitmap[7] = histoSvc()->book( "hm7", "2d hitmap8", 80, -80., 80., 80, -80., 80. );
   g_hitmap[8] = histoSvc()->book( "hm8", "2d hitmap9", 80, -80., 80., 80, -80., 80. );
   g_hitmap[9] = histoSvc()->book( "hm9", "2d hitmap10", 80, -80., 80., 80, -80., 80. );
   g_hitmap[10] = histoSvc()->book( "hm10", "2d hitmap11", 80, -80., 80., 80, -80., 80. );

#endif

  /*NTuplePtr nth(ntupleSvc(),"FILE1/Hit/2");
  if ( nth ) m_htuple = nth;
  else {
      m_htuple=ntupleSvc()->book("FILE1/Hit/2",CLID_ColumnWiseTuple,"MdcRecHit");
      if( m_htuple ) {
         m_htuple->addItem ("sigmaxy", g_sigmaxy);
         m_htuple->addItem ("sigmaz", g_sigmaz);
         m_htuple->addItem ("chi2xy", g_chi2xy);
         m_htuple->addItem ("chi2sz", g_chi2sz);
         m_htuple->addItem ("ncellMC", g_ncellMC);
         m_htuple->addItem ("ncell", g_ncell);
      }
      else    {   // did not manage to book the N tuple....
         log << MSG::ERROR <<"Cannot book N-tuple:" << long(m_tuple) << endmsg;
         return StatusCode::FAILURE;
      }
   }*/
#ifndef OnlineMode   
   num_2Dtrk=0;
   num_3Dtrk=0; 
   num_finaltrk=0;
#endif

   //m_ftFinder->init();
   m_ftFinder = new FTFinder();
   m_ftFinder->init(); 
   m_ftFinder->setBunchtime(m_bunchtime_MC);
   m_ftFinder->setT0cal(m_mdc_tcal);
   m_ftFinder->setAlgorithmPointer(this);


#ifndef OnlineMode
  StatusCode sc = service( "BesTimerSvc", m_timersvc);
  if( sc.isFailure() ) {
     log << MSG::WARNING << name() << ": Unable to locate BesTimer Service" << endreq;
     return StatusCode::FAILURE;
  }

  m_timer[1] = m_timersvc->addItem("Execution");
  m_timer[1]->propName("nExecution");
#endif
 
  return StatusCode::SUCCESS;
}

---

Member Data Documentation

double MdcFastTrkAlg::m_bunchtime_MC [private]

Definition at line 28 of file MdcFastTrkAlg.h.

Referenced by initialize(), and MdcFastTrkAlg().

FTFinder* MdcFastTrkAlg::m_ftFinder [private]

Definition at line 27 of file MdcFastTrkAlg.h.

Referenced by beginRun(), execute(), finalize(), initialize(), and MdcFastTrkAlg().

double MdcFastTrkAlg::m_mdc_tcal [private]

Definition at line 29 of file MdcFastTrkAlg.h.

Referenced by initialize(), and MdcFastTrkAlg().

BesTimer* MdcFastTrkAlg::m_timer[NTIMERS] [private]

Definition at line 36 of file MdcFastTrkAlg.h.

Referenced by execute(), and initialize().

IBesTimerSvc* MdcFastTrkAlg::m_timersvc [private]

Definition at line 34 of file MdcFastTrkAlg.h.

Referenced by finalize(), and initialize().

NTuple::Tuple* MdcFastTrkAlg::m_tuple [private]

Definition at line 32 of file MdcFastTrkAlg.h.

Referenced by execute(), and initialize().

---