{ gSystem->Load("libHAFT.so"); Char_t fname[]="HadesAcceptanceFilter_NOV02_v1.acc"; setfilename_(fname,sizeof(fname)); Int_t pid = 2; // select a particle ID Float_t p0 = 1.010; // select cut planes sitting on the grid points Float_t th0 = 45.0; // (e.g. 50 p bins, 45 theta bins, 15 phi bins centered Float_t ph0 = 31.0; // at 25,75,125,... & 1,3,5,7,9,... & 2,6,10,14,18,...) Int_t m1=0, m2=-4, m3=3, m4=4; Int_t xdim, ydim, zdim; getdimensions_(&pid, &xdim, &ydim, &zdim); Float_t plo, phi, dp, thlo, thhi, dth, phlo, phhi, dph; getlimits_(&pid, &plo, &phi, &dp, &thlo, &thhi, &dth, &phlo, &phhi, &dph); Float_t pgrid[xdim], thgrid[ydim], phgrid[zdim] ; for (Int_t i=0; iSetTitle("Momentum"); p1.SetLineColor(1); p1.SetLineWidth(2); p1->GetXaxis()->SetLabelSize(0.06); p1->GetXaxis()->SetRangeUser(plo,phi); p1->GetYaxis()->SetLabelSize(0.06); TGraph p2(np,p,acc2p); p2.SetLineColor(4); p2.SetLineWidth(2); TGraph p3(np,p,acc3p); p3.SetLineColor(3); p3.SetLineWidth(2); TGraph p4(np,p,acc4p); p4.SetLineColor(2); p4.SetLineWidth(2); c1.cd(1); gPad->SetGrid(); p1->Draw("al"); p2->Draw("l"); p3->Draw("l"); p4->Draw("l"); TGraph pg(xdim,pgrid,accgp); pg->Draw("*"); TGraph th1(nth,th,acc1th); // theta cut th1->SetTitle("Theta"); th1.SetLineColor(1); th1.SetLineWidth(2); th1->GetXaxis()->SetLabelSize(0.06); th1->GetXaxis()->SetRangeUser(thlo,thhi); th1->GetYaxis()->SetLabelSize(0.06); TGraph th2(nth,th,acc2th); th2.SetLineColor(4); th2.SetLineWidth(2); TGraph th3(nth,th,acc3th); th3.SetLineColor(3); th3.SetLineWidth(2); TGraph th4(nth,th,acc4th); th4.SetLineColor(2); th4.SetLineWidth(2); c1.cd(2); gPad->SetGrid(); th1->Draw("al"); th2->Draw("l"); th3->Draw("l"); th4->Draw("l"); TGraph thg(ydim,thgrid,accgth); thg->Draw("*"); TGraph ph1(nph,ph,acc1ph); // phi cut ph1->SetTitle("Phi"); ph1.SetLineColor(1); ph1.SetLineWidth(2); ph1->GetXaxis()->SetLabelSize(0.06); ph1->GetXaxis()->SetRangeUser(phlo,phhi); ph1->GetYaxis()->SetLabelSize(0.06); TGraph ph2(nph,ph,acc2ph); ph2.SetLineColor(4); ph2.SetLineWidth(2); TGraph ph3(nph,ph,acc3ph); ph3.SetLineColor(3); ph3.SetLineWidth(2); TGraph ph4(nph,ph,acc4ph); ph4.SetLineColor(2); ph4.SetLineWidth(2); c1.cd(3); gPad->SetGrid(); ph1->Draw("al"); ph2->Draw("l"); ph3->Draw("l"); ph4->Draw("l"); TGraph phg(zdim,phgrid,accgph); phg->Draw("*"); }