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 **************************************************************************/
16 ////////////////////////////////////////////////////////////////////////////////
18 // This class contains all code which is used to compute any of the values
19 // which can be of interest within a resonance analysis. Besides the obvious
20 // invariant mass, it allows to compute other utility values on all possible
21 // targets, in order to allow a wide spectrum of binning and checks.
22 // When needed, this object can also define a binning in the variable which
23 // it is required to compute, which is used for initializing axes of output
24 // histograms (see AliRsnFunction).
25 // The value computation requires this object to be passed the object whose
26 // informations will be used. This object can be of any allowed input type
27 // (track, pair, event), then this class must inherit from AliRsnTarget.
28 // Then, when value computation is attempted, a check on target type is done
29 // and computation is successful only if expected target matches that of the
31 // In some cases, the value computation can require a support external object,
32 // which must then be passed to this class. It can be of any type inheriting
35 // authors: A. Pulvirenti (alberto.pulvirenti@ct.infn.it)
36 // M. Vala (martin.vala@cern.ch)
38 ////////////////////////////////////////////////////////////////////////////////
40 #include "AliESDtrackCuts.h"
41 #include "AliESDpid.h"
42 #include "AliAODPid.h"
43 #include "AliCentrality.h"
45 #include "AliRsnEvent.h"
46 #include "AliRsnDaughter.h"
47 #include "AliRsnMother.h"
48 #include "AliRsnPairDef.h"
49 #include "AliRsnDaughterDef.h"
51 #include "AliRsnValue.h"
55 //_____________________________________________________________________________
56 AliRsnValue::AliRsnValue(const char *name, AliRsnTarget::ETargetType type) :
57 AliRsnTarget(name, type),
64 // Initializes the binning to an empty array.
68 //_____________________________________________________________________________
69 AliRsnValue::AliRsnValue(const AliRsnValue ©) :
71 fUseMCInfo(copy.fUseMCInfo),
72 fComputedValue(copy.fComputedValue),
73 fBinArray(copy.fBinArray)
77 // Duplicates the binning array and copies all settings.
81 //_____________________________________________________________________________
82 AliRsnValue &AliRsnValue::operator=(const AliRsnValue ©)
85 // Assignment operator.
86 // Works like copy constructor.
89 AliRsnTarget::operator=(copy);
92 fUseMCInfo = copy.fUseMCInfo;
93 fComputedValue = copy.fComputedValue;
94 fBinArray = copy.fBinArray;
99 //_____________________________________________________________________________
100 void AliRsnValue::SetBins(Int_t nbins, Double_t min, Double_t max)
103 // Set binning for the axis in equally spaced bins
104 // where the number of bins, minimum and maximum are given.
112 fBinArray.Set(nbins + 1);
114 Double_t mymax = TMath::Max(min, max);
115 Double_t mymin = TMath::Min(min, max);
118 Double_t binSize = (mymax - mymin) / ((Double_t)nbins);
120 fBinArray[0] = mymin;
121 for (k = 1; k <= nbins; k++) fBinArray[k] = fBinArray[k - 1] + binSize;
124 //_____________________________________________________________________________
125 void AliRsnValue::SetBins(Double_t min, Double_t max, Double_t step)
128 // Set binning for the axis in equally spaced bins
129 // where the bin size, minimum and maximum are given.
132 Double_t dblNbins = TMath::Abs(max - min) / step;
133 Int_t intNbins = (Int_t) (dblNbins + 0.5);
135 SetBins(intNbins, min, max);
138 //_____________________________________________________________________________
139 void AliRsnValue::SetBins(Int_t nbins, Double_t *array)
142 // Set binning for the axis in unequally spaced bins
143 // using the same way it is done in TAxis
152 fBinArray.Set(nbins);
153 for (i = 0; i < nbins; i++) fBinArray[i] = array[i];
156 //_____________________________________________________________________________
157 Bool_t AliRsnValue::Eval(TObject *)
160 // Evaluation of the required value.
161 // Checks that the passed object is of the right type
162 // and if this check is successful, computes the required value.
163 // The output of the function tells if computing was successful,
164 // and the values must be taken with GetValue().
167 AliWarning("This method must be overridden by derived classes");
171 //_____________________________________________________________________________
172 void AliRsnValue::Print(Option_t *option) const
175 // Print informations about this object.
176 // If one specifies option "BINS" all bin limits are also printed.
179 AliInfo("=== VALUE INFO =================================================");
180 AliInfo(Form(" Name : %s", GetName()));
181 AliInfo(Form(" Current computed value: %f", fComputedValue));
182 if (!strcmp(option, "BINS")) {
184 for (i = 0; i < fBinArray.GetSize(); i++) {
185 AliInfo(Form(" Bin limit #%03d = %f", i, fBinArray[i]));