- TString name1(Form("NjUEvsNj_%d",i));
- fHistNjUEvsNj[i] = new TH2F(name1, name1, 150, -0.5, 149.5, 150, -0.5, 149.5);
- fHistNjUEvsNj[i]->GetXaxis()->SetTitle("N_{jet}");
- fHistNjUEvsNj[i]->GetYaxis()->SetTitle("N_{jet_{UE}}");
- fOutput->Add(fHistNjUEvsNj[i]);
-
- TString name2(Form("NjUEoverNjVsNj_%d",i));
- fHistNjUEoverNjVsNj[i] = new TH2F(name2, name2, 150, -0.5, 149.5, 120, -0.01, 1.19);
- fHistNjUEoverNjVsNj[i]->GetXaxis()->SetTitle("N_{jet_{UE}} / N_{jet}");
- fHistNjUEoverNjVsNj[i]->GetYaxis()->SetTitle("N_{jet}");
- fOutput->Add(fHistNjUEoverNjVsNj[i]);
+ name = Form("fHistJetNconstVsPt_%d",i);
+ fHistJetNconstVsPt[i] = new TH2F(name, name, 150, -0.5, 149.5, fNbins, fMinBinPt, fMaxBinPt);
+ fHistJetNconstVsPt[i]->GetXaxis()->SetTitle("# of constituents");
+ fHistJetNconstVsPt[i]->GetYaxis()->SetTitle("p_{T,jet} (GeV/c)");
+ fOutput->Add(fHistJetNconstVsPt[i]);
+
+ for (Int_t j = 0; j < 3; j++) {
+ name = Form("NjUEoverNjVsNj_%d_%d",i,j+1);
+ fHistNjUEoverNjVsNj[i*3+j] = new TH2F(name, name, 150, -0.5, 149.5, 120, 0.01, 1.21);
+ fHistNjUEoverNjVsNj[i*3+j]->GetXaxis()->SetTitle("N_{jet}");
+ fHistNjUEoverNjVsNj[i*3+j]->GetYaxis()->SetTitle("N_{jet_{UE}} / N_{jet}");
+ fOutput->Add(fHistNjUEoverNjVsNj[i*3+j]);
+ }