const Int_t nEps = 2;
Int_t eps[nEps] = {20,10};
Int_t col[nEps] = {kBlack,kRed};
+ Int_t sty[nEps] = {1,2};
const Int_t nOmegaTau = 1;
Double_t omegaTau[nOmegaTau] = {0.32};
TString tGas[nOmegaTau] = {"NeCO2_2"};
- TString sGas[nOmegaTau] = {"Ne/CO_{2}/N_{2} (90-10-5)"};
+ TString sGas[nOmegaTau] = {"Ne-CO_{2}-N_{2} (90-10-5)"};
TF1 * fdiffR[nEps];
TF1 * fdiffPhiR[nEps];
TH1F * hdiffR[nEps];
//use always the integrate option here
Double_t integrateStep = 1.;
- TCanvas *cMap = new TCanvas("cMap","cMap",1200,500);
- cMap->Divide(2,1);
+ TCanvas *cMap[nEps];
TCanvas *cOmegaTau = new TCanvas("cOmegaTau","cOmegaTau",1200,500);
cOmegaTau->Divide(2,1);
//loop over epsilons
for(Int_t iEps = 0; iEps < nEps; ++iEps){
- cMap->cd(iEps+1)->SetPhi(150);
+ cMap[iEps] = new TCanvas(Form("cMap%d",iEps),Form("cMap%d",iEps),600,500);
+ cMap[iEps]->SetPhi(150);
// select gas
// 0 = NeCO2N2
hMap[iEps]->SetTitleSize(0.05,"XYZ");
hMap[iEps]->SetTitleOffset(1.5,"XY");
hMap[iEps]->SetTitleOffset(0.9,"Z");
- hMap[iEps]->SetTitle(Form("%s: 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
+ hMap[iEps]->SetTitle(Form("%s: 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
hMap[iEps]->DrawCopy("surf2fb");
hdiffR[iEps]->SetLineColor(col[iEps]);
hdiffPhiR[iEps]->SetLineColor(col[iEps]);
+
+ hdiffR[iEps]->SetLineStyle(sty[iEps]);
+ hdiffPhiR[iEps]->SetLineStyle(sty[iEps]);
hdiffR[iEps]->SetFillColor(col[iEps]);
hdiffPhiR[iEps]->SetFillColor(col[iEps]);
hdiffR[iEps]->SetTitleSize(0.05,"XYZ");
hdiffPhiR[iEps]->SetTitleSize(0.05,"XYZ");
- legend->AddEntry(hdiffR[iEps],Form("#epsilon = %d",eps[iEps]),"lp");
+ legend->AddEntry(hdiffR[iEps],Form("#varepsilon = %d",eps[iEps]),"lp");
cOmegaTau->cd(1);
if(iEps==0)
}
}
-
- cMap->SaveAs(Form("%s_SC_performanceNote.eps",outfilename.Data()));
- cMap->SaveAs(Form("%s_SC_performanceNote.pdf",outfilename.Data()));
+ for(Int_t iEps = 0; iEps < nEps; iEps++){
+ cMap[iEps]->SaveAs(Form("%s_epsilon%d_SC_performanceNote.eps",outfilename.Data(),eps[iEps]));
+ cMap[iEps]->SaveAs(Form("%s_epsilon%d_SC_performanceNote.pdf",outfilename.Data(),eps[iEps]));
+ }
cOmegaTau->SaveAs(Form("%s_Distortions_performanceNote.eps",outfilename.Data()));
cOmegaTau->SaveAs(Form("%s_Distortions_performanceNote.pdf",outfilename.Data()));
}
const Int_t nOmegaTau = 5;
Double_t omegaTau[nOmegaTau] = {0.34,0.32,0.43,1.77,1.84};
TString tGas[nOmegaTau] = {"NeCO2","NeCO2_2","ArCO2","NeCF4","NeCF4_2"}; // CF4 is the same as CO2 here, but different omegaTau
- TString sGas[nOmegaTau] = {"NeCO_{2} (90-10)","Ne/CO_{2}/N_{2} (90-10-5)","ArCO_{2} (90-10)","NeCF_{4} (90-10)","NeCF_{4} (80-20)"};
+ TString sGas[nOmegaTau] = {"Ne-CO_{2} (90-10)","Ne-CO_{2}-N_{2} (90-10-5)","Ar-CO_{2} (90-10)","Ne-CF_{4} (90-10)","Ne-CF_{4} (80-20)"};
TF2 * fdiffR[nEps];
TF2 * fdiffPhiR[nEps];
TH2F * hdiffR[nEps];
hMap[iEps]->SetTitleSize(0.05,"XYZ");
//hMap[iEps]->SetTitleOffset(1.5,"XY");
//hMap[iEps]->SetTitleOffset(0.9,"Z");
- hMap[iEps]->SetTitle(Form("#rho_{SC} (fC/cm^{3}) for %s, 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
+ hMap[iEps]->SetTitle(Form("#rho_{SC} (fC/cm^{3}) for %s, 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
hMap[iEps]->DrawCopy("colz");
// draw the distortion maps
hDistRMap[iEps]->SetTitleSize(0.05,"XYZ");
//hDistRMap[iEps]->SetTitleOffset(1.5,"XY");
//hDistRMap[iEps]->SetTitleOffset(0.9,"Z");
- hDistRMap[iEps]->SetTitle(Form("dr (cm) for %s, 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
+ hDistRMap[iEps]->SetTitle(Form("dr (cm) for %s, 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
hDistRMap[iEps]->DrawCopy("colz");
cDistRPMap->cd(iEps+1);
hDistRPMap[iEps]->SetTitleSize(0.05,"XYZ");
//hDistRPMap[iEps]->SetTitleOffset(1.5,"XY");
//hDistRPMap[iEps]->SetTitleOffset(0.9,"Z");
- hDistRPMap[iEps]->SetTitle(Form("d(r#varphi) (cm) for %s, 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
+ hDistRPMap[iEps]->SetTitle(Form("d(r#varphi) (cm) for %s, 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
hDistRPMap[iEps]->DrawCopy("colz");
hdiffR[iEps]->SetTitleSize(0.05,"XYZ");
//hdiffR[iEps]->SetTitleOffset(1.5,"XY");
//hdiffR[iEps]->SetTitleOffset(0.9,"Z");
- hdiffR[iEps]->SetTitle(Form("dr (cm) for %s, 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
+ hdiffR[iEps]->SetTitle(Form("dr (cm) for %s, 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
hdiffR[iEps]->DrawCopy("colz");
hdiffPhiR[iEps]->SetTitleSize(0.05,"XYZ");
//hdiffPhiR[iEps]->SetTitleOffset(1.5,"XY");
//hdiffPhiR[iEps]->SetTitleOffset(0.9,"Z");
- hdiffPhiR[iEps]->SetTitle(Form("d(r#varphi) (cm) for %s, 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
+ hdiffPhiR[iEps]->SetTitle(Form("d(r#varphi) (cm) for %s, 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
hdiffPhiR[iEps]->DrawCopy("colz");
cDistRIntMap->cd(iEps+1);
hdiffIntR[iEps]->SetTitleSize(0.05,"XYZ");
//hdiffIntR[iEps]->SetTitleOffset(1.5,"XY");
//hdiffIntR[iEps]->SetTitleOffset(0.9,"Z");
- hdiffIntR[iEps]->SetTitle(Form("dr (cm) for %s, 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
+ hdiffIntR[iEps]->SetTitle(Form("dr (cm) for %s, 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
hdiffIntR[iEps]->DrawCopy("colz");
hdiffIntPhiR[iEps]->SetTitleSize(0.05,"XYZ");
//hdiffIntPhiR[iEps]->SetTitleOffset(1.5,"XY");
//hdiffIntPhiR[iEps]->SetTitleOffset(0.9,"Z");
- hdiffIntPhiR[iEps]->SetTitle(Form("d(r#varphi) (cm) for %s, 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
+ hdiffIntPhiR[iEps]->SetTitle(Form("d(r#varphi) (cm) for %s, 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
hdiffIntPhiR[iEps]->DrawCopy("colz");
cDistZIntMap->cd(iEps+1);
hdiffIntZ[iEps]->SetTitleSize(0.05,"XYZ");
//hdiffIntZ[iEps]->SetTitleOffset(1.5,"XY");
//hdiffIntZ[iEps]->SetTitleOffset(0.9,"Z");
- hdiffIntZ[iEps]->SetTitle(Form("dz (cm) for %s, 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
+ hdiffIntZ[iEps]->SetTitle(Form("dz (cm) for %s, 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iEps]));
hdiffIntZ[iEps]->DrawCopy("colz");
}
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)","ArCO_{2} (90-10)","NeCF_{4} (90-10)","NeCF_{4} (80-20)"};
+ 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];
hMap[iOmegaTau]->SetTitleSize(0.05,"XYZ");
hMap[iOmegaTau]->SetTitleOffset(1.5,"XY");
hMap[iOmegaTau]->SetTitleOffset(0.9,"Z");
- hMap[iOmegaTau]->SetTitle(Form("%s: 50 kHz, #epsilon = %d",sGas[iOmegaTau].Data(),eps[iOmegaTau]));
+ hMap[iOmegaTau]->SetTitle(Form("%s: 50 kHz, #varepsilon = %d",sGas[iOmegaTau].Data(),eps[iOmegaTau]));
hMap[iOmegaTau]->DrawCopy("surf2fb");
hdiffR[iOmegaTau]->SetTitleSize(0.05,"XYZ");
hdiffPhiR[iOmegaTau]->SetTitleSize(0.05,"XYZ");
- legend->AddEntry(hdiffR[iOmegaTau],Form("%s (#epsilon = %d)",sGas[iOmegaTau].Data(),eps[iOmegaTau]),"lp");
+ legend->AddEntry(hdiffR[iOmegaTau],Form("%s (#varepsilon = %d)",sGas[iOmegaTau].Data(),eps[iOmegaTau]),"lp");
cOmegaTau->cd(1);
if(iOmegaTau==0)