2 /**************************************************************************
3 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
5 * Author: The ALICE Off-line Project. *
6 * Contributors are mentioned in the code where appropriate. *
8 * Permission to use, copy, modify and distribute this software and its *
9 * documentation strictly for non-commercial purposes is hereby granted *
10 * without fee, provided that the above copyright notice appears in all *
11 * copies and that both the copyright notice and this permission notice *
12 * appear in the supporting documentation. The authors make no claims *
13 * about the suitability of this software for any purpose. It is *
14 * provided "as is" without express or implied warranty. *
15 **************************************************************************/
16 //--------------------------------------------------------------------//
18 // AliCFContainer Class //
19 // Class to accumulate data on an N-dimensional grids, at different //
20 // selection stages. To be used as an input to get corrections for //
21 // Reconstruction & Trigger efficiency //
23 // -- Author : S.Arcelli //
24 //--------------------------------------------------------------------//
28 #include "AliCFGridSparse.h"
29 #include "AliCFContainer.h"
31 //____________________________________________________________________
32 ClassImp(AliCFContainer)
34 //____________________________________________________________________
35 AliCFContainer::AliCFContainer() :
41 // default constructor
45 //____________________________________________________________________
46 AliCFContainer::AliCFContainer(const Char_t* name, const Char_t* title, const Int_t nSelSteps, const Int_t nVarIn, const Int_t* nBinIn) :
47 AliCFFrame(name,title),
56 fGrid = new AliCFGridSparse*[fNStep]; //the grids at the various selection steps
57 for (Int_t istep=0; istep<fNStep; istep++) {
58 fGrid[istep] = new AliCFGridSparse(Form("%s_SelStep%d",name,istep),Form("step%d",istep),nVarIn,nBinIn);
59 fGrid[istep]->SumW2();
61 for (Int_t iVar=0; iVar<nVarIn; iVar++) SetVarTitle(iVar,Form("var%d",iVar));
62 AliInfo(Form("Grids created for %d steps required \n => Don't forget to set the bin limits !!",fNStep));
64 //____________________________________________________________________
65 AliCFContainer::AliCFContainer(const AliCFContainer& c) :
66 AliCFFrame(c.fName,c.fTitle),
75 //____________________________________________________________________
76 AliCFContainer::~AliCFContainer()
82 for ( Int_t istep=0; istep<fNStep; istep++ )
87 //____________________________________________________________________
88 AliCFContainer &AliCFContainer::operator=(const AliCFContainer &c)
93 if (this != &c) c.Copy(*this);
97 //____________________________________________________________________
98 void AliCFContainer::Copy(TObject& c) const
104 AliCFContainer& target = (AliCFContainer &) c;
105 target.fNStep = fNStep;
106 target.fGrid = new AliCFGridSparse*[fNStep];
107 for (Int_t iStep=0; iStep<fNStep; iStep++) {
108 if (fGrid[iStep]) target.fGrid[iStep] = new AliCFGridSparse(*(fGrid[iStep]));
112 //____________________________________________________________________
113 void AliCFContainer::Fill(const Double_t *var, Int_t istep, Double_t weight)
116 // Fills the grid at selection step istep for a set of values of the
117 // input variables, with a given weight (by default w=1)
119 if(istep >= fNStep || istep < 0){
120 AliError("Non-existent selection step, grid was not filled");
123 fGrid[istep]->Fill(var,weight);
126 //____________________________________________________________________
127 TH1* AliCFContainer::Project(Int_t istep, Int_t ivar1, Int_t ivar2, Int_t ivar3) const
130 // returns a projection along variables ivar1 (and ivar2 (and ivar3))
131 // at selection step istep
133 if (istep >= fNStep || istep < 0){
134 AliError("Non-existent selection step, return NULL");
137 return fGrid[istep]->Project(ivar1,ivar2,ivar3);
140 //____________________________________________________________________
141 AliCFContainer* AliCFContainer::MakeSlice(Int_t nVars, const Int_t* vars, const Double_t* varMin, const Double_t* varMax, Bool_t useBins) const
144 // Makes a slice along the "nVars" variables defined in the array "vars[nVars]" for all the container steps.
145 // The ranges of ALL the container variables must be defined in the array varMin[GetNVar()] and varMax[GetNVar()]
146 // The function returns a new container of nVars variables.
147 // If useBins=true, varMin and varMax are taken as bin numbers
149 Int_t* steps = new Int_t[fNStep];
150 for (Int_t iStep=0;iStep<fNStep;iStep++) steps[iStep]=iStep;
151 AliCFContainer* out = MakeSlice(fNStep,steps,nVars,vars,varMin,varMax,useBins);
156 //____________________________________________________________________
157 AliCFContainer* AliCFContainer::MakeSlice(Int_t nSteps, const Int_t* steps,
158 Int_t nVars, const Int_t* vars,
159 const Double_t* varMin, const Double_t* varMax, Bool_t useBins) const
162 // Makes a slice along the "nVars" variables defined in the array "vars[nVars]" for the given "nSteps" defined in "steps[nSteps]".
163 // The ranges of ALL the container variables must be defined in the array varMin[GetNVar()] and varMax[GetNVar()]
164 // The function returns a new container of nVars variables.
165 // If useBins=true, varMin and varMax are taken as bin numbers
168 if (nVars < 1 || nVars > GetNVar()) AliError("Bad number of dimensions required for the slice");
169 if (nSteps< 1 || nSteps> fNStep) AliError("Bad number of steps required for the slice");
171 AliInfo(Form("Making a slice in %d dimension(s)",nVars));
173 // create the output grids
174 AliCFGridSparse** grids = new AliCFGridSparse*[nSteps] ;
175 for (Int_t iStep=0; iStep<nSteps; iStep++) grids[iStep] = fGrid[steps[iStep]]->MakeSlice(nVars,vars,varMin,varMax,useBins);
177 TAxis ** axis = new TAxis*[nVars];
178 for (Int_t iVar=0; iVar<nVars; iVar++) axis[iVar] = ((AliCFGridSparse*)grids[0])->GetGrid()->GetAxis(iVar); //same axis for every grid
180 //define new binning for new container
181 Int_t* bins=new Int_t[nVars];
182 for (Int_t iVar=0; iVar<nVars; iVar++) bins[iVar] = axis[iVar]->GetNbins();
184 AliCFContainer* out = new AliCFContainer(fName,fTitle,nSteps,nVars,bins);
187 for (Int_t iVar=0; iVar<nVars; iVar++) {
188 Int_t nBins = bins[iVar];
189 Double_t *array = new Double_t[nBins+1];
190 for (Int_t iBin=1; iBin<=nBins+1; iBin++) {
191 array[iBin-1] = axis[iVar]->GetBinLowEdge(iBin);
193 out->SetBinLimits(iVar,array);
197 //set grid for the given steps
198 for (Int_t iStep=0; iStep<nSteps; iStep++) out->SetGrid(iStep,grids[iStep]);
205 //____________________________________________________________________
206 Long64_t AliCFContainer::Merge(TCollection* list)
208 // Merge a list of AliCorrection objects with this (needed for
210 // Returns the number of merged objects (including this).
222 while ((obj = iter())) {
223 AliCFContainer* entry = dynamic_cast<AliCFContainer*> (obj);
233 //____________________________________________________________________
234 void AliCFContainer::Add(const AliCFContainer* aContainerToAdd, Double_t c)
237 //add the content of container aContainerToAdd to the current one
239 if ((aContainerToAdd->GetNStep() != fNStep) ||
240 (aContainerToAdd->GetNVar() != GetNVar()) ||
241 (aContainerToAdd->GetNBinsTotal() != GetNBinsTotal()))
243 AliError("Different number of steps/sensitive variables/grid elements: cannot add the containers");
246 for (Int_t istep=0; istep<fNStep; istep++) {
247 fGrid[istep]->Add(aContainerToAdd->GetGrid(istep),c);
250 //____________________________________________________________________
251 Float_t AliCFContainer::GetOverFlows( Int_t ivar, Int_t istep, Bool_t exclusive) const {
253 // Get overflows in variable var at selection level istep
254 // Set 'exclusive' to true for an exclusive check on variable ivar
256 if(istep >= fNStep || istep < 0){
257 AliError("Non-existent selection step, return -1");
260 return fGrid[istep]->GetOverFlows(ivar,exclusive);
262 //____________________________________________________________________
263 Float_t AliCFContainer::GetUnderFlows( Int_t ivar, Int_t istep, Bool_t exclusive) const {
265 // Get underflows in variable var at selection level istep
266 // Set 'exclusive' to true for an exclusive check on variable ivar
268 if(istep >= fNStep || istep < 0){
269 AliError("Non-existent selection step, return -1");
272 return fGrid[istep]->GetUnderFlows(ivar,exclusive);
274 //____________________________________________________________________
275 Float_t AliCFContainer::GetEntries(Int_t istep) const {
277 // Get total entries in variable var at selection level istep
279 if(istep >= fNStep || istep < 0){
280 AliError("Non-existent selection step, return -1");
283 return fGrid[istep]->GetEntries();
285 //_____________________________________________________________________
286 Double_t AliCFContainer::GetIntegral( Int_t istep) const
289 // Get Integral over the grid at selection level istep
291 if(istep >= fNStep || istep < 0){
292 AliError("Non-existent selection step, return -1");
295 return fGrid[istep]->GetIntegral();
298 //_____________________________________________________________________
299 void AliCFContainer::SetRangeUser(Int_t ivar, Double_t varMin, Double_t varMax, Bool_t useBins) const
302 // set axis range for variable ivar
303 // put useBins=kTRUE if you want to pass bin numbers instead of values
305 if (ivar >= GetNVar() || ivar < 0){
306 AliError("Non-existent selection var");
309 for (Int_t iStep=0; iStep<GetNStep(); iStep++) fGrid[iStep]->SetRangeUser(ivar,varMin,varMax,useBins);
312 //_____________________________________________________________________
313 void AliCFContainer::SetRangeUser(const Double_t* varMin, const Double_t* varMax, Bool_t useBins) const
316 // set all axis ranges according to arrays varMin and varMax
317 // put useBins=kTRUE if you want to pass bin numbers instead of values
319 for (Int_t iStep=0; iStep<GetNStep(); iStep++) fGrid[iStep]->SetRangeUser(varMin,varMax,useBins);
322 //_____________________________________________________________________
323 void AliCFContainer::Print(const Option_t*) const {
324 AliInfo("====================================================================================");
325 AliInfo(Form("AliCFContainer : name = %s title = %s",GetName(),GetTitle()));
326 AliInfo(Form("number of steps \t %d",GetNStep()));
327 for (Int_t iStep=0;iStep<GetNStep();iStep++) AliInfo(Form("step %d \t -> %s",iStep,GetStepTitle(iStep)));
328 AliInfo(Form("number of variables \t %d",GetNVar()));
329 for (Int_t iVar=0;iVar<GetNVar();iVar++) {
330 Double_t *binLimits = new Double_t[GetNBins(iVar)+1];
331 GetBinLimits(iVar,binLimits);
332 AliInfo(Form("variable %d \t -> %s : %d bins in [%f,%f]",iVar,GetVarTitle(iVar),GetNBins(iVar),binLimits[0],binLimits[GetNBins(iVar)]));
335 AliInfo("====================================================================================");