Added macro for drawing module specific peak parameterisation using output of MakeRaw...

diff --git a/PWGGA/PHOSTasks/PHOS_PbPb/macros/production/DrawMNN.C b/PWGGA/PHOSTasks/PHOS_PbPb/macros/production/DrawMNN.C
new file mode 100644 (file)
index 0000000..1de349b
--- /dev/null
+++ b/PWGGA/PHOSTasks/PHOS_PbPb/macros/production/DrawMNN.C
@@ -0,0 +1,142 @@
+void DrawMNN()
+{
+  gROOT->LoadMacro("MakeRawProduction.C+g");
+  
+  RawProduction::Output output;
+
+  const char* example1 = "kMB/pt003_hPi0M11";
+  TCanvas* canv = new TCanvas(example1, example1);
+  gStyle->SetOptFit(1);
+  TH1* hist = output.GetHistogram(example1);
+  hist->SetAxisRange(0.05, 0.3);
+  hist->Draw();
+  canv->SaveAs("imgs/kMB_pt003_hPi0M11");
+
+  const char* example2 = "kMB/pt003_hPi0M12";
+  TCanvas* canv = new TCanvas(example2, example2);
+  gStyle->SetOptFit(1);
+  TH1* hist = output.GetHistogram(example2);
+  hist->SetAxisRange(0.05, 0.3);
+  hist->Draw();
+  canv->SaveAs(Form("imgs/kMB_pt003_hPi0M11", example2));
+  
+ 
+  
+  TStringToken triggers("kMB kCentral kSemiCentral kPHOSPb", " ");
+  while(triggers.NextToken()) {
+    canv = new TCanvas(Form("cMean_%s", triggers.Data()), Form("cMean_%s", triggers.Data()));
+
+    RawProduction::TriggerBin triggerBin(triggers);
+    TH1* hM11 = output.GetHistogram("M_M11", triggerBin);
+    TH1* hM12 = output.GetHistogram("M_M12", triggerBin);
+    TH1* hM13 = output.GetHistogram("M_M13", triggerBin);
+    TH1* hM22 = output.GetHistogram("M_M22", triggerBin);
+    TH1* hM23 = output.GetHistogram("M_M23", triggerBin);
+    TH1* hM33 = output.GetHistogram("M_M33", triggerBin);
+
+
+    hM12->SetMarkerStyle(24);
+    hM12->SetMarkerColor(kCyan+1);
+    hM12->SetLineColor(kCyan+1);
+    hM12->GetYaxis()->SetRangeUser(0.12, 0.15);
+    hM12->SetTitle("#pi^{0} Peak Mean");
+    gStyle->SetOptStat(0);
+    hM12->Draw();
+
+    hM13->SetMarkerStyle(25);
+    hM13->SetMarkerColor(kMagenta);
+    hM13->SetLineColor(kMagenta);
+    hM13->Draw("same");
+
+    hM23->SetMarkerStyle(26);
+    hM23->SetMarkerColor(kYellow+1);
+    hM23->SetLineColor(kYellow+1);
+    hM23->Draw("same");
+    
+
+    hM11->SetMarkerStyle(21);
+    hM11->SetMarkerColor(1);
+    hM11->SetLineColor(1);
+    hM11->GetYaxis()->SetRangeUser(0.12, 0.15);
+    hM11->Draw("same");
+
+    hM22->SetMarkerStyle(22);
+    hM22->SetMarkerColor(kRed);
+    hM22->SetLineColor(kRed);
+    hM22->Draw("same");
+
+    hM33->SetMarkerStyle(23);
+    hM33->SetMarkerColor(kGreen);
+    hM33->SetLineColor(kGreen);
+    hM33->Draw("same");
+
+    TLegend* leg = new TLegend(0.8,0.6,0.95,0.95);
+    leg->AddEntry(hM11, "M11", "lep");
+    leg->AddEntry(hM22, "M22", "lep");
+    leg->AddEntry(hM33, "M33", "lep");
+    leg->AddEntry(hM12, "M12", "lep");
+    leg->AddEntry(hM13, "M13", "lep");
+    leg->AddEntry(hM23, "M23", "lep");
+    leg->Draw();
+    
+    canv->SaveAs(Form("imgs/MNN_Mean_%s.pdf", triggers.Data()));
+
+    
+    
+    canv = new TCanvas(Form("cWidth_%s", triggers.Data()), Form("cWidth_%s", triggers.Data()));
+
+    TH1* hW11 = output.GetHistogram("W_M11", triggerBin);
+    TH1* hW12 = output.GetHistogram("W_M12", triggerBin);
+    TH1* hW13 = output.GetHistogram("W_M13", triggerBin);
+    TH1* hW22 = output.GetHistogram("W_M22", triggerBin);
+    TH1* hW23 = output.GetHistogram("W_M23", triggerBin);
+    TH1* hW33 = output.GetHistogram("W_M33", triggerBin);
+
+
+    hW12->SetMarkerStyle(24);
+    hW12->SetMarkerColor(kCyan+1);
+    hW12->SetLineColor(kCyan+1);
+    hW12->GetYaxis()->SetRangeUser(0., 0.02);
+    hW12->SetTitle("#pi^{0} Peak Width");
+    gStyle->SetOptStat(0);
+    hW12->Draw();
+
+    hW13->SetMarkerStyle(25);
+    hW13->SetMarkerColor(kMagenta);
+    hW13->SetLineColor(kMagenta);
+    hW13->Draw("same");
+
+    hW23->SetMarkerStyle(26);
+    hW23->SetMarkerColor(kYellow+1);
+    hW23->SetLineColor(kYellow+1);
+    hW23->Draw("same");
+    
+
+    hW11->SetMarkerStyle(21);
+    hW11->SetMarkerColor(1);
+    hW11->SetLineColor(1);
+    hW11->GetYaxis()->SetRangeUser(0.12, 0.15);
+    hW11->Draw("same");
+
+    hW22->SetMarkerStyle(22);
+    hW22->SetMarkerColor(kRed);
+    hW22->SetLineColor(kRed);
+    hW22->Draw("same");
+
+    hW33->SetMarkerStyle(23);
+    hW33->SetMarkerColor(kGreen);
+    hW33->SetLineColor(kGreen);
+    hW33->Draw("same");
+
+    leg = new TLegend(0.8,0.6,0.95,0.95);
+    leg->AddEntry(hW11, "W11", "lep");
+    leg->AddEntry(hW22, "W22", "lep");
+    leg->AddEntry(hW33, "W33", "lep");
+    leg->AddEntry(hW12, "W12", "lep");
+    leg->AddEntry(hW13, "W13", "lep");
+    leg->AddEntry(hW23, "W23", "lep");
+    leg->Draw();
+    
+    canv->SaveAs(Form("imgs/MNN_Width_%s.pdf", triggers.Data()));
+  }
+}