7 uint64_t
N()
const {
return n[0]+
n[1];}
8 uint64_t
Nx()
const {
return n[0];}
9 uint64_t
Ny()
const {
return n[1];}
14 uint64_t
x()
const {
return crd[0];}
15 uint64_t
y()
const {
return crd[1];}
16 uint64_t
z()
const {
return crd[2];}
19 uint64_t
crd[3] = {0,0,0};
20 uint64_t
n[2] = {0,0};
21 uint64_t
adc[2]= {0,0};
29 c.
t=slice.front().TDC().val;
30 c.
crd[2] = slice.front().Plane().val;
32 for(
auto const&
hit : slice){
34 const uint16_t
view =
hit.View().val-1;
40 for(
auto const&
hit : slice){
42 const uint16_t
view =
hit.View().val-1;
47 const double atten_len = 400./3.93;
51 *
exp((384.-pos )/atten_len)
52 /
exp((384./ 2. )/atten_len);
60 s<<c.
ADC()<<
" "<<c.
n[0]<<
" "<<c.
n[1]<<
" "<<c.
crd[0]<<
" "<<c.
crd[1]<<
" "<<c.
crd[2]<<
" "<<c.
t;
std::vector< DAQHit > HitList
clusterInfo CalcCluster(const novaddt::HitList &slice)
fvar< T > exp(const fvar< T > &x)
std::ostream & operator<<(std::ostream &s, const clusterInfo &c)