/* History of cvs commits:
*
* $Log$
+ * Revision 1.59 2007/03/06 06:57:46 kharlov
+ * DP:calculation of distance to CPV done in TSM
+ *
+ * Revision 1.58 2006/04/12 11:32:03 alibrary
+ * Simplification of Makefile and some small corrections
+ *
+ * Revision 1.57 2006/01/23 17:51:48 hristov
+ * Using the recommended way of forward declarations for TVector and TMatrix (see v5-08-00 release notes). Additional clean-up
+ *
+ * Revision 1.56 2005/05/28 14:19:04 schutz
+ * Compilation warnings fixed by T.P.
+ *
*/
//_________________________________________________________________________
// Oh yeah
//*-- Author: Yves Schutz (SUBATECH), Gustavo Conesa.
+// --- Standard library ---
+
// --- ROOT system ---
-class TVector3 ;
-class TMatrix ;
-class TPrincipal ;
-class TROOT ;
-class TTree ;
class TCanvas ;
class TFolder ;
-class TMatrixD ;
class TFormula;
-// --- Standard library ---
+class TPrincipal ;
+class TROOT ;
+class TTree ;
+class TVector3 ;
+#include <TMatrixDfwd.h>
+
// --- AliRoot header files ---
-class AliPHOSEmcRecPoint ;
-class AliPHOSCpvRecPoint ;
class AliPHOSClusterizerv1 ;
+class AliPHOSCpvRecPoint ;
+class AliPHOSEmcRecPoint ;
+class AliPHOSTrackSegment ;
class AliPHOSTrackSegmentMakerv1 ;
-
#include "AliPHOSPID.h"
#include "AliPID.h"
+
class AliPHOSPIDv1 : public AliPHOSPID {
public:
void SetWriting(Bool_t toWrite = kFALSE){fWrite = toWrite;}
void Print(const Option_t * = "") const ;
+ void GetVertex(void) ; //Extracts vertex in current event
+
virtual const char * Version() const { return "pid-v1" ; }
AliPHOSPIDv1 & operator = (const AliPHOSPIDv1 & /*pid*/) { return *this ;}
Double_t LandauF(Double_t x, Double_t y, Double_t *par) ; //gaussian probability, parameter dependence a+b/(x*x)+c/x
Double_t LandauPol2(Double_t x, Double_t y, Double_t *par) ; //gaussian probability, parameter dependence a+b*x+c*x*x
// Relative Distance CPV-EMC
- Float_t GetDistance (AliPHOSEmcRecPoint * emc, AliPHOSCpvRecPoint * cpv, Option_t * axis)const ;
- Int_t GetCPVBit (AliPHOSEmcRecPoint * emc, AliPHOSCpvRecPoint * cpv, Int_t EffPur, Float_t e) const;
+ Int_t GetCPVBit (AliPHOSTrackSegment * ts, Int_t EffPur, Float_t e) const;
Int_t GetPrincipalBit (TString particle, const Double_t* P, Int_t EffPur, Float_t e)const ; //Principal cut
Int_t GetHardPhotonBit(AliPHOSEmcRecPoint * emc) const;
Int_t GetHardPi0Bit (AliPHOSEmcRecPoint * emc) const;
Double_t *fPPhoton ; //! Principal photon eigenvalues
Double_t *fPPi0 ; //! Principal pi0 eigenvalues
Int_t fRecParticlesInRun ; //! Total number of recparticles in one run
- TMatrix *fParameters; //! Matrix of identification Parameters
+ TMatrixF *fParameters; //! Matrix of identification Parameters
+
+ TVector3 fVtx ; //! Vertex position in current event
//Initial pid population
Double_t fInitPID[AliPID::kSPECIESN] ; // Initial population to do bayesian PID