#include "AliMUON2DStoreValidator.h"
#include "AliMUONCalibParamNI.h"
#include "AliMUONCalibrationData.h"
+#include "AliMUONLogger.h"
+#include "AliMUONRecoParam.h"
#include "AliMUONStringIntMap.h"
-#include "AliMUONVCalibParam.h"
#include "AliMUONTrackerData.h"
+#include "AliMUONVCalibParam.h"
#include "AliMpArea.h"
#include "AliMpArrayI.h"
+#include "AliMpCDB.h"
#include "AliMpConstants.h"
#include "AliMpDDLStore.h"
#include "AliMpDEManager.h"
#include "AliMpDetElement.h"
#include "AliMpDCSNamer.h"
+#include "AliMpManuIterator.h"
#include "AliMpManuUID.h"
#include "AliCDBEntry.h"
#include <TString.h>
#include <TSystem.h>
+using std::cout;
+using std::endl;
/// \cond CLASSIMP
ClassImp(AliMUONPadStatusMaker)
/// \endcond
fGainA1Limits(0,1E30),
fGainA2Limits(-1E-30,1E30),
fGainThresLimits(0,4095),
-fHVSt12Limits(0,5000),
-fHVSt345Limits(0,5000),
fPedMeanLimits(0,4095),
fPedSigmaLimits(0,4095),
fManuOccupancyLimits(0,1.0),
-fBusPatchOccupancyLimits(0,1.0),
+fBuspatchOccupancyLimits(0,1.0),
fDEOccupancyLimits(0,1.0),
fStatus(new AliMUON2DMap(true)),
-fHV(new TExMap),
+fHV(0x0),
fPedestals(calibData.Pedestals()),
fGains(calibData.Gains()),
fTrackerData(0x0)
{
/// ctor
- AliDebug(1,Form("ped store %s gain store %s",
- fPedestals->ClassName(),
- fGains->ClassName()));
-
if ( calibData.OccupancyMap() )
{
- AliInfo("Will use occupancy map to cut, if so required in AliMUONRecoParam");
/// create a tracker data from the occupancy map
fTrackerData = new AliMUONTrackerData("OCC","OCC",*(calibData.OccupancyMap()));
- }
+ }
+ if ( calibData.HV() )
+ {
+ /// Only create the fHV internal store if there are some HV values available
+ fHV = new TExMap;
+ }
+
+ SetHVLimit(-1,0.0);
}
//_____________________________________________________________________________
AliMUONPadStatusMaker::~AliMUONPadStatusMaker()
{
/// dtor.
+
delete fStatus;
delete fHV;
delete fTrackerData;
{
/// return a human readable version of the integer status
+ if ( status == 0 )
+ {
+ return "Brave New World";
+ }
+
Int_t pedStatus;
Int_t gainStatus;
Int_t hvStatus;
/// and the switch).
/// Returns false if hv switch changed during the run.
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
+
+ if (!fHV) return kFALSE;
Bool_t error = kFALSE;
hvChannelTooLow = kFALSE;
hvChannelTooHigh = kFALSE;
hvChannelON = kTRUE;
+ Int_t chamberId = AliMpDEManager::GetChamberId(detElemId);
+
AliMpDCSNamer hvNamer("TRACKER");
- TString hvChannel(hvNamer.DCSChannelName(detElemId,sector));
+ TString hvChannel(hvNamer.DCSAliasName(detElemId,sector));
TMap* hvMap = fkCalibrationData.HV();
TPair* hvPair = static_cast<TPair*>(hvMap->FindObject(hvChannel.Data()));
}
else
{
- // find out min and max value, and makes a cut
+ // find out min value, and makes a cut
Float_t hvMin(1E9);
- Float_t hvMax(0);
TIter next(values);
AliDCSValue* val;
{
Float_t hv = val->GetFloat();
hvMin = TMath::Min(hv,hvMin);
- hvMax = TMath::Max(hv,hvMax);
}
- float lowThreshold = fHVSt12Limits.X();
- float highThreshold = fHVSt12Limits.Y();
+ float lowThreshold = HVLimit(chamberId);
if ( hvMin < lowThreshold ) hvChannelTooLow = kTRUE;
- if ( hvMax > highThreshold ) hvChannelTooHigh = kTRUE;
- if ( hvMin < 1 ) hvChannelON = kFALSE;
+ if ( hvMin < hvNamer.TrackerHVOFF() ) hvChannelON = kFALSE;
}
}
return error;
}
+//_____________________________________________________________________________
+Float_t
+AliMUONPadStatusMaker::SwitchValue(const TObjArray& dcsArray)
+{
+ /// Loop over the dcs value for a single switch to decide whether
+ /// we should consider it on or off
+
+ // we'll count the number of ON/OFF for this pad, to insure
+ // consistency (i.e. if status changed during the run, we should
+ // at least notify this fact ;-) and hope it's not the norm)
+ Int_t nTrue(0);
+ Int_t nFalse(0);
+ TIter next(&dcsArray);
+ AliDCSValue* val;
+
+ while ( ( val = static_cast<AliDCSValue*>(next()) ) )
+ {
+ if ( val->GetBool() )
+ {
+ ++nTrue;
+ }
+ else
+ {
+ ++nFalse;
+ }
+ }
+
+ if ( (nTrue>0 && nFalse>0) )
+ {
+ // change of state during the run, consider it off
+ return 0.0;
+ }
+
+ if ( nFalse )
+ {
+ /// switch = FALSE means the HV was flowding up to the PCB.
+ /// i.e. switch = FALSE = ON
+ return 1.0;
+ }
+
+ return 0.0;
+}
+
//_____________________________________________________________________________
Bool_t
AliMUONPadStatusMaker::HVSt345Status(Int_t detElemId, Int_t pcbIndex,
/// Returns false if something goes wrong (in particular if
/// hv switch changed during the run).
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
+
+ if (!fHV) return kFALSE;
Bool_t error = kFALSE;
hvChannelTooLow = kFALSE;
AliMpDCSNamer hvNamer("TRACKER");
- TString hvChannel(hvNamer.DCSChannelName(detElemId));
+ Int_t chamberId = AliMpDEManager::GetChamberId(detElemId);
+
+ TString hvChannel(hvNamer.DCSAliasName(detElemId));
TMap* hvMap = fkCalibrationData.HV();
}
else
{
- // find out min and max value, and makes a cut
+ // find out min value, and makes a cut
Float_t hvMin(1E9);
- Float_t hvMax(0);
TIter next(values);
AliDCSValue* val;
{
Float_t hv = val->GetFloat();
hvMin = TMath::Min(hv,hvMin);
- hvMax = TMath::Max(hv,hvMax);
}
- float lowThreshold = fHVSt345Limits.X();
- float highThreshold = fHVSt345Limits.Y();
+ float lowThreshold = HVLimit(chamberId);
if ( hvMin < lowThreshold ) hvChannelTooLow = kTRUE;
- else if ( hvMax > highThreshold ) hvChannelTooHigh = kTRUE;
- if ( hvMin < 1 ) hvChannelON = kFALSE;
+ if ( hvMin < hvNamer.TrackerHVOFF() ) hvChannelON = kFALSE;
}
}
- TString hvSwitch(hvNamer.DCSSwitchName(detElemId,pcbIndex));
+ TString hvSwitch(hvNamer.DCSSwitchAliasName(detElemId,pcbIndex));
TPair* switchPair = static_cast<TPair*>(hvMap->FindObject(hvSwitch.Data()));
if (!switchPair)
{
}
else
{
- // we'll count the number of ON/OFF for this pad, to insure
- // consistency (i.e. if status changed during the run, we should
- // at least notify this fact ;-) and hope it's not the norm)
- Int_t nTrue(0);
- Int_t nFalse(0);
- TIter next(values);
- AliDCSValue* val;
-
- while ( ( val = static_cast<AliDCSValue*>(next()) ) )
- {
- if ( val->GetBool() )
- {
- ++nTrue;
- }
- else
- {
- ++nFalse;
- }
- }
-
- if ( (nTrue>0 && nFalse>0) )
- {
- AliWarning(Form("Status of HV Switch %s changed during this run nTrue=%d nFalse=%d! Will consider it OFF",
- hvSwitch.Data(),nTrue,nFalse));
- error = kTRUE;
- }
-
- if ( nFalse ) hvSwitchON = kFALSE;
+ Float_t sv = SwitchValue(*values);
+ if ( sv < 0.99 ) hvSwitchON = kFALSE;
}
}
return error;
{
/// Get HV status of one manu
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
- if ( !fkCalibrationData.HV() ) return kMissing;
+ if ( !fHV ) return kMissing;
Long_t lint = fHV->GetValue(AliMpManuUID::BuildUniqueID(detElemId,manuId));
occ = fTrackerData->BusPatch(busPatchId,occIndex);
- if ( occ <= fBusPatchOccupancyLimits.X() )
+ if ( occ <= fBuspatchOccupancyLimits.X() )
{
rv |= kBusPatchOccupancyTooLow;
}
- else if ( occ > fBusPatchOccupancyLimits.Y() )
+ else if ( occ > fBuspatchOccupancyLimits.Y() )
{
rv |= kBusPatchOccupancyTooHigh;
}
if (!param)
{
// not already there, so compute it now
- AliCodeTimerAuto("ComputeStatus");
+ AliCodeTimerAuto("ComputeStatus",0);
param = ComputeStatus(detElemId,manuId);
}
return param;
/// Assign status to all manus in a given HV "zone" (defined by index, meaning
/// is different thing from St12 and St345)
- AliCodeTimerAuto("")
+ AliCodeTimerAuto("",0)
AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
fHV->Add(AliMpManuUID::BuildUniqueID(detElemId,manuId),status + 1);
}
}
+
+//_____________________________________________________________________________
+Double_t
+AliMUONPadStatusMaker::HVLimit(Int_t chamberId) const
+{
+ /// Get HV limit for a given chamber
+ if ( chamberId >=0 && chamberId < 10 )
+ {
+ return fHVLimit[chamberId];
+ }
+ return 0.0;
+}
+
+//_____________________________________________________________________________
+void
+AliMUONPadStatusMaker::SetHVLimit(Int_t chamberId, Double_t hv)
+{
+ /// Set hv limit for a given chamber (or all if chamberId==-1)
+
+ if ( chamberId == -1 )
+ {
+ for ( Int_t i = 0; i < 10; ++i )
+ {
+ fHVLimit[i] = hv;
+ }
+ }
+ else if ( chamberId >= 0 && chamberId < 10 )
+ {
+ fHVLimit[chamberId]=hv;
+ }
+ else
+ {
+ AliError(Form("chamberId=%d is invalid",chamberId));
+ }
+}
+
+//_____________________________________________________________________________
+void
+AliMUONPadStatusMaker::SetLimits(const AliMUONRecoParam& recoParams)
+{
+ /// Set the limits from the recoparam
+
+ for ( int i = 0; i < 10; ++i )
+ {
+ SetHVLimit(i,recoParams.HVLimit(i));
+ }
+
+ SetPedMeanLimits(recoParams.PedMeanLowLimit(),recoParams.PedMeanHighLimit());
+ SetPedSigmaLimits(recoParams.PedSigmaLowLimit(),recoParams.PedSigmaHighLimit());
+
+ SetGainA1Limits(recoParams.GainA1LowLimit(),recoParams.GainA1HighLimit());
+ SetGainA2Limits(recoParams.GainA2LowLimit(),recoParams.GainA2HighLimit());
+ SetGainThresLimits(recoParams.GainThresLowLimit(),recoParams.GainThresHighLimit());
+
+ SetManuOccupancyLimits(recoParams.ManuOccupancyLowLimit(),recoParams.ManuOccupancyHighLimit());
+ SetBuspatchOccupancyLimits(recoParams.BuspatchOccupancyLowLimit(),recoParams.BuspatchOccupancyHighLimit());
+ SetDEOccupancyLimits(recoParams.DEOccupancyLowLimit(),recoParams.DEOccupancyHighLimit());
+}
+
+//_____________________________________________________________________________
+void
+AliMUONPadStatusMaker::Report(UInt_t mask)
+{
+ /// Report the number of bad pads, according to the mask,
+ /// and the various reasons why they are bad (with occurence rates)
+
+ AliInfo("");
+ AliCodeTimerAuto("",0);
+
+ AliMUONLogger log(1064008);
+
+ Int_t nBadPads(0);
+ Int_t nPads(0);
+
+ AliMpManuIterator it;
+
+ Int_t detElemId, manuId;
+
+ while ( it.Next(detElemId,manuId) )
+ {
+ AliMpDetElement* de = AliMpDDLStore::Instance()->GetDetElement(detElemId);
+
+ for ( Int_t i = 0; i < AliMpConstants::ManuNofChannels(); ++i )
+ {
+ if ( de->IsConnectedChannel(manuId,i) )
+ {
+ ++nPads;
+
+ Int_t status = PadStatus(detElemId,manuId,i);
+
+ if ( mask && ( status & mask) ) // note that if mask == 0, all pads are good...
+ {
+ ++nBadPads;
+ log.Log(AsString(status));
+ }
+ }
+ }
+ }
+
+ if (!nPads)
+ {
+ AliError("Got no pad from the iterator ?! That's not normal. Please check !");
+ return;
+ }
+
+ TString msg;
+ Int_t ntimes;
+
+ cout << Form("According to mask %x (human readable form below) %6d pads are bad (over a total of %6d, i.e. %7.2f %%)",
+ mask,nBadPads,nPads,nBadPads*100.0/nPads) << endl;
+ cout << AliMUONPadStatusMaker::AsCondition(mask) << endl;
+ cout << "--------" << endl;
+
+ while ( log.Next(msg,ntimes) )
+ {
+ cout << Form("The message (%120s) occured %15d times (%7.4f %%)",msg.Data(),ntimes,ntimes*100.0/nPads) << endl;
+ }
+
+}
+
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff