200 |
|
if(!jzbregion[ij]) continue; |
201 |
|
int GlobalBin=histos[ip][il][im][ij]->FindBin(M0,M12); |
202 |
|
float currcontent=histos[ip][il][im][ij]->GetBinContent(GlobalBin); |
203 |
< |
histos[ip][il][im][ij]->SetBinContent(GlobalBin,currcontent+rpredobssign); |
203 |
> |
if(ip!=predobs.size()-1) { |
204 |
> |
histos[ip][il][im][ij]->SetBinContent(GlobalBin,currcontent+rpredobssign); |
205 |
|
if(beverbose) cout << " Written " << rpredobssign << " to " << histos[ip][il][im][ij]->GetName() << endl; |
206 |
+ |
} |
207 |
+ |
else { |
208 |
+ |
if(rpredobs==obs) { |
209 |
+ |
histos[ip][il][im][ij]->SetBinContent(GlobalBin,currcontent+rpredobssign); |
210 |
+ |
if(beverbose) cout << " Written " << rpredobssign << " to " << histos[ip][il][im][ij]->GetName() << endl; |
211 |
+ |
|
212 |
+ |
} |
213 |
+ |
else { |
214 |
+ |
histos[ip][il][im][ij]->SetBinContent(GlobalBin,currcontent-rpredobssign); |
215 |
+ |
if(beverbose) cout << " Written " << -rpredobssign << " to " << histos[ip][il][im][ij]->GetName() << endl; |
216 |
+ |
} |
217 |
+ |
} |
218 |
|
} |
219 |
|
} |
220 |
|
} |
241 |
|
f->Close(); |
242 |
|
|
243 |
|
return 0; |
244 |
< |
} |
244 |
> |
} |