void Histogram1D::multiply_with_ratio_exponented(const Histogram1D & nominator, const Histogram1D & denominator, double exponent){
check_compatibility(nominator);
check_compatibility(denominator);
const double * n_data = nominator.get_data();
const double * d_data = denominator.get_data();
double * data = get_data();
const size_t n = size();
for(size_t i=0; i<n; i++){
if(d_data[i]>0.0)
data[i] *= pow(n_data[i] / d_data[i], exponent);
}
}
void Histogram1D::operator*=(const Histogram1D & h) {
check_compatibility(h);
const double * hdata = h.get_data();
double * data = get_data();
const size_t n = size();
for (size_t i = 0; i < n; ++i) {
data[i] *= hdata[i];
}
}
