8 TDirectory *
dir = (TDirectory*)gDirectory->FindObjectAny(
"fcal");
11 TH2F* digOcc2D = (TH2F*)gDirectory->FindObjectAny(
"digOcc2D");
12 TH1I* digN = (TH1I*)gDirectory->FindObjectAny(
"digN");
13 TH2F* digPeakV2D = (TH2F*)gDirectory->FindObjectAny(
"digPeakV2D");
14 TH1I* digPeakV = (TH1I*)gDirectory->FindObjectAny(
"digPeakV");
16 double nEvents = ( digN ? digN->GetEntries() : 0 );
20 TCanvas *
c1 =
new TCanvas(
"c1",
"FCAL Monitor", 800, 800 );
27 TCanvas*
c1 = gPad->GetCanvas();
33 digN->SetFillColor( kRed );
40 TH2F* digOcc2DAvg = (TH2F*)digOcc2D->Clone(
"digOcc2DAvg" );
42 digOcc2DAvg->SetTitle(
"FCAL Pulse Occupancy per Event" );
44 for(
int x = 1;
x <= digOcc2DAvg->GetNbinsX(); ++
x ){
45 for(
int y = 1;
y <= digOcc2DAvg->GetNbinsY(); ++
y ){
47 digOcc2DAvg->SetBinContent(
x,
y, digOcc2DAvg->GetBinContent(
x,
y ) / nEvents );
51 digOcc2DAvg->SetStats( 0 );
53 digOcc2DAvg->Draw(
"colz" );
58 digPeakV->SetStats( 0 );
59 digPeakV->SetFillColor( kRed );
64 if( digPeakV2D && digOcc2D ){
66 TH2F* digPeakV2DAvg = (TH2F*)digPeakV2D->Clone(
"digPeakV2DAvg" );
67 digPeakV2DAvg->Divide( digOcc2D );
68 digPeakV2DAvg->SetMinimum( 0 );
69 digPeakV2DAvg->SetMaximum( 2 );
70 digPeakV2DAvg->SetStats( 0 );
72 digPeakV2DAvg->Draw(
"colz" );