prod/acrcaf/qa_pp/plot_macros/plotSectorQPtRatio.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[18];
  13
  14 for(Int_t sector=0;sector<18;sector++){
  15  fHpos[sector]=  (TH1F*) l->FindObject(Form("fRecQPtTpcSector%02d", sector));
  16 }
  17
  18 TH1F* fHposref =(TH1F*) l->FindObject("fQPt");
  19
  20 TLegend *legp;
  21 legp= new TLegend(0.9,0.65,0.65,0.9);
  22 legp->SetFillColor(kWhite);
  23
  24 TF1 *fun0, *fun1;
  25 fun0= new TF1("fun0","gaus",-5.0,5.0);
  26 fun1= new TF1("fun1","gaus",-5.0,5.0);
  27 fun0->SetLineColor(kBlack);
  28 fun1->SetLineColor(kRed);
  29
  30 legp->AddEntry(fun0,"positive","l");
  31 legp->AddEntry(fun1,"negative","l");
  32
  33
  34
  35 for(Int_t i=0; i<18;i++){
  36   fHpos[i]->SetLineColor(kRed);
  37   // fHpos[i]->SetLineWidth(2);
  38 }
  39
  40
  41 TCanvas * dca = new TCanvas("pt", "pt", 100, 100, 1020, 820);
  42 dca->Divide(6,3);
  43 for(Int_t i=0; i<18;i++){
  44   fHpos[i]->Divide(fHposref);
  45   dca->cd(i+1);
  46
  47   fHpos[i]->Draw();
  48
  49   fHpos[i]->SetMaximum(0.3);
  50   fHpos[i]->SetMinimum(-0.1);
  51 //   fHpos[i]->SetMaximum(100);
  52 //   fHpos[i]->SetMinimum(1)
  53     ;
  54
  55   fHpos[i]->SetTitle("");
  56   fHpos[i]->SetXTitle(Form("(Q/p_{T})_{%d}/(Q/p_{T})_{all}",i));
  57   //fHpos[i]->SetXTitle(Form("Q/p_{T} sector %2d",i));
  58
  59   gPad->SetLeftMargin(0.20);
  60   gPad->SetBottomMargin(0.17);
  61
  62   fHpos[i]->SetTitleOffset(1., "X");
  63   fHpos[i]->SetTitleOffset(.8, "Y");
  64   fHpos[i]->SetTitleSize(0.08, "X");
  65   fHpos[i]->SetTitleSize(0.05, "Y");
  66   fHpos[i]->SetLabelSize(0.1, "X");
  67   fHpos[i]->SetLabelSize(0.1, "Y");
  68   fHpos[i]->GetYaxis()->SetNdivisions(5);
  69   fHpos[i]->GetXaxis()->SetNdivisions(5);
  70
  71
  72   //gPad->SetLogy();
  73
  74 }
  75
  76
  77
  78 }