SetTPCPidResponseMaster();
SetTPCParametrisation();
+
+ SetTRDPidResponseMaster();
+ InitializeTRDResponse();
fTOFResponse.SetTimeResolution(fTOFres);
}
fBeamType="PP";
- TPRegexp reg(".*(LHC11[a-z]+[0-9]+[a-z]*)/.*");
+ TPRegexp reg(".*(LHC11[a-z]+[0-9]+[a-z_]*)/.*");
//find the period by run number (UGLY, but not stored in ESD and AOD... )
if (fRun>=114737&&fRun<=117223) { fLHCperiod="LHC10B"; fMCperiodTPC="LHC10D1"; }
else if (fRun>=118503&&fRun<=121040) { fLHCperiod="LHC10C"; fMCperiodTPC="LHC10D1"; }
// Load the TRD pid params and references from the OADB
//
if(fTRDPIDParams) return;
- AliOADBContainer contParams;
+ AliOADBContainer contParams("contParams");
+
contParams.InitFromFile(Form("%s/COMMON/PID/data/TRDPIDParams.root", fOADBPath.Data()), "AliTRDPIDParams");
fTRDPIDParams = (TObjArray *)contParams.GetObject(fRun);
- AliOADBContainer contRefs;
+ AliOADBContainer contRefs("contRefs");
contRefs.InitFromFile(Form("%s/COMMON/PID/dReferencesLQ1D.root", fOADBPath.Data()), "AliTRDPIDReference");
fTRDPIDReference = (AliTRDPIDReference *)contRefs.GetObject(fRun);
}
// returns the threshold for a given number of tracklets and a given efficiency level
//tby definition the lower of step is given.
//
- if(ntracklets > 6 || ntracklets) return NULL;
+ if(ntracklets > 6 || ntracklets <=0) return NULL;
TObjArray * entry = dynamic_cast<TObjArray *>(fkPIDParams->At(ntracklets - 1));
if(!entry) return NULL;
TIter cutIter(entry);
while((cut = dynamic_cast<TObjArray *>(cutIter()))){
effLevel = static_cast<TVectorF *>(cut->At(0));
- if(effLevel[0] > currentLower && effLevel[0] <= level){
+ if((*effLevel)[0] > currentLower && (*effLevel)[0] <= level){
// New Lower entry found
parameters = static_cast<const TVectorD *>(cut->At(1));
+ currentLower = (*effLevel)[0];
}
}