draw_hadron_check_plots.cxx

Go to the documentation of this file.

#include <iostream>
#include <fstream>
#include <iomanip>
#include <cstring>
#include <string>
#include <vector>
#include <TFile.h>
#include <TTree.h>
#include <TString.h>
#include <stdlib.h>
#include <stdio.h>
#include <TH2F.h>
#include <TCanvas.h>
#include <TColor.h>
 
using namespace std;
 
int main(int argc, char* argv[]){
  if(argc<2){
    cout << "please append the root file name" << endl;
    cout << "usage: draw_hadron_check_plots.exe filename" << endl;
    return 1;
  }
 
  TString str_filename(argv[1]);
  TFile *f = new TFile(str_filename);
  TTree *tree = (TTree*)f->Get("n103");
  TH2F *h = new TH2F("h", "dE/dx vs momentum", 48, 0.05, 2.45, 160, 0, 2000); 
  TCanvas *c1 = new TCanvas("c1","c1",10,10,800,642);
  tree->Draw("dEdx_meas:ptrk>>h");
  const UInt_t Number = 2;
  Double_t Red[Number]    = { 1.0,  0.2};
  Double_t Green[Number]  = { 1.0,  0.0};
  Double_t Blue[Number]   = { 1.0,  0.0};
  Double_t Length[Number] = { 0.0,  1.0};
  Int_t nb=200;
  TColor::CreateGradientColorTable(Number,Length,Red,Green,Blue,nb);
  h->Draw();
  c1->SaveAs("temp.eps");
  f->Close();
 
  return 0;
}