removed obsolete implementation files from package
[u/mrichter/AliRoot.git] / PWG2 / RESONANCES / AliRsnCutBetheBloch.cxx
1 //
2 // Class AliRsnCutBetheBloch
3 //
4 // General implementation of a single cut strategy, which can be:
5 // - a value contained in a given interval  [--> IsBetween()   ]
6 // - a value equal to a given reference     [--> MatchesValue()]
7 //
8 // In all cases, the reference value(s) is (are) given as data members
9 // and each kind of cut requires a given value type (Int, UInt, Double),
10 // but the cut check procedure is then automatized and chosen thanks to
11 // an enumeration of the implemented cut types.
12 // At the end, the user (or any other point which uses this object) has
13 // to use the method IsSelected() to check if this cut has been passed.
14 //
15 // authors: Martin Vala (martin.vala@cern.ch)
16 //          Alberto Pulvirenti (alberto.pulvirenti@ct.infn.it)
17 //
18 #include "TMath.h"
19
20 #include "AliLog.h"
21 #include "AliESDtrack.h"
22
23 #include "AliRsnDaughter.h"
24 #include "AliRsnCutBetheBloch.h"
25
26 ClassImp(AliRsnCutBetheBloch)
27
28 //_________________________________________________________________________________________________
29 AliRsnCutBetheBloch::AliRsnCutBetheBloch() :
30   AliRsnCut(),
31   fCorrect(kTRUE),
32   fMass(1.0),
33   fMIP(50.0)
34 {
35 //
36 // Default constructor.
37 //
38
39   fConst[0] = fConst[1] = fConst[2] = fConst[3] = fConst[4] = 0.0;
40 }
41
42 //_________________________________________________________________________________________________
43 AliRsnCutBetheBloch::AliRsnCutBetheBloch
44 (const char *name, Double_t fractionRange, Double_t mass, Double_t mip, Bool_t correct) :
45   AliRsnCut(name, -fractionRange, fractionRange),
46   fCorrect(correct),
47   fMass(mass),
48   fMIP(mip)
49 {
50 //
51 // Main constructor.
52 // the cut range is the relative fraction of the value:
53 // BB*(1-fraction) < TPC < BB*(1+fraction)
54 // which means:
55 // -fraction < (TPC - BB)/BB < fraction
56 //
57
58   fConst[0] = fConst[1] = fConst[2] = fConst[3] = fConst[4] = 0.0;
59 }
60
61 //_________________________________________________________________________________________________
62 AliRsnCutBetheBloch::AliRsnCutBetheBloch
63 (const char *name, Double_t fractionRange, AliPID::EParticleType type, Double_t mip, Bool_t correct) :
64   AliRsnCut(name, -fractionRange, fractionRange),
65   fCorrect(correct),
66   fMass(0.0),
67   fMIP(mip)
68 {
69 //
70 // Main constructor.
71 // the cut range is the relative fraction of the value:
72 // BB*(1-fraction) < TPC < BB*(1+fraction)
73 // which means:
74 // -fraction < (TPC - BB)/BB < fraction
75 //
76
77   SetMass(type);
78   fConst[0] = fConst[1] = fConst[2] = fConst[3] = fConst[4] = 0.0;
79 }
80
81 //_____________________________________________________________________________
82 Double_t AliRsnCutBetheBloch::BetheBloch(AliRsnDaughter *track)
83 {
84 //
85 // Computes the theoretical dE/dx according to
86 // a given mass hypothesis, from which betaGamma is computed
87 //
88 // This is the empirical ALEPH parameterization of the Bethe-Bloch formula.
89 // It is normalized to 1 at the minimum.
90 //
91 // The default values for the kp* parameters are for ALICE TPC.
92 // The value is computed in MIP units, multiplied by 50 to have it in energy.
93 //
94
95   Double_t betaGamma = track->P() / fMass;
96   Double_t beta = betaGamma / TMath::Sqrt(1.0 + betaGamma * betaGamma);
97   Double_t aa = TMath::Power(beta, fConst[3]);
98   Double_t bb = TMath::Power(1.0 / betaGamma, fConst[4]);
99
100   bb = TMath::Log(fConst[2] + bb);
101
102   Double_t out = (fConst[1] - aa - bb) * fConst[0] / aa;
103
104   if (fCorrect)
105   {
106     Double_t kMeanCorr = 0.1;
107     Double_t meanCorr = (1 + (out - 1) * kMeanCorr);
108     out *= meanCorr;
109   }
110
111   return out * fMIP;
112 }
113
114 //_________________________________________________________________________________________________
115 Bool_t AliRsnCutBetheBloch::IsSelected(ETarget tgt, AliRsnDaughter *track)
116 {
117 //
118 // Cut checker.
119 //
120
121   // coherence check
122   if (tgt != AliRsnCut::kParticle)
123   {
124     AliError(Form("Wrong target. Skipping cut", GetName()));
125     return kTRUE;
126   }
127
128   // retrieve the TPC signal
129   AliVParticle *vpart = track->GetRef();
130   AliESDtrack *esd = dynamic_cast<AliESDtrack*>(vpart);
131   if (!esd) {
132     AliError("ESD information unavailable");
133     return kTRUE;
134   }
135
136   // compute Bethe-Bloch with the given mass hypothesis
137   Double_t bb = BetheBloch(track);
138
139   // the cut range is the relative fraction of the value:
140   // BB*(1-fraction) < TPC < BB*(1+fraction)
141   // which means:
142   // -fraction < (TPC - BB)/BB < fraction
143   // so we must compute the cut value accordingly
144   fCutValueD = (esd->GetTPCsignal() - bb) / bb;
145
146   // then, this cut is checked inside the range
147   return OkRange();
148 }
149
150 //_________________________________________________________________________________________________
151 Bool_t AliRsnCutBetheBloch::IsSelected(ETarget /*tgt*/, AliRsnPairParticle* /*pair*/)
152 {
153 //
154 // Cut checker
155 //
156
157   AliWarning("Cannot apply this cut to pairs");
158   return kTRUE;
159 }
160
161 //_________________________________________________________________________________________________
162 Bool_t AliRsnCutBetheBloch::IsSelected(ETarget /*tgt*/, AliRsnEvent* /*event*/)
163 {
164 //
165 // Cut checker
166 //
167
168   AliWarning("Cannot apply this cut to events");
169   return kTRUE;
170 }
171
172 //_________________________________________________________________________________________________
173 Bool_t AliRsnCutBetheBloch::IsSelected(ETarget /*tgt*/, AliRsnEvent* /*ev1*/, AliRsnEvent* /*ev2*/)
174 {
175 //
176 // Cut checker
177 //
178
179   AliWarning("Cannot apply this cut to event mixing");
180   return kTRUE;
181 }