{

 double khi[10][25];
 double khi11p[25];
 double segexp[10]; 
 double esegmaexp[10];
 double qcm1[10][25], epsilon1[10][25], qprimcm1[10][25], thetacm1[10][25], phideg1[10][25], segmaDR1[10][25],costheta1[10][25];
 double alpha[25], beta[25];
ofstream fcoef("val.txt");//txt sortie "OUT"
ifstream fexp(Form("../../cross-section-q2-05-oct2013/Results/section1.txt"));//txt segmaexp
ifstream f(Form("../1-code-dr/alphabeta.txt"));//txt alpha beta
 
 for(int j=0; j<10; j++){fexp>>segexp[j]>>esegmaexp[j];} 
 
 for(int i=0; i<25; i++){
      f>>alpha[i]>>beta[i];
      //cout<<alpha[i]<<"  "<<beta[i]<<endl;
      ifstream file(Form("../1-code-dr/fort%d_%d.txt",i,i));
      //cout<<Form("../1-code-dr/fort%d_%d.txt",i,i)<<endl; //pour etre sur que je lis mes fichiers
    
 for(int j=0; j<10; j++){
     //lecture des fichiers txt
    file>>qcm1[j][i]>>epsilon1[j][i]>>qprimcm1[j][i]>>thetacm1[j][i]>>phideg1[j][i]>>segmaDR1[j][i];
    //cout<<qcm1[j][i]<<epsilon1[j][i]<<qprimcm1[j][i]<<thetacm1[j][i]<<phideg1[j][i]<<segmaDR1[j][i]<<endl;
   //calcul de Khi2 pour chaque points de sigma
    khi[j][i]=((segexp[j]-segmaDR1[j][i])*(segexp[j]-segmaDR1[j][i]))/(esegmaexp[j]*esegmaexp[j]);
    //cout<<khi[j][i]<<endl;
    }}

 // calcul de Khi2 somme
 for(int j=0; j<10; j++){for(int i=0; i<25; i++){khi11p[i] +=khi[j][i];}}
 for(int i=0; i<25; i++){cout<<khi11p[i]<<endl;}
//remplissage de fechier txt "out"
 for(int i=0; i<25; i++){fcoef<<khi11p[i]<<endl;}


}