X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;ds=sidebyside;f=EMCAL%2FAliEMCALPID.h;h=d984c4d46aee0c7879516e6b894e8c8acd96d57d;hb=79263b42879e5ceaaf382be562efc25dbaa257a2;hp=6e29f95d1d130e19af1f8b8aef304182f467a307;hpb=af885e0fc45eeab1177d3d34f858ccd8836787fb;p=u%2Fmrichter%2FAliRoot.git

diff --git a/EMCAL/AliEMCALPID.h b/EMCAL/AliEMCALPID.h
index 6e29f95d1d1..d984c4d46ae 100644
--- a/EMCAL/AliEMCALPID.h
+++ b/EMCAL/AliEMCALPID.h
@@ -1,69 +1,42 @@
-#ifndef AliEMCALPID_H
-#define AliEMCALPID_H
+#ifndef ALIEMCALPID_H
+#define ALIEMCALPID_H
 
 /* $Id$ */
-/* History of cvs commits:
- *
- * $Log$
- * Revision 1.12  2007/02/20 20:17:43  hristov
- * Corrected array size, removed warnings (icc)
- *
- * Revision 1.11  2006/12/19 08:49:35  gustavo
- * New PID class for EMCAL, bayesian analysis done with ESD data, PID information filled when calling AliEMCALPID in AliEMCALReconstructor::FillESD()
- *
- *
- */
 
 ///////////////////////////////////////////////////////////////////////////////
 // Class AliEMCALPID
+// Compute PID weights for all the clusters
 ///////////////////////////////////////////////////////////////////////////////
 
+//Root includes
 #include "TTask.h"
-#include "TArrayD.h"
-#include "AliESDEvent.h"
-#include "AliPID.h" 
+class TArrayD ;
 
-class AliEMCALPID : public TTask {
+//AliRoot includes
+class AliESDEvent ;
+#include "AliEMCALPIDUtils.h" 
+
+class AliEMCALPID : public AliEMCALPIDUtils {
 
 public:
   
   AliEMCALPID();
-  virtual ~AliEMCALPID() { }
-  
-  void     RunPID(AliESDEvent *esd);
-  void     ComputePID(Double_t energy, Double_t lambda0); // give the PID of a cluster
-  TArrayD  DistLambda0(Double_t energy, Int_t nature); // compute lambda0 distributions
-  
-  Double_t GetPID(Int_t idx) const {if (idx>=0&&idx<3) return fPID[idx]; else return 0.;}
-  Double_t GetPIDFinal(Int_t idx) const {if (idx>=0&&idx<AliPID::kSPECIESN) return fPIDFinal[idx]; else return 0.;}
-  Double_t GetPIDWeight(Int_t idx) const {if (idx>=0&&idx<3) return fPIDWeight[idx]; else return 0.;}
+  AliEMCALPID(Bool_t reconstructor);
+  //virtual ~AliEMCALPID() { }
   
-  void     SetPID(Double_t val, Int_t idx) {if (idx>=0&&idx<3) fPID[idx] = val;}
-  void     SetPIDFinal(Double_t val, Int_t idx) {if (idx>=0&&idx<AliPID::kSPECIESN) fPIDFinal[idx] = val;}
-  void     SetPIDWeight(Double_t val, Int_t idx) {if (idx>=0&&idx<3) fPIDWeight[idx] = val;}
-  void     SetPrintInfo(Bool_t yesno) {fPrintInfo = yesno;}
-   void     SetReconstructor(Bool_t yesno) {fReconstructor = yesno;}
+  void    RunPID(AliESDEvent *esd);
+  void    InitParameters();
+  void    SetReconstructor(Bool_t yesno) {fReconstructor = yesno;}
+	
  private:
   
-  Double_t Polynomial(Double_t x, Double_t *params);
-  
-  Bool_t   fPrintInfo;          // flag to decide if details about PID must be printed
-  
-  Double_t fGamma[6][6];        // Parameter to Compute PID
-  Double_t fHadron[6][6];		  // Parameter to Compute PID
-  Double_t fPiZero5to10[6][6];  // Parameter to Compute PID
-  Double_t fPiZero10to60[6][6]; // Parameter to Compute PID
-  
-  Float_t fPID[3];
+  Bool_t   fReconstructor;                // Fill esdcalocluster when called from EMCALReconstructor
   
-  Float_t fPIDFinal[AliPID::kSPECIESN+1];  // final PID format
-  Float_t fPIDWeight[3];                 // order: gamma, pi0, hadrons,
-  Double_t fProbGamma;	                // probility to be a Gamma
-  Double_t fProbPiZero;	                // probility to be a PiO
-  Double_t fProbHadron;	                // probility to be a Hadron
-  Bool_t    fReconstructor;               //Fill esdcalocluster when called from EMCALReconstructor
-  
-  ClassDef(AliEMCALPID, 0)
+  ClassDef(AliEMCALPID, 5)
+
 };
 
+
 #endif // ALIEMCALPID_H
+
+