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