4#include "AsciiDmp/AsciiData.h"
5#include "AsciiDmp/dmplib.h"
12#include "TrkReco/Helix.h"
16#include "MdcRecGeo/MdcRecGeo.h"
21 TrkReco* _trkreco = TrkReco::getTrkReco();
31 Hep3Vector fmomentum(1,1,1);
32 Helix fHelix(fPivot,fmomentum,1.0);
42 ifstream is(
"AsciiDmp_out.dat");
45 std::cout<<
" The input file is wrong! " <<std::endl;
51 vector<EVENT> evt_list;
54 for (
int ix = 0;
ix < 100;
ix++) {
62 std::cerr <<
"wrong tag, got " << ex.
got()
68 std::cerr <<
"AsciiDumpException was caught!" << std::endl;
73 kine mctrk = mckn.kineList[0];
74 cout<<
"==============New Event, MC info ================"<<endl;
75 cout<<
"nTrk:"<<mckn.nKine
76 <<
" idPart"<<mctrk.idPart
79 <<
" theta0:"<<mctrk.theta0
80 <<
" phi0:"<<mctrk.phi0<<endl;
89 MdcRecGeo* fGeom=MdcRecGeo::getMdcRecGeo();
92 for(
unsigned i=0; i<mdigi.
nDigi; i++){
98 mhit.
geo = fGeom->Wire(mdigi.digi[i].layerNo-1,mdigi.digi[i].cellNo-1);
99 mhit.
tdc = mdigi.digi[i].driftTime;
100 mhit.
adc = mdigi.digi[i].energyDeposit;
101 mhit.
ddl = mhit.
tdc * 40/10000;
110 mhit.
stat = mhit.
stat |= 1073741824;
120 evt_list.push_back(evt);
125 for(
unsigned i=0; i<nt; i++){
std::string expected() const
static Bfield * getBfield(int)
returns Bfield object.
static vector< MdcRec_wirhit > * getMdcRecWirhitCol(void)
A class to handle update timing of static objects of tracking.
void dump(const std::string &message=std::string(""), const std::string &prefix=std::string("")) const
dumps debug information.
virtual void update(void)
updates an object.