add task to manager

[u/mrichter/AliRoot.git] / PWGJE / EMCALJetTasks / AliLocalRhoParameter.h

diff --git a/PWGJE/EMCALJetTasks/AliLocalRhoParameter.h b/PWGJE/EMCALJetTasks/AliLocalRhoParameter.h
index 54a120dfe029dd52cc9a062b56025b66aca9095f..f5d925bbb38b3198c2ff7cb16e81c2a514b3a255 100644 (file)
--- a/PWGJE/EMCALJetTasks/AliLocalRhoParameter.h
+++ b/PWGJE/EMCALJetTasks/AliLocalRhoParameter.h
@@ -18,6 +18,18 @@ class AliLocalRhoParameter : public AliRhoParameter {
     Double_t denom(2*r*fLocalRho->GetParameter(0));
     return  (denom <= 0.) ? GetVal() : n*(fLocalRho->Integral(phi-r, phi+r)/denom); 
   }
+  Double_t GetLocalVal(Double_t phi, Double_t r) const {
+    return GetLocalVal(phi, r, GetVal());
+  }
+  Double_t GetLocalUncertainty(Double_t phi, Double_t r, Double_t n) const {
+      if(!fLocalRho) return 999.;
+      Double_t intError(fLocalRho->IntegralError(phi-r,phi+r));
+      Double_t absConst(TMath::Abs(n/(2*r*fLocalRho->GetParameter(0))));
+      return  intError*absConst;        // absolute error on local rho
+  }
+  Double_t GetLocalUncertainty(Double_t phi, Double_t r) const {
+      return GetLocalUncertainty(phi, r, GetVal());
+  }
  private:
   TF1*     fLocalRho;      // ! rho as function of phi