+
+ if (!fHistList) fHistList = new TList();
+ Float_t pi=TMath::Pi();
+ gStyle->SetPalette(1);
+
+ fH1DeltaPt = new TH1F("DeltaPt","Difference between the jets' Pt;#Deltap_{T} (GeV/c);Entries",150,0.,150.);
+ fH1DeltaPt->SetMarkerSize(0.6);
+ fH1DeltaPt->SetMarkerColor(4);
+ fH1DeltaPt->SetMarkerStyle(21);
+ fH1DeltaPt->SetOption("E");
+
+ fH1DeltaPhi = new TH1F("DeltaPhi","Difference in the azimuthal angle;#Delta#phi;Entries",100,0.,pi);
+ fH1DeltaPhi->SetMarkerSize(0.6);
+ fH1DeltaPhi->SetMarkerColor(4);
+ fH1DeltaPhi->SetMarkerStyle(21);
+ fH1DeltaPhi->SetOption("E");
+
+ fH1PhiImbal = new TH1F("PhiImb","Phi imbalance;#phi;Entries",100,-pi,pi);
+ fH1PhiImbal->SetMarkerSize(0.6);
+ fH1PhiImbal->SetMarkerColor(4);
+ fH1PhiImbal->SetMarkerStyle(21);
+ fH1PhiImbal->SetOption("E");
+
+ fH1Asym = new TH1F("Asym","Pt asymmetry;#Deltap_{T}/(p_{T,1}+p_{T,2});Entries",50,0.,1.);
+ fH1Asym->SetMarkerSize(0.6);
+ fH1Asym->SetMarkerColor(4);
+ fH1Asym->SetMarkerStyle(21);
+ fH1Asym->SetOption("E");
+
+ fH2Pt2vsPt1 = new TH2F("Pt2vsPt1","Pt2 versus Pt1;p_{T,1} (GeV/c);p_{T,2} (GeV/c)",250,0.,250.,250,0.,250.);
+ fH2Pt2vsPt1->SetOption("cont0");
+
+ fH2DifvsSum = new TH2F("DifvsSum","Pt difference vs Pt sum;p_{T,1}+p_{T,2} (GeV/c);#Deltap_{T} (GeV/c)",400,0.,400.,150,0.,150.);
+ fH2DifvsSum->SetOption("cont0");
+
+ fHistList->Add(fH1DeltaPt);
+ fHistList->Add(fH1DeltaPhi);
+ fHistList->Add(fH1PhiImbal);
+ fHistList->Add(fH1Asym);
+ fHistList->Add(fH2Pt2vsPt1);
+ fHistList->Add(fH2DifvsSum);