add new histogram

diff --git a/PWG0/dNdPt/AliPtResolAnalysis.cxx b/PWG0/dNdPt/AliPtResolAnalysis.cxx
index dcac70b56cb3d98abbb47402f527c6f4d75cf91c..df98479d42c6f1c42a03685ac2295ae981aa5e14 100644 (file)
--- a/PWG0/dNdPt/AliPtResolAnalysis.cxx
+++ b/PWG0/dNdPt/AliPtResolAnalysis.cxx
@@ -57,7 +57,8 @@ ClassImp(AliPtResolAnalysis)
 //_____________________________________________________________________________\r
   AliPtResolAnalysis::AliPtResolAnalysis(): AlidNdPt(),\r
   fAnalysisFolder(0),\r
-  fTrackParamHist(0)\r
+  fTrackParamHist(0),\r
+  fTrackParamHist2(0)\r
 {\r
   // default constructor\r
   Init();\r
@@ -66,7 +67,8 @@ ClassImp(AliPtResolAnalysis)
 //_____________________________________________________________________________\r
 AliPtResolAnalysis::AliPtResolAnalysis(Char_t* name, Char_t* title): AlidNdPt(name,title),\r
   fAnalysisFolder(0),\r
-  fTrackParamHist(0)\r
+  fTrackParamHist(0),\r
+  fTrackParamHist2(0)\r
 {\r
   Init();\r
 }\r
@@ -78,6 +80,7 @@ AliPtResolAnalysis::~AliPtResolAnalysis() {
   //\r
   if(fAnalysisFolder) delete fAnalysisFolder; fAnalysisFolder=0;\r
   if(fTrackParamHist) delete fTrackParamHist; fTrackParamHist=0;\r
+  if(fTrackParamHist2) delete fTrackParamHist2; fTrackParamHist2=0;\r
 }\r
 \r
 //_____________________________________________________________________________\r
@@ -92,12 +95,25 @@ void AliPtResolAnalysis::Init(){
 \r
   fTrackParamHist = new THnSparseF("fTrackParamHist","1/pT:#sigma(1/pT)",2,binsTrackParamHist,minTrackParamHist,maxTrackParamHist);\r
   fTrackParamHist->GetAxis(0)->SetTitle("1/pT (GeV/c)^{-1}");\r
-  fTrackParamHist->GetAxis(1)->SetTitle("#sigma(1/pT)" );\r
+  fTrackParamHist->GetAxis(1)->SetTitle("#sigma(1/pT)");\r
   fTrackParamHist->Sumw2();\r
   \r
+  //pt:sigma(1/pT)*pT\r
+  const Int_t ptNbins = 73;\r
+  Double_t bins[74] = {0.0, 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 0.85, 0.9, 0.95, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 1.6, 1.7, 1.8, 1.9, 2.0, 2.2, 2.4, 2.6, 2.8, 3.0, 3.2, 3.4, 3.6, 3.8, 4.0, 4.5, 5.0, 5.5, 6.0, 6.5, 7.0, 8.0, 9.0, 10.0, 11.0, 12.0, 13.0, 14.0, 15.0, 16.0, 18.0, 20.0, 22.0, 24.0, 26.0, 28.0, 30.0, 32.0, 34.0, 36.0, 40.0, 45.0, 50.0, 60.0, 70.0, 80.0, 90.0, 100.0 };\r
+\r
+  Int_t binsTrackParamHist2[2]={ptNbins,200};\r
+  Double_t minTrackParamHist2[2]={0,0}; \r
+  Double_t maxTrackParamHist2[2]={100,0.2};\r
+\r
+  fTrackParamHist2 = new THnSparseF("fTrackParamHist2","pT:#sigma(1/pT)*pT",2,binsTrackParamHist2,minTrackParamHist2,maxTrackParamHist2);\r
+  fTrackParamHist2->SetBinEdges(0,bins);\r
+  fTrackParamHist2->GetAxis(0)->SetTitle("pT (GeV/c)");\r
+  fTrackParamHist2->GetAxis(1)->SetTitle("#sigma(1/pT)*pT");\r
+  fTrackParamHist2->Sumw2();\r
+\r
   // init folder\r
   fAnalysisFolder = CreateFolder("folderdNdPt","Analysis dNdPt Folder");\r
-  \r
 }\r
 \r
 //_____________________________________________________________________________\r
@@ -200,8 +216,11 @@ void AliPtResolAnalysis::Process(AliESDEvent *const esdEvent, AliMCEvent *const
          //track->GetExternalParameters(x, p);\r
          //track->GetExternalCovariance(cov);\r
 \r
-         Double_t v[2] = {track->OneOverPt(),TMath::Sqrt(track->GetSigma1Pt2())};\r
-         fTrackParamHist->Fill(v);\r
+         Double_t v1[2] = {track->OneOverPt(),TMath::Sqrt(track->GetSigma1Pt2())};\r
+         fTrackParamHist->Fill(v1);\r
+\r
+         Double_t v2[2] = {track->Pt(),track->Pt()*TMath::Sqrt(track->GetSigma1Pt2())};\r
+         fTrackParamHist2->Fill(v2);\r
         }\r
       }  \r
     }\r
@@ -235,6 +254,7 @@ Long64_t AliPtResolAnalysis::Merge(TCollection* const list)
     \r
     //\r
     fTrackParamHist->Add(entry->fTrackParamHist);\r
+    fTrackParamHist2->Add(entry->fTrackParamHist2);\r
   }\r
 \r
 return count;\r

diff --git a/PWG0/dNdPt/AliPtResolAnalysis.h b/PWG0/dNdPt/AliPtResolAnalysis.h
index 94181496ef518b05c7f255cb6346c84d013305d6..6dc50eaf02e82b059fc17834d64e9f25cac43dd1 100644 (file)
--- a/PWG0/dNdPt/AliPtResolAnalysis.h
+++ b/PWG0/dNdPt/AliPtResolAnalysis.h
@@ -51,17 +51,19 @@ public :
   // Get analysis folder
   TFolder* GetAnalysisFolder() const {return fAnalysisFolder;}
   THnSparseF *GetTrackParamHist() const {return fTrackParamHist;} 
+  THnSparseF *GetTrackParamHist2() const {return fTrackParamHist2;} 
 
 private:
 
   // analysis folder 
   TFolder *fAnalysisFolder; // folder for analysed histograms
   THnSparseF *fTrackParamHist;  //-> sigma(1/pT):1/pT
+  THnSparseF *fTrackParamHist2;  //-> sigma(1/pT)*pT:pT
 
   AliPtResolAnalysis(const AliPtResolAnalysis&); // not implemented
   AliPtResolAnalysis& operator=(const AliPtResolAnalysis&); // not implemented
 
-  ClassDef(AliPtResolAnalysis,1);
+  ClassDef(AliPtResolAnalysis,2);
 };
 
 #endif