void plot_gain_nim(){ gROOT->Reset(); gROOT->SetStyle("Plain"); gROOT->ForceStyle(); gStyle->SetOptStat(0); gStyle->SetOptFit(0); gStyle->SetTitleX(0.1f); gStyle->SetTitleW(0.8f); gStyle->SetTextSize(1.5); gStyle->SetTitleFont(132,"xy"); gStyle->SetLabelFont(132,"xy"); gStyle->SetLabelSize(0.055,"xy"); gStyle->SetTitleSize(0.062,"xy"); gStyle->SetTitleOffset(1.1,"x"); gStyle->SetTitleOffset(0.9,"y"); gStyle->SetTitleBorderSize(0); gStyle->SetOptFit(11); FILE *file1 = fopen("my_test.txt","r"); int id[200]; float x[200]; float ex[200]; float int_led1[200]; float eint_led1[200]; float int_alpha1[200]; float eint_alpha1[200]; float ratio_ref[200]; float eratio_ref[200]; float ch[145][200]; float ech[145][200]; float ch_norm[145]; int ii = 0; float a[145]; float ea[145]; memset(ex,0,sizeof(ex)); memset(ratio_ref,0,sizeof(ratio_ref)); memset(eratio_ref,0,sizeof(eratio_ref)); char s1[] = "%d %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f %f \n"; while (!feof(file1)){ // fscanf(file1,"%d %f %f %f %f \n",&id[ii], &int_led1[ii], &eint_led1[ii], &int_alpha1[ii], &eint_alpha1[ii]); int aa = fscanf(file1,s1,&id[ii], &int_led1[ii], &eint_led1[ii], &int_alpha1[ii], &eint_alpha1[ii], &a[0], &ea[0], &a[1], &ea[1], &a[2], &ea[2],&a[3],&ea[3],&a[4],&ea[4], &a[5], &ea[5], &a[6], &ea[6], &a[7], &ea[7], &a[8], &ea[8], &a[9], &ea[9], &a[10], &ea[10], &a[11], &ea[11], &a[12],&ea[12],&a[13],&ea[13],&a[14],&ea[14],&a[15],&ea[15],&a[16],&ea[16],&a[17],&ea[17],&a[18],&ea[18],&a[19],&ea[19],&a[20],&ea[20],&a[21],&ea[21],&a[22],&ea[22],&a[23],&ea[23], &a[24],&ea[24],&a[25],&ea[25],&a[26],&ea[26],&a[27],&ea[27],&a[28],&ea[28],&a[29],&ea[29],&a[30],&ea[30],&a[31],&ea[31],&a[32],&ea[32],&a[33],&ea[33],&a[34],&ea[34],&a[35],&ea[35], &a[36],&ea[36],&a[37],&ea[37],&a[38],&ea[38],&a[39],&ea[39],&a[40],&ea[40],&a[41],&ea[41],&a[42],&ea[42],&a[43],&ea[43],&a[44],&ea[44],&a[45],&ea[45],&a[46],&ea[46],&a[47],&ea[47], &a[48],&ea[48],&a[49],&ea[49],&a[50],&ea[50],&a[51],&ea[51],&a[52],&ea[52],&a[53],&ea[53],&a[54],&ea[54],&a[55],&ea[55],&a[56],&ea[56],&a[57],&ea[57],&a[58],&ea[58],&a[59],&ea[59], &a[60],&ea[60],&a[61],&ea[61],&a[62],&ea[62],&a[63],&ea[63],&a[64],&ea[64],&a[65],&ea[65],&a[66],&ea[66],&a[67],&ea[67],&a[68],&ea[68],&a[69],&ea[69], &a[70],&ea[70],&a[71],&ea[71],&a[72],&ea[72],&a[73],&ea[73],&a[74],&ea[74],&a[75],&ea[75],&a[76],&ea[76],&a[77],&ea[77],&a[78],&ea[78],&a[79],&ea[79], &a[70],&ea[80],&a[81],&ea[81],&a[82],&ea[82],&a[83],&ea[83],&a[84],&ea[84],&a[85],&ea[85],&a[86],&ea[86],&a[87],&ea[87],&a[88],&ea[88],&a[89],&ea[89], &a[90],&ea[90],&a[91],&ea[91],&a[92],&ea[92],&a[93],&ea[93],&a[94],&ea[94],&a[95],&ea[95],&a[96],&ea[96],&a[97],&ea[97],&a[98],&ea[98],&a[99],&ea[99], &a[100],&ea[100],&a[101],&ea[101],&a[102],&ea[102],&a[103],&ea[103],&a[104],&ea[104],&a[105],&ea[105],&a[106],&ea[106],&a[107],&ea[107],&a[108],&ea[108],&a[109],&ea[109], &a[110],&ea[110],&a[111],&ea[111],&a[112],&ea[112],&a[113],&ea[113],&a[114],&ea[114],&a[115],&ea[115],&a[116],&ea[116],&a[117],&ea[117],&a[118],&ea[118],&a[119],&ea[119], &a[120],&ea[120],&a[121],&ea[121],&a[122],&ea[122],&a[123],&ea[123],&a[124],&ea[124],&a[125],&ea[125],&a[126],&ea[126],&a[127],&ea[127],&a[128],&ea[128],&a[129],&ea[129], &a[130],&ea[130],&a[131],&ea[131],&a[132],&ea[132],&a[133],&ea[133],&a[134],&ea[134],&a[135],&ea[135],&a[136],&ea[136],&a[137],&ea[137],&a[138],&ea[138],&a[139],&ea[139], &a[140],&ea[140],&a[141],&ea[141],&a[142],&ea[142],&a[143],&ea[143]); double tmp = int_led1[ii]/int_alpha1[ii]; double tmp1 = sqrt(pow(eint_led1[ii]/int_led1[ii],2) + pow(eint_alpha1[ii]/int_alpha1[ii],2)); if(ii >= 21) tmp *= 1.002; ratio_ref[ii] = tmp; eratio_ref[ii] = tmp*tmp1; cout << endl; cout << id[ii] << " Amp = " << int_led1[ii] << " " << eint_led1[ii] << " Alpha " << int_alpha1[ii] << endl; cout << " Ratio = " << ratio_ref[ii] << " Eratio = " << eratio_ref[ii] << " " << ii << endl; // x[ii] = ii + 1; x[ii] = 47 - 47*(23 - ii)/23 + 1; cout << x[ii] << endl; ii++; } double norm = ratio_ref[0]; for(int jj = 0; jj < ii; jj++){ ratio_ref[jj] = ratio_ref[jj]/norm; eratio_ref[jj] = eratio_ref[jj]/norm; } c1 = new TCanvas("Reference","Reference", 200, 10, 600, 400); gPad->SetBottomMargin(0.15); gPad->SetTopMargin(0.05); gPad->SetLeftMargin(0.11); gPad->SetRightMargin(0.05); TGraphErrors *gr = new TGraphErrors(ii+1, x, ratio_ref, ex, eratio_ref ); gr->SetMarkerStyle(20); gr->SetMarkerSize(0.9); gr->SetMarkerColor(4); // gr->GetXaxis()->SetTitle("Run period"); gr->GetXaxis()->SetTitle("Days"); gr->GetYaxis()->SetTitle("Ratio"); gr->SetTitle(""); gr->SetMinimum(0.97); gr->SetMaximum(1.03); gr->Draw("AP"); // TAxis *aa = gr->GetXaxis(); // aa->SetNdivisions(-501); // aa->ChangeLabel(1,-1,-1,-1,-1,-1,"1"); // aa->ChangeLabel(2,-1,-1,-1,-1,-1,"48"); gr->Draw("AP"); }