1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
15 #ifndef ALIHFEPIDTRD_H
16 #define ALIHFEPIDTRD_H
18 #ifndef ALIHFEPIDBASE_H
19 #include "AliHFEpidBase.h"
23 class AliAODMCParticle;
30 class AliHFEpidTRD : public AliHFEpidBase{
40 kHistTRDlikeBefore = 0,
41 kHistTRDlikeAfter = 1,
42 kHistTRDthresholds = 2,
45 kHistOverallSpecies = 5
47 AliHFEpidTRD(const Char_t *name);
48 AliHFEpidTRD(const AliHFEpidTRD &ref);
49 AliHFEpidTRD& operator=(const AliHFEpidTRD &ref);
50 virtual ~AliHFEpidTRD();
52 virtual Bool_t InitializePID();
53 virtual Int_t IsSelected(AliHFEpidObject *track);
54 virtual Bool_t HasQAhistos() const { return kTRUE; };
56 Double_t GetTRDSignalV1(AliESDtrack *track, Int_t mcPID);
57 Double_t GetTRDSignalV2(AliESDtrack *track, Int_t mcPID);
59 void SetPIDMethod(PIDMethodTRD_t method) { fPIDMethod = method; };
61 void Copy(TObject &ref) const;
62 Int_t MakePIDesd(AliESDtrack *esdTrack, AliMCParticle *mcTrack);
63 Int_t MakePIDaod(AliAODTrack *aofTrack, AliAODMCParticle *aodMC);
64 Double_t GetTRDthresholds(Double_t electronEff, Double_t p);
65 Int_t GetMCpid(AliESDtrack *track);
66 void InitParameters();
67 virtual void AddQAhistograms(TList *l);
68 void GetParameters(Double_t electronEff, Double_t *parameters);
70 void FillHistogramsLikelihood(Int_t whenFilled, Float_t p, Float_t elProb);
71 void FillHistogramsTRDSignalV1(Double_t signal, Double_t p, Int_t species);
72 void FillHistogramsTRDSignalV2(Double_t signal, Double_t p, Int_t species);
74 static const Double_t fgkVerySmall; // Check for 0
75 PIDMethodTRD_t fPIDMethod; // PID Method: 2D Likelihood or Neural Network
76 Double_t fThreshParams[kThreshParams]; // Threshold parametrisation
77 TList *fContainer; // QA Histogram Container
78 ClassDef(AliHFEpidTRD, 1) // TRD electron ID class