prod/acrcaf/qa_pp/plot_macros/plotQPtVz.C

   1 {
   2 gROOT->SetStyle("Plain");
   3 gStyle->SetOptStat(0);
   4 gStyle->SetPalette(1);
   5
   6 //open files and get lists
   7 TList *l;
   8 TFile f0("QAsym.proof.root");
   9 l = (TList *)f0.Get("QAsymHists");
  10
  11
  12 TH2F* fHpos;
  13 TH2F* fHpneg;
  14
  15 fHpos=  (TH2F*) l->FindObject("fRecQPtPosEtaVz");
  16 fHneg=  (TH2F*) l->FindObject("fRecQPtNegEtaVz");
  17
  18
  19
  20
  21
  22
  23 TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 520);
  24 dca->Divide(2,1);
  25 dca->cd(1);
  26 //fHpos->Divide(fHneg);
  27 gPad->SetRightMargin(0.15);
  28 gPad->SetLeftMargin(0.15);
  29 fHpos->SetTitleOffset(1.8,"Y");
  30 //fHpos->SetXTitle("(Charge/p_{T})_{pos} / (Charge/p_{T})_{neg}");
  31 fHpos->SetXTitle("Charge/p_{T}");
  32 fHpos->Draw("colz");
  33 fHpos->SetTitle("positive #eta");
  34 fHpos->SetYTitle("esdtrack->Zv()");
  35
  36 dca->cd(2);
  37 fHneg->Draw("colz");
  38 gPad->SetRightMargin(0.15);
  39 gPad->SetLeftMargin(0.15);
  40 fHneg->SetTitleOffset(1.8,"Y");
  41 fHneg->SetXTitle("Charge/p_{T}");
  42 fHneg->SetMaximum(fHpos->GetMaximum());
  43 fHneg->SetTitle("negative #eta");
  44 fHneg->SetYTitle("esdtrack->Zv()");
  45
  46
  47
  48 }
  49
  50
  51
  52 }