75 {
76
77 MsgStream log(
msgSvc(), name() );
78 log << MSG::INFO << "in initialize()" << endmsg;
79
80#ifndef OnlineMode
81 NTuplePtr nt(
ntupleSvc(),
"FILE102/n1" );
82 if ( nt ) m_tuple = nt;
83 else
84 {
85 m_tuple =
ntupleSvc()->book(
"FILE102/n1", CLID_ColumnWiseTuple,
"MdcRecEvent" );
86 if ( m_tuple )
87 {
89 m_tuple->addItem(
"NtrackMC",
g_ntrkMC, 0, 50 );
96 m_tuple->addItem(
"Ntrack",
g_ntrk, 0, 50 );
113 m_tuple->addItem(
"Testime",
g_estime );
114 }
115 else
116 {
117 log << MSG::ERROR << "Cannot book N-tuple:" << long( m_tuple ) << endmsg;
118 return StatusCode::FAILURE;
119 }
120 }
121
131 g_hitmap[0] =
histoSvc()->book(
"hm0",
"2d hitmap1", 80, -80., 80., 80, -80., 80. );
132 g_hitmap[1] =
histoSvc()->book(
"hm1",
"2d hitmap2", 80, -80., 80., 80, -80., 80. );
133 g_hitmap[2] =
histoSvc()->book(
"hm2",
"2d hitmap3", 80, -80., 80., 80, -80., 80. );
134 g_hitmap[3] =
histoSvc()->book(
"hm3",
"2d hitmap4", 80, -80., 80., 80, -80., 80. );
135 g_hitmap[4] =
histoSvc()->book(
"hm4",
"2d hitmap5", 80, -80., 80., 80, -80., 80. );
136 g_hitmap[5] =
histoSvc()->book(
"hm5",
"2d hitmap6", 80, -80., 80., 80, -80., 80. );
137 g_hitmap[6] =
histoSvc()->book(
"hm6",
"2d hitmap7", 80, -80., 80., 80, -80., 80. );
138 g_hitmap[7] =
histoSvc()->book(
"hm7",
"2d hitmap8", 80, -80., 80., 80, -80., 80. );
139 g_hitmap[8] =
histoSvc()->book(
"hm8",
"2d hitmap9", 80, -80., 80., 80, -80., 80. );
140 g_hitmap[9] =
histoSvc()->book(
"hm9",
"2d hitmap10", 80, -80., 80., 80, -80., 80. );
141 g_hitmap[10] =
histoSvc()->book(
"hm10",
"2d hitmap11", 80, -80., 80., 80, -80., 80. );
142
146#endif
147
148 m_ftFinder = new FTFinder();
149 m_ftFinder->init();
150 m_ftFinder->setBunchtime( m_bunchtime_MC );
151 m_ftFinder->setT0cal( m_mdc_tcal );
152 m_ftFinder->setAlgorithmPointer( this );
153
154#ifndef OnlineMode
155 StatusCode sc = service( "BesTimerSvc", m_timersvc );
156 if ( sc.isFailure() )
157 {
158 log << MSG::WARNING << name() << ": Unable to locate BesTimer Service" << endmsg;
159 return StatusCode::FAILURE;
160 }
161
162 m_timer[1] = m_timersvc->addItem( "Execution" );
163 m_timer[1]->propName( "nExecution" );
164#endif
165
166 return StatusCode::SUCCESS;
167}
IHistogram1D * g_nstereohit
NTuple::Array< float > g_pz
NTuple::Array< float > g_theta
NTuple::Array< float > g_pxMC
NTuple::Array< float > g_pzMC
NTuple::Array< float > g_pyMC
NTuple::Array< float > g_dr
IHistogram2D * g_hitmap[20]
NTuple::Array< float > g_kappa
NTuple::Array< float > g_ptMC
NTuple::Item< float > g_estime
NTuple::Array< float > g_dz
NTuple::Array< float > g_tanl
NTuple::Item< long > g_ntrk
NTuple::Item< long > g_ntrkMC
NTuple::Array< float > g_phi0MC
NTuple::Array< float > g_py
NTuple::Array< float > g_px
NTuple::Array< float > g_pt
NTuple::Array< float > g_p
NTuple::Array< float > g_phi0
NTuple::Array< float > g_vx
IHistogram1D * g_naxialhit
NTuple::Array< float > g_vy
NTuple::Array< float > g_theta0MC
NTuple::Array< float > g_vz
NTuple::Array< float > g_phi
AIDA::IHistogram1D * g_sigmaz
AIDA::IHistogram1D * g_chi2xy
AIDA::IHistogram1D * g_sigmaxy
AIDA::IHistogram1D * g_chi2sz
IHistogramSvc * histoSvc()