BadChannelMap.cxx

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
#include "CPVBase/Geometry.h"
#include "DataFormatsCPV/BadChannelMap.h"
 
#include <fairlogger/Logger.h>
 
#include <TH2.h>
 
#include <iostream>
 
using namespace o2::cpv;
 
BadChannelMap::BadChannelMap(short test)
{
  // reset all channels to be good
  mBadCells.reset();
 
  if (test == 2) {
    // Mark few channels as bad for test peurposes
    for (short i = 0; i < 60; i++) {
      // module 2
      unsigned short channelID = 3584 + i * 57;
      mBadCells.set(channelID);
      channelID = 3640 + i * 55;
      mBadCells.set(channelID);
    }
 
    for (short i = 0; i < 16; i++) {
      // module 3
      unsigned short channelID = 8972 + i * 57;
      mBadCells.set(channelID);
      channelID = 8092 + i * 57;
      mBadCells.set(channelID);
      channelID = 8147 + i * 55;
      mBadCells.set(channelID);
      channelID = 9059 + i * 55;
      mBadCells.set(channelID);
    }
  }
}
 
void BadChannelMap::getHistogramRepresentation(short module, TH2* h) const
{
  if (!h) {
    LOG(error) << "provide histogram to be filled";
  }
 
  const short MAXX = 128,
              MAXZ = 60;
  if (h->GetNbinsX() != MAXX || h->GetNbinsY() != MAXZ) {
    LOG(error) << "Wrong dimentions of input histogram:" << h->GetNbinsX() << "," << h->GetNbinsY() << " instead of " << MAXX << "," << MAXZ;
    return;
  }
 
  h->Reset();
 
  short relid[3] = {module, 1, 1};
  unsigned short absId;
  for (short ix = 1; ix <= MAXX; ix++) {
    relid[1] = ix;
    for (short iz = 1; iz <= MAXZ; iz++) {
      relid[2] = iz;
      if (Geometry::relToAbsNumbering(relid, absId)) {
        if (!isChannelGood(absId)) {
          h->SetBinContent(ix, iz, 1);
        }
      }
    }
  }
}
 
void BadChannelMap::PrintStream(std::ostream& stream) const
{
  // first sort bad channel IDs
  stream << "Number of bad cells:  " << mBadCells.count() << "\n";
  for (unsigned int cellID = 0; cellID < mBadCells.size(); cellID++) {
    if (mBadCells.test(cellID)) {
      stream << cellID << "\n";
    }
  }
}
 
std::ostream& o2::cpv::operator<<(std::ostream& stream, const BadChannelMap& bcm)
{
  bcm.PrintStream(stream);
  return stream;
}