prod/acrcaf/qa_pp/plot_macros/plotD.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 TH1F* fHpos;
  13 TH1F* fHneg;
  14
  15 fHpos=  (TH1F*) l->FindObject("fRecDcaPos");
  16 fHpos->SetLineWidth(2);
  17
  18 fHneg=  (TH1F*) l->FindObject("fRecDcaNeg");
  19 fHneg->SetLineColor(kRed);
  20 fHneg->SetLineWidth(2);
  21
  22 TH1F* fHposD;
  23 TH1F* fHnegD;
  24
  25 fHposD=  (TH1F*) l->FindObject("fRecDPos");
  26 fHposD->SetLineWidth(2);
  27
  28 fHnegD=  (TH1F*) l->FindObject("fRecDNeg");
  29 fHnegD->SetLineColor(kRed);
  30 fHnegD->SetLineWidth(2);
  31
  32
  33 fH=  (TH2F*) l->FindObject("fDiffDcaD");
  34
  35
  36 TLegend *legp;
  37 legp= new TLegend(0.9,0.65,0.65,0.9);
  38 legp->SetFillColor(kWhite);
  39
  40 TF1 *fun0, *fun1;
  41 fun0= new TF1("fun0","gaus",-5.0,5.0);
  42 fun1= new TF1("fun1","gaus",-5.0,5.0);
  43 fun0->SetLineColor(kBlack);
  44 fun1->SetLineColor(kRed);
  45
  46 legp->AddEntry(fun0,"positive","l");
  47 legp->AddEntry(fun1,"negative","l");
  48
  49 TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 480);
  50 dca->Divide(2,2);
  51 dca->cd(1);
  52 fHpos->DrawClone("");
  53 fHpos->SetTitle("");
  54 fHpos->SetYTitle("");
  55 fHneg->DrawClone("same");
  56 gPad->SetLogy();
  57 legp->Draw();
  58
  59 dca->cd(2);
  60 fHposD->DrawClone();
  61 fHnegD->DrawClone("same");
  62 legp->Draw();
  63 gPad->SetLogy();
  64
  65
  66 dca->cd(3);
  67 fH->Draw("colz");
  68
  69 }
  70
  71
  72
  73 }