Coverity fix
[u/mrichter/AliRoot.git] / TRD / AliTRDdEdxCalibHistArray.cxx
CommitLineData
6951a056 1/**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
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
16//
17// xx
18// xx
19// xx
20// xx
21// xx
22//
23// Xianguo Lu
24// lu@physi.uni-heidelberg.de
25// Xianguo.Lu@cern.ch
26//
27//
28#include "THnBase.h"
6856b6a8 29#include "THn.h"
6951a056 30#include "TCollection.h"
31
32#include "AliTRDdEdxBaseUtils.h"
33#include "AliTRDdEdxCalibHistArray.h"
34
35ClassImp(AliTRDdEdxCalibHistArray);
36
37AliTRDdEdxCalibHistArray::AliTRDdEdxCalibHistArray(const Bool_t kNoInv):
38 TObjArray(kNoInv ? 4: 8)
39{
40 //
41 //constructor
42 //
43 SetName(GetArrayName());
44 SetOwner(kTRUE);
45
687aa844 46 const Int_t nbin[2]={AliTRDdEdxBaseUtils::NTRDtimebin(), 200};
47 const Double_t xmin[2]={0, 0.01};
48 const Double_t xmax[2]={nbin[0], 10};
6856b6a8 49 const TString aname[2]={"globalTimeBin", "trdqovertpc"};
50 const TString atitle[2]={"det * AliTRDseedV1::kNtb + itb", "TRD-Cluster-Q / TPC-Signal"};
6951a056 51
52 for(Int_t iter=0; iter<GetSize(); iter++){
6856b6a8 53 THnBase *hi = new THnF(GetNameAt(iter), "", 2, nbin, xmin, xmax);
54 for(Int_t iaxis=0; iaxis<2; iaxis++){
55 TAxis *xi = hi->GetAxis(iaxis);
56 xi->SetName(aname[iaxis]);
57 xi->SetTitle(atitle[iaxis]);
396eba11 58 //only log for y-axis
59 if(iaxis==1){
60 AliTRDdEdxBaseUtils::BinLogX(xi);
61 }
6856b6a8 62 }
6951a056 63 AddAt(hi, iter);
64 }
65}
66
67AliTRDdEdxCalibHistArray::AliTRDdEdxCalibHistArray(const AliTRDdEdxCalibHistArray &obj):
68 TObjArray(obj)
69{
70 //
71 //copy constructor
72 //
73}
74
75AliTRDdEdxCalibHistArray & AliTRDdEdxCalibHistArray::operator=(const AliTRDdEdxCalibHistArray &obj)
76{
77 //
78 //assignment operator
79 //
80
81 if(&obj == this) return *this;
82
83 TObjArray::operator=(obj);
84
85 return *this;
86}
87
88Long64_t AliTRDdEdxCalibHistArray::Merge(const TCollection* list)
89{
90 //
91 // Merge list of objects (needed by PROOF)
92 //
93
94 if(!list)
95 return 0;
96
97 if(list->IsEmpty())
98 return 1;
99
100 TIterator* iter = list->MakeIterator();
101 TObject* obj = 0;
102
103 Int_t count=0;
104 while((obj = iter->Next()) != 0)
105 {
106 AliTRDdEdxCalibHistArray * entry = dynamic_cast<AliTRDdEdxCalibHistArray*>(obj);
107 if (entry == 0) continue;
108
109 if(GetSize()!= entry->GetSize()){
110 printf("AliTRDdEdxCalibHistArray::Merge GetSize()!= entry->GetSize() %d %d\n", GetSize(), entry->GetSize()); exit(1);
111 }
112
113 for(Int_t ii=0; ii<GetSize(); ii++){
114 THnBase *h0 = (THnBase*) At(ii);
115 THnBase *h1 = (THnBase*) entry->At(ii);
116 h0->Add(h1);
117 }
118
119 count++;
120 }
121
122 return count;
123
124}