#include "Cal/AliTRDCalPad.h"
#include "Cal/AliTRDCalDet.h"
#include "Cal/AliTRDCalDCS.h"
+#include "Cal/AliTRDCalDCSv2.h"
#include "Cal/AliTRDCalPID.h"
#include "Cal/AliTRDCalMonitoring.h"
#include "Cal/AliTRDCalChamberStatus.h"
#include "Cal/AliTRDCalPadStatus.h"
#include "Cal/AliTRDCalSingleChamberStatus.h"
+#include "Cal/AliTRDCalTrkAttach.h"
ClassImp(AliTRDcalibDB)
,fPRFhi(0)
,fPRFwid(0)
,fPRFpad(0)
+ ,fPIDResponse(NULL)
{
//
// Default constructor
,fPRFhi(0)
,fPRFwid(0)
,fPRFpad(0)
+ ,fPIDResponse(NULL)
{
//
// Copy constructor (not that it make any sense for a singleton...)
delete [] fPRFsmp;
fPRFsmp = 0;
}
+ if(fPIDResponse) delete fPIDResponse;
Invalidate();
case kIDPIDLQ :
return CacheCDBEntry(kIDPIDLQ ,"TRD/Calib/PIDLQ");
break;
+ case kIDPIDLQ1D:
+ return CacheCDBEntry(kIDPIDLQ1D ,"TRD/Calib/PIDLQ1D");
+ break;
case kIDRecoParam :
- return CacheCDBEntry(kIDRecoParam ,"TRD/Calib/RecoParam");
+ return CacheCDBEntry(kIDRecoParam ,"TRD/Calib/RecoParam");
+ break;
+ case kIDAttach :
+ return CacheCDBEntry(kIDAttach ,"TRD/Calib/TrkAttach");
break;
-
}
return 0;
//
// Caches the entry <id> with cdb path <cdbPath>
//
-
+
if (!fCDBCache[id]) {
fCDBEntries[id] = GetCDBEntry(cdbPath);
if (fCDBEntries[id]) {
fCDBCache[id] = fCDBEntries[id]->GetObject();
}
- }
-
+ }
+
return fCDBCache[id];
-
}
//_____________________________________________________________________________
return -1;
}
- Double_t mean = 0.0;
+ Double_t sum = 0.0;
for (Int_t channel = 0; channel < roc->GetNchannels(); ++channel) {
- mean += (calDet->GetValue(det) + roc->GetValue(channel)) / roc->GetNchannels();
+ sum += roc->GetValue(channel);
}
-
- return mean;
+ sum /= roc->GetNchannels();
+ sum += calDet->GetValue(det);
+ return sum;
}
// Returns Number of time bins from the DCS
//
- const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
+ Int_t nMixed = -2; // not the same number for all chambers
+ Int_t nUndef = -1; // default value - has not been set!
+ Int_t nTbSor = nUndef;
+ Int_t nTbEor = nUndef;
+ Int_t calver = 0; // Check CalDCS version
- if(!dcsArr){
- printf("No DCS Object found\n");
- //return -1;
- return 30;
- }
- const AliTRDCalDCS *calDCSsor = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(0)); // Take SOR
- const AliTRDCalDCS *calDCSeor = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(1));
-
- // prefer SOR
- if(!calDCSsor){
- if(!calDCSeor){
- printf("No calDCSeor found\n");
- //return -1;
- return 30;
+ const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
+ if (!dcsArr) {
+ AliError("No DCS object found!");
+ return nUndef;
+ }
+
+ if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCS")) calver = 1;
+ if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
+
+ if (calver == 1) {
+ // DCS object
+ const AliTRDCalDCS *calDCSsor = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(0));
+ const AliTRDCalDCS *calDCSeor = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(1));
+ if (!calDCSsor) {
+ // the SOR file is mandatory
+ AliError("NO SOR AliTRDCalDCS object found in CDB file!");
+ return nUndef;
}
- if(calDCSeor->GetGlobalNumberOfTimeBins() > 0.0) return calDCSeor->GetGlobalNumberOfTimeBins();
- else return 30;
- }
- // if SOR is available and the number of timebins is > -1, take this, otherwise check EOR
- Int_t nTimeSOR = calDCSsor->GetGlobalNumberOfTimeBins();
- if(nTimeSOR > -1){
- // Make a consistency check
- if(calDCSeor){
- Int_t nTimeEOR = calDCSeor->GetGlobalNumberOfTimeBins();
- if((nTimeEOR > -1) && (nTimeSOR != nTimeEOR)){
- // Parameter inconsistency found, return -2 to be able to catch the error
- //return -2;
- printf("Inconsistency\n");
- return 30;
- }
+ if (!calDCSeor) {
+ // this can happen if the run is shorter than a couple of seconds.
+ AliWarning("NO EOR AliTRDCalDCS object found in CDB file.");
+ }
+
+ // get the numbers
+ nTbSor = calDCSsor->GetGlobalNumberOfTimeBins();
+ if (calDCSeor) nTbEor = calDCSeor->GetGlobalNumberOfTimeBins();
+
+ } else if (calver == 2) {
+ // DCSv2 object
+ const AliTRDCalDCSv2 *calDCSsorv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(0));
+ const AliTRDCalDCSv2 *calDCSeorv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(1));
+ if (!calDCSsorv2) {
+ // the SOR file is mandatory
+ AliError("NO SOR AliTRDCalDCSv2 object found in CDB file!");
+ return nUndef;
}
- // Consisency check passed or not done
- if(nTimeSOR > 0.0) return nTimeSOR;
- else return 30;
- } else {
- // SOR has unphysical time parameter, take EOR
- if(calDCSeor) {
- if(calDCSeor->GetGlobalNumberOfTimeBins() > 0.0) return calDCSeor->GetGlobalNumberOfTimeBins();
- else return 30;
+ if (!calDCSeorv2) {
+ // this can happen if the run is shorter than a couple of seconds.
+ AliWarning("NO EOR AliTRDCalDCSv2 object found in CDB file.");
}
- return 30; // Both SOR and EOR not available
+
+ // get the numbers
+ nTbSor = calDCSsorv2->GetGlobalNumberOfTimeBins();
+ if (calDCSeorv2) nTbEor = calDCSeorv2->GetGlobalNumberOfTimeBins();
+
+ } else AliError("NO DCS/DCSv2 OCDB entry found!");
+
+ // if they're the same return the value
+ // -2 means mixed, -1: no data, >= 0: good number of time bins
+ if (nTbSor == nTbEor) return nTbSor;
+
+ // if they're differing:
+ if (nTbSor == nMixed || nTbEor == nMixed) {
+ AliWarning("Inconsistent number of time bins found!");
+ return nMixed;
}
+
+ // one is undefined, the other ok -> return that one
+ if (nTbSor == nUndef) return nTbEor;
+ if (nTbEor == nUndef) return nTbSor;
+
+ // only remains: two different numbers >= 0
+ return nMixed;
}
//_____________________________________________________________________________
filterType = "";
return;
}
- const AliTRDCalDCS *calDCS = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(1)); // Take EOR
-
- if(!calDCS){
- filterType = "";
+
+ Int_t esor = 0; // Take SOR
+ Int_t calver = 0; // Check CalDCS version
+ if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCS")) calver = 1;
+ if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
+
+ if (calver == 1) {
+ // DCS object
+ const AliTRDCalDCS *calDCS = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(esor));
+ if(!calDCS){
+ filterType = "";
+ return;
+ }
+ filterType = calDCS->GetGlobalFilterType();
+ } else if (calver == 2) {
+ // DCSv2 object
+ const AliTRDCalDCSv2 *calDCSv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(esor));
+ if(!calDCSv2){
+ filterType = "";
+ return;
+ }
+ filterType = calDCSv2->GetGlobalFilterType();
+ } else AliError("NO DCS/DCSv2 OCDB entry found!");
+
+}
+
+//_____________________________________________________________________________
+void AliTRDcalibDB::GetGlobalConfiguration(TString &config){
+ //
+ // Get Configuration from the DCS
+ //
+
+ const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
+ if(!dcsArr){
+ config = "";
return;
- }
- filterType = calDCS->GetGlobalFilterType();
+ }
+
+ Int_t esor = 0; // Take SOR
+ Int_t calver = 0; // Check CalDCS version
+ if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCS")) calver = 1;
+ if (!strcmp(dcsArr->At(0)->ClassName(),"AliTRDCalDCSv2")) calver = 2;
+
+ if (calver == 1) {
+ // DCS object
+ const AliTRDCalDCS *calDCS = dynamic_cast<const AliTRDCalDCS *>(dcsArr->At(esor));
+ if(!calDCS){
+ config = "";
+ return;
+ }
+ config = calDCS->GetGlobalConfigName();
+ } else if (calver == 2) {
+ // DCSv2 object
+ const AliTRDCalDCSv2 *calDCSv2 = dynamic_cast<const AliTRDCalDCSv2 *>(dcsArr->At(esor));
+ if(!calDCSv2){
+ config = "";
+ return;
+ }
+ config = calDCSv2->GetGlobalConfigName();
+ } else AliError("NO DCS/DCSv2 OCDB entry found!");
+
+}
+
+//_____________________________________________________________________________
+Bool_t AliTRDcalibDB::HasOnlineFilterPedestal()
+{
+ //
+ // Checks whether pedestal filter was applied online
+ //
+ TString cname;
+ // Temporary: Get the filter config from the configuration name
+ GetGlobalConfiguration(cname);
+ TString filterconfig = cname(cname.First("_") + 1, cname.First("-") - cname.First("_") - 1);
+ // TString filterconfig;
+ //GetFilterType(filterconfig);
+ return filterconfig.Contains("p");
+}
+
+//_____________________________________________________________________________
+Bool_t AliTRDcalibDB::HasOnlineFilterGain(){
+ //
+ // Checks whether online gain filter was applied
+ //
+ TString cname;
+ // Temporary: Get the filter config from the configuration name
+ GetGlobalConfiguration(cname);
+ TString filterconfig = cname(cname.First("_") + 1, cname.First("-") - cname.First("_") - 1);
+ //TString filterconfig;
+ //GetFilterType(filterconfig);
+ return filterconfig.Contains("g");
+}
+//_____________________________________________________________________________
+Bool_t AliTRDcalibDB::HasOnlineTailCancellation(){
+ //
+ // Checks whether online tail cancellation was applied
+ //
+ TString cname;
+ // Temporary: Get the filter config from the configuration name
+ GetGlobalConfiguration(cname);
+ TString filterconfig = cname(cname.First("_") + 1, cname.First("-") - cname.First("_") - 1);
+ //TString filterconfig;
+ //GetFilterType(filterconfig);
+ return filterconfig.Contains("t");
}
//_____________________________________________________________________________
}
+//_____________________________________________________________________________
+AliTRDPIDResponse *AliTRDcalibDB::GetPIDResponse(AliTRDPIDResponse::ETRDPIDMethod method){
+ if(!fPIDResponse){
+ fPIDResponse = new AliTRDPIDResponse;
+ // Load Reference Histos from OCDB
+ fPIDResponse->SetPIDmethod(method);
+ fPIDResponse->Load(dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDPIDLQ1D)));
+ }
+ return fPIDResponse;
+}
+
+//_____________________________________________________________________________
+const AliTRDCalTrkAttach* AliTRDcalibDB::GetAttachObject()
+{
+ //
+ // Returns the object storing likelihood distributions for cluster to track attachment
+ //
+ return dynamic_cast<const AliTRDCalTrkAttach*>(GetCachedCDBObject(kIDAttach));
+}
+
+
//_____________________________________________________________________________
const AliTRDCalMonitoring *AliTRDcalibDB::GetMonitoringObject()
{
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff