+ const Int_t kWS(AliPID::kSPECIES*AliTRDCalPID::kNMom);
+ Float_t *data[2*kWS];
+ for(Int_t i(0); i<2*kWS; i++) data[i]=new Float_t[kMaxStat];
+ Int_t ndata[kWS]; memset(ndata, 0, kWS*sizeof(Int_t));
+
+ AliDebug(1, Form("Loading data[%d]", fData->GetEntries()));
+ for(Int_t itrk=0; itrk < fData->GetEntries(); itrk++){
+ if(!(fData->GetEntry(itrk))) continue;
+ Int_t sbin(fPIDbin);
+ for(Int_t ily=fPIDdataArray->fNtracklets; ily--;){
+ Int_t pbin(fPIDdataArray->fData[ily].fPLbin & 0xf);
+
+ Double_t dedx[] = {0., 0.};
+ if(!AliTRDCalPIDLQ::CookdEdx(fPIDdataArray->fData[ily].fdEdx, dedx)) continue;
+ Int_t idx=AliTRDCalPIDLQ::GetModelID(pbin,sbin);
+ if(ndata[idx]==kMaxStat) continue;
+
+ // store data
+ data[idx][ndata[idx]] = dedx[0];
+ data[idx+kWS][ndata[idx]] = dedx[1];
+ ndata[idx]++;
+ }
+ }