9 TDirectory *
dir = (TDirectory*)gDirectory->FindObjectAny(
"fcal");
12 TH2F* digOcc2D = (TH2F*)gDirectory->FindObjectAny(
"digOcc2D");
13 TH2F* digPed2D = (TH2F*)gDirectory->FindObjectAny(
"digPed2D");
14 TH2F* digPedSq2D = (TH2F*)gDirectory->FindObjectAny(
"digPedSq2D");
15 TH1I* digPed = (TH1I*)gDirectory->FindObjectAny(
"digPed");
16 TProfile* digPedChan = (TProfile*)gDirectory->FindObjectAny(
"digPedChan");
20 TCanvas *
c1 =
new TCanvas(
"c1",
"FCAL Monitor", 800, 800 );
27 TCanvas*
c1 = gPad->GetCanvas();
32 digPed->SetStats( 0 );
33 digPed->SetFillColor( kRed );
40 digPedChan->SetStats( 0 );
45 if( digPed2D && digOcc2D && digPed ){
47 TH2F* digPed2DAvg = (TH2F*)digPed2D->Clone(
"digPed2DAvg" );
48 digPed2DAvg->Divide( digOcc2D );
49 double avgPed = digPed->GetMean();
50 digPed2DAvg->SetTitle(
"FCAL Pedestal - Average Pedestal" );
52 for(
int x = 1;
x <= digPed2DAvg->GetNbinsX(); ++
x ){
53 for(
int y = 1;
y <= digPed2DAvg->GetNbinsY(); ++
y ){
55 digPed2DAvg->SetBinContent(
x,
y, digPed2DAvg->GetBinContent(
x,
y ) - avgPed );
59 digPed2DAvg->SetStats( 0 );
60 digPed2DAvg->SetMinimum( -0.2*avgPed );
61 digPed2DAvg->SetMaximum( 0.2*avgPed );
63 digPed2DAvg->Draw(
"colz" );
66 if( digPed2D && digPedSq2D && digOcc2D && digPed ){
68 TH2F* digPed2DAvg = (TH2F*)digPed2D->Clone(
"digPed2DAvg" );
69 TH2F* digPed2DRMS = (TH2F*)digPed2D->Clone(
"digPed2DRMS" );
70 TH2F* digPedSq2DAvg = (TH2F*)digPedSq2D->Clone(
"digPedSq2DAvg" );
71 digPed2DAvg->Divide( digOcc2D );
72 digPedSq2DAvg->Divide( digOcc2D );
74 digPed2DRMS->SetTitle(
"FCAL Pedestal RMS [ADC Counts]" );
76 for(
int x = 1;
x <= digPed2DAvg->GetNbinsX(); ++
x ){
77 for(
int y = 1;
y <= digPed2DAvg->GetNbinsY(); ++
y ){
79 double var = digPedSq2DAvg->GetBinContent(
x,
y );
80 var -= ( digPed2DAvg->GetBinContent(
x,
y ) *
81 digPed2DAvg->GetBinContent(
x,
y ) );
83 if( digOcc2D->GetBinContent(
x,
y ) != 0 ){
85 digPed2DRMS->SetBinContent(
x,
y, TMath::Sqrt( var ) );
91 digPed2DRMS->SetBinContent(
x,
y, -1 );
96 digPed2DRMS->SetStats( 0 );
97 digPed2DRMS->SetMaximum( 10 );
98 digPed2DRMS->SetMinimum( 0 );
100 digPed2DRMS->Draw(
"colz" );