prod/acrcaf/qa_pp/plot_macros/plot.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[7];
  13 TH1F* fHneg[7];
  14
  15
  16 fHpos[0]=  (TH1F*) l->FindObject("fRecPtPos");
  17 fHneg[0]=  (TH1F*) l->FindObject("fRecPtNeg");
  18 fHpos[1]=  (TH1F*) l->FindObject("fRecPhiPos");
  19 fHneg[1]=  (TH1F*) l->FindObject("fRecPhiNeg");
  20 fHpos[2]=  (TH1F*) l->FindObject("fRecEtaPos");
  21 fHneg[2]=  (TH1F*) l->FindObject("fRecEtaNeg");
  22 fHpos[3]=  (TH1F*) l->FindObject("fRecEtaPtPos");
  23 fHneg[3]=  (TH1F*) l->FindObject("fRecEtaPtNeg");
  24
  25 fHpos[4]=  (TH1F*) l->FindObject("fRecPtPosEta");
  26 fHneg[4]=  (TH1F*) l->FindObject("fRecPtNegEta");
  27 fHpos[5]=  (TH1F*) l->FindObject("fRecPhiPosEta");
  28 fHneg[5]=  (TH1F*) l->FindObject("fRecPhiNegEta");
  29 fHpos[6]=  (TH1F*) l->FindObject("fRecQPtPosEta");
  30 fHneg[6]=  (TH1F*) l->FindObject("fRecQPtNegEta");
  31
  32
  33 TLegend *legp1;
  34 legp1= new TLegend(0.9,0.65,0.65,0.9);
  35 legp1->SetFillColor(kWhite);
  36
  37 TLegend *legp2;
  38 legp2= new TLegend(0.9,0.65,0.65,0.9);
  39 legp2->SetFillColor(kWhite);
  40
  41 TF1 *fun0, *fun1;
  42 fun0= new TF1("fun0","gaus",-5.0,5.0);
  43 fun1= new TF1("fun1","gaus",-5.0,5.0);
  44 fun0->SetLineColor(kBlack);
  45 fun1->SetLineColor(kRed);
  46
  47 legp1->AddEntry(fun0,"pos. charge","l");
  48 legp1->AddEntry(fun1,"neg. charge","l");
  49
  50 legp2->AddEntry(fun0,"pos. #eta","l");
  51 legp2->AddEntry(fun1,"neg. #eta","l");
  52
  53
  54
  55 for(Int_t i=0; i<7;i++){
  56   fHpos[i]->SetLineColor(kBlack);
  57   fHpos[i]->SetLineWidth(2);
  58   fHneg[i]->SetLineColor(kRed);
  59   fHneg[i]->SetLineWidth(2);
  60 }
  61
  62
  63 TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
  64 dca->Divide(3,3);
  65 for(Int_t i=0; i<7;i++){
  66   dca->cd(i+1);
  67   fHpos[i]->Draw();
  68   fHneg[i]->Draw("same");
  69   if(i==0||i==4)gPad->SetLogy();
  70   if (i < 4)
  71     legp1->Draw();
  72   else
  73     legp2->Draw();
  74 }
  75
  76
  77
  78 }