135 os <<
"PHG4CylinderCellGeom::identify - layer: " <<
layer
141 os <<
", zbins: " <<
nzbins
142 <<
", zmin: " <<
zmin
143 <<
", zstepsize: " <<
zstep;
146 os <<
", etabins: " <<
nzbins
147 <<
", etamin: " <<
zmin
148 <<
", etastepsize: " <<
zstep;
151 os <<
", etabins: " <<
nzbins
152 <<
", etamin: " <<
zmin
153 <<
", etastepsize: " <<
zstep;
156 os <<
", etabins: " <<
nzbins<<
" for Spacal";
159 os <<
"no valid binning method: " <<
binning << std::endl;
163 os <<
", phimin: " <<
phimin
170 std::pair<double, double>
173 if (ibin < 0 || ibin >
nzbins)
175 std::cout <<
PHWHERE <<
" Asking for invalid bin in z: " << ibin << std::endl;
180 double zhigh = zlow +
zstep;
181 return std::make_pair(zlow, zhigh);
184 std::pair<double, double>
187 if (ibin < 0 || ibin >
nzbins)
189 std::cout <<
PHWHERE <<
" Asking for invalid bin in z: " << ibin << std::endl;
195 double zhigh = zlow +
zstep;
196 return std::make_pair(zlow, zhigh);
200 std::pair<double, double>
205 std::cout <<
PHWHERE <<
"Asking for invalid bin in phi: " << ibin << std::endl;
210 double phihigh = philow +
phistep;
211 return std::make_pair(philow, phihigh);
224 return floor( (z-
zmin)/zstep );
236 return floor( (eta-
zmin)/zstep );
242 double norm_phi =
phi;
245 int nwraparound = -floor((phi-
phimin) * 0.5 /
M_PI);
246 norm_phi += 2*
M_PI*nwraparound;
249 return floor( (norm_phi-
phimin)/phistep );
255 if (ibin < 0 || ibin >
nzbins)
257 std::cout <<
PHWHERE <<
"Asking for invalid bin in z: " << ibin << std::endl;
267 if (ibin < 0 || ibin >
nzbins)
269 std::cout <<
PHWHERE <<
"Asking for invalid bin in eta: " << ibin << std::endl;
270 std::cout <<
"minbin: 0, maxbin " << nzbins << std::endl;
282 std::cout <<
PHWHERE <<
"Asking for invalid bin in phi: " << ibin << std::endl;
299 return "Eta/Phi bins";
302 return "Eta/numslat bins";
305 return "SPACAL Tower bins";
334 std::cout << src<<
" : ";
355 std::cout << src<<
" : ";