42 double B_NEW = (a_new+b_new)/2;
55 double CURVEG1 = fabs((1-
csi)*secondG1 + (b_new*
csi - a_new)*firstG1);
61 double firstG2 = ((2*(a_new-2*B_NEW+(
m-
shift))*sigma2)/exp((a_new-(
m-
shift))*(a_new-(
m-
shift))/(2*sigma2*sigma2)) - (2*(b_new-2*B_NEW+(
m-
shift))*sigma2)/exp((b_new-(
m-
shift))*(b_new-(
m-
shift))/(2*sigma2*sigma2))+ sqrt(2*TMath::Pi())*((B_NEW-(
m-
shift))*(B_NEW-(
m-
shift)) + sigma2*sigma2)*TMath::Erf((-a_new+(
m-
shift))/(sqrt(2)*sigma2)) - sqrt(2*TMath::Pi())*((B_NEW-(
m-
shift))*(B_NEW-(
m-
shift)) + sigma2*sigma2) * TMath::Erf((-b_new+(
m-
shift))/(sqrt(2)*sigma2)))/(2*sqrt(2*TMath::Pi()));
62 double secondG2 = (((2*sigma2*(a_new*a_new + B_NEW*B_NEW + a_new*(
m-
shift) + (
m-
shift)*(
m-
shift) - 2*B_NEW*(a_new+(
m-
shift)) + 2*(sigma2*sigma2)))/exp((a_new-(
m-
shift))*(a_new-(
m-
shift))/(2*(sigma2*sigma2))) - (2*sigma2*(b_new*b_new + B_NEW*B_NEW + b_new*(
m-
shift) + (
m-
shift)*(
m-
shift) - 2*B_NEW*(b_new + (
m-
shift)) + 2*(sigma2*sigma2)))/exp((b_new - (
m-
shift))*(b_new - (
m-
shift))/(2*(sigma2*sigma2))) - sqrt(2*TMath::Pi())*(-((B_NEW - (
m-
shift))*(B_NEW - (
m-
shift)) *(
m-
shift)) + (2*B_NEW - 3*(
m-
shift))*(sigma2*sigma2))*TMath::Erf((-a_new + (
m-
shift))/(sqrt(2)*sigma2)) + sqrt(2*TMath::Pi())* (-((B_NEW - (
m-
shift))*(B_NEW - (
m-
shift))*(
m-
shift)) + (2*B_NEW - 3*(
m-
shift))*(sigma2*sigma2)) *TMath::Erf((-b_new + (
m-
shift))/(sqrt(2)*sigma2)))/(2 *sqrt(2*TMath::Pi())));
63 double CURVEG2 = fabs((1-
csi)*secondG2 + (b_new*
csi - a_new)*firstG2);
RooRealProxy fraction_sigma