1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186
|
{
double khi1[10];
double khi2[10];
double khi3[10],khi4[10],khi5[10],khi6[10],khi7[10],khi8[10],khi9[10],khi10[10],khi11[10];
double khi24[10],khi12[10],khi13[10],khi14[10],khi15[10],khi16[10],khi17[10],khi18[10],khi19[10],khi20[10],khi21[10],khi22[10],khi23[10],khi25[10];
double khi11p;
double khi21p, khi31,khi41,khi51,khi61,khi71,khi81,khi91,khi101,khi111,khi121;
double khi131,khi141,khi151,khi161,khi171,khi181,khi191,khi201,khi211,khi221,khi231,khi241,khi251;
double segexp[10];
double esegmaexp[10];
ofstream fcoef("../3-root/val.txt");
ifstream f1(Form("../1-code-dr/fort.0303"));
ifstream f2(Form("../1-code-dr/fort.0307"));
ifstream f3(Form("../1-code-dr/fort.0311"));
ifstream f4(Form("../1-code-dr/fort.0315"));
ifstream f5(Form("../1-code-dr/fort.0319"));
ifstream f6(Form("../1-code-dr/fort.0703"));
ifstream f7(Form("../1-code-dr/fort.0707"));
ifstream f8(Form("../1-code-dr/fort.0711"));
ifstream f9(Form("../1-code-dr/fort.0715"));
ifstream f10(Form("../1-code-dr/fort.0719"));
ifstream f11(Form("../1-code-dr/fort.1103"));
ifstream f12(Form("../1-code-dr/fort.1107"));
ifstream f13(Form("../1-code-dr/fort.1111"));
ifstream f14(Form("../1-code-dr/fort.1115"));
ifstream f15(Form("../1-code-dr/fort.1119"));
ifstream f16(Form("../1-code-dr/fort.1503"));
ifstream f17(Form("../1-code-dr/fort.1507"));
ifstream f18(Form("../1-code-dr/fort.1511"));
ifstream f19(Form("../1-code-dr/fort.1515"));
ifstream f20(Form("../1-code-dr/fort.1519"));
ifstream f21(Form("../1-code-dr/fort.1903"));
ifstream f22(Form("../1-code-dr/fort.1907"));
ifstream f23(Form("../1-code-dr/fort.1911"));
ifstream f24(Form("../1-code-dr/fort.1915"));
ifstream f25(Form("../1-code-dr/fort.1919"));
ifstream fexp(Form("../../cross-section-q2-05-oct2013/Results/section1.txt"));
double qcm1[10], epsilon1[10], qprimcm1[10], thetacm1[10], phideg1[10], segmaDR1[10],costheta1[10];
double qcm2[10], epsilon2[10], qprimcm2[10], thetacm2[10], phideg2[10], segmaDR2[10],costheta2[10];
double qcm3[10], epsilon3[10], qprimcm3[10], thetacm3[10], phideg3[10], segmaDR3[10],costheta3[10];
double qcm4[10], epsilon4[10], qprimcm4[10], thetacm4[10], phideg4[10], segmaDR4[10],costheta4[10];
double qcm5[10], epsilon5[10], qprimcm5[10], thetacm5[10], phideg5[10], segmaDR5[10],costheta5[10];
double qcm6[10], epsilon6[10], qprimcm6[10], thetacm6[10], phideg6[10], segmaDR6[10],costheta6[10];
double qcm7[10], epsilon7[10], qprimcm7[10], thetacm7[10], phideg7[10], segmaDR7[10],costheta7[10];
double qcm8[10], epsilon8[10], qprimcm8[10], thetacm8[10], phideg8[10], segmaDR8[10],costheta8[10];
double qcm9[10], epsilon9[10], qprimcm9[10], thetacm9[10], phideg9[10], segmaDR9[10],costheta9[10];
double qcm10[10], epsilon10[10], qprimcm10[10], thetacm10[10], phideg10[10], segmaDR10[10],costheta10[10];
double qcm11[10], epsilon11[10], qprimcm11[10], thetacm11[10], phideg11[10], segmaDR11[10],costheta11[10];
double qcm12[10], epsilon12[10], qprimcm12[10], thetacm12[10], phideg12[10], segmaDR12[10],costheta12[10];
double qcm13[10], epsilon13[10], qprimcm13[10], thetacm13[10], phideg13[10], segmaDR13[10],costheta13[10];
double qcm14[10], epsilon14[10], qprimcm14[10], thetacm14[10], phideg14[10], segmaDR14[10],costheta14[10];
double qcm15[10], epsilon15[10], qprimcm15[10], thetacm15[10], phideg15[10], segmaDR15[10],costheta15[10];
double qcm16[10], epsilon16[10], qprimcm16[10], thetacm16[10], phideg16[10], segmaDR16[10],costheta16[10];
double qcm17[10], epsilon17[10], qprimcm17[10], thetacm17[10], phideg17[10], segmaDR17[10],costheta17[10];
double qcm18[10], epsilon18[10], qprimcm18[10], thetacm18[10], phideg18[10], segmaDR18[10],costheta18[10];
double qcm19[10], epsilon19[10], qprimcm19[10], thetacm19[10], phideg19[10], segmaDR19[10],costheta19[10];
double qcm20[10], epsilon20[10], qprimcm20[10], thetacm20[10], phideg20[10], segmaDR20[10],costheta20[10];
double qcm21[10], epsilon21[10], qprimcm21[10], thetacm21[10], phideg21[10], segmaDR21[10],costheta21[10];
double qcm22[10], epsilon22[10], qprimcm22[10], thetacm22[10], phideg22[10], segmaDR22[10],costheta22[10];
double qcm23[10], epsilon23[10], qprimcm23[10], thetacm23[10], phideg23[10], segmaDR23[10],costheta23[10];
double qcm24[10], epsilon24[10], qprimcm24[10], thetacm24[10], phideg24[10], segmaDR24[10],costheta24[10];
double qcm25[10], epsilon25[10], qprimcm25[10], thetacm25[10], phideg25[10], segmaDR25[10],costheta25[10];
for(int i=0; i<10; i++){
//lecture des fichiers txt
f1>>qcm1[i]>>epsilon1[i]>>qprimcm1[i]>>thetacm1[i]>>phideg1[i]>>segmaDR1[i];
f2>>qcm2[i]>>epsilon2[i]>>qprimcm2[i]>>thetacm2[i]>>phideg2[i]>>segmaDR2[i];
f3>>qcm3[i]>>epsilon3[i]>>qprimcm3[i]>>thetacm3[i]>>phideg3[i]>>segmaDR3[i];
f4>>qcm4[i]>>epsilon4[i]>>qprimcm4[i]>>thetacm4[i]>>phideg4[i]>>segmaDR4[i];
f5>>qcm5[i]>>epsilon5[i]>>qprimcm5[i]>>thetacm5[i]>>phideg5[i]>>segmaDR5[i];
f6>>qcm6[i]>>epsilon6[i]>>qprimcm6[i]>>thetacm6[i]>>phideg6[i]>>segmaDR6[i];
f7>>qcm7[i]>>epsilon7[i]>>qprimcm7[i]>>thetacm7[i]>>phideg7[i]>>segmaDR7[i];
f8>>qcm8[i]>>epsilon8[i]>>qprimcm8[i]>>thetacm8[i]>>phideg8[i]>>segmaDR8[i];
f9>>qcm9[i]>>epsilon9[i]>>qprimcm9[i]>>thetacm9[i]>>phideg9[i]>>segmaDR9[i];
f10>>qcm10[i]>>epsilon10[i]>>qprimcm10[i]>>thetacm10[i]>>phideg10[i]>>segmaDR10[i];
f11>>qcm11[i]>>epsilon11[i]>>qprimcm11[i]>>thetacm11[i]>>phideg11[i]>>segmaDR11[i];
f12>>qcm12[i]>>epsilon12[i]>>qprimcm12[i]>>thetacm12[i]>>phideg12[i]>>segmaDR12[i];
f13>>qcm13[i]>>epsilon13[i]>>qprimcm13[i]>>thetacm13[i]>>phideg13[i]>>segmaDR13[i];
f14>>qcm14[i]>>epsilon14[i]>>qprimcm14[i]>>thetacm14[i]>>phideg14[i]>>segmaDR14[i];
f15>>qcm15[i]>>epsilon15[i]>>qprimcm15[i]>>thetacm15[i]>>phideg15[i]>>segmaDR15[i];
f16>>qcm16[i]>>epsilon16[i]>>qprimcm16[i]>>thetacm16[i]>>phideg16[i]>>segmaDR16[i];
f17>>qcm17[i]>>epsilon17[i]>>qprimcm17[i]>>thetacm17[i]>>phideg17[i]>>segmaDR17[i];
f18>>qcm18[i]>>epsilon18[i]>>qprimcm18[i]>>thetacm18[i]>>phideg18[i]>>segmaDR18[i];
f19>>qcm19[i]>>epsilon19[i]>>qprimcm19[i]>>thetacm19[i]>>phideg19[i]>>segmaDR19[i];
f20>>qcm20[i]>>epsilon20[i]>>qprimcm20[i]>>thetacm20[i]>>phideg20[i]>>segmaDR20[i];
f21>>qcm21[i]>>epsilon21[i]>>qprimcm21[i]>>thetacm21[i]>>phideg21[i]>>segmaDR21[i];
f22>>qcm22[i]>>epsilon22[i]>>qprimcm22[i]>>thetacm22[i]>>phideg22[i]>>segmaDR22[i];
f23>>qcm23[i]>>epsilon23[i]>>qprimcm23[i]>>thetacm23[i]>>phideg23[i]>>segmaDR23[i];
f24>>qcm24[i]>>epsilon24[i]>>qprimcm24[i]>>thetacm24[i]>>phideg24[i]>>segmaDR24[i];
f25>>qcm25[i]>>epsilon25[i]>>qprimcm25[i]>>thetacm25[i]>>phideg25[i]>>segmaDR25[i];
//cout<<qcm1[i]<<epsilon1[i]<<qprimcm1[i]<<thetacm1[i]<<phideg1[i]<<segmaDR1[i]<<endl;
fexp>>segexp[i]>>esegmaexp[i];
//calcul de Khi2 pour chaque points de sigma
khi1[i]=((segexp[i]-segmaDR1[i])*(segexp[i]-segmaDR1[i]))/(esegmaexp[i]*esegmaexp[i]);
khi2[i]=((segexp[i]-segmaDR2[i])*(segexp[i]-segmaDR2[i]))/(esegmaexp[i]*esegmaexp[i]);
khi3[i]=((segexp[i]-segmaDR3[i])*(segexp[i]-segmaDR3[i]))/(esegmaexp[i]*esegmaexp[i]);
khi4[i]=((segexp[i]-segmaDR4[i])*(segexp[i]-segmaDR4[i]))/(esegmaexp[i]*esegmaexp[i]);
khi5[i]=((segexp[i]-segmaDR5[i])*(segexp[i]-segmaDR5[i]))/(esegmaexp[i]*esegmaexp[i]);
khi6[i]=((segexp[i]-segmaDR6[i])*(segexp[i]-segmaDR6[i]))/(esegmaexp[i]*esegmaexp[i]);
khi7[i]=((segexp[i]-segmaDR7[i])*(segexp[i]-segmaDR7[i]))/(esegmaexp[i]*esegmaexp[i]);
khi8[i]=((segexp[i]-segmaDR8[i])*(segexp[i]-segmaDR8[i]))/(esegmaexp[i]*esegmaexp[i]);
khi9[i]=((segexp[i]-segmaDR9[i])*(segexp[i]-segmaDR9[i]))/(esegmaexp[i]*esegmaexp[i]);
khi10[i]=((segexp[i]-segmaDR10[i])*(segexp[i]-segmaDR10[i]))/(esegmaexp[i]*esegmaexp[i]);
khi11[i]=((segexp[i]-segmaDR11[i])*(segexp[i]-segmaDR11[i]))/(esegmaexp[i]*esegmaexp[i]);
khi12[i]=((segexp[i]-segmaDR12[i])*(segexp[i]-segmaDR12[i]))/(esegmaexp[i]*esegmaexp[i]);
khi13[i]=((segexp[i]-segmaDR13[i])*(segexp[i]-segmaDR13[i]))/(esegmaexp[i]*esegmaexp[i]);
khi14[i]=((segexp[i]-segmaDR14[i])*(segexp[i]-segmaDR14[i]))/(esegmaexp[i]*esegmaexp[i]);
khi15[i]=((segexp[i]-segmaDR15[i])*(segexp[i]-segmaDR15[i]))/(esegmaexp[i]*esegmaexp[i]);
khi16[i]=((segexp[i]-segmaDR16[i])*(segexp[i]-segmaDR16[i]))/(esegmaexp[i]*esegmaexp[i]);
khi17[i]=((segexp[i]-segmaDR17[i])*(segexp[i]-segmaDR17[i]))/(esegmaexp[i]*esegmaexp[i]);
khi18[i]=((segexp[i]-segmaDR18[i])*(segexp[i]-segmaDR18[i]))/(esegmaexp[i]*esegmaexp[i]);
khi19[i]=((segexp[i]-segmaDR19[i])*(segexp[i]-segmaDR19[i]))/(esegmaexp[i]*esegmaexp[i]);
khi20[i]=((segexp[i]-segmaDR20[i])*(segexp[i]-segmaDR20[i]))/(esegmaexp[i]*esegmaexp[i]);
khi21[i]=((segexp[i]-segmaDR21[i])*(segexp[i]-segmaDR21[i]))/(esegmaexp[i]*esegmaexp[i]);
khi22[i]=((segexp[i]-segmaDR22[i])*(segexp[i]-segmaDR22[i]))/(esegmaexp[i]*esegmaexp[i]);
khi23[i]=((segexp[i]-segmaDR23[i])*(segexp[i]-segmaDR23[i]))/(esegmaexp[i]*esegmaexp[i]);
khi24[i]=((segexp[i]-segmaDR24[i])*(segexp[i]-segmaDR24[i]))/(esegmaexp[i]*esegmaexp[i]);
khi25[i]=((segexp[i]-segmaDR25[i])*(segexp[i]-segmaDR25[i]))/(esegmaexp[i]*esegmaexp[i]);
}
// calcul de Khi2 somme
for(int i=0; i<10; i++){
khi11p +=khi1[i];
khi21p +=khi2[i];
khi31 +=khi3[i];
khi41 +=khi4[i];
khi51 +=khi5[i];
khi61 +=khi6[i];
khi71 +=khi7[i];
khi81 +=khi8[i];
khi91 +=khi9[i];
khi101 +=khi10[i];
khi111 +=khi11[i];
khi121 +=khi12[i];
khi131 +=khi13[i];
khi141 +=khi14[i];
khi151 +=khi15[i];
khi161 +=khi16[i];
khi171 +=khi17[i];
khi181 +=khi18[i];
khi191 +=khi19[i];
khi201 +=khi20[i];
khi211 +=khi21[i];
khi221 +=khi22[i];
khi231 +=khi23[i];
khi241 +=khi24[i];
khi251 +=khi25[i];
//cout<<khi1[i]<<endl;
}
//remplissage de fechier txt "out"
cout<<khi11p<<endl;// j'ajoute p ici si car j'ai un tableau a le mem nom que khi11
fcoef<<khi11p<<endl;
fcoef<<khi21p<<endl;
fcoef<<khi31<<endl;
fcoef<<khi41<<endl;
fcoef<<khi51<<endl;
fcoef<<khi61<<endl;
fcoef<<khi71<<endl;
fcoef<<khi81<<endl;
fcoef<<khi91<<endl;
fcoef<<khi101<<endl;
fcoef<<khi111<<endl;
fcoef<<khi121<<endl;
fcoef<<khi131<<endl;
fcoef<<khi141<<endl;
fcoef<<khi151<<endl;
fcoef<<khi161<<endl;
fcoef<<khi171<<endl;
fcoef<<khi181<<endl;
fcoef<<khi191<<endl;
fcoef<<khi201<<endl;
fcoef<<khi211<<endl;
fcoef<<khi221<<endl;
fcoef<<khi231<<endl;
fcoef<<khi241<<endl;
fcoef<<khi251<<endl;
} |
Partager