+ case 4: // plot a 3x3 canvas with tracking related histograms
+ gPad->SetTopMargin(0.05); gPad->SetBottomMargin(0.001);
+ gPad->SetLeftMargin(0.001); gPad->SetRightMargin(0.001);
+ gPad->Divide(3,3,0.,0.);
+ l=gPad->GetListOfPrimitives();
+ // eta-phi distr. for positive TPC tracks
+ pad = ((TVirtualPad*)l->At(0)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ h3F = dynamic_cast<TH3F*>(fHistos->At(kTPCRefTracksPos));
+ h2FtpcP = (TH2F*)Proj3D((TH3F*)h3F, 0x0, 1, h3F->GetZaxis()->GetNbins(), nada)->Clone();
+ h2FtpcP->SetStats(kFALSE);
+ h2FtpcP->GetXaxis()->SetTitle("#eta");
+ h2FtpcP->GetXaxis()->CenterTitle();
+ h2FtpcP->GetXaxis()->SetTitleSize(0.07);
+ h2FtpcP->GetXaxis()->SetTitleOffset(0.8);
+ h2FtpcP->GetXaxis()->SetLabelSize(0.05);
+ h2FtpcP->GetYaxis()->SetTitle("detector #varphi");
+ h2FtpcP->GetYaxis()->CenterTitle();
+ h2FtpcP->GetYaxis()->SetTitleSize(0.07);
+ h2FtpcP->GetYaxis()->SetTitleOffset(0.8);
+ h2FtpcP->GetYaxis()->SetLabelSize(0.05);
+ h2FtpcP->SetTitle("");
+ h2FtpcP->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TPC positive ref. tracks");
+ //-----------------
+ // eta-phi distr. for negative TPC tracks
+ pad = ((TVirtualPad*)l->At(1)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ h3F = dynamic_cast<TH3F*>(fHistos->At(kTPCRefTracksNeg));
+ h2FtpcN = (TH2F*)Proj3D((TH3F*)h3F, 0x0, 1, h3F->GetZaxis()->GetNbins(), nada)->Clone();
+ h2FtpcN->SetStats(kFALSE);
+ h2FtpcN->GetXaxis()->SetTitle("#eta");
+ h2FtpcN->GetXaxis()->CenterTitle();
+ h2FtpcN->GetXaxis()->SetTitleSize(0.07);
+ h2FtpcN->GetXaxis()->SetTitleOffset(0.8);
+ h2FtpcN->GetXaxis()->SetLabelSize(0.05);
+ h2FtpcN->GetYaxis()->SetTitle("detector #varphi");
+ h2FtpcN->GetYaxis()->CenterTitle();
+ h2FtpcN->GetYaxis()->SetTitleSize(0.07);
+ h2FtpcN->GetYaxis()->SetTitleOffset(0.8);
+ h2FtpcN->GetYaxis()->SetLabelSize(0.05);
+ h2FtpcN->SetTitle("");
+ h2FtpcN->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TPC negative ref. tracks");
+ // eta-phi distr. for positive TRD tracks
+ pad = ((TVirtualPad*)l->At(3)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ h3F = dynamic_cast<TH3F*>(fHistos->At(kTRDRefTracksPos));
+ h2FtrdP = (TH2F*)Proj3D((TH3F*)h3F, 0x0, 1, h3F->GetZaxis()->GetNbins(), nada)->Clone();
+ h2FtrdP->SetStats(kFALSE);
+ h2FtrdP->GetXaxis()->SetTitle("#eta");
+ h2FtrdP->GetXaxis()->CenterTitle();
+ h2FtrdP->GetXaxis()->SetTitleSize(0.07);
+ h2FtrdP->GetXaxis()->SetTitleOffset(0.8);
+ h2FtrdP->GetXaxis()->SetLabelSize(0.05);
+ h2FtrdP->GetYaxis()->SetTitle("detector #varphi");
+ h2FtrdP->GetYaxis()->CenterTitle();
+ h2FtrdP->GetYaxis()->SetTitleSize(0.07);
+ h2FtrdP->GetYaxis()->SetTitleOffset(0.8);
+ h2FtrdP->GetYaxis()->SetLabelSize(0.05);
+ h2FtrdP->SetMaximum(h2FtpcP->GetMaximum());
+ h2FtrdP->SetTitle("");
+ h2FtrdP->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TRD positive ref. tracks");
+ //-----------------
+ // eta-phi distr. for negative TRD tracks
+ pad = ((TVirtualPad*)l->At(4)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ h3F = dynamic_cast<TH3F*>(fHistos->At(kTRDRefTracksNeg));
+ h2FtrdN = (TH2F*)Proj3D((TH3F*)h3F, 0x0, 1, h3F->GetZaxis()->GetNbins(), nada)->Clone();
+ h2FtrdN->SetStats(kFALSE);
+ h2FtrdN->GetXaxis()->SetTitle("#eta");
+ h2FtrdN->GetXaxis()->CenterTitle();
+ h2FtrdN->GetXaxis()->SetTitleSize(0.07);
+ h2FtrdN->GetXaxis()->SetTitleOffset(0.8);
+ h2FtrdN->GetXaxis()->SetLabelSize(0.05);
+ h2FtrdN->GetYaxis()->SetTitle("detector #varphi");
+ h2FtrdN->GetYaxis()->CenterTitle();
+ h2FtrdN->GetYaxis()->SetTitleSize(0.07);
+ h2FtrdN->GetYaxis()->SetTitleOffset(0.8);
+ h2FtrdN->GetYaxis()->SetLabelSize(0.05);
+ h2FtrdN->SetMaximum(h2FtpcN->GetMaximum());
+ h2FtrdN->SetTitle("");
+ h2FtrdN->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TRD negative ref. tracks");
+ // eta-phi efficiency for positive TRD tracks
+ pad = ((TVirtualPad*)l->At(6)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ h2Feff = (TH2F*)h2FtrdP->Clone();
+ h2Feff->Reset();
+ h2Feff->SetStats(kFALSE);
+ h2Feff->Divide(h2FtrdP, h2FtpcP);
+ h2Feff->GetXaxis()->SetTitle("#eta");
+ h2Feff->GetXaxis()->CenterTitle();
+ h2Feff->GetXaxis()->SetTitleSize(0.07);
+ h2Feff->GetXaxis()->SetTitleOffset(0.8);
+ h2Feff->GetXaxis()->SetLabelSize(0.05);
+ h2Feff->GetYaxis()->SetTitle("detector #varphi");
+ h2Feff->GetYaxis()->CenterTitle();
+ h2Feff->GetYaxis()->SetTitleSize(0.07);
+ h2Feff->GetYaxis()->SetTitleOffset(0.8);
+ h2Feff->GetYaxis()->SetLabelSize(0.05);
+ h2Feff->SetMaximum(1.0);
+ h2Feff->SetTitle("");
+ h2Feff->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "Efficiency positive tracks");
+ // eta-phi efficiency for negative TRD tracks
+ pad = ((TVirtualPad*)l->At(7)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ h2Feff = (TH2F*)h2FtrdN->Clone();
+ h2Feff->Reset();
+ h2Feff->SetStats(kFALSE);
+ h2Feff->Divide(h2FtrdN, h2FtpcN);
+ h2Feff->GetXaxis()->SetTitle("#eta");
+ h2Feff->GetXaxis()->CenterTitle();
+ h2Feff->GetXaxis()->SetTitleSize(0.07);
+ h2Feff->GetXaxis()->SetTitleOffset(0.8);
+ h2Feff->GetXaxis()->SetLabelSize(0.05);
+ h2Feff->GetYaxis()->SetTitle("detector #varphi");
+ h2Feff->GetYaxis()->CenterTitle();
+ h2Feff->GetYaxis()->SetTitleSize(0.07);
+ h2Feff->GetYaxis()->SetTitleOffset(0.8);
+ h2Feff->GetYaxis()->SetLabelSize(0.05);
+ h2Feff->SetMaximum(1.0);
+ h2Feff->SetTitle("");
+ h2Feff->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "Efficiency negative tracks");
+
+ // <ntracklets> vs (phi,eta)
+ pad = ((TVirtualPad*)l->At(2)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ hProf2D = dynamic_cast<TProfile2D*>(fHistos->At(kTRDEtaPhiAvNtrkl));
+ hProf2D->SetStats(kFALSE);
+ hProf2D->SetTitle("");
+ hProf2D->GetXaxis()->SetTitle("#eta");
+ hProf2D->GetXaxis()->SetTitleOffset(0.8);
+ hProf2D->GetXaxis()->SetTitleSize(0.07);
+ hProf2D->GetXaxis()->CenterTitle();
+ hProf2D->GetXaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->SetTitle("detector #varphi");
+ hProf2D->GetYaxis()->SetTitleOffset(0.8);
+ hProf2D->GetYaxis()->SetTitleSize(0.07);
+ hProf2D->GetYaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->CenterTitle();
+ hProf2D->SetMinimum(0.);
+ hProf2D->SetMaximum(6.);
+ hProf2D->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TRD <N_{tracklets}>");
+ // TPC-TRD matching efficiency vs pt
+ pad = ((TVirtualPad*)l->At(5)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.02);
+ pad->SetTopMargin(0.02); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ hFeffP = EfficiencyTRD(1);
+ hFeffN = EfficiencyTRD(-1);
+ h2F=new TH2F("rangeEffPt", "",10,0.,10.,10,0.,1.1);
+ h2F->SetStats(kFALSE);
+ h2F->GetXaxis()->SetTitle("p_{T} [GeV/c]");
+ h2F->GetXaxis()->SetTitleOffset(0.8);
+ h2F->GetXaxis()->SetTitleSize(0.07);
+ h2F->GetXaxis()->CenterTitle();
+ h2F->GetXaxis()->SetLabelSize(0.05);
+ h2F->GetYaxis()->SetTitle("TRD-TPC matching efficiency");
+ h2F->GetYaxis()->SetTitleOffset(0.8);
+ h2F->GetYaxis()->SetTitleSize(0.07);
+ h2F->GetYaxis()->SetLabelSize(0.05);
+ h2F->GetYaxis()->CenterTitle();
+ h2F->Draw();
+ line.SetLineStyle(2);
+ line.SetLineWidth(2);
+ line.DrawLine(h2F->GetXaxis()->GetXmin(), 0.7, h2F->GetXaxis()->GetXmax(), 0.7);
+ line.DrawLine(h2F->GetXaxis()->GetXmin(), 0.9, h2F->GetXaxis()->GetXmax(), 0.9);
+ hFeffP->SetMarkerStyle(20);
+ hFeffP->SetMarkerColor(2);
+ hFeffN->SetMarkerStyle(22);
+ hFeffN->SetMarkerColor(4);
+ hFeffP->Draw("same");
+ hFeffN->Draw("same");
+ leg=new TLegend(0.65, 0.2, 0.95, 0.4);
+ leg->SetFillColor(0);
+ leg->AddEntry(hFeffP, "positives", "p");
+ leg->AddEntry(hFeffN, "negatives", "p");
+ leg->Draw();
+ // create trending values for the TPC-TRD matching efficiency
+ // fit the efficiency histos with a constant in the range [1.0,1.5] GeV/c
+ fitFunc = new TF1("constantFunc","[0]",1.0,1.5);
+ hFeffP->Fit(fitFunc,"Q0","",1.0,1.5);
+ PutTrendValue("TrackingEffPos1GeV", fitFunc->GetParameter(0));
+ PutTrendValue("TrackingEffPos1GeVErr", fitFunc->GetParError(0));
+ hFeffN->Fit(fitFunc,"Q0","",1.0,1.5);
+ PutTrendValue("TrackingEffNeg1GeV", fitFunc->GetParameter(0));
+ PutTrendValue("TrackingEffNeg1GeVErr", fitFunc->GetParError(0));
+ // Nclusters per TRD track
+ pad = ((TVirtualPad*)l->At(8)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.12);
+ pad->SetTopMargin(0.02); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ pad->SetLogz();
+ h2F = dynamic_cast<TH2F*>(fHistos->At(kNClsTrackTRD));
+ h2F->SetStats(kFALSE);
+ h2F->SetTitle("");
+ h2F->GetXaxis()->SetTitle("p [GeV/c]");
+ h2F->GetXaxis()->SetTitleOffset(0.8);
+ h2F->GetXaxis()->SetTitleSize(0.07);
+ h2F->GetXaxis()->CenterTitle();
+ h2F->GetXaxis()->SetLabelSize(0.05);
+ h2F->GetYaxis()->SetTitle("#clusters per TRD track");
+ h2F->GetYaxis()->SetTitleOffset(0.8);
+ h2F->GetYaxis()->SetTitleSize(0.07);
+ h2F->GetYaxis()->CenterTitle();
+ h2F->GetYaxis()->SetLabelSize(0.05);
+ h2F->Draw("colz");
+ break;
+ case 5: // plot a 3x3 canvas with PID related histograms
+ gPad->SetTopMargin(0.05); gPad->SetBottomMargin(0.001);
+ gPad->SetLeftMargin(0.001); gPad->SetRightMargin(0.001);
+ gPad->Divide(3,3,0.,0.);
+ l=gPad->GetListOfPrimitives();
+ // eta-phi distr. for <Qtot> in layer 0
+ pad = ((TVirtualPad*)l->At(0)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ hProf2D = dynamic_cast<TProfile2D*>(fHistos->At(kTRDEtaPhiAvQtot+0));
+ hProf2D->SetStats(kFALSE);
+ hProf2D->SetTitle("");
+ hProf2D->GetXaxis()->SetTitle("#eta");
+ hProf2D->GetXaxis()->SetTitleOffset(0.8);
+ hProf2D->GetXaxis()->SetTitleSize(0.07);
+ hProf2D->GetXaxis()->CenterTitle();
+ hProf2D->GetXaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->SetTitle("detector #varphi");
+ hProf2D->GetYaxis()->SetTitleOffset(0.8);
+ hProf2D->GetYaxis()->SetTitleSize(0.07);
+ hProf2D->GetYaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->CenterTitle();
+ hProf2D->SetMinimum(0.);
+ hProf2D->SetMaximum(25.);
+ hProf2D->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TRD <Q_{tot}> Layer 0");
+ // eta-phi distr. for <Qtot> in layer 1
+ pad = ((TVirtualPad*)l->At(3)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ hProf2D = dynamic_cast<TProfile2D*>(fHistos->At(kTRDEtaPhiAvQtot+1));
+ hProf2D->SetStats(kFALSE);
+ hProf2D->SetTitle("");
+ hProf2D->GetXaxis()->SetTitle("#eta");
+ hProf2D->GetXaxis()->SetTitleOffset(0.8);
+ hProf2D->GetXaxis()->SetTitleSize(0.07);
+ hProf2D->GetXaxis()->CenterTitle();
+ hProf2D->GetXaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->SetTitle("detector #varphi");
+ hProf2D->GetYaxis()->SetTitleOffset(0.8);
+ hProf2D->GetYaxis()->SetTitleSize(0.07);
+ hProf2D->GetYaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->CenterTitle();
+ hProf2D->SetMinimum(0.);
+ hProf2D->SetMaximum(25.);
+ hProf2D->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TRD <Q_{tot}> Layer 1");
+ // eta-phi distr. for <Qtot> in layer 2
+ pad = ((TVirtualPad*)l->At(6)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ hProf2D = dynamic_cast<TProfile2D*>(fHistos->At(kTRDEtaPhiAvQtot+2));
+ hProf2D->SetStats(kFALSE);
+ hProf2D->SetTitle("");
+ hProf2D->GetXaxis()->SetTitle("#eta");
+ hProf2D->GetXaxis()->SetTitleOffset(0.8);
+ hProf2D->GetXaxis()->SetTitleSize(0.07);
+ hProf2D->GetXaxis()->CenterTitle();
+ hProf2D->GetXaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->SetTitle("detector #varphi");
+ hProf2D->GetYaxis()->SetTitleOffset(0.8);
+ hProf2D->GetYaxis()->SetTitleSize(0.07);
+ hProf2D->GetYaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->CenterTitle();
+ hProf2D->SetMinimum(0.);
+ hProf2D->SetMaximum(25.);
+ hProf2D->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TRD <Q_{tot}> Layer 2");
+ // eta-phi distr. for <Qtot> in layer 3
+ pad = ((TVirtualPad*)l->At(1)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ hProf2D = dynamic_cast<TProfile2D*>(fHistos->At(kTRDEtaPhiAvQtot+3));
+ hProf2D->SetStats(kFALSE);
+ hProf2D->SetTitle("");
+ hProf2D->GetXaxis()->SetTitle("#eta");
+ hProf2D->GetXaxis()->SetTitleOffset(0.8);
+ hProf2D->GetXaxis()->SetTitleSize(0.07);
+ hProf2D->GetXaxis()->CenterTitle();
+ hProf2D->GetXaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->SetTitle("detector #varphi");
+ hProf2D->GetYaxis()->SetTitleOffset(0.8);
+ hProf2D->GetYaxis()->SetTitleSize(0.07);
+ hProf2D->GetYaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->CenterTitle();
+ hProf2D->SetMinimum(0.);
+ hProf2D->SetMaximum(25.);
+ hProf2D->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TRD <Q_{tot}> Layer 3");
+ // eta-phi distr. for <Qtot> in layer 4
+ pad = ((TVirtualPad*)l->At(4)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ hProf2D = dynamic_cast<TProfile2D*>(fHistos->At(kTRDEtaPhiAvQtot+4));
+ hProf2D->SetStats(kFALSE);
+ hProf2D->SetTitle("");
+ hProf2D->GetXaxis()->SetTitle("#eta");
+ hProf2D->GetXaxis()->SetTitleOffset(0.8);
+ hProf2D->GetXaxis()->SetTitleSize(0.07);
+ hProf2D->GetXaxis()->CenterTitle();
+ hProf2D->GetXaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->SetTitle("detector #varphi");
+ hProf2D->GetYaxis()->SetTitleOffset(0.8);
+ hProf2D->GetYaxis()->SetTitleSize(0.07);
+ hProf2D->GetYaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->CenterTitle();
+ hProf2D->SetMinimum(0.);
+ hProf2D->SetMaximum(25.);
+ hProf2D->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TRD <Q_{tot}> Layer 4");
+ // eta-phi distr. for <Qtot> in layer 5
+ pad = ((TVirtualPad*)l->At(7)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ hProf2D = dynamic_cast<TProfile2D*>(fHistos->At(kTRDEtaPhiAvQtot+5));
+ hProf2D->SetStats(kFALSE);
+ hProf2D->SetTitle("");
+ hProf2D->GetXaxis()->SetTitle("#eta");
+ hProf2D->GetXaxis()->SetTitleOffset(0.8);
+ hProf2D->GetXaxis()->SetTitleSize(0.07);
+ hProf2D->GetXaxis()->CenterTitle();
+ hProf2D->GetXaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->SetTitle("detector #varphi");
+ hProf2D->GetYaxis()->SetTitleOffset(0.8);
+ hProf2D->GetYaxis()->SetTitleSize(0.07);
+ hProf2D->GetYaxis()->SetLabelSize(0.05);
+ hProf2D->GetYaxis()->CenterTitle();
+ hProf2D->SetMinimum(0.);
+ hProf2D->SetMaximum(25.);
+ hProf2D->Draw("colz");
+ lat->DrawLatex(-0.9, 3.6, "TRD <Q_{tot}> Layer 5");
+ // PH versus slice number
+ pad = ((TVirtualPad*)l->At(2)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ h2F = dynamic_cast<TH2F*>(fHistos->At(kPHSlice));
+ hF = Proj2D((TH2F*)h2F);
+ h2F->SetStats(kFALSE);
+ h2F->SetTitle("");
+ h2F->GetXaxis()->SetTitle("slice");
+ h2F->GetXaxis()->SetTitleOffset(0.8);
+ h2F->GetXaxis()->SetTitleSize(0.07);
+ h2F->GetXaxis()->CenterTitle();
+ h2F->GetXaxis()->SetLabelSize(0.05);
+ h2F->GetYaxis()->SetTitle("PH");
+ h2F->GetYaxis()->SetTitleOffset(0.8);
+ h2F->GetYaxis()->SetTitleSize(0.07);
+ h2F->GetYaxis()->SetLabelSize(0.05);
+ h2F->GetYaxis()->CenterTitle();
+ h2F->Draw("colz");
+ hF->SetLineWidth(2);
+ hF->Draw("same");
+ // Qtot vs P
+ pad = ((TVirtualPad*)l->At(5)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ pad->SetLogz();
+ h2F = dynamic_cast<TH2F*>(fHistos->At(kQtotP));
+ h2F->SetStats(kFALSE);
+ h2F->SetTitle("");
+ h2F->GetXaxis()->SetTitle("P [GeV/c]");
+ h2F->GetXaxis()->SetTitleOffset(0.8);
+ h2F->GetXaxis()->SetTitleSize(0.07);
+ h2F->GetXaxis()->CenterTitle();
+ h2F->GetXaxis()->SetLabelSize(0.05);
+ h2F->GetYaxis()->SetRangeUser(0.0,100.0);
+ h2F->GetYaxis()->SetTitle("Q_{tot}");
+ h2F->GetYaxis()->SetTitleOffset(0.8);
+ h2F->GetYaxis()->SetTitleSize(0.07);
+ h2F->GetYaxis()->SetLabelSize(0.05);
+ h2F->GetYaxis()->CenterTitle();
+ h2F->Draw("colz");
+ // create trending value for the average Qtot at 1 GeV/c
+ hProf = h2F->ProfileX("profileQtot",1,h2F->GetYaxis()->FindBin(40.));
+ PutTrendValue("AvQtot1GeV", hProf->GetBinContent(hProf->GetXaxis()->FindBin(1.)));
+ PutTrendValue("AvQtot1GeVErr", hProf->GetBinError(hProf->GetXaxis()->FindBin(1.)));
+ // PH versus slice number for TPC pions and electrons
+ pad = ((TVirtualPad*)l->At(8)); pad->cd();
+ pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
+ pad->SetTopMargin(0.1); pad->SetBottomMargin(0.15);
+ pad->SetGridx(kFALSE); pad->SetGridy(kFALSE);
+ h2FtrdP = dynamic_cast<TH2F*>(fHistos->At(kPHSliceTPCpions));
+ h2FtrdN = dynamic_cast<TH2F*>(fHistos->At(kPHSliceTPCelectrons));
+ hFeffP = Proj2D((TH2F*)h2FtrdP);
+ hFeffN = Proj2D((TH2F*)h2FtrdN);
+ h2F = new TH2F("PHvsSlice","",10,h2FtrdN->GetXaxis()->GetXmin(),h2FtrdN->GetXaxis()->GetXmax(),
+ 10,h2FtrdN->GetYaxis()->GetXmin(),h2FtrdN->GetYaxis()->GetXmax());
+ h2F->SetStats(kFALSE);
+ h2F->SetTitle("");
+ h2F->GetXaxis()->SetTitle("slice");
+ h2F->GetXaxis()->SetTitleOffset(0.8);
+ h2F->GetXaxis()->SetTitleSize(0.07);
+ h2F->GetXaxis()->CenterTitle();
+ h2F->GetXaxis()->SetLabelSize(0.05);
+ h2F->GetYaxis()->SetTitle("PH");
+ h2F->GetYaxis()->SetTitleOffset(0.8);
+ h2F->GetYaxis()->SetTitleSize(0.07);
+ h2F->GetYaxis()->SetLabelSize(0.05);
+ h2F->GetYaxis()->CenterTitle();
+ h2F->Draw();
+ hFeffN->SetLineWidth(2);
+ hFeffN->SetLineColor(2);
+ hFeffP->SetLineWidth(2);
+ hFeffP->SetLineColor(4);
+ hFeffN->Draw("same");
+ hFeffP->Draw("same");
+ leg=new TLegend(0.65, 0.8, 0.95, 0.95);
+ leg->SetFillColor(0);
+ leg->AddEntry(hFeffP, "TPC pions", "l");
+ leg->AddEntry(hFeffN, "TPC electrons", "l");
+ leg->Draw();
+ break;