//dEdx Truncation Study //Test of different Truncations #include "TGraph.h" #include "TGraphErrors.h" #include "TString.h" #include "TMultiGraph.h" #include "TLegend.h" #include "TCanvas.h" #include "TLatex.h" #include "TAxis.h" #include "TROOT.h" #include "TFile.h" #include #include #include #include struct FitResults { unsigned int run; unsigned int data; unsigned int lowtrunc; unsigned int hightrunc; double meanel; double meanelerr; double sigmael; double sigmaelerr; double meanpi; double meanpierr; double sigmapi; double sigmapierr; double separation1; double separation1err; double separation2; double separation2err; }; std::string inttostring(Int_t input) { std::string s; std::stringstream out; out << input; s = out.str(); return s; } //temporary solution: a map - better with root trees std::map fitresultsmap; void LoadData(TString datafile, unsigned int datasetnumber){ ifstream instr; instr.open(datafile.Data()); if (instr.is_open()) std::cout<<"opened file: "<>run>>data>>lowtrunc>>hightrunc>>meanel>>meanelerr>>sigmael>>sigmaelerr>>meanpi>>meanpierr>>sigmapi>>sigmapierr>>separation1>>separation1err>>separation2>>separation2err; if(instr.eof()) break; cout << "reading data: " << run <<" "<< data << " "<< lowtrunc <<" "<< hightrunc<ProcessLine(".x ~/rootlogon_Bernhard.C"); gROOT->SetStyle("col"); gROOT->ForceStyle(); } LoadData("/nfs/hicran/project/panda/SIM/pgadow/fopiROOT_trunk/macro/ALICE/plots/dEdx/FITDATA.dat", datasetnumber); TString dirname = "/nfs/mds/data/tpc/alice/ps_test_beam/dEdx/fitresults/"; TString filename = "dEdxFitResults"; filename+= "_"+inttostring(datasetnumber)+"_TruncationStudy"; TString savefilename = dirname + filename + ".root"; TFile *savefile = new TFile(savefilename.Data(),"RECREATE"); savefile->cd(); //################################################################################################################################################ //1 GeV - Beam ################################################################################################################################### //################################################################################################################################################ //################################################################################################################################################ //Standard Gain: Mean ############################################################################################################################ //################################################################################################################################################ TCanvas *cstdmeanel = new TCanvas("cstdmeanel","dEdx standard gain mean electrons",1000,700); cstdmeanel->SetFillColor(00); cstdmeanel->SetGrid(); cstdmeanel->cd(); TMultiGraph *mgstdmeanel = new TMultiGraph(); TMultiGraph *mgstdmeanpi = new TMultiGraph(); const int nstdmean = 5; double gainstdmean[nstdmean] = {69., 70., 71., 72., 73.}; double meanel0025[nstdmean]; double meanel0525[nstdmean]; double meanel1025[nstdmean]; double meanel0030[nstdmean]; double meanel0530[nstdmean]; double meanel1030[nstdmean]; double meanpi0025[nstdmean]; double meanpi0525[nstdmean]; double meanpi1025[nstdmean]; double meanpi0030[nstdmean]; double meanpi0530[nstdmean]; double meanpi1030[nstdmean]; double dmeanel0025[nstdmean]; double dmeanel0525[nstdmean]; double dmeanel1025[nstdmean]; double dmeanel0030[nstdmean]; double dmeanel0530[nstdmean]; double dmeanel1030[nstdmean]; double dmeanpi0025[nstdmean]; double dmeanpi0525[nstdmean]; double dmeanpi1025[nstdmean]; double dmeanpi0030[nstdmean]; double dmeanpi0530[nstdmean]; double dmeanpi1030[nstdmean]; meanel0025[0] = (fitresultsmap[6810025]).meanel; meanel0025[1] = (fitresultsmap[6820025]).meanel; meanel0025[2] = (fitresultsmap[6830025]).meanel; meanel0025[3] = (fitresultsmap[6840025]).meanel; meanel0025[4] = (fitresultsmap[6850025]).meanel; dmeanel0025[0] = (fitresultsmap[6810025]).meanelerr; dmeanel0025[1] = (fitresultsmap[6820025]).meanelerr; dmeanel0025[2] = (fitresultsmap[6830025]).meanelerr; dmeanel0025[3] = (fitresultsmap[6840025]).meanelerr; dmeanel0025[4] = (fitresultsmap[6850025]).meanelerr; meanel0525[0] = (fitresultsmap[6810525]).meanel; meanel0525[1] = (fitresultsmap[6820525]).meanel; meanel0525[2] = (fitresultsmap[6830525]).meanel; meanel0525[3] = (fitresultsmap[6840525]).meanel; meanel0525[4] = (fitresultsmap[6850525]).meanel; dmeanel0525[0] = (fitresultsmap[6810525]).meanelerr; dmeanel0525[1] = (fitresultsmap[6820525]).meanelerr; dmeanel0525[2] = (fitresultsmap[6830525]).meanelerr; dmeanel0525[3] = (fitresultsmap[6840525]).meanelerr; dmeanel0525[4] = (fitresultsmap[6850525]).meanelerr; meanel1025[0] = (fitresultsmap[6811025]).meanel; meanel1025[1] = (fitresultsmap[6821025]).meanel; meanel1025[2] = (fitresultsmap[6831025]).meanel; meanel1025[3] = (fitresultsmap[6841025]).meanel; meanel1025[4] = (fitresultsmap[6851025]).meanel; dmeanel1025[0] = (fitresultsmap[6811025]).meanelerr; dmeanel1025[1] = (fitresultsmap[6821025]).meanelerr; dmeanel1025[2] = (fitresultsmap[6831025]).meanelerr; dmeanel1025[3] = (fitresultsmap[6841025]).meanelerr; dmeanel1025[4] = (fitresultsmap[6851025]).meanelerr; meanel0030[0] = (fitresultsmap[6810030]).meanel; meanel0030[1] = (fitresultsmap[6820030]).meanel; meanel0030[2] = (fitresultsmap[6830030]).meanel; meanel0030[3] = (fitresultsmap[6840030]).meanel; meanel0030[4] = (fitresultsmap[6850030]).meanel; dmeanel0030[0] = (fitresultsmap[6810030]).meanelerr; dmeanel0030[1] = (fitresultsmap[6820030]).meanelerr; dmeanel0030[2] = (fitresultsmap[6830030]).meanelerr; dmeanel0030[3] = (fitresultsmap[6840030]).meanelerr; dmeanel0030[4] = (fitresultsmap[6850030]).meanelerr; meanel0530[0] = (fitresultsmap[6810530]).meanel; meanel0530[1] = (fitresultsmap[6820530]).meanel; meanel0530[2] = (fitresultsmap[6830530]).meanel; meanel0530[3] = (fitresultsmap[6840530]).meanel; meanel0530[4] = (fitresultsmap[6850530]).meanel; dmeanel0530[0] = (fitresultsmap[6810530]).meanelerr; dmeanel0530[1] = (fitresultsmap[6820530]).meanelerr; dmeanel0530[2] = (fitresultsmap[6830530]).meanelerr; dmeanel0530[3] = (fitresultsmap[6840530]).meanelerr; dmeanel0530[4] = (fitresultsmap[6850530]).meanelerr; meanel1030[0] = (fitresultsmap[6811030]).meanel; meanel1030[1] = (fitresultsmap[6821030]).meanel; meanel1030[2] = (fitresultsmap[6831030]).meanel; meanel1030[3] = (fitresultsmap[6841030]).meanel; meanel1030[4] = (fitresultsmap[6851030]).meanel; dmeanel1030[0] = (fitresultsmap[6811030]).meanelerr; dmeanel1030[1] = (fitresultsmap[6821030]).meanelerr; dmeanel1030[2] = (fitresultsmap[6831030]).meanelerr; dmeanel1030[3] = (fitresultsmap[6841030]).meanelerr; dmeanel1030[4] = (fitresultsmap[6851030]).meanelerr; meanpi0025[0] = (fitresultsmap[6810025]).meanpi; meanpi0025[1] = (fitresultsmap[6820025]).meanpi; meanpi0025[2] = (fitresultsmap[6830025]).meanpi; meanpi0025[3] = (fitresultsmap[6840025]).meanpi; meanpi0025[4] = (fitresultsmap[6850025]).meanpi; dmeanpi0025[0] = (fitresultsmap[6810025]).meanpierr; dmeanpi0025[1] = (fitresultsmap[6820025]).meanpierr; dmeanpi0025[2] = (fitresultsmap[6830025]).meanpierr; dmeanpi0025[3] = (fitresultsmap[6840025]).meanpierr; dmeanpi0025[4] = (fitresultsmap[6850025]).meanpierr; meanpi0525[0] = (fitresultsmap[6810525]).meanpi; meanpi0525[1] = (fitresultsmap[6820525]).meanpi; meanpi0525[2] = (fitresultsmap[6830525]).meanpi; meanpi0525[3] = (fitresultsmap[6840525]).meanpi; meanpi0525[4] = (fitresultsmap[6850525]).meanpi; dmeanpi0525[0] = (fitresultsmap[6810525]).meanpierr; dmeanpi0525[1] = (fitresultsmap[6820525]).meanpierr; dmeanpi0525[2] = (fitresultsmap[6830525]).meanpierr; dmeanpi0525[3] = (fitresultsmap[6840525]).meanpierr; dmeanpi0525[4] = (fitresultsmap[6850525]).meanpierr; meanpi1025[0] = (fitresultsmap[6811025]).meanpi; meanpi1025[1] = (fitresultsmap[6821025]).meanpi; meanpi1025[2] = (fitresultsmap[6831025]).meanpi; meanpi1025[3] = (fitresultsmap[6841025]).meanpi; meanpi1025[4] = (fitresultsmap[6851025]).meanpi; dmeanpi1025[0] = (fitresultsmap[6811025]).meanpierr; dmeanpi1025[1] = (fitresultsmap[6821025]).meanpierr; dmeanpi1025[2] = (fitresultsmap[6831025]).meanpierr; dmeanpi1025[3] = (fitresultsmap[6841025]).meanpierr; dmeanpi1025[4] = (fitresultsmap[6851025]).meanpierr; meanpi0030[0] = (fitresultsmap[6810030]).meanpi; meanpi0030[1] = (fitresultsmap[6820030]).meanpi; meanpi0030[2] = (fitresultsmap[6830030]).meanpi; meanpi0030[3] = (fitresultsmap[6840030]).meanpi; meanpi0030[4] = (fitresultsmap[6850030]).meanpi; dmeanpi0030[0] = (fitresultsmap[6810030]).meanpierr; dmeanpi0030[1] = (fitresultsmap[6820030]).meanpierr; dmeanpi0030[2] = (fitresultsmap[6830030]).meanpierr; dmeanpi0030[3] = (fitresultsmap[6840030]).meanpierr; dmeanpi0030[4] = (fitresultsmap[6850030]).meanpierr; meanpi0530[0] = (fitresultsmap[6810530]).meanpi; meanpi0530[1] = (fitresultsmap[6820530]).meanpi; meanpi0530[2] = (fitresultsmap[6830530]).meanpi; meanpi0530[3] = (fitresultsmap[6840530]).meanpi; meanpi0530[4] = (fitresultsmap[6850530]).meanpi; dmeanpi0530[0] = (fitresultsmap[6810530]).meanpierr; dmeanpi0530[1] = (fitresultsmap[6820530]).meanpierr; dmeanpi0530[2] = (fitresultsmap[6830530]).meanpierr; dmeanpi0530[3] = (fitresultsmap[6840530]).meanpierr; dmeanpi0530[4] = (fitresultsmap[6850530]).meanpierr; meanpi1030[0] = (fitresultsmap[6811030]).meanpi; meanpi1030[1] = (fitresultsmap[6821030]).meanpi; meanpi1030[2] = (fitresultsmap[6831030]).meanpi; meanpi1030[3] = (fitresultsmap[6841030]).meanpi; meanpi1030[4] = (fitresultsmap[6851030]).meanpi; dmeanpi1030[0] = (fitresultsmap[6811030]).meanpierr; dmeanpi1030[1] = (fitresultsmap[6821030]).meanpierr; dmeanpi1030[2] = (fitresultsmap[6831030]).meanpierr; dmeanpi1030[3] = (fitresultsmap[6841030]).meanpierr; dmeanpi1030[4] = (fitresultsmap[6851030]).meanpierr; //0025 TGraphErrors* grmeanel0025 = new TGraphErrors(nstdmean,gainstdmean,meanel0025,0,dmeanel0025); grmeanel0025->SetLineColor(1); grmeanel0025->SetLineWidth(1); grmeanel0025->SetMarkerColor(1); grmeanel0025->SetMarkerStyle(20); grmeanel0025->SetTitle("electrons 0025"); mgstdmeanel->Add(grmeanel0025); TGraphErrors* grmeanpi0025 = new TGraphErrors(nstdmean,gainstdmean,meanpi0025,0,dmeanpi0025); grmeanpi0025->SetLineColor(1); grmeanpi0025->SetLineWidth(1); grmeanpi0025->SetMarkerColor(1); grmeanpi0025->SetMarkerStyle(20); grmeanpi0025->SetTitle("pions 0025"); mgstdmeanpi->Add(grmeanpi0025); //0525 TGraphErrors* grmeanel0525 = new TGraphErrors(nstdmean,gainstdmean,meanel0525,0,dmeanel0525); grmeanel0525->SetLineColor(2); grmeanel0525->SetLineWidth(1); grmeanel0525->SetMarkerColor(2); grmeanel0525->SetMarkerStyle(20); grmeanel0525->SetTitle("electrons 0525"); mgstdmeanel->Add(grmeanel0525); TGraphErrors* grmeanpi0525 = new TGraphErrors(nstdmean,gainstdmean,meanpi0525,0,dmeanpi0525); grmeanpi0525->SetLineColor(2); grmeanpi0525->SetLineWidth(1); grmeanpi0525->SetMarkerColor(2); grmeanpi0525->SetMarkerStyle(20); grmeanpi0525->SetTitle("pions 0525"); mgstdmeanpi->Add(grmeanpi0525); //1025 TGraphErrors* grmeanel1025 = new TGraphErrors(nstdmean,gainstdmean,meanel1025,0,dmeanel1025); grmeanel1025->SetLineColor(3); grmeanel1025->SetLineWidth(1); grmeanel1025->SetMarkerColor(3); grmeanel1025->SetMarkerStyle(24); grmeanel1025->SetTitle("electrons 1025"); mgstdmeanel->Add(grmeanel1025); TGraphErrors* grmeanpi1025 = new TGraphErrors(nstdmean,gainstdmean,meanpi1025,0,dmeanpi1025); grmeanpi1025->SetLineColor(3); grmeanpi1025->SetLineWidth(1); grmeanpi1025->SetMarkerColor(3); grmeanpi1025->SetMarkerStyle(20); grmeanpi1025->SetTitle("pions 1025"); mgstdmeanpi->Add(grmeanpi1025); //0030 TGraphErrors* grmeanel0030 = new TGraphErrors(nstdmean,gainstdmean,meanel0030,0,dmeanel0030); grmeanel0030->SetLineColor(4); grmeanel0030->SetLineWidth(1); grmeanel0030->SetMarkerColor(4); grmeanel0030->SetMarkerStyle(20); grmeanel0030->SetTitle("electrons 0030"); mgstdmeanel->Add(grmeanel0030); TGraphErrors* grmeanpi0030 = new TGraphErrors(nstdmean,gainstdmean,meanpi0030,0,dmeanpi0030); grmeanpi0030->SetLineColor(4); grmeanpi0030->SetLineWidth(1); grmeanpi0030->SetMarkerColor(4); grmeanpi0030->SetMarkerStyle(20); grmeanpi0030->SetTitle("pions 0030"); mgstdmeanpi->Add(grmeanpi0030); //0530 TGraphErrors* grmeanel0530 = new TGraphErrors(nstdmean,gainstdmean,meanel0530,0,dmeanel0530); grmeanel0530->SetLineColor(7); grmeanel0530->SetLineWidth(1); grmeanel0530->SetMarkerColor(7); grmeanel0530->SetMarkerStyle(20); grmeanel0530->SetTitle("electrons 0530"); mgstdmeanel->Add(grmeanel0530); TGraphErrors* grmeanpi0530 = new TGraphErrors(nstdmean,gainstdmean,meanpi0530,0,dmeanpi0530); grmeanpi0530->SetLineColor(7); grmeanpi0530->SetLineWidth(1); grmeanpi0530->SetMarkerColor(7); grmeanpi0530->SetMarkerStyle(20); grmeanpi0530->SetTitle("pions 0530"); mgstdmeanpi->Add(grmeanpi0530); //1030 TGraphErrors* grmeanel1030 = new TGraphErrors(nstdmean,gainstdmean,meanel1030,0,dmeanel1030); grmeanel1030->SetLineColor(6); grmeanel1030->SetLineWidth(1); grmeanel1030->SetMarkerColor(6); grmeanel1030->SetMarkerStyle(20); grmeanel1030->SetTitle("electrons 1030"); mgstdmeanel->Add(grmeanel1030); TGraphErrors* grmeanpi1030 = new TGraphErrors(nstdmean,gainstdmean,meanpi1030,0,dmeanpi1030); grmeanpi1030->SetLineColor(6); grmeanpi1030->SetLineWidth(1); grmeanpi1030->SetMarkerColor(6); grmeanpi1030->SetMarkerStyle(20); grmeanpi1030->SetTitle("pions 1030"); mgstdmeanpi->Add(grmeanpi1030); mgstdmeanel->Draw("ap"); mgstdmeanel->GetXaxis()->SetTitle("HV scaling factor (%)"); mgstdmeanel->GetYaxis()->SetTitle(" electrons (a.u.)"); mgstdmeanel->GetXaxis()->SetLimits(68,74); mgstdmeanel->GetYaxis()->SetRangeUser(0.001,330); TLegend *legstdmeanel=new TLegend(0.2,0.6,0.5,.9); legstdmeanel->SetBorderSize(1); legstdmeanel->SetFillColor(10); legstdmeanel->AddEntry(grmeanel0025,"truncations [0,75]","p"); legstdmeanel->AddEntry(grmeanel0525,"truncations [5,75]","p"); legstdmeanel->AddEntry(grmeanel1025,"truncations [10,75]","p"); legstdmeanel->AddEntry(grmeanel0030,"truncations [0,70]","p"); legstdmeanel->AddEntry(grmeanel0530,"truncations [5,70]","p"); legstdmeanel->AddEntry(grmeanel1030,"truncations [10,70]","p"); legstdmeanel->Draw("same"); cstdmeanel->Update(); cstdmeanel->Modified(); cstdmeanel->Write(); cstdmeanel->SaveAs(dirname+filename+"StdMeanEl.pdf"); TCanvas *cstdmeanpi = new TCanvas("cstdmeanpi","dEdx standard gain mean pions",1000,700); cstdmeanpi->SetFillColor(00); cstdmeanpi->SetGrid(); cstdmeanpi->cd(); mgstdmeanpi->Draw("ap"); mgstdmeanpi->GetXaxis()->SetTitle("HV scaling factor (%)"); mgstdmeanpi->GetYaxis()->SetTitle(" pions (a.u.)"); mgstdmeanpi->GetXaxis()->SetLimits(68,74); mgstdmeanpi->GetYaxis()->SetRangeUser(0.001,220); TLegend *legstdmeanpi=new TLegend(0.2,0.6,0.5,.9); legstdmeanpi->SetBorderSize(1); legstdmeanpi->SetFillColor(10); legstdmeanpi->AddEntry(grmeanpi0025,"truncations [0,75]","p"); legstdmeanpi->AddEntry(grmeanpi0525,"truncations [5,75]","p"); legstdmeanpi->AddEntry(grmeanpi1025,"truncations [10,75]","p"); legstdmeanpi->AddEntry(grmeanpi0030,"truncations [0,70]","p"); legstdmeanpi->AddEntry(grmeanpi0530,"truncations [5,70]","p"); legstdmeanpi->AddEntry(grmeanpi1030,"truncations [10,70]","p"); legstdmeanpi->Draw("same"); cstdmeanpi->Update(); cstdmeanpi->Modified(); cstdmeanpi->Write(); cstdmeanpi->SaveAs(dirname+filename+"StdMeanPi.pdf"); //################################################################################################################################################ //Standard Gain: Sigma ########################################################################################################################### //################################################################################################################################################ TCanvas *cstdsigmael = new TCanvas("cstdsigmael","dEdx standard gain sigma electrons",1000,700); cstdsigmael->SetFillColor(00); cstdsigmael->SetGrid(); cstdsigmael->cd(); TMultiGraph *mgstdsigmael = new TMultiGraph(); TMultiGraph *mgstdsigmapi = new TMultiGraph(); const int nstdsigma = 5; double gainstdsigma[nstdsigma] = {69., 70., 71., 72., 73.}; double sigmael0025[nstdsigma]; double sigmael0525[nstdsigma]; double sigmael1025[nstdsigma]; double sigmael0030[nstdsigma]; double sigmael0530[nstdsigma]; double sigmael1030[nstdsigma]; double sigmapi0025[nstdsigma]; double sigmapi0525[nstdsigma]; double sigmapi1025[nstdsigma]; double sigmapi0030[nstdsigma]; double sigmapi0530[nstdsigma]; double sigmapi1030[nstdsigma]; double dsigmael0025[nstdsigma]; double dsigmael0525[nstdsigma]; double dsigmael1025[nstdsigma]; double dsigmael0030[nstdsigma]; double dsigmael0530[nstdsigma]; double dsigmael1030[nstdsigma]; double dsigmapi0025[nstdsigma]; double dsigmapi0525[nstdsigma]; double dsigmapi1025[nstdsigma]; double dsigmapi0030[nstdsigma]; double dsigmapi0530[nstdsigma]; double dsigmapi1030[nstdsigma]; sigmael0025[0] = (fitresultsmap[6810025]).sigmael*100; sigmael0025[1] = (fitresultsmap[6820025]).sigmael*100; sigmael0025[2] = (fitresultsmap[6830025]).sigmael*100; sigmael0025[3] = (fitresultsmap[6840025]).sigmael*100; sigmael0025[4] = (fitresultsmap[6850025]).sigmael*100; dsigmael0025[0] = (fitresultsmap[6810025]).sigmaelerr*100; dsigmael0025[1] = (fitresultsmap[6820025]).sigmaelerr*100; dsigmael0025[2] = (fitresultsmap[6830025]).sigmaelerr*100; dsigmael0025[3] = (fitresultsmap[6840025]).sigmaelerr*100; dsigmael0025[4] = (fitresultsmap[6850025]).sigmaelerr*100; sigmael0525[0] = (fitresultsmap[6810525]).sigmael*100; sigmael0525[1] = (fitresultsmap[6820525]).sigmael*100; sigmael0525[2] = (fitresultsmap[6830525]).sigmael*100; sigmael0525[3] = (fitresultsmap[6840525]).sigmael*100; sigmael0525[4] = (fitresultsmap[6850525]).sigmael*100; dsigmael0525[0] = (fitresultsmap[6810525]).sigmaelerr*100; dsigmael0525[1] = (fitresultsmap[6820525]).sigmaelerr*100; dsigmael0525[2] = (fitresultsmap[6830525]).sigmaelerr*100; dsigmael0525[3] = (fitresultsmap[6840525]).sigmaelerr*100; dsigmael0525[4] = (fitresultsmap[6850525]).sigmaelerr*100; sigmael1025[0] = (fitresultsmap[6811025]).sigmael*100; sigmael1025[1] = (fitresultsmap[6821025]).sigmael*100; sigmael1025[2] = (fitresultsmap[6831025]).sigmael*100; sigmael1025[3] = (fitresultsmap[6841025]).sigmael*100; sigmael1025[4] = (fitresultsmap[6851025]).sigmael*100; dsigmael1025[0] = (fitresultsmap[6811025]).sigmaelerr*100; dsigmael1025[1] = (fitresultsmap[6821025]).sigmaelerr*100; dsigmael1025[2] = (fitresultsmap[6831025]).sigmaelerr*100; dsigmael1025[3] = (fitresultsmap[6841025]).sigmaelerr*100; dsigmael1025[4] = (fitresultsmap[6851025]).sigmaelerr*100; sigmael0030[0] = (fitresultsmap[6810030]).sigmael*100; sigmael0030[1] = (fitresultsmap[6820030]).sigmael*100; sigmael0030[2] = (fitresultsmap[6830030]).sigmael*100; sigmael0030[3] = (fitresultsmap[6840030]).sigmael*100; sigmael0030[4] = (fitresultsmap[6850030]).sigmael*100; dsigmael0030[0] = (fitresultsmap[6810030]).sigmaelerr*100; dsigmael0030[1] = (fitresultsmap[6820030]).sigmaelerr*100; dsigmael0030[2] = (fitresultsmap[6830030]).sigmaelerr*100; dsigmael0030[3] = (fitresultsmap[6840030]).sigmaelerr*100; dsigmael0030[4] = (fitresultsmap[6850030]).sigmaelerr*100; sigmael0530[0] = (fitresultsmap[6810530]).sigmael*100; sigmael0530[1] = (fitresultsmap[6820530]).sigmael*100; sigmael0530[2] = (fitresultsmap[6830530]).sigmael*100; sigmael0530[3] = (fitresultsmap[6840530]).sigmael*100; sigmael0530[4] = (fitresultsmap[6850530]).sigmael*100; dsigmael0530[0] = (fitresultsmap[6810530]).sigmaelerr*100; dsigmael0530[1] = (fitresultsmap[6820530]).sigmaelerr*100; dsigmael0530[2] = (fitresultsmap[6830530]).sigmaelerr*100; dsigmael0530[3] = (fitresultsmap[6840530]).sigmaelerr*100; dsigmael0530[4] = (fitresultsmap[6850530]).sigmaelerr*100; sigmael1030[0] = (fitresultsmap[6811030]).sigmael*100; sigmael1030[1] = (fitresultsmap[6821030]).sigmael*100; sigmael1030[2] = (fitresultsmap[6831030]).sigmael*100; sigmael1030[3] = (fitresultsmap[6841030]).sigmael*100; sigmael1030[4] = (fitresultsmap[6851030]).sigmael*100; dsigmael1030[0] = (fitresultsmap[6811030]).sigmaelerr*100; dsigmael1030[1] = (fitresultsmap[6821030]).sigmaelerr*100; dsigmael1030[2] = (fitresultsmap[6831030]).sigmaelerr*100; dsigmael1030[3] = (fitresultsmap[6841030]).sigmaelerr*100; dsigmael1030[4] = (fitresultsmap[6851030]).sigmaelerr*100; sigmapi0025[0] = (fitresultsmap[6810025]).sigmapi*100; sigmapi0025[1] = (fitresultsmap[6820025]).sigmapi*100; sigmapi0025[2] = (fitresultsmap[6830025]).sigmapi*100; sigmapi0025[3] = (fitresultsmap[6840025]).sigmapi*100; sigmapi0025[4] = (fitresultsmap[6850025]).sigmapi*100; dsigmapi0025[0] = (fitresultsmap[6810025]).sigmapierr*100; dsigmapi0025[1] = (fitresultsmap[6820025]).sigmapierr*100; dsigmapi0025[2] = (fitresultsmap[6830025]).sigmapierr*100; dsigmapi0025[3] = (fitresultsmap[6840025]).sigmapierr*100; dsigmapi0025[4] = (fitresultsmap[6850025]).sigmapierr*100; sigmapi0525[0] = (fitresultsmap[6810525]).sigmapi*100; sigmapi0525[1] = (fitresultsmap[6820525]).sigmapi*100; sigmapi0525[2] = (fitresultsmap[6830525]).sigmapi*100; sigmapi0525[3] = (fitresultsmap[6840525]).sigmapi*100; sigmapi0525[4] = (fitresultsmap[6850525]).sigmapi*100; dsigmapi0525[0] = (fitresultsmap[6810525]).sigmapierr*100; dsigmapi0525[1] = (fitresultsmap[6820525]).sigmapierr*100; dsigmapi0525[2] = (fitresultsmap[6830525]).sigmapierr*100; dsigmapi0525[3] = (fitresultsmap[6840525]).sigmapierr*100; dsigmapi0525[4] = (fitresultsmap[6850525]).sigmapierr*100; sigmapi1025[0] = (fitresultsmap[6811025]).sigmapi*100; sigmapi1025[1] = (fitresultsmap[6821025]).sigmapi*100; sigmapi1025[2] = (fitresultsmap[6831025]).sigmapi*100; sigmapi1025[3] = (fitresultsmap[6841025]).sigmapi*100; sigmapi1025[4] = (fitresultsmap[6851025]).sigmapi*100; dsigmapi1025[0] = (fitresultsmap[6811025]).sigmapierr*100; dsigmapi1025[1] = (fitresultsmap[6821025]).sigmapierr*100; dsigmapi1025[2] = (fitresultsmap[6831025]).sigmapierr*100; dsigmapi1025[3] = (fitresultsmap[6841025]).sigmapierr*100; dsigmapi1025[4] = (fitresultsmap[6851025]).sigmapierr*100; sigmapi0030[0] = (fitresultsmap[6810030]).sigmapi*100; sigmapi0030[1] = (fitresultsmap[6820030]).sigmapi*100; sigmapi0030[2] = (fitresultsmap[6830030]).sigmapi*100; sigmapi0030[3] = (fitresultsmap[6840030]).sigmapi*100; sigmapi0030[4] = (fitresultsmap[6850030]).sigmapi*100; dsigmapi0030[0] = (fitresultsmap[6810030]).sigmapierr*100; dsigmapi0030[1] = (fitresultsmap[6820030]).sigmapierr*100; dsigmapi0030[2] = (fitresultsmap[6830030]).sigmapierr*100; dsigmapi0030[3] = (fitresultsmap[6840030]).sigmapierr*100; dsigmapi0030[4] = (fitresultsmap[6850030]).sigmapierr*100; sigmapi0530[0] = (fitresultsmap[6810530]).sigmapi*100; sigmapi0530[1] = (fitresultsmap[6820530]).sigmapi*100; sigmapi0530[2] = (fitresultsmap[6830530]).sigmapi*100; sigmapi0530[3] = (fitresultsmap[6840530]).sigmapi*100; sigmapi0530[4] = (fitresultsmap[6850530]).sigmapi*100; dsigmapi0530[0] = (fitresultsmap[6810530]).sigmapierr*100; dsigmapi0530[1] = (fitresultsmap[6820530]).sigmapierr*100; dsigmapi0530[2] = (fitresultsmap[6830530]).sigmapierr*100; dsigmapi0530[3] = (fitresultsmap[6840530]).sigmapierr*100; dsigmapi0530[4] = (fitresultsmap[6850530]).sigmapierr*100; sigmapi1030[0] = (fitresultsmap[6811030]).sigmapi*100; sigmapi1030[1] = (fitresultsmap[6821030]).sigmapi*100; sigmapi1030[2] = (fitresultsmap[6831030]).sigmapi*100; sigmapi1030[3] = (fitresultsmap[6841030]).sigmapi*100; sigmapi1030[4] = (fitresultsmap[6851030]).sigmapi*100; dsigmapi1030[0] = (fitresultsmap[6811030]).sigmapierr*100; dsigmapi1030[1] = (fitresultsmap[6821030]).sigmapierr*100; dsigmapi1030[2] = (fitresultsmap[6831030]).sigmapierr*100; dsigmapi1030[3] = (fitresultsmap[6841030]).sigmapierr*100; dsigmapi1030[4] = (fitresultsmap[6851030]).sigmapierr*100; //0025 TGraphErrors* grsigmael0025 = new TGraphErrors(nstdsigma,gainstdsigma,sigmael0025,0,dsigmael0025); grsigmael0025->SetLineColor(1); grsigmael0025->SetLineWidth(1); grsigmael0025->SetMarkerColor(1); grsigmael0025->SetMarkerStyle(20); grsigmael0025->SetTitle("electrons 0025"); mgstdsigmael->Add(grsigmael0025); TGraphErrors* grsigmapi0025 = new TGraphErrors(nstdsigma,gainstdsigma,sigmapi0025,0,dsigmapi0025); grsigmapi0025->SetLineColor(1); grsigmapi0025->SetLineWidth(1); grsigmapi0025->SetMarkerColor(1); grsigmapi0025->SetMarkerStyle(20); grsigmapi0025->SetTitle("pions 0025"); mgstdsigmapi->Add(grsigmapi0025); //0525 TGraphErrors* grsigmael0525 = new TGraphErrors(nstdsigma,gainstdsigma,sigmael0525,0,dsigmael0525); grsigmael0525->SetLineColor(2); grsigmael0525->SetLineWidth(1); grsigmael0525->SetMarkerColor(2); grsigmael0525->SetMarkerStyle(20); grsigmael0525->SetTitle("electrons 0525"); mgstdsigmael->Add(grsigmael0525); TGraphErrors* grsigmapi0525 = new TGraphErrors(nstdsigma,gainstdsigma,sigmapi0525,0,dsigmapi0525); grsigmapi0525->SetLineColor(2); grsigmapi0525->SetLineWidth(1); grsigmapi0525->SetMarkerColor(2); grsigmapi0525->SetMarkerStyle(20); grsigmapi0525->SetTitle("pions 0525"); mgstdsigmapi->Add(grsigmapi0525); //1025 TGraphErrors* grsigmael1025 = new TGraphErrors(nstdsigma,gainstdsigma,sigmael1025,0,dsigmael1025); grsigmael1025->SetLineColor(3); grsigmael1025->SetLineWidth(1); grsigmael1025->SetMarkerColor(3); grsigmael1025->SetMarkerStyle(20); grsigmael1025->SetTitle("electrons 1025"); mgstdsigmael->Add(grsigmael1025); TGraphErrors* grsigmapi1025 = new TGraphErrors(nstdsigma,gainstdsigma,sigmapi1025,0,dsigmapi1025); grsigmapi1025->SetLineColor(3); grsigmapi1025->SetLineWidth(1); grsigmapi1025->SetMarkerColor(3); grsigmapi1025->SetMarkerStyle(20); grsigmapi1025->SetTitle("pions 1025"); mgstdsigmapi->Add(grsigmapi1025); //0030 TGraphErrors* grsigmael0030 = new TGraphErrors(nstdsigma,gainstdsigma,sigmael0030,0,dsigmael0030); grsigmael0030->SetLineColor(4); grsigmael0030->SetLineWidth(1); grsigmael0030->SetMarkerColor(4); grsigmael0030->SetMarkerStyle(20); grsigmael0030->SetTitle("electrons 0030"); mgstdsigmael->Add(grsigmael0030); TGraphErrors* grsigmapi0030 = new TGraphErrors(nstdsigma,gainstdsigma,sigmapi0030,0,dsigmapi0030); grsigmapi0030->SetLineColor(4); grsigmapi0030->SetLineWidth(1); grsigmapi0030->SetMarkerColor(4); grsigmapi0030->SetMarkerStyle(20); grsigmapi0030->SetTitle("pions 0030"); mgstdsigmapi->Add(grsigmapi0030); //0530 TGraphErrors* grsigmael0530 = new TGraphErrors(nstdsigma,gainstdsigma,sigmael0530,0,dsigmael0530); grsigmael0530->SetLineColor(7); grsigmael0530->SetLineWidth(1); grsigmael0530->SetMarkerColor(7); grsigmael0530->SetMarkerStyle(20); grsigmael0530->SetTitle("electrons 0530"); mgstdsigmael->Add(grsigmael0530); TGraphErrors* grsigmapi0530 = new TGraphErrors(nstdsigma,gainstdsigma,sigmapi0530,0,dsigmapi0530); grsigmapi0530->SetLineColor(7); grsigmapi0530->SetLineWidth(1); grsigmapi0530->SetMarkerColor(7); grsigmapi0530->SetMarkerStyle(20); grsigmapi0530->SetTitle("pions 0530"); mgstdsigmapi->Add(grsigmapi0530); //1030 TGraphErrors* grsigmael1030 = new TGraphErrors(nstdsigma,gainstdsigma,sigmael1030,0,dsigmael1030); grsigmael1030->SetLineColor(6); grsigmael1030->SetLineWidth(1); grsigmael1030->SetMarkerColor(6); grsigmael1030->SetMarkerStyle(20); grsigmael1030->SetTitle("electrons 1030"); mgstdsigmael->Add(grsigmael1030); TGraphErrors* grsigmapi1030 = new TGraphErrors(nstdsigma,gainstdsigma,sigmapi1030,0,dsigmapi1030); grsigmapi1030->SetLineColor(6); grsigmapi1030->SetLineWidth(1); grsigmapi1030->SetMarkerColor(6); grsigmapi1030->SetMarkerStyle(20); grsigmapi1030->SetTitle("pions 1030"); mgstdsigmapi->Add(grsigmapi1030); mgstdsigmael->Draw("ap"); mgstdsigmael->GetXaxis()->SetTitle("HV scaling factor (%)"); mgstdsigmael->GetYaxis()->SetTitle("#sigma(dE/dx) electrons (%)"); mgstdsigmael->GetXaxis()->SetLimits(68,74); mgstdsigmael->GetYaxis()->SetRangeUser(9,12); TLegend *legstdsigmael=new TLegend(0.6,0.6,0.95,.9); legstdsigmael->SetBorderSize(1); legstdsigmael->SetFillColor(10); legstdsigmael->AddEntry(grsigmael0025,"truncations [ 0,75]","p"); legstdsigmael->AddEntry(grsigmael0525,"truncations [ 5,75]","p"); legstdsigmael->AddEntry(grsigmael1025,"truncations [10,75]","p"); legstdsigmael->AddEntry(grsigmael0030,"truncations [ 0,70]","p"); legstdsigmael->AddEntry(grsigmael0530,"truncations [ 5,70]","p"); legstdsigmael->AddEntry(grsigmael1030,"truncations [10,70]","p"); legstdsigmael->Draw("same"); cstdsigmael->Update(); cstdsigmael->Modified(); cstdsigmael->Write(); cstdsigmael->SaveAs(dirname+filename+"StdSigmaEl.pdf"); TCanvas *cstdsigmapi = new TCanvas("cstdsigmapi","dEdx standard gain sigma pions",1000,700); cstdsigmapi->SetFillColor(00); cstdsigmapi->SetGrid(); cstdsigmapi->cd(); mgstdsigmapi->Draw("ap"); mgstdsigmapi->GetXaxis()->SetTitle("HV scaling factor (%)"); mgstdsigmapi->GetYaxis()->SetTitle("#sigma(dE/dx) pions (%)"); mgstdsigmapi->GetXaxis()->SetLimits(68,74); mgstdsigmapi->GetYaxis()->SetRangeUser(10,14); TLegend *legstdsigmapi=new TLegend(0.6,0.6,0.95,.9); legstdsigmapi->SetBorderSize(1); legstdsigmapi->SetFillColor(10); legstdsigmapi->AddEntry(grsigmapi0025,"truncations [0,75]","p"); legstdsigmapi->AddEntry(grsigmapi0525,"truncations [5,75]","p"); legstdsigmapi->AddEntry(grsigmapi1025,"truncations [10,75]","p"); legstdsigmapi->AddEntry(grsigmapi0030,"truncations [0,70]","p"); legstdsigmapi->AddEntry(grsigmapi0530,"truncations [5,70]","p"); legstdsigmapi->AddEntry(grsigmapi1030,"truncations [10,70]","p"); legstdsigmapi->Draw("same"); cstdsigmapi->Update(); cstdsigmapi->Modified(); cstdsigmapi->Write(); cstdsigmapi->SaveAs(dirname+filename+"StdSigmaPi.pdf"); //################################################################################################################################################ //IBF Gain: Mean TF 2 400############################################################################################################################ //################################################################################################################################################ TCanvas *cibfmeanel = new TCanvas("cibfmeanel","dEdx ibf gain mean electrons",1000,700); cibfmeanel->SetFillColor(00); cibfmeanel->SetGrid(); cibfmeanel->cd(); TMultiGraph *mgibfmeanel = new TMultiGraph(); TMultiGraph *mgibfmeanpi = new TMultiGraph(); const int nibfmean = 4; double gainibfmean[nibfmean] = {100., 103., 105., 107.}; double meanibfel0025[nibfmean]; double meanibfel0525[nibfmean]; double meanibfel1025[nibfmean]; double meanibfel0030[nibfmean]; double meanibfel0530[nibfmean]; double meanibfel1030[nibfmean]; double meanibfpi0025[nibfmean]; double meanibfpi0525[nibfmean]; double meanibfpi1025[nibfmean]; double meanibfpi0030[nibfmean]; double meanibfpi0530[nibfmean]; double meanibfpi1030[nibfmean]; double dmeanibfel0025[nibfmean]; double dmeanibfel0525[nibfmean]; double dmeanibfel1025[nibfmean]; double dmeanibfel0030[nibfmean]; double dmeanibfel0530[nibfmean]; double dmeanibfel1030[nibfmean]; double dmeanibfpi0025[nibfmean]; double dmeanibfpi0525[nibfmean]; double dmeanibfpi1025[nibfmean]; double dmeanibfpi0030[nibfmean]; double dmeanibfpi0530[nibfmean]; double dmeanibfpi1030[nibfmean]; meanibfel0025[0] = (fitresultsmap[6870025]).meanel; meanibfel0025[1] = (fitresultsmap[6920025]).meanel; meanibfel0025[2] = (fitresultsmap[6990025]).meanel; meanibfel0025[3] = (fitresultsmap[7040025]).meanel; dmeanibfel0025[0] = (fitresultsmap[6870025]).meanelerr; dmeanibfel0025[1] = (fitresultsmap[6920025]).meanelerr; dmeanibfel0025[2] = (fitresultsmap[6990025]).meanelerr; dmeanibfel0025[3] = (fitresultsmap[7040025]).meanelerr; meanibfel0525[0] = (fitresultsmap[6870525]).meanel; meanibfel0525[1] = (fitresultsmap[6920525]).meanel; meanibfel0525[2] = (fitresultsmap[6990525]).meanel; meanibfel0525[3] = (fitresultsmap[7040525]).meanel; dmeanibfel0525[0] = (fitresultsmap[6870525]).meanelerr; dmeanibfel0525[1] = (fitresultsmap[6920525]).meanelerr; dmeanibfel0525[2] = (fitresultsmap[6990525]).meanelerr; dmeanibfel0525[3] = (fitresultsmap[7040525]).meanelerr; meanibfel1025[0] = (fitresultsmap[6871025]).meanel; meanibfel1025[1] = (fitresultsmap[6921025]).meanel; meanibfel1025[2] = (fitresultsmap[6991025]).meanel; meanibfel1025[3] = (fitresultsmap[7041025]).meanel; dmeanibfel1025[0] = (fitresultsmap[6871025]).meanelerr; dmeanibfel1025[1] = (fitresultsmap[6921025]).meanelerr; dmeanibfel1025[2] = (fitresultsmap[6991025]).meanelerr; dmeanibfel1025[3] = (fitresultsmap[7041025]).meanelerr; meanibfel0030[0] = (fitresultsmap[6870030]).meanel; meanibfel0030[1] = (fitresultsmap[6920030]).meanel; meanibfel0030[2] = (fitresultsmap[6990030]).meanel; meanibfel0030[3] = (fitresultsmap[7040030]).meanel; dmeanibfel0030[0] = (fitresultsmap[6870030]).meanelerr; dmeanibfel0030[1] = (fitresultsmap[6920030]).meanelerr; dmeanibfel0030[2] = (fitresultsmap[6990030]).meanelerr; dmeanibfel0030[3] = (fitresultsmap[7040030]).meanelerr; meanibfel0530[0] = (fitresultsmap[6870530]).meanel; meanibfel0530[1] = (fitresultsmap[6920530]).meanel; meanibfel0530[2] = (fitresultsmap[6990530]).meanel; meanibfel0530[3] = (fitresultsmap[7040530]).meanel; dmeanibfel0530[0] = (fitresultsmap[6870530]).meanelerr; dmeanibfel0530[1] = (fitresultsmap[6920530]).meanelerr; dmeanibfel0530[2] = (fitresultsmap[6990530]).meanelerr; dmeanibfel0530[3] = (fitresultsmap[7040530]).meanelerr; meanibfel1030[0] = (fitresultsmap[6871030]).meanel; meanibfel1030[1] = (fitresultsmap[6921030]).meanel; meanibfel1030[2] = (fitresultsmap[6991030]).meanel; meanibfel1030[3] = (fitresultsmap[7041030]).meanel; dmeanibfel1030[0] = (fitresultsmap[6871030]).meanelerr; dmeanibfel1030[1] = (fitresultsmap[6921030]).meanelerr; dmeanibfel1030[2] = (fitresultsmap[6991030]).meanelerr; dmeanibfel1030[3] = (fitresultsmap[7041030]).meanelerr; meanibfpi0025[0] = (fitresultsmap[6870025]).meanpi; meanibfpi0025[1] = (fitresultsmap[6920025]).meanpi; meanibfpi0025[2] = (fitresultsmap[6990025]).meanpi; meanibfpi0025[3] = (fitresultsmap[7040025]).meanpi; dmeanibfpi0025[0] = (fitresultsmap[6870025]).meanpierr; dmeanibfpi0025[1] = (fitresultsmap[6920025]).meanpierr; dmeanibfpi0025[2] = (fitresultsmap[6990025]).meanpierr; dmeanibfpi0025[3] = (fitresultsmap[7040025]).meanpierr; meanibfpi0525[0] = (fitresultsmap[6870525]).meanpi; meanibfpi0525[1] = (fitresultsmap[6920525]).meanpi; meanibfpi0525[2] = (fitresultsmap[6990525]).meanpi; meanibfpi0525[3] = (fitresultsmap[7040525]).meanpi; dmeanibfpi0525[0] = (fitresultsmap[6870525]).meanpierr; dmeanibfpi0525[1] = (fitresultsmap[6920525]).meanpierr; dmeanibfpi0525[2] = (fitresultsmap[6990525]).meanpierr; dmeanibfpi0525[3] = (fitresultsmap[7040525]).meanpierr; meanibfpi1025[0] = (fitresultsmap[6871025]).meanpi; meanibfpi1025[1] = (fitresultsmap[6921025]).meanpi; meanibfpi1025[2] = (fitresultsmap[6991025]).meanpi; meanibfpi1025[3] = (fitresultsmap[7041025]).meanpi; dmeanibfpi1025[0] = (fitresultsmap[6871025]).meanpierr; dmeanibfpi1025[1] = (fitresultsmap[6921025]).meanpierr; dmeanibfpi1025[2] = (fitresultsmap[6991025]).meanpierr; dmeanibfpi1025[3] = (fitresultsmap[7041025]).meanpierr; meanibfpi0030[0] = (fitresultsmap[6870030]).meanpi; meanibfpi0030[1] = (fitresultsmap[6920030]).meanpi; meanibfpi0030[2] = (fitresultsmap[6990030]).meanpi; meanibfpi0030[3] = (fitresultsmap[7040030]).meanpi; dmeanibfpi0030[0] = (fitresultsmap[6870030]).meanpierr; dmeanibfpi0030[1] = (fitresultsmap[6920030]).meanpierr; dmeanibfpi0030[2] = (fitresultsmap[6990030]).meanpierr; dmeanibfpi0030[3] = (fitresultsmap[7040030]).meanpierr; meanibfpi0530[0] = (fitresultsmap[6870530]).meanpi; meanibfpi0530[1] = (fitresultsmap[6920530]).meanpi; meanibfpi0530[2] = (fitresultsmap[6990530]).meanpi; meanibfpi0530[3] = (fitresultsmap[7040530]).meanpi; dmeanibfpi0530[0] = (fitresultsmap[6870530]).meanpierr; dmeanibfpi0530[1] = (fitresultsmap[6920530]).meanpierr; dmeanibfpi0530[2] = (fitresultsmap[6990530]).meanpierr; dmeanibfpi0530[3] = (fitresultsmap[7040530]).meanpierr; meanibfpi1030[0] = (fitresultsmap[6871030]).meanpi; meanibfpi1030[1] = (fitresultsmap[6921030]).meanpi; meanibfpi1030[2] = (fitresultsmap[6991030]).meanpi; meanibfpi1030[3] = (fitresultsmap[7041030]).meanpi; dmeanibfpi1030[0] = (fitresultsmap[6871030]).meanpierr; dmeanibfpi1030[1] = (fitresultsmap[6921030]).meanpierr; dmeanibfpi1030[2] = (fitresultsmap[6991030]).meanpierr; dmeanibfpi1030[3] = (fitresultsmap[7041030]).meanpierr; //0025 TGraphErrors* grmeanibfel0025 = new TGraphErrors(nibfmean,gainibfmean,meanibfel0025,0,dmeanibfel0025); grmeanibfel0025->SetLineColor(1); grmeanibfel0025->SetLineWidth(1); grmeanibfel0025->SetMarkerColor(1); grmeanibfel0025->SetMarkerStyle(20); grmeanibfel0025->SetTitle("electrons 0025"); mgibfmeanel->Add(grmeanibfel0025); TGraphErrors* grmeanibfpi0025 = new TGraphErrors(nibfmean,gainibfmean,meanibfpi0025,0,dmeanibfpi0025); grmeanibfpi0025->SetLineColor(1); grmeanibfpi0025->SetLineWidth(1); grmeanibfpi0025->SetMarkerColor(1); grmeanibfpi0025->SetMarkerStyle(20); grmeanibfpi0025->SetTitle("pions 0025"); mgibfmeanpi->Add(grmeanibfpi0025); //0525 TGraphErrors* grmeanibfel0525 = new TGraphErrors(nibfmean,gainibfmean,meanibfel0525,0,dmeanibfel0525); grmeanibfel0525->SetLineColor(2); grmeanibfel0525->SetLineWidth(1); grmeanibfel0525->SetMarkerColor(2); grmeanibfel0525->SetMarkerStyle(20); grmeanibfel0525->SetTitle("electrons 0525"); mgibfmeanel->Add(grmeanibfel0525); TGraphErrors* grmeanibfpi0525 = new TGraphErrors(nibfmean,gainibfmean,meanibfpi0525,0,dmeanibfpi0525); grmeanibfpi0525->SetLineColor(2); grmeanibfpi0525->SetLineWidth(1); grmeanibfpi0525->SetMarkerColor(2); grmeanibfpi0525->SetMarkerStyle(20); grmeanibfpi0525->SetTitle("pions 0525"); mgibfmeanpi->Add(grmeanibfpi0525); //1025 TGraphErrors* grmeanibfel1025 = new TGraphErrors(nibfmean,gainibfmean,meanibfel1025,0,dmeanibfel1025); grmeanibfel1025->SetLineColor(3); grmeanibfel1025->SetLineWidth(1); grmeanibfel1025->SetMarkerColor(3); grmeanibfel1025->SetMarkerStyle(20); grmeanibfel1025->SetTitle("electrons 1025"); mgibfmeanel->Add(grmeanibfel1025); TGraphErrors* grmeanibfpi1025 = new TGraphErrors(nibfmean,gainibfmean,meanibfpi1025,0,dmeanibfpi1025); grmeanibfpi1025->SetLineColor(3); grmeanibfpi1025->SetLineWidth(1); grmeanibfpi1025->SetMarkerColor(3); grmeanibfpi1025->SetMarkerStyle(20); grmeanibfpi1025->SetTitle("pions 1025"); mgibfmeanpi->Add(grmeanibfpi1025); //0030 TGraphErrors* grmeanibfel0030 = new TGraphErrors(nibfmean,gainibfmean,meanibfel0030,0,dmeanibfel0030); grmeanibfel0030->SetLineColor(4); grmeanibfel0030->SetLineWidth(1); grmeanibfel0030->SetMarkerColor(4); grmeanibfel0030->SetMarkerStyle(20); grmeanibfel0030->SetTitle("electrons 0030"); mgibfmeanel->Add(grmeanibfel0030); TGraphErrors* grmeanibfpi0030 = new TGraphErrors(nibfmean,gainibfmean,meanibfpi0030,0,dmeanibfpi0030); grmeanibfpi0030->SetLineColor(4); grmeanibfpi0030->SetLineWidth(1); grmeanibfpi0030->SetMarkerColor(4); grmeanibfpi0030->SetMarkerStyle(20); grmeanibfpi0030->SetTitle("pions 0030"); mgibfmeanpi->Add(grmeanibfpi0030); //0530 TGraphErrors* grmeanibfel0530 = new TGraphErrors(nibfmean,gainibfmean,meanibfel0530,0,dmeanibfel0530); grmeanibfel0530->SetLineColor(7); grmeanibfel0530->SetLineWidth(1); grmeanibfel0530->SetMarkerColor(7); grmeanibfel0530->SetMarkerStyle(20); grmeanibfel0530->SetTitle("electrons 0530"); mgibfmeanel->Add(grmeanibfel0530); TGraphErrors* grmeanibfpi0530 = new TGraphErrors(nibfmean,gainibfmean,meanibfpi0530,0,dmeanibfpi0530); grmeanibfpi0530->SetLineColor(7); grmeanibfpi0530->SetLineWidth(1); grmeanibfpi0530->SetMarkerColor(7); grmeanibfpi0530->SetMarkerStyle(20); grmeanibfpi0530->SetTitle("pions 0530"); mgibfmeanpi->Add(grmeanibfpi0530); //1030 TGraphErrors* grmeanibfel1030 = new TGraphErrors(nibfmean,gainibfmean,meanibfel1030,0,dmeanibfel1030); grmeanibfel1030->SetLineColor(6); grmeanibfel1030->SetLineWidth(1); grmeanibfel1030->SetMarkerColor(6); grmeanibfel1030->SetMarkerStyle(20); grmeanibfel1030->SetTitle("electrons 1030"); mgibfmeanel->Add(grmeanibfel1030); TGraphErrors* grmeanibfpi1030 = new TGraphErrors(nibfmean,gainibfmean,meanibfpi1030,0,dmeanibfpi1030); grmeanibfpi1030->SetLineColor(6); grmeanibfpi1030->SetLineWidth(1); grmeanibfpi1030->SetMarkerColor(6); grmeanibfpi1030->SetMarkerStyle(20); grmeanibfpi1030->SetTitle("pions 1030"); mgibfmeanpi->Add(grmeanibfpi1030); mgibfmeanel->Draw("ap"); mgibfmeanel->GetXaxis()->SetTitle("HV scaling factor (%)"); mgibfmeanel->GetYaxis()->SetTitle(" electrons (a.u.)"); mgibfmeanel->GetXaxis()->SetLimits(99,108); mgibfmeanel->GetYaxis()->SetRangeUser(0.001,330); TLegend *legibfmeanel=new TLegend(0.2,0.6,0.5,.9); legibfmeanel->SetBorderSize(1); legibfmeanel->SetFillColor(10); legibfmeanel->AddEntry(grmeanibfel0025,"truncations [0,75]","p"); legibfmeanel->AddEntry(grmeanibfel0525,"truncations [5,75]","p"); legibfmeanel->AddEntry(grmeanibfel1025,"truncations [10,75]","p"); legibfmeanel->AddEntry(grmeanibfel0030,"truncations [0,70]","p"); legibfmeanel->AddEntry(grmeanibfel0530,"truncations [5,70]","p"); legibfmeanel->AddEntry(grmeanibfel1030,"truncations [10,70]","p"); legibfmeanel->Draw("same"); cibfmeanel->Update(); cibfmeanel->Modified(); cibfmeanel->Write(); cibfmeanel->SaveAs(dirname+filename+"ibfMeanEl.pdf"); TCanvas *cibfmeanpi = new TCanvas("cibfmeanpi","dEdx ibf gain mean pions",1000,700); cibfmeanpi->SetFillColor(00); cibfmeanpi->SetGrid(); cibfmeanpi->cd(); mgibfmeanpi->Draw("ap"); mgibfmeanpi->GetXaxis()->SetTitle("HV scaling factor (%)"); mgibfmeanpi->GetYaxis()->SetTitle(" pions (a.u.)"); mgibfmeanpi->GetXaxis()->SetLimits(99,108); mgibfmeanpi->GetYaxis()->SetRangeUser(0.001,220); TLegend *legibfmeanpi=new TLegend(0.2,0.6,0.5,.9); legibfmeanpi->SetBorderSize(1); legibfmeanpi->SetFillColor(10); legibfmeanpi->AddEntry(grmeanibfpi0025,"truncations [0,75]","p"); legibfmeanpi->AddEntry(grmeanibfpi0525,"truncations [5,75]","p"); legibfmeanpi->AddEntry(grmeanibfpi1025,"truncations [10,75]","p"); legibfmeanpi->AddEntry(grmeanibfpi0030,"truncations [0,70]","p"); legibfmeanpi->AddEntry(grmeanibfpi0530,"truncations [5,70]","p"); legibfmeanpi->AddEntry(grmeanibfpi1030,"truncations [10,70]","p"); legibfmeanpi->Draw("same"); cibfmeanpi->Update(); cibfmeanpi->Modified(); cibfmeanpi->Write(); cibfmeanpi->SaveAs(dirname+filename+"ibfMeanPi.pdf"); //################################################################################################################################################ //ibf Gain: Sigma ########################################################################################################################### //################################################################################################################################################ TCanvas *cibfsigmael = new TCanvas("cibfsigmael","dEdx ibf gain sigma electrons",1000,700); cibfsigmael->SetFillColor(00); cibfsigmael->SetGrid(); cibfsigmael->cd(); TMultiGraph *mgibfsigmael = new TMultiGraph(); TMultiGraph *mgibfsigmapi = new TMultiGraph(); const int nibfsigma = 4; double gainibfsigma[nibfsigma] = {100., 103., 105., 107.}; double sigmaelibf0025[nibfsigma]; double sigmaelibf0525[nibfsigma]; double sigmaelibf1025[nibfsigma]; double sigmaelibf0030[nibfsigma]; double sigmaelibf0530[nibfsigma]; double sigmaelibf1030[nibfsigma]; double sigmapiibf0025[nibfsigma]; double sigmapiibf0525[nibfsigma]; double sigmapiibf1025[nibfsigma]; double sigmapiibf0030[nibfsigma]; double sigmapiibf0530[nibfsigma]; double sigmapiibf1030[nibfsigma]; double dsigmaelibf0025[nibfsigma]; double dsigmaelibf0525[nibfsigma]; double dsigmaelibf1025[nibfsigma]; double dsigmaelibf0030[nibfsigma]; double dsigmaelibf0530[nibfsigma]; double dsigmaelibf1030[nibfsigma]; double dsigmapiibf0025[nibfsigma]; double dsigmapiibf0525[nibfsigma]; double dsigmapiibf1025[nibfsigma]; double dsigmapiibf0030[nibfsigma]; double dsigmapiibf0530[nibfsigma]; double dsigmapiibf1030[nibfsigma]; sigmaelibf0025[0] = (fitresultsmap[6870025]).sigmael*100; sigmaelibf0025[1] = (fitresultsmap[6920025]).sigmael*100; sigmaelibf0025[2] = (fitresultsmap[6990025]).sigmael*100; sigmaelibf0025[3] = (fitresultsmap[7040025]).sigmael*100; dsigmaelibf0025[0] = (fitresultsmap[6870025]).sigmaelerr*100; dsigmaelibf0025[1] = (fitresultsmap[6920025]).sigmaelerr*100; dsigmaelibf0025[2] = (fitresultsmap[6990025]).sigmaelerr*100; dsigmaelibf0025[3] = (fitresultsmap[7040025]).sigmaelerr*100; sigmaelibf0525[0] = (fitresultsmap[6870525]).sigmael*100; sigmaelibf0525[1] = (fitresultsmap[6920525]).sigmael*100; sigmaelibf0525[2] = (fitresultsmap[6990525]).sigmael*100; sigmaelibf0525[3] = (fitresultsmap[7040525]).sigmael*100; dsigmaelibf0525[0] = (fitresultsmap[6870525]).sigmaelerr*100; dsigmaelibf0525[1] = (fitresultsmap[6920525]).sigmaelerr*100; dsigmaelibf0525[2] = (fitresultsmap[6990525]).sigmaelerr*100; dsigmaelibf0525[3] = (fitresultsmap[7040525]).sigmaelerr*100; sigmaelibf1025[0] = (fitresultsmap[6871025]).sigmael*100; sigmaelibf1025[1] = (fitresultsmap[6921025]).sigmael*100; sigmaelibf1025[2] = (fitresultsmap[6991025]).sigmael*100; sigmaelibf1025[3] = (fitresultsmap[7041025]).sigmael*100; dsigmaelibf1025[0] = (fitresultsmap[6871025]).sigmaelerr*100; dsigmaelibf1025[1] = (fitresultsmap[6921025]).sigmaelerr*100; dsigmaelibf1025[2] = (fitresultsmap[6991025]).sigmaelerr*100; dsigmaelibf1025[3] = (fitresultsmap[7041025]).sigmaelerr*100; sigmaelibf0030[0] = (fitresultsmap[6870030]).sigmael*100; sigmaelibf0030[1] = (fitresultsmap[6920030]).sigmael*100; sigmaelibf0030[2] = (fitresultsmap[6990030]).sigmael*100; sigmaelibf0030[3] = (fitresultsmap[7040030]).sigmael*100; dsigmaelibf0030[0] = (fitresultsmap[6870030]).sigmaelerr*100; dsigmaelibf0030[1] = (fitresultsmap[6920030]).sigmaelerr*100; dsigmaelibf0030[2] = (fitresultsmap[6990030]).sigmaelerr*100; dsigmaelibf0030[3] = (fitresultsmap[7040030]).sigmaelerr*100; sigmaelibf0530[0] = (fitresultsmap[6870530]).sigmael*100; sigmaelibf0530[1] = (fitresultsmap[6920530]).sigmael*100; sigmaelibf0530[2] = (fitresultsmap[6990530]).sigmael*100; sigmaelibf0530[3] = (fitresultsmap[7040530]).sigmael*100; dsigmaelibf0530[0] = (fitresultsmap[6870530]).sigmaelerr*100; dsigmaelibf0530[1] = (fitresultsmap[6920530]).sigmaelerr*100; dsigmaelibf0530[2] = (fitresultsmap[6990530]).sigmaelerr*100; dsigmaelibf0530[3] = (fitresultsmap[7040530]).sigmaelerr*100; sigmaelibf1030[0] = (fitresultsmap[6871030]).sigmael*100; sigmaelibf1030[1] = (fitresultsmap[6921030]).sigmael*100; sigmaelibf1030[2] = (fitresultsmap[6991030]).sigmael*100; sigmaelibf1030[3] = (fitresultsmap[7041030]).sigmael*100; dsigmaelibf1030[0] = (fitresultsmap[6871030]).sigmaelerr*100; dsigmaelibf1030[1] = (fitresultsmap[6921030]).sigmaelerr*100; dsigmaelibf1030[2] = (fitresultsmap[6991030]).sigmaelerr*100; dsigmaelibf1030[3] = (fitresultsmap[7041030]).sigmaelerr*100; sigmapiibf0025[0] = (fitresultsmap[6870025]).sigmapi*100; sigmapiibf0025[1] = (fitresultsmap[6920025]).sigmapi*100; sigmapiibf0025[2] = (fitresultsmap[6990025]).sigmapi*100; sigmapiibf0025[3] = (fitresultsmap[7040025]).sigmapi*100; dsigmapiibf0025[0] = (fitresultsmap[6870025]).sigmapierr*100; dsigmapiibf0025[1] = (fitresultsmap[6920025]).sigmapierr*100; dsigmapiibf0025[2] = (fitresultsmap[6990025]).sigmapierr*100; dsigmapiibf0025[3] = (fitresultsmap[7040025]).sigmapierr*100; sigmapiibf0525[0] = (fitresultsmap[6870525]).sigmapi*100; sigmapiibf0525[1] = (fitresultsmap[6920525]).sigmapi*100; sigmapiibf0525[2] = (fitresultsmap[6990525]).sigmapi*100; sigmapiibf0525[3] = (fitresultsmap[7040525]).sigmapi*100; dsigmapiibf0525[0] = (fitresultsmap[6870525]).sigmapierr*100; dsigmapiibf0525[1] = (fitresultsmap[6920525]).sigmapierr*100; dsigmapiibf0525[2] = (fitresultsmap[6990525]).sigmapierr*100; dsigmapiibf0525[3] = (fitresultsmap[7040525]).sigmapierr*100; sigmapiibf1025[0] = (fitresultsmap[6871025]).sigmapi*100; sigmapiibf1025[1] = (fitresultsmap[6921025]).sigmapi*100; sigmapiibf1025[2] = (fitresultsmap[6991025]).sigmapi*100; sigmapiibf1025[3] = (fitresultsmap[7041025]).sigmapi*100; dsigmapiibf1025[0] = (fitresultsmap[6871025]).sigmapierr*100; dsigmapiibf1025[1] = (fitresultsmap[6921025]).sigmapierr*100; dsigmapiibf1025[2] = (fitresultsmap[6991025]).sigmapierr*100; dsigmapiibf1025[3] = (fitresultsmap[7041025]).sigmapierr*100; sigmapiibf0030[0] = (fitresultsmap[6870030]).sigmapi*100; sigmapiibf0030[1] = (fitresultsmap[6920030]).sigmapi*100; sigmapiibf0030[2] = (fitresultsmap[6990030]).sigmapi*100; sigmapiibf0030[3] = (fitresultsmap[7040030]).sigmapi*100; dsigmapiibf0030[0] = (fitresultsmap[6870030]).sigmapierr*100; dsigmapiibf0030[1] = (fitresultsmap[6920030]).sigmapierr*100; dsigmapiibf0030[2] = (fitresultsmap[6990030]).sigmapierr*100; dsigmapiibf0030[3] = (fitresultsmap[7040030]).sigmapierr*100; sigmapiibf0530[0] = (fitresultsmap[6870530]).sigmapi*100; sigmapiibf0530[1] = (fitresultsmap[6920530]).sigmapi*100; sigmapiibf0530[2] = (fitresultsmap[6990530]).sigmapi*100; sigmapiibf0530[3] = (fitresultsmap[7040530]).sigmapi*100; dsigmapiibf0530[0] = (fitresultsmap[6870530]).sigmapierr*100; dsigmapiibf0530[1] = (fitresultsmap[6920530]).sigmapierr*100; dsigmapiibf0530[2] = (fitresultsmap[6990530]).sigmapierr*100; dsigmapiibf0530[3] = (fitresultsmap[7040530]).sigmapierr*100; sigmapiibf1030[0] = (fitresultsmap[6871030]).sigmapi*100; sigmapiibf1030[1] = (fitresultsmap[6921030]).sigmapi*100; sigmapiibf1030[2] = (fitresultsmap[6991030]).sigmapi*100; sigmapiibf1030[3] = (fitresultsmap[7041030]).sigmapi*100; dsigmapiibf1030[0] = (fitresultsmap[6871030]).sigmapierr*100; dsigmapiibf1030[1] = (fitresultsmap[6921030]).sigmapierr*100; dsigmapiibf1030[2] = (fitresultsmap[6991030]).sigmapierr*100; dsigmapiibf1030[3] = (fitresultsmap[7041030]).sigmapierr*100; //0025 TGraphErrors* grsigmaelibf0025 = new TGraphErrors(nibfsigma,gainibfsigma,sigmaelibf0025,0,dsigmaelibf0025); grsigmaelibf0025->SetLineColor(1); grsigmaelibf0025->SetLineWidth(1); grsigmaelibf0025->SetMarkerColor(1); grsigmaelibf0025->SetMarkerStyle(20); grsigmaelibf0025->SetTitle("electrons 0025"); mgibfsigmael->Add(grsigmaelibf0025); TGraphErrors* grsigmapiibf0025 = new TGraphErrors(nibfsigma,gainibfsigma,sigmapiibf0025,0,dsigmapiibf0025); grsigmapiibf0025->SetLineColor(1); grsigmapiibf0025->SetLineWidth(1); grsigmapiibf0025->SetMarkerColor(1); grsigmapiibf0025->SetMarkerStyle(20); grsigmapiibf0025->SetTitle("pions 0025"); mgibfsigmapi->Add(grsigmapiibf0025); //0525 TGraphErrors* grsigmaelibf0525 = new TGraphErrors(nibfsigma,gainibfsigma,sigmaelibf0525,0,dsigmaelibf0525); grsigmaelibf0525->SetLineColor(2); grsigmaelibf0525->SetLineWidth(1); grsigmaelibf0525->SetMarkerColor(2); grsigmaelibf0525->SetMarkerStyle(20); grsigmaelibf0525->SetTitle("electrons 0525"); mgibfsigmael->Add(grsigmaelibf0525); TGraphErrors* grsigmapiibf0525 = new TGraphErrors(nibfsigma,gainibfsigma,sigmapiibf0525,0,dsigmapiibf0525); grsigmapiibf0525->SetLineColor(2); grsigmapiibf0525->SetLineWidth(1); grsigmapiibf0525->SetMarkerColor(2); grsigmapiibf0525->SetMarkerStyle(20); grsigmapiibf0525->SetTitle("pions 0525"); mgibfsigmapi->Add(grsigmapiibf0525); //1025 TGraphErrors* grsigmaelibf1025 = new TGraphErrors(nibfsigma,gainibfsigma,sigmaelibf1025,0,dsigmaelibf1025); grsigmaelibf1025->SetLineColor(3); grsigmaelibf1025->SetLineWidth(1); grsigmaelibf1025->SetMarkerColor(3); grsigmaelibf1025->SetMarkerStyle(20); grsigmaelibf1025->SetTitle("electrons 1025"); mgibfsigmael->Add(grsigmaelibf1025); TGraphErrors* grsigmapiibf1025 = new TGraphErrors(nibfsigma,gainibfsigma,sigmapiibf1025,0,dsigmapiibf1025); grsigmapiibf1025->SetLineColor(3); grsigmapiibf1025->SetLineWidth(1); grsigmapiibf1025->SetMarkerColor(3); grsigmapiibf1025->SetMarkerStyle(20); grsigmapiibf1025->SetTitle("pions 1025"); mgibfsigmapi->Add(grsigmapiibf1025); //0030 TGraphErrors* grsigmaelibf0030 = new TGraphErrors(nibfsigma,gainibfsigma,sigmaelibf0030,0,dsigmaelibf0030); grsigmaelibf0030->SetLineColor(4); grsigmaelibf0030->SetLineWidth(1); grsigmaelibf0030->SetMarkerColor(4); grsigmaelibf0030->SetMarkerStyle(20); grsigmaelibf0030->SetTitle("electrons 0030"); mgibfsigmael->Add(grsigmaelibf0030); TGraphErrors* grsigmapiibf0030 = new TGraphErrors(nibfsigma,gainibfsigma,sigmapiibf0030,0,dsigmapiibf0030); grsigmapiibf0030->SetLineColor(4); grsigmapiibf0030->SetLineWidth(1); grsigmapiibf0030->SetMarkerColor(4); grsigmapiibf0030->SetMarkerStyle(20); grsigmapiibf0030->SetTitle("pions 0030"); mgibfsigmapi->Add(grsigmapiibf0030); //0530 TGraphErrors* grsigmaelibf0530 = new TGraphErrors(nibfsigma,gainibfsigma,sigmaelibf0530,0,dsigmaelibf0530); grsigmaelibf0530->SetLineColor(7); grsigmaelibf0530->SetLineWidth(1); grsigmaelibf0530->SetMarkerColor(7); grsigmaelibf0530->SetMarkerStyle(20); grsigmaelibf0530->SetTitle("electrons 0530"); mgibfsigmael->Add(grsigmaelibf0530); TGraphErrors* grsigmapiibf0530 = new TGraphErrors(nibfsigma,gainibfsigma,sigmapiibf0530,0,dsigmapiibf0530); grsigmapiibf0530->SetLineColor(7); grsigmapiibf0530->SetLineWidth(1); grsigmapiibf0530->SetMarkerColor(7); grsigmapiibf0530->SetMarkerStyle(20); grsigmapiibf0530->SetTitle("pions 0530"); mgibfsigmapi->Add(grsigmapiibf0530); //1030 TGraphErrors* grsigmaelibf1030 = new TGraphErrors(nibfsigma,gainibfsigma,sigmaelibf1030,0,dsigmaelibf1030); grsigmaelibf1030->SetLineColor(6); grsigmaelibf1030->SetLineWidth(1); grsigmaelibf1030->SetMarkerColor(6); grsigmaelibf1030->SetMarkerStyle(20); grsigmaelibf1030->SetTitle("electrons 1030"); mgibfsigmael->Add(grsigmaelibf1030); TGraphErrors* grsigmapiibf1030 = new TGraphErrors(nibfsigma,gainibfsigma,sigmapiibf1030,0,dsigmapiibf1030); grsigmapiibf1030->SetLineColor(6); grsigmapiibf1030->SetLineWidth(1); grsigmapiibf1030->SetMarkerColor(6); grsigmapiibf1030->SetMarkerStyle(20); grsigmapiibf1030->SetTitle("pions 1030"); mgibfsigmapi->Add(grsigmapiibf1030); mgibfsigmael->Draw("ap"); mgibfsigmael->GetXaxis()->SetTitle("HV scaling factor (%)"); mgibfsigmael->GetYaxis()->SetTitle("#sigma(dE/dx) electrons (%)"); mgibfsigmael->GetXaxis()->SetLimits(99,108); mgibfsigmael->GetYaxis()->SetRangeUser(9,12); TLegend *legibfsigmael=new TLegend(0.6,0.6,0.95,.9); legibfsigmael->SetBorderSize(1); legibfsigmael->SetFillColor(10); legibfsigmael->AddEntry(grsigmaelibf0025,"truncations [ 0,75]","p"); legibfsigmael->AddEntry(grsigmaelibf0525,"truncations [ 5,75]","p"); legibfsigmael->AddEntry(grsigmaelibf1025,"truncations [10,75]","p"); legibfsigmael->AddEntry(grsigmaelibf0030,"truncations [ 0,70]","p"); legibfsigmael->AddEntry(grsigmaelibf0530,"truncations [ 5,70]","p"); legibfsigmael->AddEntry(grsigmaelibf1030,"truncations [10,70]","p"); legibfsigmael->Draw("same"); cibfsigmael->Update(); cibfsigmael->Modified(); cibfsigmael->Write(); cibfsigmael->SaveAs(dirname+filename+"ibfSigmaEl.pdf"); TCanvas *cibfsigmapi = new TCanvas("cibfsigmapi","dEdx ibf gain sigma pions",1000,700); cibfsigmapi->SetFillColor(00); cibfsigmapi->SetGrid(); cibfsigmapi->cd(); mgibfsigmapi->Draw("ap"); mgibfsigmapi->GetXaxis()->SetTitle("HV scaling factor (%)"); mgibfsigmapi->GetYaxis()->SetTitle("#sigma(dE/dx) pions (%)"); mgibfsigmapi->GetXaxis()->SetLimits(99,108); mgibfsigmapi->GetYaxis()->SetRangeUser(10,14); TLegend *legibfsigmapi=new TLegend(0.6,0.6,0.95,.9); legibfsigmapi->SetBorderSize(1); legibfsigmapi->SetFillColor(10); legibfsigmapi->AddEntry(grsigmapiibf0025,"truncations [0,75]","p"); legibfsigmapi->AddEntry(grsigmapiibf0525,"truncations [5,75]","p"); legibfsigmapi->AddEntry(grsigmapiibf1025,"truncations [10,75]","p"); legibfsigmapi->AddEntry(grsigmapiibf0030,"truncations [0,70]","p"); legibfsigmapi->AddEntry(grsigmapiibf0530,"truncations [5,70]","p"); legibfsigmapi->AddEntry(grsigmapiibf1030,"truncations [10,70]","p"); legibfsigmapi->Draw("same"); cibfsigmapi->Update(); cibfsigmapi->Modified(); cibfsigmapi->Write(); cibfsigmapi->SaveAs(dirname+filename+"ibfSigmaPi.pdf"); //----------------------------------------------------------------------------------- savefile->Close(); }