///////////////////////////////////////////////////////////////////////////////
#include <TClonesArray.h>
+#include <TObjArray.h>
#include "AliCDBManager.h"
#include "AliCDBEntry.h"
#include "Cal/AliTRDCalROC.h"
#include "Cal/AliTRDCalPad.h"
#include "Cal/AliTRDCalDet.h"
-#include "Cal/AliTRDCalFEE.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 kIDFEE :
return CacheCDBEntry(kIDFEE ,"TRD/Calib/FEE");
break;
+ case kIDDCS :
+ return CacheCDBEntry(kIDDCS ,"TRD/Calib/DCS");
+ break;
case kIDPIDNN :
return CacheCDBEntry(kIDPIDNN ,"TRD/Calib/PIDNN");
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;
}
}
//_____________________________________________________________________________
-Int_t AliTRDcalibDB::GetNumberOfTimeBins()
+Int_t AliTRDcalibDB::GetNumberOfTimeBinsDCS()
{
//
- // Returns the number of time bins which are read-out.
+ // Returns Number of time bins from the DCS
//
- const AliTRDCalFEE *calFEE = dynamic_cast<const AliTRDCalFEE *>
- (GetCachedCDBObject(kIDFEE));
- if (!calFEE) {
- return -1;
+ 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
+
+ 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) {
+ // 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;
+ }
+ if (!calDCSeorv2) {
+ // this can happen if the run is shorter than a couple of seconds.
+ AliWarning("NO EOR AliTRDCalDCSv2 object found in CDB file.");
+ }
+
+ // 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;
+}
+
+//_____________________________________________________________________________
+void AliTRDcalibDB::GetFilterType(TString &filterType)
+{
+ //
+ // Returns the filter type
+ //
+
+ const TObjArray *dcsArr = dynamic_cast<const TObjArray *>(GetCachedCDBObject(kIDDCS));
+ if(!dcsArr){
+ filterType = "";
+ return;
}
- return calFEE->GetNumberOfTimeBins();
+ 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;
+ }
+
+ 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");
+}
+
//_____________________________________________________________________________
Char_t AliTRDcalibDB::GetPadStatus(Int_t det, Int_t col, Int_t row)
{
//
// Returns status, see name of functions for details ;-)
//
+
const AliTRDCalPadStatus * cal = dynamic_cast<const AliTRDCalPadStatus *>
(GetCachedCDBObject(kIDPadStatus));
if (!cal) {
return cal->IsMasked(det);
}
+//_____________________________________________________________________________
+Bool_t AliTRDcalibDB::IsHalfChamberMasked(Int_t det, Int_t side){
+ //
+ // Returns status, see name of functions for details ;-)
+ //
+ const AliTRDCalChamberStatus * cal = dynamic_cast<const AliTRDCalChamberStatus *>
+ (GetCachedCDBObject(kIDChamberStatus));
+ if (!cal) {
+ return -1;
+ }
+
+ return side > 0 ? cal->IsHalfChamberSideBMasked(det) : cal->IsHalfChamberSideAMasked(det);
+
+}
//_____________________________________________________________________________
const AliTRDCalPID *AliTRDcalibDB::GetPIDObject(AliTRDpidUtil::ETRDPIDMethod method)
{
}
+//_____________________________________________________________________________
+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()
{
}
+