plotTurnOn_multiMethods.C

#include <TFile.h>
#include <TH1D.h>
#include <TCanvas.h>
#include <TGraphAsymmErrors.h>
#include <TLegend.h>
#include <TLine.h>
#include <iostream>

void plotTurnOn()
{
  const Int_t INPUTFILES = 4;


  // TString outFileTag = "MBData_500Hz";
  // TString inFileName[INPUTFILES] = {"hist_MBData_nominal.root",
  // 				    "hist_MBData_sigmaSubtraction.root",
  // 				    "hist_MBData_oneByOne.root",
  // 				    "hist_MBData_twoByTwo.root"
  // };
  // // TString labels[INPUTFILES] = {"nominal, L1_52, 158Hz",
  // // 				"sigma subtracted, L1_24, 187Hz",
  // // 				"1x1 jets, L1_24, 180Hz",
  // // 				"2x2 jets, L1_52, 210Hz"};
  // // const Double_t L1_THRESHOLD[INPUTFILES] = {52, 24, 24, 40};
  // TString labels[INPUTFILES] = {"nominal, L1_44, 550Hz",
  // 				"sigma subtracted, L1_20, 505Hz",
  // 				"1x1 jets, L1_20, 441Hz",
  // 				"2x2 jets, L1_36, 444Hz"};
  // const Double_t L1_THRESHOLD[INPUTFILES] = {44, 20, 20, 36};

  TString outFileTag = "Hydjet276_150Hz";
  TString inFileName[INPUTFILES] = {"hist_Hydjet276_nominal.root",
				    "hist_Hydjet276_sigmaSubtraction.root",
				    "hist_Hydjet276_oneByOne.root",
				    "hist_Hydjet276_twoByTwo.root"};
  //Double_t 
  

  // TString outFileTag = "Hydjet502_800Hz";
  // TString inFileName[INPUTFILES] = {"hist_Hydjet502_nominal.root",
  // 				    "hist_Hydjet502_sigmaSubtraction.root",
  // 				    "hist_Hydjet502_oneByOne.root",
  // 				    "hist_Hydjet502_twoByTwo.root"
  // };
  // // // TString labels[INPUTFILES] = {"nominal, L1_80, 150Hz",
  // // // 				"sigma subtracted, L1_36, 163Hz",
  // // // 				"1x1 jets, L1_36, 167Hz",
  // // // 				"2x2 jets, L1_88, 155Hz"};
  // // // const Double_t L1_THRESHOLD[INPUTFILES] = {80, 36, 36, 64};
  // TString labels[INPUTFILES] = {"nominal, L1_68, 602Hz",
  // 				"sigma subtracted, L1_28, 756Hz",
  // 				"1x1 jets, L1_28, 801Hz",
  // 				"2x2 jets, L1_52, 814Hz"};
  // const Double_t L1_THRESHOLD[INPUTFILES] = {68, 28, 28, 52};

  TFile *inFile[INPUTFILES];
  const Int_t COLORS[INPUTFILES] = {kBlack, kRed, kBlue, kGreen+3};//, kMagenta+3};
  TGraphAsymmErrors *asymm[INPUTFILES];//[2];

  for(int i = 0; i < INPUTFILES; i++)
  {
    //for(int j = 0; j < 2; j++)
    {
      inFile[i] = TFile::Open(inFileName[i]);
      asymm[i] = (TGraphAsymmErrors*)inFile[i]->Get(Form("asymm_pt_%d_0",(int)L1_THRESHOLD[i]));
      asymm[i]->SetMarkerColor(COLORS[i]);
      asymm[i]->SetLineColor(COLORS[i]);
    }
    //asymm[i][1]->SetMarkerStyle(25);
  }

  // these values MUST MATCH those used in makeTurnOn.C
  const int nBins = 75;
  const double maxPt = 300;

  TH1D *hEmpty = new TH1D("hEmpty",Form(";Jet p_{T} (GeV);Efficiency"),nBins,0,maxPt);

  TCanvas *c1 = new TCanvas();
  hEmpty->SetMinimum(0);
  hEmpty->SetMaximum(1.2);
  hEmpty->Draw();
  //c1->SetLogy();

  TLine *line = new TLine(0,1,maxPt,1);
  line->Draw();

  for(int i = 0; i < INPUTFILES; i++)
  {
    //for(int j = 0; j < 2; j++)
    {
      asymm[i]->Draw("p");
    }
  }

  TLegend *leg = new TLegend(0.55,0.2,0.9,0.5,"|#eta| < 2");
  leg->SetFillColor(0);
  leg->SetTextFont(42);
  leg->SetTextSizePixels(18);

  for(int i = 0; i < INPUTFILES; i++)
  {
    //for(int j = 0; j < 2; j++)
    {
      leg->AddEntry(asymm[i], Form("%s", labels[i].Data()), "lp");
    }
  }

  leg->Draw();

  c1->SaveAs(Form("%s_turnon.pdf",outFileTag.Data()));

}

int main(int argc, char **argv)
{
  if(argc == 1)
  {
    plotTurnOn();
    return 0;
  }
  else
  {
    std::cout << "Usage:\nplotTurnOn_multiMethods.exe <output_tag>" << std::cout;
    std::cout << "An output pdf will be named <output_tag>_turnon.pdf" << std::cout;
    return 1;
  }
}