4 /* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * See cxx source for full Copyright notice */
7 ///////////////////////////////////////////////////////////////////////////////
9 /// particle id probability densities //
11 ///////////////////////////////////////////////////////////////////////////////
19 class AliPID : public TObject {
22 kSPECIES = 5, // Number of particle species recognized by the PID
23 kSPECIESN = 10 // Number of charged+neutral particle species recognized by the PHOS/EMCAL PID
38 static Float_t ParticleMass(Int_t iType) {
39 if(!fgkParticleMass[0]) Init();
40 return fgkParticleMass[iType];
42 static const char* ParticleName(Int_t iType)
43 {return fgkParticleName[iType];};
44 static const char* ParticleShortName(Int_t iType)
45 {return fgkParticleShortName[iType];};
46 static const char* ParticleLatexName(Int_t iType)
47 {return fgkParticleLatexName[iType];};
48 static Int_t ParticleCode(Int_t iType)
49 {return fgkParticleCode[iType];};
52 AliPID(const Double_t* probDensity, Bool_t charged = kTRUE);
53 AliPID(const Float_t* probDensity, Bool_t charged = kTRUE);
54 AliPID(const AliPID& pid);
55 AliPID& operator = (const AliPID& pid);
57 Double_t GetProbability(EParticleType iType,
58 const Double_t* prior) const;
59 Double_t GetProbability(EParticleType iType) const;
60 void GetProbabilities(Double_t* probabilities,
61 const Double_t* prior) const;
62 void GetProbabilities(Double_t* probabilities) const;
63 EParticleType GetMostProbable(const Double_t* prior) const;
64 EParticleType GetMostProbable() const;
66 void SetProbabilities(const Double_t* probabilities,
67 Bool_t charged = kTRUE);
69 static void SetPriors(const Double_t* prior,
70 Bool_t charged = kTRUE);
71 static void SetPrior(EParticleType iType, Double_t prior);
73 AliPID& operator *= (const AliPID& pid);
79 Bool_t fCharged; // flag for charged/neutral
80 Double_t fProbDensity[kSPECIESN]; // probability densities
81 static Double_t fgPrior[kSPECIESN]; // a priori probabilities
83 static /*const*/ Float_t fgkParticleMass[kSPECIESN+1]; // particle masses
84 static const char* fgkParticleName[kSPECIESN+1]; // particle names
85 static const char* fgkParticleShortName[kSPECIESN+1]; // particle names
86 static const char* fgkParticleLatexName[kSPECIESN+1]; // particle names
87 static const Int_t fgkParticleCode[kSPECIESN+1]; // particle codes
89 ClassDef(AliPID, 2) // particle id probability densities
93 AliPID operator * (const AliPID& pid1, const AliPID& pid2);