code1.cpp

void DRTileMPV(int a, float cc) {
    TFile *file;
    switch (a) {
          case 22:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz22mm_e200.root");
          goto X;
          }
          case 24:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz24mm_e200.root");
          goto X;
          }
          case 26:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz26mm_e200.root");
          goto X;
          }
          case 28:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz28mm_e200.root");
          goto X;
          }
          case 210:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz210mm_e200.root");
          goto X;
          }
          case 212:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz212mm_e200.root");
          goto X;
          }
          case 214:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz214mm_e200.root");
          goto X;
          }
          case 42:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz42mm_e200.root");
          goto X;
          }
          case 44:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz44mm_e200.root");
          goto X;
          }
          case 46:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz46mm_e200.root");
          goto X;
          }
          case 48:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz48mm_e200.root");
          goto X;
          }
          case 410:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz410mm_e200.root");
          goto X;
          }
          case 412:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz412mm_e200.root");
          goto X;
          }
          case 414:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz414mm_e200.root");
          goto X;
          }
          case 62:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz62mm_e200.root");
          goto X;
          }
          case 64:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz64mm_e200.root");
          goto X;
          }
          case 66:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz66mm_e200.root");
          goto X;
          }
          case 68:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz68mm_e200.root");
          goto X;
          }
          case 610:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz610mm_e200.root");
          goto X;
          }
          case 612:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz612mm_e200.root");
          goto X;
          }
          case 614:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz614mm_e200.root");
          goto X;
          }
          case 82:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz82mm_e200.root");
          goto X;
          }
          case 84:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz84mm_e200.root");
          goto X;
          }
          case 86:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz86mm_e200.root");
          goto X;
          }
          case 88:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz88mm_e200.root");
          goto X;
          }
          case 810:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz810mm_e200.root");
          goto X;
          }
          case 812:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz812mm_e200.root");
          goto X;
          }
          case 814:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz814mm_e200.root");
          goto X;
          }
          case 102:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz102mm_e200.root");
          goto X;
          }
          case 104:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz104mm_e200.root");
          goto X;
          }
          case 106:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz106mm_e200.root");
          goto X;
          }
          case 108:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz108mm_e200.root");
          goto X;
          }
          case 1010:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz1010mm_e200.root");
          goto X;
          }
          case 1012:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz1012mm_e200.root");
          goto X;
          }
          case 1014:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz1014mm_e200.root");
          goto X;
          }
          case 122:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz122mm_e200.root");
          goto X;
          }
          case 124:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz124mm_e200.root");
          goto X;
          }
          case 126:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz126mm_e200.root");
          goto X;
          }
          case 128:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz128mm_e200.root");
          goto X;
          }
          case 1210:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz1210mm_e200.root");
          goto X;
          }
          case 1212:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz1212mm_e200.root");
          goto X;
          }
          case 1214:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz1214mm_e200.root");
          goto X;
          }
          case 142:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz142mm_e200.root");
          goto X;
          }
          case 144:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz144mm_e200.root");
          goto X;
          }
          case 146:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz146mm_e200.root");
          goto X;
          }
          case 148:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz148mm_e200.root");
          goto X;
          }
          case 1410:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz1410mm_e200.root");
          goto X;
          }
          case 1412:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz1412mm_e200.root");
          goto X;
          }
          case 1414:
          {
          file = new TFile("dimple_scint3_box15_gap01_xz1414mm_e200.root");
          goto X;
          }
          default:
          {
          std::cout << "Error in \"a\" " << std::endl;
          goto X;
          }

    }
    X : 
    TH1F *hist;
    TF1 *fitFunc1;
    TF1 *fitFunc2;

    TCanvas *c1 = new TCanvas("c1", "c1", 1);

    printf("File open... ");
    if (!file) {
        printf("[FAIL]\n");
        return 0;
    } else printf("[OK]\n");

    printf("Loading histogram... ");
    hist = (TH1F *)file->Get("fCounterSipm");
    if (!hist) {
        printf("[FAIL]\n");
        return 0;
    } else printf("[OK]\n");

      fitFunc1 = new TF1("fitFunc1", "gaus");

    hist->GetYaxis()->SetTitle("Number of events");
    hist->GetXaxis()->SetTitle("Number of photons");
    hist->Fit(fitFunc1, "0");
    float mean1 = fitFunc1->GetParameter(1);
    float sigma1 = fitFunc1->GetParameter(2);
    std::cout << "Mean1 : " << mean1 << std::endl;
    std::cout << "Sigma1 : " << sigma1 << std::endl;


   fitFunc2 = new TF1("fitFunc2", "gaus", mean1 - cc*sigma1, mean1 + cc*sigma1);
   hist->Fit(fitFunc2, "R");
   float mean2 = fitFunc2->GetParameter(1);
   float sigma2 = fitFunc2->GetParameter(2);
   std::cout << "Mean2 : " << mean2 << std::endl;
   std::cout << "Sigma2 : " << sigma2 << std::endl;

   int flag = 0;
   double chi2 = 0;
   chi2 = fitFunc2->GetChisquare() / fitFunc2->GetNDF();
   if (chi2 > 5) flag = 1;
    
   c1->SaveAs("pic.jpg");

//   hist->Draw();

   FILE *out = fopen("out.txt", "a+");
       printf("File open... ");
       if (!out) {
          printf("[FAIL]\n");
          return 0;
       } else printf("[OK]\n");
   fprintf(out,"%i %E %E %E %i\n", a, mean2, sigma2, chi2, flag);
  // out << mean2 << " " << sigma2 << " " << chi2 << " " << flag << std::endl;
  //  fclose(out);

}