void draw_qa(){ TFile *file = new TFile("reco.richqa.root"); TDirectory *current = gDirectory; current->cd("RichRingQaHist"); //rich->cd(); //gROOT->SetStyle("Plain"); // gStyle->SetOptStat(0000); gStyle->SetPalette(1,0); gStyle->SetHistLineWidth(3); gROOT->ForceStyle(); //gStyle->SetOptStat(0); //gStyle->SetOptFit(1); TCanvas *c1 = new TCanvas("CbmRichQa1","c1",1200,1000); c1->Divide(3,3); c1->cd(1); fh_FakeNofHits->SetLineColor(kRed); fh_TrueElNofHits->Draw(); fh_FakeNofHits->Draw("SAME"); c1->cd(2); fh_FakeDistance->SetLineColor(kRed); fh_TrueElDistance->Draw(); fh_FakeDistance->Draw("SAME"); c1->cd(3); fh_FakeAngle->SetLineColor(kRed); fh_TrueElAngle->Draw(); fh_FakeAngle->Draw("SAME"); c1->cd(4); fh_FakeNofHitsOnRing->SetLineColor(kRed); fh_TrueElNofHitsOnRing->Draw(); fh_FakeNofHitsOnRing->Draw("SAME"); c1->cd(5); fh_FakeRadPos->SetLineColor(kRed); fh_TrueElRadPos->Draw(); fh_FakeRadPos->Draw("SAME"); c1->cd(6); fh_FakeChi2->SetLineColor(kRed); fh_TrueElChi2->Draw(); fh_FakeChi2->Draw("SAME"); c1->cd(7); fh_FakeRadius->SetLineColor(kRed); fh_TrueElRadius->Draw(); fh_FakeRadius->Draw("SAME"); TCanvas *c2 = new TCanvas("CbmRichQa2","c2",1200,1000); fh_TrueFoundElRingsProjHitCutMom->Sumw2(); fh_MCElRingsProjHitCutMom->Sumw2(); TH1D* th1 = new TH1D("th1","Efficiency electrons vs momentum",40,0,10); th1->GetXaxis()->SetTitle("momentum, GeV"); th1->GetYaxis()->SetTitle("efficiency"); th1->Divide(fh_TrueFoundElRingsProjHitCutMom,fh_MCElRingsProjHitCutMom,100); th1->Draw(); TCanvas *c3 = new TCanvas("CbmRichQa3","c3",1200,1000); fh_TrueFoundElRingsProjHitCutRadPos->Sumw2(); fh_MCElRingsProjHitCutRadPos->Sumw2(); TH1D* th2 = new TH1D("th1","Efficiency electrons vs RadialPosition",40,0,150); th2->GetXaxis()->SetTitle("radial position, cm"); th2->GetYaxis()->SetTitle("efficiency"); th2->Divide(fh_TrueFoundElRingsProjHitCutRadPos,fh_MCElRingsProjHitCutRadPos,100); th2->Draw(); }