//---------------------------------------------------------------------- /// show the jets contained in filename (as produced by /// ClusterSequence::print_jets_for_root()), with an optional label TCanvas * showjets (const char * filename, const char * label = 0) { // display the various 2-d drawing options gROOT->Reset(); // set up canvas TCanvas * lego = new TCanvas("lego","lego options",400,50,800,600); lego->SetTheta(30.0); lego->SetPhi(20.0); // orientation used for plots in subtraction paper //lego->SetTheta(62.15); //lego->SetPhi(9.15); ////vector<double> col set_default_colours(lego); TPaveLabel pl; JetHist * jets = new JetHist(filename); jets->stack.Draw("lego1"); // cyl does not work with 5.16 if (label != 0) { Float_t x1=0.63, y1=0.875, x2=0.95, y2=0.925; pl.DrawPaveLabel(x1,y1,x2,y2,label,"brNDC"); } else if (jets->comment() != "") { Float_t x1=0.15, y1=0.875, x2=0.95, y2=0.925; pl.DrawPaveLabel(x1,y1,x2,y2,jets->comment().c_str(),"brNDC"); } // normal histogram labels not working, so draw them by hand TLatex l; l.SetTextAlign(22); l.SetTextSize(0.05); //l.DrawLatex(0.0,0.85,"anti-k_{t}, R=1"); l.SetTextSize(0.04); l.DrawLatex(0.20,-0.98,"y"); l.SetTextAlign(32); l.DrawLatex(-0.7,0.8,"p_{t} [GeV]"); l.DrawLatex(-0.6,-0.78,"#phi"); // do not delete jets -- otherwise you lose everything!; return lego; /// lego->Update(); }