diff --git a/macros/CCEScan/pulseFit.C b/macros/CCEScan/pulseFit.C index 75ea27b..11d7320 100755 --- a/macros/CCEScan/pulseFit.C +++ b/macros/CCEScan/pulseFit.C @@ -286,7 +286,8 @@ - //MIchele - STart + //Evaluation of the systematic error on the integral of the pulse shape using a polinomial + //Define the interval were most likely the zeros should be double low_limit = -50; double high_limit = 50; @@ -299,7 +300,7 @@ ge_pulse->Fit("poly3","M"); Info("pulseFit","#chi^{2}/ndf: %6.2f/%d",poly3->GetChisquare(),poly3->GetNDF()); - double x1_poly3 = poly3->GetX(0,low_limit,0); + double x1_poly3 = poly3->GetX(0,low_limit,0); //GetX will return x* such that f(x*)=0 only if this exists, will (probably) return the argmin(f(x)) in the interval otherwis double x2_poly3 = poly3->GetX(0,0, high_limit); double z1_poly3 = -99; double z2_poly3 = -99; @@ -453,7 +454,9 @@ //Michele 3/9/2017 - //Estimate the systematics + //Estimate the systematics: # if a polynomial with zeros is found take as sys the semidifference + // between its integral and the integral obtained from the fitted theoretical model + // # if no polynomial behaved properly set the systematic error to -999 double int_sys = -999.0; double b_int_sys = true;