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 **************************************************************************/
24 // lu@physi.uni-heidelberg.de
30 #include "TCollection.h"
33 #include "AliTRDdEdxBaseUtils.h"
34 #include "AliTRDdEdxCalibHistArray.h"
36 ClassImp(AliTRDdEdxCalibHistArray);
38 AliTRDdEdxCalibHistArray::AliTRDdEdxCalibHistArray(const Bool_t kNoInv):
39 TObjArray(kNoInv ? 4: 8)
44 SetName(GetArrayName());
47 const Int_t nbin[2]={AliTRDdEdxBaseUtils::NTRDtimebin(), 200};
48 const Double_t xmin[2]={0, 0.01};
49 const Double_t xmax[2]={nbin[0], 10};
50 const TString aname[2]={"globalTimeBin", "trdqovertpc"};
51 const TString atitle[2]={"det * AliTRDseedV1::kNtb + itb", "TRD-Cluster-Q / TPC-Signal"};
53 for(Int_t iter=0; iter<GetSize(); iter++){
54 THnBase *hi = new THnF(GetNameAt(iter), "", 2, nbin, xmin, xmax);
55 for(Int_t iaxis=0; iaxis<2; iaxis++){
56 TAxis *xi = hi->GetAxis(iaxis);
57 xi->SetName(aname[iaxis]);
58 xi->SetTitle(atitle[iaxis]);
61 AliTRDdEdxBaseUtils::BinLogX(xi);
68 AliTRDdEdxCalibHistArray::AliTRDdEdxCalibHistArray(const AliTRDdEdxCalibHistArray &obj):
76 AliTRDdEdxCalibHistArray & AliTRDdEdxCalibHistArray::operator=(const AliTRDdEdxCalibHistArray &obj)
82 if(&obj == this) return *this;
84 TObjArray::operator=(obj);
89 Long64_t AliTRDdEdxCalibHistArray::Merge(const TCollection* list)
92 // Merge list of objects (needed by PROOF)
101 TIterator* iter = list->MakeIterator();
105 while((obj = iter->Next()) != 0)
107 AliTRDdEdxCalibHistArray * entry = dynamic_cast<AliTRDdEdxCalibHistArray*>(obj);
108 if (entry == 0) continue;
110 if(GetSize()!= entry->GetSize()){
111 AliFatal(Form("GetSize()!= entry->GetSize() %d %d\n", GetSize(), entry->GetSize()));
114 for(Int_t ii=0; ii<GetSize(); ii++){
115 THnBase *h0 = (THnBase*) At(ii);
116 THnBase *h1 = (THnBase*) entry->At(ii);