75 {
77 Gaudi::svcLocator()->service(
"MessageSvc",
msgSvc);
78 MsgStream log(
msgSvc,
"MdcCalRecHit");
79 log << MSG::DEBUG << "MdcCalRecHit::setKalHit()" << endreq;
80
83
84 identifier = (*it_hit) -> getMdcId();
85 m_layid = mdcid.
layer(identifier);
86 m_cellid = mdcid.
wire(identifier);
87 m_stat = 1;
88
89 m_lr = (*it_hit) -> getFlagLR();
90 if(-1 == m_lr) m_lr = 0;
91
92 m_docaInc = (*it_hit) -> getDocaIncl();
93 m_docaExc = (*it_hit) -> getDocaExcl();
94 if(1 == m_lr){
95 m_docaInc *= -1.0;
96 m_docaExc *= -1.0;
97 }
98
99 m_dmeas = (*it_hit) -> getDD();
100 m_errdmeas = 0.013;
101
102 m_tdc = (*it_hit) ->
getTdc();
103 m_tdrift = (*it_hit) -> getDT();
104
105 m_qhit = (*it_hit) -> getAdc();
106 m_zhit = (*it_hit) ->
getZhit();
108 m_chisqadd = 0.0;
109 m_tof = (*it_hit) ->
getTof();
110
111 m_Dr = (*it_hit) -> getDrExcl();
112 m_Fi0 = (*it_hit) -> getFi0Excl();
113 m_Cpa = (*it_hit) -> getCpaExcl();
114 m_Dz = (*it_hit) -> getDzExcl();
115 m_Tanl = (*it_hit) -> getTanlExcl();
116 m_helix = (*it_hit)->getHelixExcl();
117 m_helixerr = (*it_hit)->getErrorExcl();
118
119
120 m_docaInc *= m_df;
121 m_docaExc *= m_df;
122 m_dmeas *= m_df;
123 m_errdmeas *= m_df;
124 m_zhit *= m_df;
125 m_qhit *= m_qf;
126
127 m_resiInc = fabs(m_dmeas) - fabs(m_docaInc);
128 if( 0 == m_lr ) m_resilrInc = -1.0 * m_resiInc;
129 else m_resilrInc = m_resiInc;
130
131 m_resiExc = fabs(m_dmeas) - fabs(m_docaExc);
132 if( 0 == m_lr ) m_resilrExc = -1.0 * m_resiExc;
133 else m_resilrExc = m_resiExc;
134
135}
static int layer(const Identifier &id)
Values of different levels (failure returns 0)
static int wire(const Identifier &id)
1.12.0