16 {
18 Gaudi::svcLocator()->service(
"MessageSvc",
msgSvc);
19 MsgStream log(
msgSvc,
"MdcAliRecHit");
20 log << MSG::DEBUG << "MdcAliRecHit::setRecHit()" << endreq;
21
24
25 log << MSG::DEBUG << "hit Id: " << (*it_hit)->getId()
26 << " hits drift time " << (*it_hit)->getDriftT()
27 << " hits doca " << (*it_hit)->getDoca()
28 << endreq;
29
30 identifier = (*it_hit)->getMdcId();
31 m_layid = mdcid.
layer(identifier);
32 m_cellid = mdcid.
wire(identifier);
33 m_lr = (*it_hit) -> getFlagLR();
34 m_stat = (*it_hit) ->
getStat();
35
36 m_docaInc = (*it_hit) ->
getDoca();
37 m_docaExc = m_docaInc;
38
39 if( 0 == m_lr ){
40 m_dmeas = (*it_hit) -> getDriftDistLeft();
41 m_errdmeas = (*it_hit) -> getErrDriftDistLeft();
42 }else{
43 m_dmeas = (*it_hit) -> getDriftDistRight();
44 m_errdmeas = (*it_hit) -> getErrDriftDistRight();
45 }
46
47 m_rawtime = (*it_hit) -> getTdc();
48 m_tdrift = (*it_hit) -> getDriftT();
49
50 m_qhit = (*it_hit) -> getAdc();
51 m_zhit = (*it_hit) ->
getZhit();
54
55 double flen = (*it_hit) -> getFltLen();
56 m_tof = flen / 30.0;
57
58
59 m_docaInc *= m_df;
60 m_docaExc *= m_df;
61 m_dmeas *= m_df;
62 m_errdmeas *= m_df;
63 m_zhit *= m_df;
64 m_qhit *= m_qf;
65
66 m_resiInc = fabs(m_dmeas) - fabs(m_docaInc);
67 if( 0 == m_lr ) m_resilrInc = -1.0 * m_resiInc;
68 else m_resilrInc = m_resiInc;
69
70 m_resiExc = fabs(m_dmeas) - fabs(m_docaExc);
71 if( 0 == m_lr ) m_resilrExc = -1.0 * m_resiExc;
72 else m_resilrExc = m_resiExc;
73}
double getChisqAdd() const
bool getDoca(double trkpar[], double wpos[], double &doca, double whitPos[], double zini)