const Int_t numberOfCentralityBins = 8; TString centralityArray[numberOfCentralityBins] = {"0-10","10-20","20-30","30-40","40-50","50-60","60-70","70-80"}; void drawCorrelationFunctionPsiSummarySummary(const char* lhcPeriod = "LHC11h", Int_t gTrainID = 250, Double_t psiMin = -0.5, Double_t psiMax = 3.5){ TFile *fPar[3][4]; TGraphErrors *gPar[3][4][18]; Int_t iCentrality[3] = {1,3,5}; TString sType[4] = {"PP","NN","PN","NP"}; for(Int_t iCent = 0 ; iCent < 3; iCent++){ for(Int_t iType = 0 ; iType < 4; iType++){ // open file fPar[iCent][iType] = TFile::Open(Form("PbPb/%s/Train%d/figs/correlationFunctionFit_Cent%d_%s_FitParameters.root",lhcPeriod,gTrainID,iCentrality[iCent],sType[iType].Data())); if(!fPar[iCent][iType]){ cerr<<"FILE "<