PWGGA/PHOSTasks/PHOS_PbPb/macros/DrawPi0Flow.C

   1 DrawPi0Flow(const TString filename = "Pi0Flow_000167920.root")
   2 {
   3   gStyle->SetOptStat(0);
   4   TFile * f = new TFile(filename) ;
   5   TList *histoList = (TList*)f->Get("PHOSPi0Flow");
   6
   7   TH1F *hev        = (TH1F*)histoList->FindObject("hTotSelEvents") ;
   8   TH2F* hZvertex   = (TH2F*)histoList->FindObject("hZvertex");
   9   TH2F* hCenTrack  = (TH2F*)histoList->FindObject("hCenTrack");
  10   TH2F *hCenPHOS   = (TH2F*)histoList->FindObject("hCenPHOS") ;
  11   TH2F* phiRPV0A   = (TH2F*)histoList->FindObject("phiRPV0A");
  12   TH2F* phiRPV0C   = (TH2F*)histoList->FindObject("phiRPV0C");
  13   TH2F* hCellNXZM1 = (TH2F*)histoList->FindObject("hCellNXZM1");
  14   TH2F* hCellNXZM2 = (TH2F*)histoList->FindObject("hCellNXZM2");
  15   TH2F* hCellNXZM3 = (TH2F*)histoList->FindObject("hCellNXZM3");
  16   TH2F *hPi0All_cen0     = (TH2F*)histoList->FindObject("hPi0All_cen0");
  17   TH2F *hPi0Allcore_cen0 = (TH2F*)histoList->FindObject("hPi0Allcore_cen0");
  18
  19   //-----------------------------------------------------------------------------
  20   TCanvas *c1 = new TCanvas("c1","Event selection");
  21   hev->GetXaxis()->SetBinLabel(1,"Total");
  22   hev->GetXaxis()->SetBinLabel(2,"Total");
  23   hev->GetXaxis()->SetBinLabel(3,"Goog vertex");
  24   hev->GetXaxis()->SetBinLabel(4,"No pileup");
  25   hev->GetXaxis()->SetBinLabel(5,"Good centrality");
  26   hev->GetXaxis()->SetBinLabel(6,"Good RP");
  27   hev->GetXaxis()->SetBinLabel(7,"MC");
  28   hev->GetXaxis()->SetBinLabel(8,"QA filled");
  29   hev->GetXaxis()->SetBinLabel(9,"Clusters filled");
  30   hev->GetXaxis()->SetBinLabel(10,"Pi0 filled");
  31   hev->GetXaxis()->SetBinLabel(11,"Mixed filled");
  32   hev->GetXaxis()->SetBinLabel(12,"Lists updated");
  33   hev->SetYTitle("N_{events}");
  34   hev->Draw();
  35   c1->Print("PHOS_EvSel.eps");
  36
  37   //-----------------------------------------------------------------------------
  38   TCanvas *c2 = new TCanvas("c2","Track multiplicity vs centrality");
  39   hCenTrack->SetXTitle("centrality (%)");
  40   hCenTrack->SetYTitle("Number of tracks");
  41   hCenTrack->Draw("colz");
  42   c2->Print("TrackMultCentrality.eps");
  43
  44   //-----------------------------------------------------------------------------
  45   TCanvas *c3 = new TCanvas("c3","PHOS multiplicity vs centrality");
  46   hCenPHOS->SetXTitle("centrality (%)");
  47   hCenPHOS->SetYTitle("Number of PHOS clusters");
  48   hCenPHOS->Draw("colz");
  49   c3->Print("PHOSMultCentrality.eps");
  50
  51   //-----------------------------------------------------------------------------
  52   TCanvas *c4 = new TCanvas("c4","RP phi");
  53   phiRPV0A1 = phiRPV0A->ProjectionX();
  54   phiRPV0C1 = phiRPV0C->ProjectionX();
  55   phiRPV0A1->SetLineColor(kRed);
  56   phiRPV0C1->SetLineColor(kBlue);
  57   phiRPV0A1->SetXTitle("RP #phi");
  58   phiRPV0C1->SetXTitle("RP #phi");
  59   phiRPV0A1->Draw();
  60   phiRPV0C1->Draw("same");
  61   leg = new TLegend(0.7,0.8,0.89,0.89);
  62   leg->SetFillColor(kWhite);
  63   leg->SetBorderSize(1);
  64   leg->AddEntry(phiRPV0A1,"V0A","l");
  65   leg->AddEntry(phiRPV0C1,"V0C","l");
  66   leg->Draw();
  67   c4->Print("V0RPphi.eps");
  68
  69   //-----------------------------------------------------------------------------
  70   TCanvas *c5 = new TCanvas("c5","XZ in M1");
  71   c5->SetLogz();
  72   hCellNXZM1->SetTitle("Cell occupancy in module 4");
  73   hCellNXZM1->SetXTitle("X (cells)");
  74   hCellNXZM1->SetYTitle("X (cells)");
  75   hCellNXZM1->Draw("colz");
  76   c5->Print("PHOS_XYM4.eps");
  77
  78   //-----------------------------------------------------------------------------
  79   TCanvas *c6 = new TCanvas("c6","XZ in M2");
  80   c6->SetLogz();
  81   hCellNXZM2->SetTitle("Cell occupancy in module 3");
  82   hCellNXZM2->SetXTitle("X (cells)");
  83   hCellNXZM2->SetYTitle("X (cells)");
  84   hCellNXZM2->Draw("colz");
  85   c6->Print("PHOS_XYM3.eps");
  86
  87   //-----------------------------------------------------------------------------
  88   TCanvas *c7 = new TCanvas("c7","XZ in M3");
  89   c7->SetLogz();
  90   hCellNXZM3->SetTitle("Cell occupancy in module 2");
  91   hCellNXZM3->SetXTitle("X (cells)");
  92   hCellNXZM3->SetYTitle("X (cells)");
  93   hCellNXZM3->Draw("colz");
  94   c7->Print("PHOS_XYM2.eps");
  95
  96   //-----------------------------------------------------------------------------
  97   TCanvas *c8 = new TCanvas("c8","gg mass vs pt, PID=All");
  98   c8->SetLogz();
  99   hPi0All_cen0->SetTitle("PID: All");
 100   hPi0All_cen0->SetXTitle("M_{#gamma#gamma} (GeV/c^{2})");
 101   hPi0All_cen0->SetYTitle("p_{T} (GeV/c)");
 102   hPi0All_cen0->Draw("colz");
 103   c3->Print("PHOS_MggAll.eps");
 104
 105   //-----------------------------------------------------------------------------
 106   TCanvas *c9 = new TCanvas("c9","gg mass vs pt, PID=Allcore");
 107   c9->SetLogz();
 108   hPi0Allcore_cen0->SetTitle("PID: Allcore");
 109   hPi0Allcore_cen0->SetXTitle("M_{#gamma#gamma} (GeV/c^{2})");
 110   hPi0Allcore_cen0->SetYTitle("p_{T} (GeV/c)");
 111   hPi0Allcore_cen0->Draw("colz");
 112   c9->Print("PHOS_MggAllcore.eps");
 113
 114 }