void epsilonScanFineBinningForPerformanceNote(Int_t radiusScale=1, Int_t epsScale=1,Int_t iOmegaTau = 0, Double_t x = 85., Double_t y = 0., Double_t z = 10.,Double_t integrateStep = 0.,Bool_t bLowE = kFALSE){ // // do for one gas (with omegaTau values accordingly) the correction // use the integrate along drift line option if integrateStep > 0 // else old method // // 1. Initialzation form space charge maps // AliTPCSpaceCharge3D *spaceCharge = new AliTPCSpaceCharge3D; // omega tau parameters and TF1 const Int_t nOmegaTau = 1; Double_t omegaTau[nOmegaTau] = {0.32}; TString tGas[nOmegaTau] = {"NeCO2_2"}; TString tGasOut[nOmegaTau] = {"NeCO2_2"}; TString sGas[nOmegaTau] = {"Ne-CO_{2}-N_{2} (90-10-5)"}; TString sLowE = ""; //if(bLowE) sLowE = "_lowE"; if(bLowE) sLowE = "_lowE2"; //second iteration const Int_t nEps = 51; Double_t epsilon[nEps]; for(Int_t iEps = 0; iEps < nEps; iEps++){ epsilon[iEps] = (Double_t)iEps; } TH1F *hR = new TH1F("hR",Form("distortion at x = %.0f, y = %.0f, z = %.0f",x,y,z),nEps,epsilon[0],epsilon[nEps-1]+1); TH1F *hPhiR = new TH1F("hPhiR",Form("distortion at x = %.0f, y = %.0f, z = %.0f",x,y,z),nEps,epsilon[0],epsilon[nEps-1]+1); //use always the integrate option here Double_t rate = 50.; for(Int_t iEps = 0; iEps < nEps; iEps++){ // select gas // 0 = NeCO2N2 cout<<"Open file = "<