make_eop_maps.cxx File Reference

#include "TH2.h"
#include "TCanvas.h"
#include "TFile.h"
#include "TString.h"

Include dependency graph for make_eop_maps.cxx:

Go to the source code of this file.

Functions
void	make_eop_maps (TString path, Float_t range=0.2)
void	make_sagitta_maps (TString path, Float_t range=0.004)

Function Documentation

make_eop_maps()

void make_eop_maps	(	TString	path,
		Float_t	range = 0.2 )

Definition at line 21 of file make_eop_maps.cxx.

                                                 { //zmin=0.8, Float_t zmax=1.2){
 
  // Make the ratio histogram from the input file
  TFile *f = new TFile(path);
  TH2D *h1 = (TH2D*)f->Get("meanPlotPos");
  TH2D *h2 = (TH2D*)f->Get("meanPlotNeg");
  TH2D *outhist = h1->Clone();
  outhist->Divide(h2);
 
  // Set up the output canvas
  TCanvas *c1 = new TCanvas("c1","c1",800,600);
  c1->SetRightMargin(0.13);
 
  // Set up ranges
  Float_t zmin = 1.-range;
  Float_t zmax = 1.+range;
 
  // Plot aesthetics
  outhist->SetTitle("E/p ratio <+>/<->");
  outhist->GetZaxis()->SetRangeUser(zmin,zmax);
  outhist->GetXaxis()->SetTitle("#eta");
  outhist->GetXaxis()->SetTitleOffset(1);
  outhist->GetYaxis()->SetTitle("#phi [rad]");
  outhist->GetYaxis()->SetTitleOffset(.7);
  TPaletteAxis *palette = (TPaletteAxis*)outhist->GetListOfFunctions()->FindObject("palette");
  palette->SetX1NDC(0.87);
  palette->SetX2NDC(0.92);
  palette->SetY1NDC(0.1);
  palette->SetY2NDC(0.9);
 
  // Draw and print
  outhist->Draw("colz");
  c1->Print("EoverP.pdf","Portrait pdf");
 
  return;
}

make_sagitta_maps()

void make_sagitta_maps	(	TString	path,
		Float_t	range = 0.004 )

Definition at line 58 of file make_eop_maps.cxx.

                                                       {
 
  // Make the ratio histogram from the input file
  TFile *f = new TFile(path);
  TH2D *outhist = (TH2D*)f->Get("FinalCorrections");
 
  // Set up the output canvas
  TCanvas *c1 = new TCanvas("c1","c1",800,600);
  c1->SetRightMargin(0.16);
 
  // Set up ranges
  range = fabs(range);
  Float_t zmin = (-1.)*range;
  Float_t zmax = range;
 
  // Plot aesthetics
  outhist->SetTitle("Sagitta Bias");
  outhist->GetZaxis()->SetRangeUser(zmin,zmax);
  outhist->GetXaxis()->SetTitle("#eta");
  outhist->GetXaxis()->SetTitleOffset(1);
  outhist->GetYaxis()->SetTitle("#phi [rad]");
  outhist->GetYaxis()->SetTitleOffset(.7);
  TPaletteAxis *palette = (TPaletteAxis*)outhist->GetListOfFunctions()->FindObject("palette");
  palette->SetX1NDC(0.84);
  palette->SetX2NDC(0.89);
  palette->SetY1NDC(0.1);
  palette->SetY2NDC(0.9);
  TGaxis::SetMaxDigits(2);
 
  // Draw and print
  outhist->Draw("colz");
  c1->Print("Sagitta.pdf","Portrait pdf");
 
  return;
}