From e16bcc1febb1d2359ebfe83bac3ef982d1f2f09c Mon Sep 17 00:00:00 2001
From: hdalsgaa <hdalsgaa@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Wed, 15 Sep 2010 12:55:00 +0000
Subject: [PATCH] script to draw dN/deta

---
 PWG2/FORWARD/analysis/drawdNdeta.C | 73 ++++++++++++++++++++++++++++++
 1 file changed, 73 insertions(+)
 create mode 100644 PWG2/FORWARD/analysis/drawdNdeta.C

diff --git a/PWG2/FORWARD/analysis/drawdNdeta.C b/PWG2/FORWARD/analysis/drawdNdeta.C
new file mode 100644
index 00000000000..89211517764
--- /dev/null
+++ b/PWG2/FORWARD/analysis/drawdNdeta.C
@@ -0,0 +1,73 @@
+void drawdNdeta(const Char_t* filename="fmdana.root", 
+		Int_t sample = 0, 
+		Int_t rebin =1,
+		Float_t v1 = -10, 
+		Float_t v2 = 10, 
+		Bool_t realdata=kTRUE, 
+		Float_t energy = 900, 
+		Int_t magfield=0){
+  
+  gSystem->Load("libANALYSIS"); 
+  gSystem->Load("libANALYSISalice"); 
+  gSystem->Load("libPWG0base"); 
+  gSystem->Load("libPWG0dep"); 
+  gSystem->Load("libPWG2forward"); 
+  gStyle->SetTextFont(132);
+  gStyle->SetLabelFont(132,"X");
+  gStyle->SetLabelFont(132,"Y"); 
+  gStyle->SetLabelFont(132,"Z"); 
+  gStyle->SetTitleFont(132,"X");
+  gStyle->SetTitleFont(132,"Y");
+  gStyle->SetTitleFont(132,"Z");
+  
+  AliFMDAnaParameters* pars = AliFMDAnaParameters::Instance();
+  if(energy == 900)
+    pars->SetEnergy(AliFMDAnaParameters::k900);
+  else if(energy == 10000)
+    pars->SetEnergy(AliFMDAnaParameters::k10000);
+  else if(energy == 7000)
+    pars->SetEnergy(AliFMDAnaParameters::k7000);
+  else {
+    std::cout<<"invalid energy - quitting"<<std::endl;
+    return; }
+  
+  if(magfield == 0)
+    pars->SetMagField(AliFMDAnaParameters::k0G); 
+  else if(magfield == 1)
+    pars->SetMagField(AliFMDAnaParameters::k5G); 
+  else
+    return;
+  pars->SetRealData(realdata);
+  pars->Init();
+  pars->Print();
+  AliFMDDndeta t;
+  t.SetNbinsToCut(2);
+  t.SetVtxCut1(v1);
+  t.SetVtxCut2(v2);
+  //  t.SetVtxCut(6);
+  t.Init(filename);
+  if(sample == 0) {
+    t.GenerateMult(AliFMDDndeta::kMult);
+    t.DrawDndeta(AliFMDDndeta::kMult,rebin,realdata);
+  }
+  if(sample == 1) {
+    t.GenerateMult(AliFMDDndeta::kMultTrVtx);
+    t.DrawDndeta(AliFMDDndeta::kMultTrVtx,rebin,realdata);
+  }
+  if(sample == 2) {
+    t.GenerateMult(AliFMDDndeta::kHits);
+    t.DrawDndeta(AliFMDDndeta::kHits,rebin);
+  }
+  if(sample == 3) {
+    t.GenerateMult(AliFMDDndeta::kHitsTrVtx);
+    t.DrawDndeta(AliFMDDndeta::kHitsTrVtx,rebin);
+  }
+  if(sample == 4) {
+    t.GenerateMult(AliFMDDndeta::kMultNSD);
+    t.DrawDndeta(AliFMDDndeta::kMultNSD,rebin,realdata);
+  }
+  
+  
+  if(sample == 10) 
+    t.CreateSharingEfficiency(filename,kTRUE);
+}
-- 
2.39.3