#include #include #include #include #include #include #include #include using namespace std; Float_t lowWplot = 1.0; Float_t highWplot = 3.0; //for 5G //Float_t highWplot = 1.2; //to isolate elastic //Float_t highWplot = 2.5; //for 4G Float_t divideWbins = 5.0; Int_t Wbplot = int(((float(numWbins)*(highWplot-lowWplot)/(highW - lowW))/divideWbins)+0.5); Float_t lowxplot = 0.0; Float_t highxplot = 1.0; Int_t xbplot = Wbplot; int NBX, NBY; double LOX, HIX, LOY, HIY; void createHistos(); TH1D* hi_gcpart = new TH1D("hi_gcpart", "gcpart", 80, -1, 5); TH1D* hi_gpart = new TH1D("hi_gpart", "gpart", 80, -1, 12); TH1D* hi_clas_z = new TH1D("hi_clas_z", "Z of vertex as of CLAS, in BONUS's coordinate system.", 300, -200, 200); TH1D* hi_beta = new TH1D("hi_beta", "beta", 300, -0.1, 2.1); TH1D* hi_pe = new TH1D("hi_pe", "electron Momentum", 100, 0.0, 2.5); TH1D* hi_ppro = new TH1D("hi_ppro", "CLAS proton Momentum", 100, 0.0, 2.5); TH1D* hi_ppim = new TH1D("hi_ppim", "CLAS pi minus Momentum", 100, 0.0, 2.5); TH1D* hi_m = new TH1D("hi_m", "CLAS Mass", 250, -1, 1); TH1D* hi_emiss = new TH1D("hi_emiss", "missing energy", 250, 0.5, 2.5); TH1D* hwb4corr[6]; TH1D* hwafcorr[6]; TH2F* h2_dpVphisec[12]; TProfile* prof[12]; TH2F* h2_ratVx = new TH2F("h2_ratVx", "", 40, 0.1, 1, 40, 0, 1.2); TProfile* prof_ratVx; TH2D* h2_pimCorrVw = new TH2D("h2_pimCorrVw", "", Wbplot, lowWplot, highWplot, 500, 0.95, 1.01); TProfile* prof_pimCorrVw; TH2D* h2_epemCorrVw = new TH2D("h2_epemCorrVw", "", Wbplot, lowWplot, highWplot, 500, 0.95, 1.01); TProfile* prof_epemCorrVw; TH2D* h2_jxAcc_w = new TH2D("h2_jxAcc_w", "", Wbplot, lowWplot, highWplot, 100, 0.0, 1.0); TProfile* prof_jxAcc_w; TH2D* h2_jxAcc_x = new TH2D("h2_jxAcc_x", "", xbplot, lowxplot, highxplot, 100, 0.0, 1.0); TProfile* prof_jxAcc_x; TH2D* h2_jxEff = new TH2D("h2_jxEff", "", Wbplot, lowWplot, highWplot, 100, 0.0, 1.0); TProfile* prof_jxEff; TH1D* dqdx_rat_pbinned[20]; TH1D* hi_w = new TH1D("hi_w", "", Wbplot, lowWplot, highWplot); TH1D* hi_wdelta0 = new TH1D("hi_wdelta0", "invariant mass of final state", Wbplot, lowWplot, highWplot); TH1D* hi_wstar = new TH1D("hi_wstar", "untag(black) tag(red)", Wbplot, lowWplot, highWplot); TH1D* hi_wbkg = new TH1D("hi_wbkg", "W bg", Wbplot, lowWplot, highWplot); TH1D* hi_wstarbkg = new TH1D("hi_wstarbkg", "W* bg", Wbplot, lowWplot, highWplot); TH1D* hi_wstarcorr = new TH1D("hi_wstarcorr", "W* bg corrected", Wbplot, lowWplot, highWplot); TH1D* hi_wstarcorr2 = new TH1D("hi_wstarcorr2", "", Wbplot, lowWplot, highWplot); TH1D* hi_wrat = new TH1D("hi_wrat", "", Wbplot, lowWplot, highWplot); TH1D* hi_w_pro = new TH1D("hi_w_pro", "W for 'protons'", Wbplot, lowWplot, highWplot); TH1D* hi_w_other = new TH1D("hi_w_other", "W for 'others'", Wbplot, lowWplot, highWplot); TH1D* hi_x_bj_other = new TH1D("hi_x_bj_other", "x_bj for 'others'", 200, 0.0,10.0); TH1D* hi_x_bj_pro = new TH1D("hi_x_bj_pro", "x_bj for 'protons'", 200, 0.0,10.0); TH1D* hi_Q2 = new TH1D("hi_Q2", "Q2", 120, 0.2, 2.8); TH1D* hi_mmiss = new TH1D("hi_mmiss", "Missing Mass", 200, 0.5, 1.5); TH1D* hi_x = new TH1D("hi_x", "", xbplot, lowxplot, highxplot); TH1D* hi_xstar = new TH1D("hi_xstar", "untag(black) tag(red)", xbplot, lowxplot, highxplot); TH1D* hi_xrat = new TH1D("hi_xrat", "", xbplot, lowxplot, highxplot); TH1D* hi_wmod = new TH1D("hi_wmod", "", Wbplot, lowWplot, highWplot); TH1D* hi_wstarmod = new TH1D("hi_wstarmod", "untag(black) tag(red)", Wbplot, lowWplot, highWplot); TH1D* hi_wratmod = new TH1D("hi_wratmod", "", Wbplot, lowWplot, highWplot); TH1D* hi_xmod = new TH1D("hi_xmod", "", xbplot, lowxplot, highxplot); TH1D* hi_xstarmod = new TH1D("hi_xstarmod", "untag(black) tag(red)", xbplot, lowxplot, highxplot); TH1D* hi_xratmod = new TH1D("hi_xratmod", "", xbplot, lowxplot, highxplot); TH1D* hi_wratstate = new TH1D("hi_wratstate", "", Wbplot, lowWplot, highWplot); TH1D* hi_wratsyse = new TH1D("hi_wratsyse", "", Wbplot, lowWplot, highWplot); TH1D* hi_wratmodstate = new TH1D("hi_wratmodstate", "", Wbplot, lowWplot, highWplot); TH1D* hi_wratmodsyse = new TH1D("hi_wratmodsyse", "", Wbplot, lowWplot, highWplot); TH1D* hi_xratstate = new TH1D("hi_xratstate", "", xbplot, lowxplot, highxplot); TH1D* hi_xratsyse = new TH1D("hi_xratsyse", "", xbplot, lowxplot, highxplot); TH1D* hi_xratmodstate = new TH1D("hi_xratmodstate", "", xbplot, lowxplot, highxplot); TH1D* hi_xratmodsyse = new TH1D("hi_xratmodsyse", "", xbplot, lowxplot, highxplot); TH1D* hi_du_w = new TH1D("hi_du_w", "", Wbplot, lowWplot, highWplot); TH1D* hi_du_x = new TH1D("hi_du_x", "", xbplot, lowxplot, highxplot); TH1D* hi_den_w = new TH1D("hi_den_w", "", Wbplot, lowWplot, highWplot); TH1D* hi_den_x = new TH1D("hi_den_x", "", xbplot, lowxplot, highxplot); TH1D* hi_D2Data = new TH1D("hi_D2Data", "Inclusive Data", Wbplot, lowWplot, highWplot); TH1D* hi_sigrD2 = new TH1D("hi_sigrD2", "D2 radiated Model", Wbplot, lowWplot, highWplot); TH1D* hi_accCorr_w = new TH1D("hi_accCorr", "(data/model)", Wbplot, lowWplot, highWplot); TH1D* hi_accCorr_x = new TH1D("hi_accCorr", "(data/model)", xbplot, lowxplot, highxplot); TH1D* hi_pimCorr = new TH1D("hi_pimCorr", "corr to tag/untag", Wbplot, lowWplot, highWplot); // -------------- BONUS --------------------------------------- TH1D* hi_bo_px = new TH1D("hi_bo_px", "hi_bo_px", 120, -150, 150); TH1D* hi_bo_py = new TH1D("hi_bo_py", "hi_bo_py", 120, -150, 150); TH1D* hi_bo_pz = new TH1D("hi_bo_pz", "hi_bo_pz", 50, -150, 150); TH1D* hi_pre_px = new TH1D("hi_pre_px", "Bonus(red) Clas(Blue)", 120, -150, 150); TH1D* hi_pre_py = new TH1D("hi_pre_py", "Bonus(red) Clas(Blue)", 120, -150, 150); TH1D* hi_pre_pz = new TH1D("hi_pre_pz", "Bonus(red) Clas(Blue)", 50, -150, 150); TH1D* hi_pre_th = new TH1D("hi_pre_th", "Predicted Spec Theta", 120, -180, 180); TH1D* hi_pre_ph = new TH1D("hi_pre_ph", "Predicted Spec Phi", 120, 0, 360); TH1D* hi_dpx = new TH1D("hi_dpx", "hi_dpx", 120, -100, 100); TH1D* hi_dpy = new TH1D("hi_dpy", "hi_dpy", 120, -100, 100); TH1D* hi_dpz = new TH1D("hi_dpz", "hi_dpz", 120, -100, 100); TH1D* bo_z_back = new TH1D("bo_z_back", "bo_z_back", 300, -100, 150); TH1D* hi_bo_theta = new TH1D("hi_bo_theta", "Scattering angle CLAS(blue), BoNuS(red)", 100, 0, 180); TH1D* hi_theta_clas = new TH1D("hi_theta_clas", "Scattering angle CLAS(blue), BoNuS(red)", 100,0, 180); TH1D* hi_theta_diff = new TH1D("hi_theta_diff", "Difference of calculated and reconstructed angles", 200, -0.4, 0.4); TH1D* hi_bo_phi = new TH1D("hi_bo_phi", "Azimuthal angle in CLAS(blue),BONUS(red)", 200,0,360); TH1D* hi_phi_clas = new TH1D("hi_phi_clas", "Azimuthal angle in CLAS(blue),BONUS(red)", 200,0,360); TH1D* hi_phi_diff = new TH1D("hi_phi_diff", "d#phi", 120,-20,20); TH1D* hi_dphi = new TH1D("hi_dphi", "(elastic only) phi_clas_proton - phi_clas_electron", 200,-360, 360); TH1D* hi_dphi2 = new TH1D("hi_dphi2", "phi_bo_proton - phi_bonus_electron", 200, 150,210); TH1D* hi_theta_meas = new TH1D("hi_theta_meas", "d#theta", 120,-20,20); TH1D* hi_mom_diff = new TH1D("hi_mom_diff", "Momentum Resolution",100, -1,1); TH1D* hi_pt_diff = new TH1D("hi_pt_diff", "P_{perp} Resolution",100, -1,1); TH1D* hi_bo_z = new TH1D("hi_bo_z", "Z of vertex as of BONUS", 300,-100,150); TH1D* hi_diff_z = new TH1D("hi_diff_z", "Z(CLAS uncorrected blue)-Z(BoNuS)",120,-50,50); TH1D* hi_cutdiff_z = new TH1D("hi_cutdiff_z", "Z(CLAS corrected)-Z(BoNuS)",120,-50,50); TH1D* hi_clasdiff_z = new TH1D("hi_clasdiff_z" ,"Z(CLAS)-Z(CLAS_AVERAGE)",120,-1,1); TH1F* hi_r = new TH1F("hi_r", "track radius of curvature", 200, -100, 100); TH1F* hi_x2 = new TH1F("hi_x2", "chi^2 of helix fit to track", 200, 0, 10); TH1F* hi_npt = new TH1F("hi_npt", "number of pads hit in track",40, 0, 40); TH1D* hi_edist = new TH1D("hi_edist", "edist",120,-20,20); TH1D* hi_sdist = new TH1D("hi_sdist", "sdist",120,-20,20); TH1D* hi_p_corr = new TH1D("hi_p_corr", "p corrected for energy loss - All cuts except VIP", 100,0,300); TH1D* hi_p_rtpc = new TH1D("hi_p_rtpc", "BoNuS Proton Momentum(red) CLAS(blue)", 50,0,200); TH1D* hi_pt_rtpc = new TH1D("hi_pt_rtpc", "BoNuS Proton p_{t}(red) CLAS(blue)", 50,0,200); TH1D* hi_p_raw = new TH1D("hi_p_raw", "raw(green) protons(red) others(blue) pro_corrected(orange)", 100,0,300); TH1D* hi_p_pro = new TH1D("hi_p_pro", "", 100,0,300); TH1D* hi_p_others = new TH1D("hi_p_others", "", 100,0,300); TH1D* hi_p_pre = new TH1D("hi_p_pre", "BoNuS Proton Momentum(red) CLAS(blue)", 50,0,200); TH1D* hi_pt_m = new TH1D("hi_pt_m", "BoNuS Proton p_{t}(red) CLAS(blue)", 50,0,200); TH1D* hi_betheprof = new TH1D("hi_betheprof", "Profile of dqdx vs p", 150,0,400); TH1D* hi_diff_dqdx = new TH1D("hi_diff_dqdx", "", 150,-100,100); TH1D* hi_rat_dqdx = new TH1D("hi_rat_dqdx", "", 150,0,4); TH1D* hi_edistL = new TH1D("hi_edistL", "edist Left",120,-20,20); TH1D* hi_sdistL = new TH1D("hi_sdistL", "sdist Left",120,-20,20); TH1D* hi_dzL = new TH1D("hi_dzL", "Z(CLAS corrected)-Z(BoNuS)",120,-50,50); TH1D* hi_dthetaL = new TH1D("hi_dthetaL","d#theta", 120,-20,20); TH1D* hi_dphiL = new TH1D("hi_dphiL", "d#phi", 120,-20,20); TH1D* hi_edistR = new TH1D("hi_edistR", "edist Right",120,-20,20); TH1D* hi_sdistR = new TH1D("hi_sdistR", "sdist Right",120,-20,20); TH1D* hi_dzR = new TH1D("hi_dzR", "Z(CLAS corrected)-Z(BoNuS)",120,-50,50); TH1D* hi_dthetaR = new TH1D("hi_dthetaR","d#theta", 120,-20,20); TH1D* hi_dphiR = new TH1D("hi_dphiR", "d#phi", 120,-20,20); /////2 dimensional histograms//////////// TH2D* h2_ed_sd = new TH2D("h2_ed_sd","edist vs sdist",120,-20,20,120,-20,20); TH2D* h2_ed_sd_dz = new TH2D("h2_ed_sd_dz","edist vs sdist vs dz",120,-20,20,120,-20,20); TH2D* h2_phi = new TH2D("h2_phi","bonus_phi vs clas_phi",100,-1,7,100,-1,7); TH2D* h2_phiVdth = new TH2D("h2_phiVdth","phi vs dtheta",50,-10,10,50,0,360); TH2D* h2_bethe = new TH2D("h2_bethe","dQ/dx vs momentum ALL",500,0,250,500,0,200); TH2D* h2_bethepro = new TH2D("h2_bethepro","dQ/dx vs momentum (pass RTPC proton cut)",500,0,400,500,0,600); TProfile* prof_bethe; TH2D* h2_avgadc = new TH2D("h2_avgadc","Avg ADC per hit vs raw momentum",500,0,400,500,0,100); TH2D* h2_avgadc_st = new TH2D("h2_avgadc_st","Avg ADC per hit * sin(theta) vs raw momentum",500,0,400,500,0,100); TH2D* h2_nhits = new TH2D("h2_nhits","Number of hits vs raw momentum",500,0,400,100,0,100); TH2D* h2_mom = new TH2D("h2_mom","BoNuS momentum vs CLAS mom",500,0,200,500,0,200); TH2D* h2_sd = new TH2D("h2_sd","sdist vs phi_track vs z_track",100,-120,120,100,-0.5,2*PI+0.5); TH2D* h2_chiVr = new TH2D("h2_chiVr","chi^2 of line fit vs 1/(fit radius)",250,-0.04,0.04,120,0,1e-07); TH2D* h2_dphiVmom = new TH2D("h2_dphiVmom","dphi vs momentum ",50,100,1000,50,-10,10); TH2D* h2_dphiVdz = new TH2D("h2_dphiVdz","dphi(proton phi - electron phi) vs proton vertex",200,-100,100,60,150,210); TH2D* h2_dphiVpphi = new TH2D("h2_dphiVpphi","dphi(proton phi - electron phi-180) vs phi proton",360,0,360,60,-30,30); TH2D* h2_sdVz = new TH2D("h2_sdVz","sdist vs z",50,-100,100,50,-20,30); TH2D* h2_dpVphi = new TH2D("h2_dpVphi","dp vs phi proton",500,0,360,500,-1,1); TH2D* h2_dpVth = new TH2D("h2_dpVth","dp vs theta proton",500,0,180,500,-1,1); TH2D* h2_dthVdphi = new TH2D("h2_dthVdphi","dtheta vs dphi proton",50,-10,10,50,-10,10); TH2D* h2_dpVchi2 = new TH2D("h2_dpVchi2","dp vs chi^2",100,0,1.6,2000,-1000,1000); TH2D* h2_dpVp = new TH2D("h2_dpVp","dp/p vs measured momentum",100,0,250,100,-1,1); TProfile* prof_dpVp; TH2D* h2_bethewstar = new TH2D("h2_bethewstar", "Profile dqdx vs p vs W*", Wbplot, lowWplot, highWplot,150,0,400); TH2D* h2_ps_w_sig = new TH2D("h2_ps_w_sig", "p_{spec} vs W vs #sigma_{n}/#sigma_{d}", Wbplot, lowWplot, highWplot,25,0,250); TH2D* h2_theta = new TH2D("h2_theta","bonus_theta vs clas_theta (electrons)",250,0,60,250,0,60); TH2D* h2_zbVzc = new TH2D("h2_zVz","z_BoNuS vs z_CLAS (electrons)",300, -60,100,300,-60,100); TH2D* h2_zVw = new TH2D("h2_zVw","z vertex vs corrected W",300, lowWplot, highWplot,300,-100,175); TH2D* h2_zVphi_nocorr = new TH2D("h2_zVphi_nocorr","z vertex vs phi (nocorr)",300, 0,360,300,120,160); TH2D* h2_zVphi_corr = new TH2D("h2_zVphi_corr","z vertex vs phi (corr)",300, 0, 360,300,120,160); TH2D* h2_wVphib4 = new TH2D("h2_wVphib4","w vs phi before corr",300,0,360 ,300, 0.9,1.0); TH2D* h2_wVphiaf = new TH2D("h2_wVphiaf","w vs phi after corr",300,0,360 ,300, 0.9,1.0); TH2D* h2_dpVdcb4 = new TH2D("h2_dpVdcb4","dp/p vs phi (W<1.0) before corr",300,0,360 ,300, -0.1,0.1); TH2D* h2_dpVdcaf = new TH2D("h2_dpVdcaf","dp/p vs phi (W<1.0) after corr",300,0,360 ,300, -0.1,0.1); TH2D* h2_wVQ2 = new TH2D("h2_wVQ2","w vs Q^2 ",300,0 ,3.5 ,300, 0.5,3.5); TH2D* h2_psVcostpq = new TH2D("h2_psVcostpq","W^{2} < 3.5",300,-1 ,1 ,300, 40,200); TH2D* h2_psValpha = new TH2D("h2_psValpha","W^{2} < 3.5",300,0.7 ,1.3 ,300, 50,200); TH2D* h2_psVcostpqgt = new TH2D("h2_psVcostpqgt","W^{2} > 3.5",300,-1 ,1 ,300, 40,200); TH2D* h2_psValphagt = new TH2D("h2_psValphagt","W^{2} > 3.5",300,0.7 ,1.3 ,300, 50,200); TH2D* hi_pmzVpmx = new TH2D("hi_pmzVpmx","pmz VS pmx",120,-150 ,150 ,120,-150,150); TH2D* hi_pmzVpmy = new TH2D("hi_pmzVpmy","pmz VS pmy",120,-150 ,150 ,120,-150,150); TH2D* h2_wVwstar = new TH2D("h2_wVwstar","w vs w* ",Wbplot, lowWplot, highWplot,Wbplot, lowWplot, highWplot); TH2D* h2_deltaVps = new TH2D("h2_deltaVps","delta vs ps ",300,40 ,200 ,300, -0.1,0.5); void createHistos() { hi_gcpart->SetLineWidth(2); hi_gcpart->SetLineColor(8); hi_gcpart->GetXaxis()->SetTitle("gcpart"); hi_gcpart->GetYaxis()->SetTitle("Entries"); hi_gpart->SetLineWidth(2); hi_gpart->SetLineColor(4); hi_gpart->GetXaxis()->SetTitle("gpart"); hi_gpart->GetYaxis()->SetTitle("Entries"); hi_clas_z->SetLineWidth(2); hi_clas_z->SetLineColor(6); hi_clas_z->GetXaxis()->SetTitle("Z in CLAS"); hi_clas_z->GetYaxis()->SetTitle("Entries"); hi_beta->SetLineWidth(2); hi_beta->SetLineColor(9); hi_beta->GetXaxis()->SetTitle("beta"); hi_beta->GetYaxis()->SetTitle("Entries"); hi_pe->SetLineWidth(2); hi_pe->SetLineColor(13); hi_pe->GetXaxis()->SetTitle("Electron Momentum, GeV/c"); hi_pe->GetYaxis()->SetTitle("Entries"); hi_ppro->SetLineWidth(2); hi_ppro->SetLineColor(13); hi_ppro->GetXaxis()->SetTitle("CLAS Proton Momentum, GeV/c"); hi_ppro->GetYaxis()->SetTitle("Entries"); hi_ppim->SetLineWidth(2); hi_ppim->SetLineColor(13); hi_ppim->GetXaxis()->SetTitle("CLAS pi minus Momentum, GeV/c"); hi_ppim->GetYaxis()->SetTitle("Entries"); hi_m->SetLineWidth(2); hi_m->SetLineColor(13); hi_m->GetXaxis()->SetTitle("CLAS m, GeV"); hi_m->GetYaxis()->SetTitle("Entries"); hi_emiss->SetLineWidth(2); hi_emiss->SetLineColor(13); hi_emiss->GetXaxis()->SetTitle("missing energy, GeV/c"); hi_emiss->GetYaxis()->SetTitle("Entries"); hi_x->SetLineWidth(1); hi_x->SetLineColor(1); hi_x->SetMarkerColor(12); hi_x->SetMarkerStyle(kFullCircle); hi_x->SetMarkerSize(0.5); hi_x->GetXaxis()->SetTitle("x_{B}"); hi_x->GetYaxis()->SetTitleOffset(0.9); hi_x->GetYaxis()->SetTitleSize(.06); hi_x->GetXaxis()->SetTitleSize(.06); hi_x->GetXaxis()->SetLabelSize(.07); hi_x->GetYaxis()->SetLabelSize(.07); hi_x->Sumw2(); hi_xstar->SetLineWidth(1); hi_xstar->SetLineColor(1); hi_xstar->SetMarkerColor(kRed); hi_xstar->SetMarkerStyle(kFullCircle); hi_xstar->SetMarkerSize(0.5); hi_xstar->GetXaxis()->SetTitle("x*"); hi_xstar->GetYaxis()->SetTitleOffset(0.9); hi_xstar->GetYaxis()->SetTitleSize(.06); hi_xstar->GetXaxis()->SetTitleSize(.06); hi_xstar->GetXaxis()->SetLabelSize(.07); hi_xstar->GetYaxis()->SetLabelSize(.07); hi_xstar->Sumw2(); hi_xrat->SetMarkerColor(12); hi_xrat->SetMarkerStyle(kFullCircle); hi_xrat->SetMarkerSize(0.5); hi_xrat->GetXaxis()->SetTitle("x_{B}"); hi_xrat->GetYaxis()->SetTitle("F_{2}^{n}/F_{2}^{d}"); hi_xrat->GetYaxis()->SetTitleOffset(0.9); hi_xrat->GetYaxis()->SetTitleSize(.06); hi_xrat->GetXaxis()->SetTitleSize(.06); hi_xrat->GetXaxis()->SetLabelSize(.07); hi_xrat->GetYaxis()->SetLabelSize(.07); hi_xrat->Sumw2(); hi_du_x->SetMarkerColor(12); hi_du_x->SetMarkerStyle(kFullCircle); hi_du_x->SetMarkerSize(0.5); hi_du_x->GetXaxis()->SetTitle("x_{B}"); hi_du_x->GetYaxis()->SetTitle("d/u"); hi_du_x->GetYaxis()->SetTitleOffset(0.9); hi_du_x->GetYaxis()->SetTitleSize(.06); hi_du_x->GetXaxis()->SetTitleSize(.06); hi_du_x->GetXaxis()->SetLabelSize(.07); hi_du_x->GetYaxis()->SetLabelSize(.07); hi_du_x->Sumw2(); hi_du_w->SetMarkerColor(12); hi_du_w->SetMarkerStyle(kFullCircle); hi_du_w->SetMarkerSize(0.5); hi_du_w->GetXaxis()->SetTitle("W(GeV)"); hi_du_w->GetYaxis()->SetTitle("d/u"); hi_du_w->GetYaxis()->SetTitleOffset(0.9); hi_du_w->GetYaxis()->SetTitleSize(.06); hi_du_w->GetXaxis()->SetTitleSize(.06); hi_du_w->GetXaxis()->SetLabelSize(.07); hi_du_w->GetYaxis()->SetLabelSize(.07); hi_du_w->Sumw2(); hi_xmod->SetLineWidth(1); hi_xmod->SetLineColor(1); hi_xmod->SetMarkerColor(12); hi_xmod->SetMarkerStyle(kFullCircle); hi_xmod->SetMarkerSize(0.5); hi_xmod->GetXaxis()->SetTitle("x_{B}"); hi_xmod->Sumw2(); hi_xstarmod->SetLineWidth(1); hi_xstarmod->SetLineColor(1); hi_xstarmod->SetMarkerColor(kRed); hi_xstarmod->SetMarkerStyle(kFullCircle); hi_xstarmod->SetMarkerSize(0.5); hi_xstarmod->GetXaxis()->SetTitle("x*"); hi_xstarmod->Sumw2(); hi_xratmod->SetMarkerColor(12); hi_xratmod->SetMarkerStyle(kFullCircle); hi_xratmod->SetMarkerSize(0.5); hi_xratmod->GetXaxis()->SetTitle("x_{B}"); hi_xratmod->GetYaxis()->SetTitle("F_{2}^{n}/F_{2}^{p}"); hi_xratmod->GetYaxis()->SetTitleOffset(0.9); hi_xratmod->GetYaxis()->SetTitleSize(.06); hi_xratmod->GetXaxis()->SetTitleSize(.06); hi_xratmod->GetXaxis()->SetLabelSize(.07); hi_xratmod->GetYaxis()->SetLabelSize(.07); hi_xratmod->Sumw2(); hi_wmod->SetLineWidth(1); hi_wmod->SetLineColor(1); hi_wmod->SetMarkerColor(12); hi_wmod->SetMarkerStyle(kFullCircle); hi_wmod->SetMarkerSize(0.5); hi_wmod->GetXaxis()->SetTitle("W(GeV)"); hi_wmod->Sumw2(); hi_wstarmod->SetLineWidth(1); hi_wstarmod->SetLineColor(1); hi_wstarmod->SetMarkerColor(kRed); hi_wstarmod->SetMarkerStyle(kFullCircle); hi_wstarmod->SetMarkerSize(0.5); hi_wstarmod->GetXaxis()->SetTitle("W*(GeV)"); hi_wstarmod->Sumw2(); hi_wratmod->SetMarkerColor(12); hi_wratmod->SetMarkerStyle(kFullCircle); hi_wratmod->SetMarkerSize(0.5); hi_wratmod->GetXaxis()->SetTitle("W(GeV)"); hi_wratmod->GetYaxis()->SetTitle("F_{2}^{n}/F_{2}^{p}"); hi_wratmod->GetYaxis()->SetTitleOffset(0.9); hi_wratmod->GetYaxis()->SetTitleSize(.06); hi_wratmod->GetXaxis()->SetTitleSize(.06); hi_wratmod->GetXaxis()->SetLabelSize(.07); hi_wratmod->GetYaxis()->SetLabelSize(.07); hi_wratmod->Sumw2(); hi_den_w->Sumw2(); hi_den_x->Sumw2(); hi_wratstate->Sumw2(); hi_wratstate->SetLineColor(kRed); //hi_wratsyse->Sumw2(); hi_wratmodstate->Sumw2(); hi_wratmodstate->SetLineColor(kRed); //hi_wratmodsyse->Sumw2(); hi_xratstate->Sumw2(); hi_xratstate->SetLineColor(kRed); //hi_xratsyse->Sumw2(); hi_xratmodstate->Sumw2(); hi_xratmodstate->SetLineColor(kRed); //hi_xratmodsyse->Sumw2(); hi_D2Data->SetLineWidth(2); hi_D2Data->SetLineColor(1); hi_D2Data->GetXaxis()->SetTitle("W"); hi_D2Data->GetYaxis()->SetTitleOffset(1.35); hi_D2Data->Sumw2(); hi_sigrD2->SetLineWidth(2); hi_sigrD2->SetLineColor(1); hi_sigrD2->GetXaxis()->SetTitle("W"); hi_sigrD2->GetYaxis()->SetTitleOffset(1.35); hi_sigrD2->Sumw2(); hi_accCorr_w->SetLineWidth(2); hi_accCorr_w->SetLineColor(kRed); hi_accCorr_w->GetXaxis()->SetTitle("W, GeV"); hi_accCorr_w->GetYaxis()->SetTitleOffset(1.35); hi_accCorr_w->Sumw2(); hi_accCorr_x->SetLineWidth(2); hi_accCorr_x->SetLineColor(kRed); hi_accCorr_x->GetXaxis()->SetTitle("x_{B}"); hi_accCorr_x->GetYaxis()->SetTitleOffset(1.35); hi_accCorr_x->Sumw2(); hi_w->SetLineWidth(1);// hi_w->SetLineColor(1); hi_w->GetXaxis()->SetTitle("W (GeV)"); hi_w->SetMarkerColor(12); hi_w->SetMarkerStyle(kFullCircle); hi_w->SetMarkerSize(0.5); hi_w->GetYaxis()->SetTitleOffset(0.9); hi_w->GetYaxis()->SetTitleSize(.06); hi_w->GetXaxis()->SetTitleSize(.06); hi_w->GetXaxis()->SetLabelSize(.07); hi_w->GetYaxis()->SetLabelSize(.07); hi_w->Sumw2(); hi_wstar->SetLineWidth(1);// hi_wstar->SetLineColor(kRed); hi_wstar->SetFillColor(kRed); hi_wstar->SetMarkerColor(kRed); hi_wstar->SetMarkerStyle(kFullCircle); hi_wstar->SetMarkerSize(0.5); hi_wstar->GetXaxis()->SetTitle("W*(GeV)"); hi_wstar->GetYaxis()->SetTitleOffset(0.9); hi_wstar->GetYaxis()->SetTitleSize(.06); hi_wstar->GetXaxis()->SetTitleSize(.06); hi_wstar->GetXaxis()->SetLabelSize(.07); hi_wstar->GetYaxis()->SetLabelSize(.07); //hi_wstar->GetXaxis()->SetTitleSize(0.6); hi_wstar->Sumw2(); hi_wbkg->SetLineWidth(2); hi_wbkg->SetLineColor(1); hi_wbkg->GetXaxis()->SetTitle("W, GeV"); //hi_w->GetYaxis()->SetTitle("Entries"); hi_wbkg->Sumw2(); hi_wstarbkg->SetLineWidth(2); hi_wstarbkg->SetLineColor(kRed); hi_wstarbkg->SetFillColor(kRed); hi_wstarbkg->GetXaxis()->SetTitle("W, GeV"); //hi_wstar->GetYaxis()->SetTitle("Entries"); hi_wstarbkg->Sumw2(); hi_wstarcorr->SetLineWidth(2); hi_wstarcorr->SetLineColor(kRed); hi_wstarcorr->SetFillColor(kRed); hi_wstarcorr->GetXaxis()->SetTitle("W, GeV"); //hi_wstar->GetYaxis()->SetTitle("Entries"); hi_wstarcorr->Sumw2(); hi_wstarcorr2->SetLineWidth(2); hi_wstarcorr2->SetLineColor(kBlue); hi_wstarcorr2->SetFillColor(kRed); hi_wstarcorr2->GetXaxis()->SetTitle("W, GeV"); //hi_wstar->GetYaxis()->SetTitle("Entries"); hi_wstarcorr2->GetYaxis()->SetTitleOffset(1.35); hi_wstarcorr2->Sumw2(); hi_wrat->SetMarkerColor(12); hi_wrat->SetMarkerStyle(kFullCircle); hi_wrat->SetMarkerSize(0.5); hi_wrat->GetXaxis()->SetTitle("W, GeV"); hi_wrat->GetYaxis()->SetTitle("F_{2}^{n}/F_{2}^{d}"); hi_wrat->GetYaxis()->SetTitleOffset(0.9); hi_wrat->GetYaxis()->SetTitleSize(.06); hi_wrat->GetXaxis()->SetTitleSize(.06); hi_wrat->GetXaxis()->SetLabelSize(.07); hi_wrat->GetYaxis()->SetLabelSize(.07); hi_w_pro->SetLineWidth(2); hi_w_pro->SetLineColor(12); hi_w_pro->GetXaxis()->SetTitle("W, GeV"); hi_w_pro->GetYaxis()->SetTitle("Entries"); hi_w_pro->GetYaxis()->SetTitleOffset(1.35); hi_w_other->SetLineWidth(2); hi_w_other->SetLineColor(kRed); hi_w_other->GetXaxis()->SetTitle("W, GeV"); hi_w_other->GetYaxis()->SetTitle("Entries"); hi_w_other->GetYaxis()->SetTitleOffset(1.35); hi_x_bj_other->SetLineWidth(2); hi_x_bj_other->SetLineColor(kRed); hi_x_bj_other->GetXaxis()->SetTitle("x_bj, GeV"); hi_x_bj_other->GetYaxis()->SetTitle("Entries"); hi_x_bj_other->GetYaxis()->SetTitleOffset(1.35); hi_x_bj_pro->SetLineWidth(2); hi_x_bj_pro->SetLineColor(12); hi_x_bj_pro->GetXaxis()->SetTitle("x_bj, GeV"); hi_x_bj_pro->GetYaxis()->SetTitle("Entries"); hi_x_bj_pro->GetYaxis()->SetTitleOffset(1.35); hi_Q2->SetFillColor(kYellow); hi_Q2->GetXaxis()->SetTitle("Q^{2}, GeV^{2}"); hi_Q2->GetYaxis()->SetTitleOffset(1.35); hi_mmiss->SetLineWidth(2); hi_mmiss->SetLineColor(1); hi_mmiss->GetXaxis()->SetTitle("Missing Mass, GeV"); hi_mmiss->GetYaxis()->SetTitle("Entries"); hi_mmiss->GetYaxis()->SetTitleOffset(1.35); hi_bo_px->SetLineWidth(1); hi_bo_px->SetLineColor(kRed); hi_bo_px->GetXaxis()->SetTitle("bonus Px, MeV"); hi_bo_px->GetYaxis()->SetTitleOffset(1.5); hi_bo_py->SetLineWidth(1); hi_bo_py->SetLineColor(kRed); hi_bo_py->GetXaxis()->SetTitle("bonus Py, MeV"); hi_bo_py->GetYaxis()->SetTitleOffset(1.5); hi_bo_pz->SetLineWidth(1); hi_bo_pz->SetLineColor(kRed); hi_bo_pz->GetXaxis()->SetTitle("bonus Pz, MeV"); hi_pre_px->SetLineWidth(1); hi_pre_px->SetLineColor(kBlue); hi_pre_px->GetXaxis()->SetTitle("missing Px, MeV"); hi_pre_px->GetYaxis()->SetTitleOffset(1.5); hi_pre_py->SetLineWidth(1); hi_pre_py->SetLineColor(kBlue); hi_pre_py->GetXaxis()->SetTitle("missing Py, MeV"); hi_pre_py->GetYaxis()->SetTitleOffset(1.5); hi_pre_pz->SetLineWidth(1); hi_pre_pz->SetLineColor(kBlue); hi_pre_pz->GetXaxis()->SetTitle("missing Pz, MeV"); hi_pre_th->SetLineWidth(2); hi_pre_th->SetLineColor(kRed); hi_pre_th->GetXaxis()->SetTitle("missing spec theta, deg"); hi_pre_th->GetYaxis()->SetTitleOffset(1.5); hi_pre_ph->SetLineWidth(2); hi_pre_ph->SetLineColor(kRed); hi_pre_ph->GetXaxis()->SetTitle("missing spec phi, deg"); hi_pre_ph->GetYaxis()->SetTitleOffset(1.5); hi_dpx->SetLineWidth(2); hi_dpx->SetLineColor(kBlue); hi_dpx->GetXaxis()->SetTitle("missing px - bonus px, MeV"); hi_dpx->GetYaxis()->SetTitleOffset(1.5); hi_dpy->SetLineWidth(2); hi_dpy->SetLineColor(kBlue); hi_dpy->GetXaxis()->SetTitle("missing py - bonus py, MeV"); hi_dpy->GetYaxis()->SetTitleOffset(1.5); hi_dpz->SetLineWidth(2); hi_dpz->SetLineColor(kBlue); hi_dpz->GetXaxis()->SetTitle("missing Pz - bonus pz, MeV"); bo_z_back->GetXaxis()->SetTitle("bonus z, mm"); hi_bo_theta->SetLineWidth(2); hi_bo_theta->SetLineColor(kRed); hi_bo_theta->GetXaxis()->SetTitle("Theta, deg"); hi_theta_clas->SetLineWidth(2); hi_theta_clas->SetLineColor(kBlue); hi_theta_clas->GetXaxis()->SetTitle("Polar angle in CLAS, deg"); hi_theta_diff->SetLineWidth(1); hi_theta_diff->SetLineColor(6); hi_theta_diff->GetXaxis()->SetTitle("Theta_CLAS(calc.) - Theta_BONUS(rec-ed), rad"); hi_bo_phi->SetLineWidth(2); hi_bo_phi->SetLineColor(kRed); hi_bo_phi->GetXaxis()->SetTitle("Phi_BONUS(rec-ed), deg"); hi_phi_clas->SetLineWidth(2); hi_phi_clas->SetLineColor(kBlue); hi_phi_clas->GetXaxis()->SetTitle("Azimuthal angle in CLAS, deg"); hi_phi_diff->SetLineWidth(3); hi_phi_diff->SetLineColor(6); hi_phi_diff->GetXaxis()->SetTitle("Phi_CLAS - Phi_BONUS,deg"); hi_dphi->SetLineWidth(3); hi_dphi->SetLineColor(6); hi_dphi->GetXaxis()->SetTitle("dphi,deg"); hi_dphi2->SetLineWidth(3); hi_dphi2->SetLineColor(6); hi_dphi2->GetXaxis()->SetTitle("dphi2,deg"); hi_theta_meas->SetLineWidth(2); hi_theta_meas->SetLineColor(2); hi_theta_meas->GetXaxis()->SetTitle("Theta_CLAS - Theta_RTPC, deg"); hi_mom_diff->SetLineWidth(1); hi_mom_diff->GetXaxis()->SetTitle("(predicted p - measured p) / predicted p"); hi_pt_diff->SetLineWidth(1); hi_pt_diff->GetXaxis()->SetTitle("#Delta p_{t}/p_{t}"); hi_bo_z->SetLineWidth(2); hi_bo_z->SetLineColor(8); hi_bo_z->GetXaxis()->SetTitle("Z in BONUS, mm"); hi_diff_z->SetLineWidth(2); hi_diff_z->SetLineColor(kBlue); hi_diff_z->GetXaxis()->SetTitle("Z(CLAS)-Z(BoNuS), mm"); hi_cutdiff_z->SetLineWidth(2); hi_cutdiff_z->SetLineColor(kRed); hi_cutdiff_z->GetXaxis()->SetTitle("Z(CLAS)-Z(BoNuS), mm"); hi_clasdiff_z->SetLineWidth(2); hi_clasdiff_z->SetLineColor(kGreen); hi_clasdiff_z->GetXaxis()->SetTitle("Z(CLAS)-Z(avg), cm"); hi_r->SetLineWidth(2); hi_r->SetLineColor(kBlue); hi_r->GetXaxis()->SetTitle("mm"); hi_x2->SetLineWidth(2); hi_x2->SetLineColor(kBlue); hi_x2->GetXaxis()->SetTitle("#Chi^2"); hi_npt->SetLineWidth(2); hi_npt->SetLineColor(kBlue); hi_npt->GetXaxis()->SetTitle("number"); hi_edist->SetLineWidth(2); hi_edist->SetLineColor(kRed); hi_edist->GetXaxis()->SetTitle("edist, mm"); hi_sdist->SetLineWidth(2); hi_sdist->SetLineColor(kRed); hi_sdist->GetXaxis()->SetTitle("sdist, mm"); hi_p_corr->SetLineWidth(2); hi_p_corr->SetLineColor(kOrange); hi_p_corr->GetXaxis()->SetTitle("p, MeV"); hi_p_pre->SetLineWidth(1); hi_p_pre->SetLineColor(kBlue); hi_p_pre->GetXaxis()->SetTitle("p_predicted (blue) measured(red), MeV"); hi_pt_m->SetLineWidth(1); hi_pt_m->SetLineColor(kBlue); hi_pt_m->GetXaxis()->SetTitle("predicted p_{perp} (blue) measured(red), MeV"); hi_p_rtpc->SetLineWidth(1); hi_p_rtpc->SetLineColor(kRed); hi_p_rtpc->GetXaxis()->SetTitle("p, MeV"); hi_pt_rtpc->SetLineWidth(1); hi_pt_rtpc->SetLineColor(kRed); hi_pt_rtpc->GetXaxis()->SetTitle("p, MeV"); hi_p_raw->SetLineWidth(2); hi_p_raw->SetLineColor(kGreen); hi_p_raw->GetXaxis()->SetTitle("p, MeV"); hi_p_pro->SetLineWidth(2); hi_p_pro->SetLineColor(kRed); hi_p_others->SetLineWidth(2); hi_p_others->SetLineColor(kBlue); hi_betheprof->SetLineWidth(2); hi_betheprof->SetLineColor(kBlue); hi_betheprof->GetXaxis()->SetTitle("p_bonus, MeV/c"); hi_diff_dqdx->SetLineWidth(2); hi_diff_dqdx->SetLineColor(kBlue); hi_diff_dqdx->GetXaxis()->SetTitle("measured - expected dqdx"); hi_rat_dqdx->SetLineWidth(2); hi_rat_dqdx->SetLineColor(kBlue); hi_rat_dqdx->GetXaxis()->SetTitle("measured/expected dqdx"); hi_edistL->SetLineWidth(2); hi_edistL->SetLineColor(kRed); hi_edistL->GetXaxis()->SetTitle("edist, mm"); hi_sdistL->SetLineWidth(2); hi_sdistL->SetLineColor(kRed); hi_sdistL->GetXaxis()->SetTitle("sdist, mm"); hi_dzL->SetLineWidth(2); hi_dzL->SetLineColor(kRed); hi_dzL->GetXaxis()->SetTitle("Z(CLAS)-Z(BoNuS), mm"); hi_dthetaL->SetLineWidth(2); hi_dthetaL->SetLineColor(kRed); hi_dthetaL->GetXaxis()->SetTitle("Theta_CLAS - Theta_RTPC, deg"); hi_dphiL->SetLineWidth(3); hi_dphiL->SetLineColor(kRed); hi_dphiL->GetXaxis()->SetTitle("Phi_CLAS - Phi_BONUS,deg"); hi_edistR->SetLineWidth(2); hi_edistR->SetLineColor(kBlue); hi_edistR->GetXaxis()->SetTitle("edist, mm"); hi_sdistR->SetLineWidth(2); hi_sdistR->SetLineColor(kBlue); hi_sdistR->GetXaxis()->SetTitle("sdist, mm"); hi_dzR->SetLineWidth(2); hi_dzR->SetLineColor(kBlue); hi_dzR->GetXaxis()->SetTitle("Z(CLAS)-Z(BoNuS), mm"); hi_dthetaR->SetLineWidth(2); hi_dthetaR->SetLineColor(kBlue); hi_dthetaR->GetXaxis()->SetTitle("Theta_CLAS - Theta_RTPC, deg"); hi_dphiR->SetLineWidth(3); hi_dphiR->SetLineColor(kBlue); hi_dphiR->GetXaxis()->SetTitle("Phi_CLAS - Phi_BONUS,deg"); ///////////////////2D histos///////////////////////////////////// h2_ed_sd->SetXTitle("sdist (mm)"); h2_ed_sd->SetYTitle("edist (mm)"); h2_ed_sd_dz->SetXTitle("sdist (mm)"); h2_ed_sd_dz->SetYTitle("edist (mm)"); h2_phi->SetXTitle("clas_phi in BONUS coords(rad)"); h2_phi->SetYTitle("bonus_phi(rad)"); h2_phiVdth->SetXTitle("theta_CLAS - theta_BONUS (deg)"); h2_phiVdth->SetYTitle("phi(deg)"); h2_bethe->SetXTitle("mom from curvature (MeV)"); h2_bethe->SetYTitle("dQdx (ADC counts/mm)"); h2_bethepro->SetXTitle("mom from curvature (MeV)"); h2_bethepro->SetYTitle("dQdx (ADC counts/mm)"); h2_avgadc_st->SetXTitle("mom from curvature (MeV)"); h2_avgadc_st->SetYTitle("avg ADC * sin(#theta) / hit"); h2_avgadc->SetXTitle("mom from curvature (MeV)"); h2_avgadc->SetYTitle("avg ADC / hit"); h2_nhits->SetXTitle("mom from curvature (MeV)"); h2_nhits->SetYTitle("number of hits"); h2_mom->SetXTitle("p_CLAS, MeV"); //h2_mom->SetMarkerStyle(kFullCircle); //h2_mom->SetMarkerColor(kRed); h2_mom->SetYTitle("p_BONUS, MeV"); h2_sd->SetXTitle("z_track, mm"); h2_sd->SetYTitle("phi_track, rad"); h2_chiVr->SetXTitle("inverse fit radius, 1/mm"); h2_chiVr->SetYTitle("chi^2 - kind of"); h2_dphiVmom->SetXTitle("CLAS measured proton momentum in MeV"); h2_dphiVmom->SetYTitle("dphi(rad)"); h2_dphiVdz->SetXTitle("proton z-vertex(mm)"); h2_dphiVdz->SetYTitle("dphi(deg)"); h2_dphiVpphi->SetXTitle("phi proton(deg)"); h2_dphiVpphi->SetYTitle("dphi(deg)"); h2_sdVz->SetXTitle("z(mm)"); h2_sdVz->SetYTitle("sdist(mm)"); h2_dpVphi->SetXTitle("phi(deg)"); h2_dpVphi->SetYTitle("dp/p"); //h2_dpVphi->SetMarkerStyle(kFullCircle); h2_dpVth->SetXTitle("theta(deg)"); h2_dpVth->SetYTitle("dp/p"); //h2_dpVth->SetMarkerStyle(kFullCircle); h2_dthVdphi->SetXTitle("dphi(deg)"); h2_dthVdphi->SetYTitle("dth(deg)"); h2_dpVchi2->SetXTitle("chi^2"); h2_dpVchi2->SetYTitle("dp(MeV)"); h2_dpVp->SetXTitle("measured momentum (MeV)"); h2_dpVp->SetYTitle("(predicted - measured momentum) / measured"); h2_bethewstar->SetXTitle("W* (GeV)"); h2_bethewstar->SetYTitle("p_bonus (MeV)"); h2_ps_w_sig->SetXTitle("W (GeV)"); h2_ps_w_sig->SetYTitle("p_{spec} (MeV/c)"); h2_theta->SetXTitle("clas_theta (deg)"); h2_theta->SetYTitle("bonus_theta(deg)"); h2_zbVzc->SetXTitle("z CLAS (mm)"); h2_zbVzc->SetYTitle("z BoNuS (mm)"); h2_zVw->SetXTitle("corrected W"); h2_zVw->SetYTitle("z (mm)"); //h2_zVw->SetMaximum(100); h2_wVphib4->SetXTitle("phi (deg)"); h2_wVphib4->SetYTitle("W"); h2_wVphiaf->SetXTitle("phi (deg)"); h2_wVphiaf->SetYTitle("W"); h2_dpVdcb4->SetXTitle("phiDC (deg)"); h2_dpVdcb4->SetYTitle("(MEASURED p - CALC p)/p"); h2_dpVdcaf->SetXTitle("phiDC (deg)"); h2_dpVdcaf->SetYTitle("(MEASURED p - CALC p)/p"); h2_wVQ2->SetXTitle("Q^{2}"); h2_wVQ2->SetXTitle("Q^{2}"); h2_wVQ2->SetYTitle("W (GeV)"); h2_wVwstar->SetXTitle("w*(GeV)"); h2_wVwstar->SetYTitle("w (GeV)"); h2_psVcostpq->GetYaxis()->SetTitleOffset(1.35); h2_psVcostpq->SetXTitle("cos(#theta_{pq})"); h2_psVcostpq->SetYTitle("p_{spec} (MeV/c)"); h2_psValpha->GetYaxis()->SetTitleOffset(1.35); h2_psValpha->SetXTitle("#alpha_{s}"); h2_psValpha->SetYTitle("p_{spec} (MeV/c)"); h2_psVcostpqgt->GetYaxis()->SetTitleOffset(1.35); h2_psVcostpqgt->SetXTitle("cos(#theta_{pq})"); h2_psVcostpqgt->SetYTitle("p_{spec} (MeV/c)"); h2_psValphagt->GetYaxis()->SetTitleOffset(1.35); h2_psValphagt->SetXTitle("#alpha_{s}"); h2_psValphagt->SetYTitle("p_{spec} (MeV/c)"); h2_deltaVps->SetXTitle("p_{spec} (MeV/c)"); h2_deltaVps->SetYTitle("delta (W - W*)"); h2_ratVx->SetXTitle("x_{bj}"); h2_ratVx->SetYTitle("#sigma_{n}/#sigma_{p}"); //prof_ratVx->SetXTitle("x_{bj}"); //prof_ratVx->SetYTitle("#sigma_{n}/#sigma_{p}"); }