- if(trdQuality) {
- // Track TRD quality vs (eta,phi) and vs TRD n-clusters
- TH3D* trdQuality3D = (TH3D*)trdQuality->Project(0, trdQuality->GetVar(fgkVarNames[kTrackEtaTRD]),
- trdQuality->GetVar(fgkVarNames[kTrackPhiTRD]),
- trdQuality->GetVar(fgkVarNames[kTrackTrdQuality]));
- trdQuality3D->GetZaxis()->SetRangeUser(0.1,100.0);
- trdQuality3D->SetName("trdQuality3D");
- TProfile2D* prof2DQuality = trdQuality3D->Project3DProfile("yx");
- prof2DQuality->SetName("prof2DQuality");
- trdQuality->SetRangeUser(trdQuality->GetVar(fgkVarNames[kTrackCharge]), -1.5,-0.5);
- TH2D* trdQualityPneg = (TH2D*)trdQuality->Project(0, trdQuality->GetVar(fgkVarNames[kTrackP]),
- trdQuality->GetVar(fgkVarNames[kTrackTrdQuality]));
- trdQualityPneg->SetName("trdQualityPneg");
- trdQuality->SetRangeUser(trdQuality->GetVar(fgkVarNames[kTrackCharge]), +0.5,1.5);
- TH2D* trdQualityPpos = (TH2D*)trdQuality->Project(0, trdQuality->GetVar(fgkVarNames[kTrackP]),
- trdQuality->GetVar(fgkVarNames[kTrackTrdQuality]));
- trdQualityPpos->SetName("trdQualityPpos");
- trdQualityPneg->GetYaxis()->SetRangeUser(0.1,100.0);
- TProfile* trdQualityPnegprof = trdQualityPneg->ProfileX();
- trdQualityPnegprof->SetName("trdQualityPnegprof");
- trdQualityPpos->GetYaxis()->SetRangeUser(0.1,100.0);
- TProfile* trdQualityPposprof = trdQualityPpos->ProfileX();
- trdQualityPposprof->SetName("trdQualityPposprof");
- pad = ((TVirtualPad*)l->At(3)); pad->cd();
- pad->SetLeftMargin(0.15); pad->SetRightMargin(0.1);
- pad->SetTopMargin(0.06); pad->SetBottomMargin(0.15);
- prof2DQuality->SetStats(kFALSE);
- prof2DQuality->SetTitle("");
- SetStyle(prof2DQuality->GetXaxis(), "#eta", 0.06, 1.0, kTRUE, 0.06);
- SetStyle(prof2DQuality->GetYaxis(), "#varphi (rad.)", 0.06, 1.0, kTRUE, 0.06);
- prof2DQuality->SetMaximum(1.6); prof2DQuality->SetMinimum(0.4);
- prof2DQuality->Draw("colz");
- lat->DrawLatex(0.2, 0.95, "TRD quality");
- pad = ((TVirtualPad*)l->At(6)); pad->cd();
- pad->SetLeftMargin(0.15); pad->SetRightMargin(0.01);
- pad->SetTopMargin(0.06); pad->SetBottomMargin(0.15);
- trdQualityPnegprof->SetStats(kFALSE);
- trdQualityPnegprof->SetTitle("");
- SetStyle(trdQualityPnegprof->GetXaxis(), "P (GeV/c)", 0.06, 1.0, kTRUE, 0.06);
- SetStyle(trdQualityPnegprof->GetYaxis(), "<TRD quality>", 0.06, 1.0, kTRUE, 0.06);
- SetStyle(trdQualityPnegprof, 1, 2, 2, 20, 2, 1);
- SetStyle(trdQualityPposprof, 1, 4, 2, 20, 4, 1);
- trdQualityPnegprof->GetYaxis()->SetRangeUser(0.5, 1.39);
- trdQualityPnegprof->Draw();
- trdQualityPposprof->Draw("same");
- lat->DrawLatex(0.2, 0.95, "TRD quality");
- }
-