//
//
#include "THnBase.h"
-#include "THnSparse.h"
+#include "THn.h"
#include "TCollection.h"
+#include "AliLog.h"
#include "AliTRDdEdxBaseUtils.h"
#include "AliTRDdEdxCalibHistArray.h"
SetName(GetArrayName());
SetOwner(kTRUE);
- const Int_t nbin[2]={AliTRDdEdxBaseUtils::NTRDtimebin(), 11250};
- const Double_t xmin[2]={0, 0};
- const Double_t xmax[2]={nbin[0], 20};
+ const Int_t nbin[2]={AliTRDdEdxBaseUtils::NTRDtimebin(), 200};
+ const Double_t xmin[2]={0, 0.01};
+ const Double_t xmax[2]={static_cast<Double_t>(nbin[0]), 10};
+ const TString aname[2]={"globalTimeBin", "trdqovertpc"};
+ const TString atitle[2]={"det * AliTRDseedV1::kNtb + itb", "TRD-Cluster-Q / TPC-Signal"};
for(Int_t iter=0; iter<GetSize(); iter++){
- THnBase *hi = new THnSparseS(GetNameAt(iter), "", 2, nbin, xmin, xmax);
+ THnBase *hi = new THnF(GetNameAt(iter), "", 2, nbin, xmin, xmax);
+ for(Int_t iaxis=0; iaxis<2; iaxis++){
+ TAxis *xi = hi->GetAxis(iaxis);
+ xi->SetName(aname[iaxis]);
+ xi->SetTitle(atitle[iaxis]);
+ //only log for y-axis
+ if(iaxis==1){
+ AliTRDdEdxBaseUtils::BinLogX(xi);
+ }
+ }
AddAt(hi, iter);
}
}
if (entry == 0) continue;
if(GetSize()!= entry->GetSize()){
- printf("AliTRDdEdxCalibHistArray::Merge GetSize()!= entry->GetSize() %d %d\n", GetSize(), entry->GetSize()); exit(1);
+ AliFatal(Form("GetSize()!= entry->GetSize() %d %d\n", GetSize(), entry->GetSize()));
}
for(Int_t ii=0; ii<GetSize(); ii++){