31 bound_t invalid_bound(NAN, NAN);
38 bound_t invalid_bound(NAN, NAN);
44 os <<
"RawTowerGeomContainer_Cylinderv1: radius: " <<
radius
49 for (bound_map_t::const_iterator iter =
eta_bound_map.begin();
52 os <<
"eta_bin[" << i <<
"](" << iter->first <<
", " << iter->second
58 for (bound_map_t::const_iterator iter =
phi_bound_map.begin();
61 os <<
"phi_bin[" << i <<
"](" << iter->first <<
", " << iter->second
76 <<
"RawTowerGeomContainer_Cylinderv1::get_etabounds - Asking for invalid bin in eta: "
90 <<
"RawTowerGeomContainer_Cylinderv1::get_phibounds - Asking for invalid bin in phi: "
104 double min_deta = 10;
106 for (bound_map_t::const_iterator iter =
eta_bound_map.begin();
109 const double mean_eta = 0.5 * (iter->first + iter->second);
111 if (eta >= iter->first && eta < iter->
second)
120 const double deta = fabs(mean_eta - eta);
134 <<
"RawTowerGeomContainer_Cylinderv1::get_etabin - ERROR - Asking for invalid bin in eta "
149 double min_dphi = 10;
151 for (bound_map_t::const_iterator iter =
phi_bound_map.begin();
154 const double mean_phi = 0.5 * (iter->first + iter->second);
156 const double phi_fold = phi -
round((phi - mean_phi) / 2. /
M_PI) * 2 *
M_PI;
158 if (phi_fold >= iter->first && phi_fold < iter->
second)
167 const double dphi = fabs(mean_phi - phi_fold);
181 <<
"RawTowerGeomContainer_Cylinderv1::get_phibin - ERROR - Asking for invalid bin in phi "
195 <<
"RawTowerGeomContainer_Cylinderv1::get_etacenter - Asking for invalid bin in eta: "
197 cout <<
"minbin: 0, maxbin " <<
get_etabins() << endl;
204 const std::pair<double, double>& bounds)
209 <<
"RawTowerGeomContainer_Cylinderv1::set_bounds - Asking for invalid bin in eta: "
211 cout <<
"minbin: 0, maxbin " <<
get_etabins() << endl;
215 std::pair<double, double> b_reg(bounds);
216 if (b_reg.first > b_reg.second)
218 b_reg.second = bounds.first;
219 b_reg.first = bounds.second;
231 <<
"RawTowerGeomContainer_Cylinderv1::get_phicenter - Asking for invalid bin in phi: "
233 cout <<
"minbin: 0, maxbin " <<
get_phibins() << endl;
240 const std::pair<double, double>& bounds)
245 <<
"RawTowerGeomContainer_Cylinderv1::set_bounds - Asking for invalid bin in phi: "
247 cout <<
"minbin: 0, maxbin " <<
get_phibins() << endl;
251 std::pair<double, double> b_reg(bounds);
252 if (b_reg.first > b_reg.second)
254 b_reg.second = bounds.first;
255 b_reg.first = bounds.second;