pderosa888/Angular_Distribution
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

fixed AD fit, chi-squared wrong

Browse Source
This commit is contained in:
Peter DeRosa 2022-07-18 20:16:40 -04:00
parent 434dd48296
commit f0214213b2
1 changed files with 7 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
AD.cxx
View File

@ -642,21 +642,22 @@ int main(int argc,char **argv){
//delta loop; //delta loop;
for(int i = 0; i< points; i++){ for(int i = 0; i< points; i++){
delta = i*step + delta_min; atan_delta = i*step + delta_min;
// delta = tan(atan_delta); delta = tan(atan_delta);
//now sum over all theta; //now sum over all theta;
for(int j = 0; j < dangler.size(); j++){ for(int j = 0; j < dangler.size(); j++){
Tangle = dangler[j]; Tangle = dangler[j];
//calculate Y_intensity_theory and and experiment and sum over theta //calculate Y_intensity_theory and and experiment and sum over theta
// rd0T = (rk01 + 2*delta*rk11 + pow(delta,2)*rk21)/(1+pow(delta,2)); // rd0T = (1 + 2*delta + pow(delta,2))/(1+pow(delta,2));
rd2T = (rk01 + 2*delta*rk11 + pow(delta,2)*rk21)/(1+pow(delta,2)); rd2T = (rk01 + 2*delta*rk11 + pow(delta,2)*rk21)/(1+pow(delta,2));
rd4T = (rk02 + 2*delta*rk12 + pow(delta,2)*rk22)/(1+pow(delta,2)); rd4T = (rk02 + 2*delta*rk12 + pow(delta,2)*rk22)/(1+pow(delta,2));
YT0 = 1; // P0(cos(theta)) = 1 YT0 = 1.; // P0(cos(theta)) = 1
YT2 = QD2*Bk11*rd2T*((1.5*pow(cos(Tangle),2)-.5)); YT2 = QD2*Bk11*rd2T*((1.5*pow(cos(Tangle),2)-.5));
@ -666,13 +667,13 @@ int main(int argc,char **argv){
YE = 1 + a2E*(1.5*pow(cos(Tangle),2)-.5) + a4E*(35./8.* pow(cos(Tangle),4) - 30./8.*pow(cos(Tangle),2) + 3./8.); YE = 1 + a2E*(1.5*pow(cos(Tangle),2)-.5) + a4E*(35./8.* pow(cos(Tangle),4) - 30./8.*pow(cos(Tangle),2) + 3./8.);
X2_total += pow((YT- YE),2)/(200.); X2_total += pow((YT- YE),2)/(dangler.size()-1)*(pow(Y_err,2));
} }
//Figure out denominator. //Figure out denominator.
chisqr.push_back(log(X2_total)); chisqr.push_back(log(X2_total));
tdelta.push_back(delta); tdelta.push_back(atan_delta);
X2_total = 0.; X2_total = 0.;