From: kleinb <kleinb@f7af4fe6-9843-0410-8265-dc069ae4e863>
Date: Sun, 15 Mar 2009 15:56:36 +0000 (+0000)
Subject: Take the correct weight for different pT hard bins from the task output
X-Git-Url: http://git.uio.no/git/?p=u%2Fmrichter%2FAliRoot.git;a=commitdiff_plain;h=1e8729f08df5cc5399533d3d9ed426e4d8d16fd2

Take the correct weight for different pT hard bins from the task output
---

diff --git a/PWG4/macros/runJetSpectrumUnfolding.C b/PWG4/macros/runJetSpectrumUnfolding.C
index ebe4b66d185..fbe4d92df8c 100644
--- a/PWG4/macros/runJetSpectrumUnfolding.C
+++ b/PWG4/macros/runJetSpectrumUnfolding.C
@@ -447,19 +447,21 @@ void correct(){
 void mergeJetAnaOutput(){
   // This is used to merge the analysis-output from different 
   // data samples/pt_hard bins
-  // in case the eventweigth was set to xsection/ntrials already this
-  // is not needed. Both methods only work in cse we do not mix different 
+  // in case the eventweigth was set to xsection/ntrials already, this
+  // is not needed. Both methods only work in case we do not mix different 
   // pt_hard bins, and do not have overlapping bins
 
   const Int_t nBins = 2;
   // LHC08q jetjet100: Mean = 1.42483e-03, RMS = 6.642e-05
   // LHC08r jetjet50: Mean = 2.44068e-02, RMS = 1.144e-03
   // LHC08v jetjet15-50: Mean = 2.168291 , RMS = 7.119e-02
-  const Float_t xsection[nBins] = {2.168291,2.44068e-02};
+  // const Float_t xsection[nBins] = {2.168291,2.44068e-02};
+  Float_t xsection[nBins] = {0,0};
   Float_t nTrials[nBins] = {0,0};
   Float_t sf[nBins] = {0,0};
 
-  const char *cFile[nBins] = {"pwg4spec_15-50_all.root","pwg4spec_50_all.root"};
+  const char *cFile[nBins] = {"pwg4spec_0000000-1000000_LHC08v_jetjet15-50.root",
+			      "pwg4spec_0000000-1000000_LHC08r_jetjet50.root"};
 
 
   TList *lIn[nBins];
@@ -468,6 +470,8 @@ void mergeJetAnaOutput(){
     fIn[ib] = TFile::Open(cFile[ib]);
     lIn[ib] = (TList*)fIn[ib]->Get("pwg4spec");
     TH1* hTrials = (TH1F*)lIn[ib]->FindObject("fh1PtHard_Trials");
+    TProfile* hXsec = (TProfile*)lIn[ib]->FindObject("fh1Xsec");
+    xsection[ib] = hXsec->GetBinContent(1);
     nTrials[ib] = hTrials->Integral();
     sf[ib] = xsection[ib]/ nTrials[ib];
   }