3225{
3226
3227 this->SetName("resplots");
3228
3230
3231 TDirectory *trt = nullptr;
3232 if (!isAr)
3233 trt = (TDirectory *)
file->FindKey(
"TRT_all")->ReadObj();
3234 else
3235 trt = (TDirectory *)
file->FindKey(
"TRT_Ar_all")->ReadObj();
3236
3237 TH2F *reshist1 =
nullptr;
3238 if (trt->FindKey("residual"))
3239 reshist1 = (
TH2F *)trt->FindKey(
"residual")->ReadObj();
3240 else
3241 throw(string("residual histograms for entire TRT not found!"));
3242
3243 TH2F *treshist1 =
nullptr;
3244 if (trt->FindKey("timeresidual"))
3245 treshist1 = (
TH2F *)trt->FindKey(
"timeresidual")->ReadObj();
3246 else
3247 throw(string("timeresidual histograms for entire TRT not found!"));
3248
3249 reshist1->SetTitle(Form("residual whole TRT"));
3250 if (isAr)
3251 reshist1->SetTitle(Form("Argon residual whole TRT"));
3252 reshist1->SetFillColor(46);
3253 treshist1->SetTitle(Form("time-residual whole TRT"));
3254 if (isAr)
3255 treshist1->SetTitle(Form("Argon time-residual whole TRT"));
3256 treshist1->SetFillColor(38);
3257
3258 this->Divide(1, 2, 0.01, 0.01);
3259 this->SetLogy();
3260
3262 reshist1->Draw();
3264 treshist1->Draw();
3265}
TH2F(name, title, nxbins, bins_par2, bins_par3, bins_par4, bins_par5=None, bins_par6=None, path='', **kwargs)