2 // *** Class AliRsnPairFunctions ***
4 // "Core" method for defining the work on a pari of particles.
5 // For one analysis, one must setup one of this for each pair he wants to analyze,
6 // adding to it all analysis which he desires to do.
7 // Here he defines the cuts, and the particle types and charges, and can add
8 // functions which do different operations on the same pair, and some binning
9 // with respect to some kinematic variables (eta, momentum)
11 // authors: A. Pulvirenti (email: alberto.pulvirenti@ct.infn.it)
12 // M. Vala (email: martin.vala@cern.ch)
19 #include "AliRsnMother.h"
20 #include "AliRsnEvent.h"
21 #include "AliRsnFunction.h"
22 #include "AliRsnCutSet.h"
23 #include "AliRsnValue.h"
25 #include "AliRsnPairFunctions.h"
27 ClassImp(AliRsnPairFunctions)
29 //_____________________________________________________________________________
30 AliRsnPairFunctions::AliRsnPairFunctions(const char *name, AliRsnPairDef *def) :
31 AliRsnPair(name, def),
32 fFunctions("AliRsnFunction", 0)
35 // Default constructor
38 AliDebug(AliLog::kDebug + 2, "<-");
39 AliDebug(AliLog::kDebug + 2, "->");
42 //_____________________________________________________________________________
43 AliRsnPairFunctions::AliRsnPairFunctions(const AliRsnPairFunctions& copy) :
45 fFunctions(copy.fFunctions)
48 // Default constructor
51 AliDebug(AliLog::kDebug + 2, "<-");
52 AliDebug(AliLog::kDebug + 2, "->");
55 //_____________________________________________________________________________
56 AliRsnPairFunctions& AliRsnPairFunctions::operator=(const AliRsnPairFunctions& copy)
58 AliRsnPair::operator=(copy);
60 Int_t i, n = copy.fFunctions.GetEntries();
61 for (i = 0; i < n; i++) {
62 AliRsnFunction *fcn = (AliRsnFunction*)copy.fFunctions[i];
63 if (fcn) AddFunction(fcn);
69 //_____________________________________________________________________________
70 AliRsnPairFunctions::~AliRsnPairFunctions()
76 AliDebug(AliLog::kDebug + 2, "<-");
77 AliDebug(AliLog::kDebug + 2, "->");
80 //_____________________________________________________________________________
81 void AliRsnPairFunctions::Compute()
84 // Makes computations using the two passed daughter objects.
85 // Checks all cuts and then computes the corresponding pair object
86 // and then fill the list of required values using it.
89 AliDebug(AliLog::kDebug + 2, "<-");
91 TObjArrayIter nextFcn(&fFunctions);
92 AliRsnFunction *fcn = 0x0;
94 while ((fcn = (AliRsnFunction*)nextFcn())) {
95 fcn->SetPairDef(fPairDef);
96 fcn->SetPair(&fMother);
100 AliDebug(AliLog::kDebug + 2, "->");
103 //_____________________________________________________________________________
104 void AliRsnPairFunctions::Init(const char *prefix, TList *list)
107 // Generates needed histograms, giving them a name based on
108 // the flags defined here, on the pair definition, and attaches
109 // a prefix to it, according to the argument.
111 // All generated histograms are stored into the output TList.
114 AliDebug(AliLog::kDebug + 2, "<-");
118 AliRsnFunction *fcn = 0;
119 for (i = 0; i < fFunctions.GetEntries(); i++) {
120 fcn = (AliRsnFunction*)fFunctions.At(i);
125 hName += fcn->GetName();
126 if (fcn->IsUsingTH1()) list->Add(fcn->CreateHistogram(hName.Data(), ""));
127 else list->Add(fcn->CreateHistogramSparse(hName.Data(), ""));
130 AliDebug(AliLog::kDebug + 2, "->");
133 //_____________________________________________________________________________
134 void AliRsnPairFunctions::AddFunction(AliRsnFunction *const fcn)
137 // Adds a new computing function
140 AliDebug(AliLog::kDebug + 2, "<-");
143 Int_t size = fFunctions.GetEntries();
144 new(fFunctions[size]) AliRsnFunction(*fcn);
146 AliDebug(AliLog::kDebug + 2, "->");