Double_t histFunction2D::evaluate(Double_t *x, Double_t *par) {
Double_t xx = x[0];
int binY=GetBinNumber(xx,lowEdgeY,highEdgeY,nbinY);
Double_t total=0;
for (int j=0;j<nbinX;j++) {
total+=fabs(par[j])*hist->GetBinContent(j+1,binY);
}
// add a constant term
total+=par[(int)nbinX];
return total;
}
void SearchBinEventCount::Fill(double HT, double MHT, int NJets, int BTags, double Weight)
{
double bin = GetBinNumber(HT,MHT,NJets,BTags);
if(bin<bins_.size()+2)
{
fullTH1D_->Fill(bin-0.01, Weight);
unsigned int splitHist=0;
// std::cout<<"bin before split: "<<bin<<std::endl;
for(int ii=0;bin>splitAfter_;ii++)
{
splitHist++;
bin = bin-splitAfter_;
}
// if(splitHist==3)std::cout<<"BinForSplit: "<<bin<<" with splitHistNumber "<<splitHist<<" and TH1DSearchBinsSplit_.size(): "<<TH1DSearchBinsSplit_.size()<<std::endl;
splitTH1D_[splitHist]->Fill(bin-0.1, Weight);
}
}
