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
24 kSPECIESLN = 4 // Number of light nuclei: deuteron, triton, helium-3 and alpha
43 static Float_t ParticleMass(Int_t iType) {
44 if(!fgkParticleMass[0]) Init();
45 return fgkParticleMass[iType];
47 static const char* ParticleName(Int_t iType)
48 {return fgkParticleName[iType];};
49 static const char* ParticleShortName(Int_t iType)
50 {return fgkParticleShortName[iType];};
51 static const char* ParticleLatexName(Int_t iType)
52 {return fgkParticleLatexName[iType];};
53 static Int_t ParticleCode(Int_t iType)
54 {return fgkParticleCode[iType];};
57 AliPID(const Double_t* probDensity, Bool_t charged = kTRUE);
58 AliPID(const Float_t* probDensity, Bool_t charged = kTRUE);
59 AliPID(const AliPID& pid);
60 AliPID& operator = (const AliPID& pid);
62 Double_t GetProbability(EParticleType iType,
63 const Double_t* prior) const;
64 Double_t GetProbability(EParticleType iType) const;
65 void GetProbabilities(Double_t* probabilities,
66 const Double_t* prior) const;
67 void GetProbabilities(Double_t* probabilities) const;
68 EParticleType GetMostProbable(const Double_t* prior) const;
69 EParticleType GetMostProbable() const;
71 void SetProbabilities(const Double_t* probabilities,
72 Bool_t charged = kTRUE);
74 static void SetPriors(const Double_t* prior,
75 Bool_t charged = kTRUE);
76 static void SetPrior(EParticleType iType, Double_t prior);
78 AliPID& operator *= (const AliPID& pid);
84 Bool_t fCharged; // flag for charged/neutral
85 Double_t fProbDensity[kSPECIESN]; // probability densities
86 static Double_t fgPrior[kSPECIESN]; // a priori probabilities
88 static /*const*/ Float_t fgkParticleMass[kSPECIESN+kSPECIESLN+1]; // particle masses
89 static const char* fgkParticleName[kSPECIESN+kSPECIESLN+1]; // particle names
90 static const char* fgkParticleShortName[kSPECIESN+kSPECIESLN+1]; // particle names
91 static const char* fgkParticleLatexName[kSPECIESN+kSPECIESLN+1]; // particle names
92 static const Int_t fgkParticleCode[kSPECIESN+kSPECIESLN+1]; // particle codes
94 ClassDef(AliPID, 2) // particle id probability densities
98 AliPID operator * (const AliPID& pid1, const AliPID& pid2);