EXP1803

{

        gSystem->Load("libeloss"); // Loading of the library for the calculating of the energy losses in the matter.

        TH2D *h_sq20 = new TH2D("h_sq20","sq20 thikness",8,1,9,15,1,16);

        char dname[32]; 

        int i,j;

        double abc[10];

        double a[8][15],d_sq20[8][15];

        double de[8][15],e[8][15];

        double sq20=18;

        TELoss eloss(1,2.33); //  TELoss 

        eloss.AddEL(14,28); //  TELoss.  -- Si

        eloss.SetZP(2,4); // (Z) (A) 

        eloss.SetEtab(1000,0.1,10); //  1000  0.1  10 

        FILE *pf1;

         pf1=fopen("sq20_calib_03.txt","r");//File with positions of the 7.63 alpha peak after passing through the 20mkm of Si

        for(j=0; j<15; j++)

          {

          fscanf(pf1,"%lf %lf %lf %lf %lf %lf %lf %lf\n",&a[0][j],&a[1][j],&a[2][j],&a[3][j],&a[4][j],&a[5][j],&a[6][j],&a[7][j]);

          printf("%lf %lf %lf %lf %lf %lf %lf %lf\n",a[0][j],a[1][j],a[2][j],a[3][j],a[4][j],a[5][j],a[6][j],a[7][j]);

          }

        fclose(pf1);  

        FILE *FW = fopen("sq20_de_04.txt","wt"); // File with dE of the 7.63 alpha after passing through the 20mkm of Si

        for(j=0; j<15; j++)

          {

              for(i=0; i<8; i++)

                  {

                  de[i][j]=7.63-a[i][j];        

                  }

              fprintf(FW,"%lf %lf %lf %lf %lf %lf %lf %lf\n",de[0][j],de[1][j],de[2][j],de[3][j],de[4][j],de[5][j],de[6][j],de[7][j]);

          }

        fclose(FW);

        FILE *FW1 = fopen("sq20_thickness.txt","wt"); // File with thicknesses of SQ20

        FILE *FW2 = fopen("sq20_list_04.txt","wt"); // List of calculations

        for(j=0; j<15; j++)

          {

              for(i=0; i<8; i++)

                  {

                        sq20=18;

                          do 

                        {        

                        e[i][j] = 7.63-eloss.GetE(7.63, sq20);

                        fprintf(FW2,"%lf %lf %lf %d %d\n",e[i][j],de[i][j],sq20,i,j);                        

                        sq20=sq20+0.01;

                        }while((e[i][j]-de[i][j])<=0.01);

                        d_sq20[i][j]=sq20;

                        h_sq20->SetBinContent(i+1,j+1,sq20);

                  }

              fprintf(FW1,"%lf %lf %lf %lf %lf %lf %lf %lf\n",d_sq20[0][j],d_sq20[1][j],d_sq20[2][j],d_sq20[3][j],d_sq20[4][j],d_sq20[5][j],d_sq20[6][j],d_sq20[7][j]);

          }

        FILE *FW3 = fopen("sq20_coef.txt","wt"); // File with coefficients from thickness

        for(j=0; j<15; j++)

          {

          fprintf(FW3,"%lf %lf %lf %lf %lf %lf %lf %lf\n",22/d_sq20[0][j],22/d_sq20[1][j],22/d_sq20[2][j],22/d_sq20[3][j],22/d_sq20[4][j],22/d_sq20[5][j],22/d_sq20[6][j],22/d_sq20[7][j]);

          }

        h_sq20->Draw("lego2z");

        fclose(FW3);

        fclose(FW2);

         fclose(FW1);

}