void SCdistortionsForTDRfinalAllGases(Double_t radiusScale=1.5, Int_t epsScale=1){ // // do for one file (given by directory) the correction (specify the gas = iOmegaTau) // use the integrate along drift line option // // 1. Initialzation form space charge maps // AliTPCSpaceCharge3D *spaceCharge = new AliTPCSpaceCharge3D; const Double_t fgke0 = 8.854187817e-12; // vacuum permittivity [A·s/(V·m)] // omega tau parameters and TF1 const Int_t nOmegaTau = 4; Double_t omegaTau[nOmegaTau] = {0.32}; Double_t omegaTau[nOmegaTau] = {0.32,0.43,1.77,1.84}; TString tGas[nOmegaTau] = {"NeCO2_2","ArCO2","NeCF4","NeCF4_2"}; // CF4 is the same as CO2 here, but different omegaTau TString sGas[nOmegaTau] = {"Ne-CO_{2}-N_{2} (90-10-5)","Ar-CO_{2} (90-10)","Ne-CF_{4} (90-10)","Ne-CF_{4} (80-20)"}; Int_t eps[nOmegaTau] = {20,10,20,20}; Int_t col[nOmegaTau] = {kBlack,kRed,kOrange-3,kGreen+2}; TF1 * fdiffR[nOmegaTau]; TF1 * fdiffPhiR[nOmegaTau]; TH1F * hdiffR[nOmegaTau]; TH1F * hdiffPhiR[nOmegaTau]; TH2F * hMap[nOmegaTau]; TLegend *legend = new TLegend(0.25,0.6,0.85,0.85,Form("#rho_{SC} ~ r^{-%.1f} for 50 kHz",radiusScale),"brNDC"); setupLegend(legend,0.05); //use always the integrate option here Double_t integrateStep = 1.; TCanvas *cMap = new TCanvas("cMap","cMap",1200,900); cMap->Divide(2,2); TCanvas *cOmegaTau = new TCanvas("cOmegaTau","cOmegaTau",1200,500); cOmegaTau->Divide(2,1); TString outfilename = Form("SCdistortions_SC_NeCO2_50kHz_radiusScaling%.0f_epsScaling%d",radiusScale,epsScale); if(radiusScale>1.1 && radiusScale < 1.9) outfilename = Form("SCdistortions_SC_NeCO2_50kHz_radiusScaling%.1f_epsScaling%d",radiusScale,epsScale); //loop over gases for(Int_t iOmegaTau = 0; iOmegaTau < nOmegaTau; ++iOmegaTau){ cMap->cd(iOmegaTau+1)->SetPhi(150); // select gas // 0 = NeCO2N2 if(radiusScale>1.1 && radiusScale < 1.9){ cout<<"Open file = "<