PWGHF/vertexingHF/macros/PlotOutputMCCheck.C

   1 #if !defined(__CINT__) || defined(__MAKECINT__)
   2 #include <TCanvas.h>
   3 #include <TH1.h>
   4 #include <TH2.h>
   5 #include <TStyle.h>
   6 #include <TFile.h>
   7 #endif
   8
   9 /* $Id$ */
  10
  11 // Macro to plot the output of AliAnalysisTaskCheckHFMCProd
  12 // Author: F. Prino, prino@to.infn.it
  13
  14
  15 void PlotOutputMCCheck(){
  16   TFile *fil=new TFile("AnalysisResultsMerged.root");
  17   TDirectoryFile* df=(TDirectoryFile*)fil->Get("HFMCCheck");
  18   TList* l=(TList*)df->Get("clistHFMCCheck");
  19   l->ls();
  20
  21   TH1F* hNEvents=(TH1F*)l->FindObject("hNEvents");
  22   Int_t nAnalEv=hNEvents->GetBinContent(1);
  23   printf("Number of events= %d\n",nAnalEv);
  24
  25
  26   TCanvas* cv=new TCanvas("cv","Vertex");
  27   cv->Divide(3,3);
  28   cv->cd(1);
  29   TH1F* hSPD3DvX=(TH1F*)l->FindObject("hSPD3DvX");
  30   hSPD3DvX->Draw();
  31   cv->cd(2);
  32   TH1F* hSPD3DvY=(TH1F*)l->FindObject("hSPD3DvY");
  33   hSPD3DvY->Draw();
  34   cv->cd(3);
  35   TH1F* hSPD3DvZ=(TH1F*)l->FindObject("hSPD3DvZ");
  36   hSPD3DvZ->Draw();
  37   cv->cd(4);
  38   TH1F* hSPDZvX=(TH1F*)l->FindObject("hSPDZvX");
  39   hSPDZvX->Draw();
  40   cv->cd(5);
  41   TH1F* hSPDZvY=(TH1F*)l->FindObject("hSPDZvY");
  42   hSPDZvY->Draw();
  43   cv->cd(6);
  44   TH1F* hSPDZvZ=(TH1F*)l->FindObject("hSPDZvZ");
  45   hSPDZvZ->Draw();
  46   cv->cd(7);
  47   TH1F* hTRKvX=(TH1F*)l->FindObject("hTRKvX");
  48   hTRKvX->Draw();
  49   cv->cd(8);
  50   TH1F* hTRKvY=(TH1F*)l->FindObject("hTRKvY");
  51   hTRKvY->Draw();
  52   cv->cd(9);
  53   TH1F* hTRKvZ=(TH1F*)l->FindObject("hTRKvZ");
  54   hTRKvZ->Draw();
  55
  56   TCanvas* c1=new TCanvas("c1","Multipl");
  57   c1->Divide(3,1);
  58   c1->cd(1);
  59   TH1F* hTracklets=(TH1F*)l->FindObject("hTracklets");
  60   hTracklets->Draw();
  61   c1->cd(2);
  62   TH1F* hTracks=(TH1F*)l->FindObject("hTracks");
  63   hTracks->Draw();
  64   c1->cd(3);
  65   TH1F* hSelTracks=(TH1F*)l->FindObject("hSelTracks");
  66   hSelTracks->Draw();
  67
  68   TH1F* hncharmed=(TH1F*)l->FindObject("hncharmed");
  69   TCanvas* cn=new TCanvas("cn","ncharm");
  70   hncharmed->Draw("box");
  71   hncharmed->GetXaxis()->SetTitle("dNch/dy");
  72   hncharmed->GetYaxis()->SetTitle("N Charm hadrons in golden channels");
  73   cn->Update();
  74
  75   TH1F* hnbvsnc=(TH1F*)l->FindObject("hnbvsnc");
  76   TCanvas* cnhf=new TCanvas("cnhf","nb/c");
  77   hnbvsnc->Draw("box");
  78   hnbvsnc->GetXaxis()->SetTitle("Nc");
  79   hnbvsnc->GetYaxis()->SetTitle("Nb");
  80   cnhf->Update();
  81
  82   TH2F* hyptd0prompt=(TH2F*)l->FindObject("hyptd0prompt");
  83   TH2F* hyptd0feeddown=(TH2F*)l->FindObject("hyptd0feeddown");
  84   TH2F* hyptD02=(TH2F*)l->FindObject("hyptD02");
  85   TH2F* hyptD04=(TH2F*)l->FindObject("hyptD04");
  86
  87   TCanvas* cd0=new TCanvas("cd0","D0");
  88   cd0->Divide(2,2);
  89   cd0->cd(1);
  90   hyptd0prompt->Draw("colz");
  91   cd0->cd(2);
  92   hyptd0feeddown->Draw("colz");
  93   cd0->cd(3);
  94   hyptD02->Draw("colz");
  95   cd0->cd(4);
  96   hyptD04->Draw("colz");
  97
  98   TH2F* hyptdplusprompt=(TH2F*)l->FindObject("hyptdplusprompt");
  99   TH2F* hyptdplusfeeddown=(TH2F*)l->FindObject("hyptdplusfeeddown");
 100   TH2F* hyptDplusnonreson=(TH2F*)l->FindObject("hyptDplusnonreson");
 101   TH2F* hyptDplusreson=(TH2F*)l->FindObject("hyptDplusreson");
 102
 103   TCanvas* cdplus=new TCanvas("cdplus","Dplus");
 104   cdplus->Divide(2,2);
 105   cdplus->cd(1);
 106   hyptdplusprompt->Draw("colz");
 107   cdplus->cd(2);
 108   hyptdplusfeeddown->Draw("colz");
 109   cdplus->cd(3);
 110   hyptDplusnonreson->Draw("colz");
 111   cdplus->cd(4);
 112   hyptDplusreson->Draw("colz");
 113
 114   TH2F* hyptdsprompt=(TH2F*)l->FindObject("hyptdsprompt");
 115   TH2F* hyptdsfeeddown=(TH2F*)l->FindObject("hyptdsfeedown");
 116   TH2F* hyptdsphi=(TH2F*)l->FindObject("hyptdsphi");
 117   TH2F* hyptdsK0st=(TH2F*)l->FindObject("hyptdsk0st");
 118
 119   TCanvas* cds=new TCanvas("cds","Ds");
 120   cds->Divide(2,2);
 121   cds->cd(1);
 122   hyptdsprompt->Draw("colz");
 123   cds->cd(2);
 124   hyptdsfeeddown->Draw("colz");
 125   cds->cd(3);
 126   hyptdsphi->Draw("colz");
 127   cds->cd(4);
 128   hyptdsK0st->Draw("colz");
 129
 130   TH2F* hyptdstarprompt=(TH2F*)l->FindObject("hyptdstarprompt");
 131   TH2F* hyptdstarfeedown=(TH2F*)l->FindObject("hyptdstarfeedown");
 132   TH2F* hyptlcprompt=(TH2F*)l->FindObject("hyptlcprompt");
 133   TH2F* hyptlcfeedown=(TH2F*)l->FindObject("hyptlcfeedown");
 134
 135   TCanvas* cdstlc=new TCanvas("cdstls","Dstar LambdaC");
 136   cdstlc->Divide(2,2);
 137   cdstlc->cd(1);
 138   hyptdstarprompt->Draw("colz");
 139   cdstlc->cd(2);
 140   hyptdstarfeedown->Draw("colz");
 141   cdstlc->cd(3);
 142   hyptlcprompt->Draw("colz");
 143   cdstlc->cd(4);
 144   hyptlcfeedown->Draw("colz");
 145
 146 }