Re: TH2D::Integral()

From: Mario Kadastik <mario.kadastik_at_cern.ch>
Date: Tue, 2 Jun 2009 16:27:41 +0300


> the integral for each box, and add those values. Contrary to my
> expectations, the value does not turn out to be same as that
> obtained by integrating full histo.

<snip>

> // Box1
> float box1=BoxIntegral(h2,0,x_cut,0,y_cut);
>
> //Box2
> float box2=BoxIntegral(h2,0,x_cut,y_cut,1);
>
> //Box3
> float box3=BoxIntegral(h2,x_cut,1,0,y_cut);
>
> //Box4
> float box4=BoxIntegral(h2,x_cut,1,0,y_cut);

Aren't you taking a small overlap in the region where x_cut and y_cut lines overlap? I'd guess that if you exclude the overlap region you could well end up with the right result. Though the 25% increase is a bit more than I'd have expected from the overlap in case of 100 x 100 matrix.

Mario Received on Tue Jun 02 2009 - 15:27:45 CEST

This archive was generated by hypermail 2.2.0 : Tue Jun 02 2009 - 17:50:02 CEST