void pl_over_CluArea() { TH3* tHisto3; TProfile2D* tProfile; TCanvas* tCanvas = new TCanvas("tCanvas", "", 0., 0., 1575., 750.); tCanvas->Divide(3, 2); tCanvas->cd(1); gPad->SetLeftMargin(0.15); gPad->SetRightMargin(0.20); gPad->SetTopMargin(0.1); gPad->SetBottomMargin(0.2); tHisto3 = dynamic_cast(gROOT->FindObjectAny("hSel2CluSizeArea")); tProfile = tHisto3->Project3DProfile("yx"); tProfile->SetTitle("Sel2 surface cluster size"); tProfile->GetXaxis()->SetTitle("local X [cm]"); tProfile->GetXaxis()->SetTitleOffset(1.2); tProfile->GetYaxis()->SetTitle("local Y [cm]"); tProfile->GetZaxis()->SetTitle("mean cluster size [cells]"); tProfile->GetZaxis()->SetTitleSize(0.08); tProfile->GetZaxis()->SetTitleOffset(0.9); tProfile->SetStats(0); tProfile->Draw("colz"); for(Int_t iCell = 1; iCell <= tProfile->GetNcells(); iCell++) { if(tProfile->GetBinEntries(iCell) < 100.) { tProfile->SetBinContent(iCell, 0.); tProfile->SetBinError(iCell, 0.); } } gPad->Modified(); gPad->Update(); tCanvas->cd(2); gPad->SetLeftMargin(0.15); gPad->SetRightMargin(0.20); gPad->SetTopMargin(0.1); gPad->SetBottomMargin(0.2); tHisto3 = dynamic_cast(gROOT->FindObjectAny("hDUTCluSizeArea")); tProfile = tHisto3->Project3DProfile("yx"); tProfile->SetTitle("DUT surface cluster size"); tProfile->GetXaxis()->SetTitle("local X [cm]"); tProfile->GetXaxis()->SetTitleOffset(1.2); tProfile->GetYaxis()->SetTitle("local Y [cm]"); tProfile->GetZaxis()->SetTitle("mean cluster size [cells]"); tProfile->GetZaxis()->SetTitleSize(0.08); tProfile->GetZaxis()->SetTitleOffset(0.9); tProfile->SetStats(0); tProfile->Draw("colz"); for(Int_t iCell = 1; iCell <= tProfile->GetNcells(); iCell++) { if(tProfile->GetBinEntries(iCell) < 100.) { tProfile->SetBinContent(iCell, 0.); tProfile->SetBinError(iCell, 0.); } } gPad->Modified(); gPad->Update(); tCanvas->cd(3); gPad->SetLeftMargin(0.15); gPad->SetRightMargin(0.20); gPad->SetTopMargin(0.1); gPad->SetBottomMargin(0.2); tHisto3 = dynamic_cast(gROOT->FindObjectAny("hMRefCluSizeArea")); tProfile = tHisto3->Project3DProfile("yx"); tProfile->SetTitle("MRef surface cluster size"); tProfile->GetXaxis()->SetTitle("local X [cm]"); tProfile->GetXaxis()->SetTitleOffset(1.2); tProfile->GetYaxis()->SetTitle("local Y [cm]"); tProfile->GetZaxis()->SetTitle("mean cluster size [cells]"); tProfile->GetZaxis()->SetTitleSize(0.08); tProfile->GetZaxis()->SetTitleOffset(0.9); tProfile->SetStats(0); tProfile->Draw("colz"); for(Int_t iCell = 1; iCell <= tProfile->GetNcells(); iCell++) { if(tProfile->GetBinEntries(iCell) < 100.) { tProfile->SetBinContent(iCell, 0.); tProfile->SetBinError(iCell, 0.); } } gPad->Modified(); gPad->Update(); tCanvas->cd(4); gPad->SetLeftMargin(0.15); gPad->SetRightMargin(0.20); gPad->SetTopMargin(0.1); gPad->SetBottomMargin(0.2); tHisto3 = dynamic_cast(gROOT->FindObjectAny("hSelHitTupleSel2SelCluSizeArea")); tProfile = tHisto3->Project3DProfile("yx"); tProfile->SetTitle("Sel2 selector surface CS"); tProfile->GetXaxis()->SetTitle("local X [cm]"); tProfile->GetXaxis()->SetTitleOffset(1.2); tProfile->GetYaxis()->SetTitle("local Y [cm]"); tProfile->GetZaxis()->SetTitle("mean cluster size [cells]"); tProfile->GetZaxis()->SetTitleSize(0.08); tProfile->GetZaxis()->SetTitleOffset(0.9); tProfile->SetStats(0); tProfile->Draw("colz"); for(Int_t iCell = 1; iCell <= tProfile->GetNcells(); iCell++) { if(tProfile->GetBinEntries(iCell) < 100.) { tProfile->SetBinContent(iCell, 0.); tProfile->SetBinError(iCell, 0.); } } gPad->Modified(); gPad->Update(); tCanvas->cd(5); gPad->SetLeftMargin(0.15); gPad->SetRightMargin(0.20); gPad->SetTopMargin(0.1); gPad->SetBottomMargin(0.2); tHisto3 = dynamic_cast(gROOT->FindObjectAny("hSelHitTupleDUTNNCluSizeArea")); tProfile = tHisto3->Project3DProfile("yx"); tProfile->SetTitle("DUT selector surface CS"); tProfile->GetXaxis()->SetTitle("local X [cm]"); tProfile->GetXaxis()->SetTitleOffset(1.2); tProfile->GetYaxis()->SetTitle("local Y [cm]"); tProfile->GetZaxis()->SetTitle("mean cluster size [cells]"); tProfile->GetZaxis()->SetTitleSize(0.08); tProfile->GetZaxis()->SetTitleOffset(0.9); tProfile->SetStats(0); tProfile->Draw("colz"); for(Int_t iCell = 1; iCell <= tProfile->GetNcells(); iCell++) { if(tProfile->GetBinEntries(iCell) < 100.) { tProfile->SetBinContent(iCell, 0.); tProfile->SetBinError(iCell, 0.); } } gPad->Modified(); gPad->Update(); tCanvas->cd(6); gPad->SetLeftMargin(0.15); gPad->SetRightMargin(0.20); gPad->SetTopMargin(0.1); gPad->SetBottomMargin(0.2); tHisto3 = dynamic_cast(gROOT->FindObjectAny("hSelHitTupleMRefSelCluSizeArea")); tProfile = tHisto3->Project3DProfile("yx"); tProfile->SetTitle("MRef selector surface CS"); tProfile->GetXaxis()->SetTitle("local X [cm]"); tProfile->GetXaxis()->SetTitleOffset(1.2); tProfile->GetYaxis()->SetTitle("local Y [cm]"); tProfile->GetZaxis()->SetTitle("mean cluster size [cells]"); tProfile->GetZaxis()->SetTitleSize(0.08); tProfile->GetZaxis()->SetTitleOffset(0.9); tProfile->SetStats(0); tProfile->Draw("colz"); for(Int_t iCell = 1; iCell <= tProfile->GetNcells(); iCell++) { if(tProfile->GetBinEntries(iCell) < 100.) { tProfile->SetBinContent(iCell, 0.); tProfile->SetBinError(iCell, 0.); } } gPad->Modified(); gPad->Update(); tCanvas->SaveAs("pl_over_CluArea.pdf"); }