CDCBadWires.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: CDC group                                                *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
#pragma once
 
#include <map>
#include <algorithm>
#include <iostream>
#include <fstream>
#include <iomanip>
#include <TObject.h>
#include <cdc/dataobjects/WireID.h>
 
namespace Belle2 {
  class CDCBadWires: public TObject {
  public:
 
    CDCBadWires() {}
 
    void setWire(const WireID& wid, double eff = 0)
    {
      m_wires.insert(std::pair(wid.getEWire(), eff));
    }
 
    void setWire(unsigned short iCLayer, unsigned short iWire, double eff = 0)
    {
      m_wires.insert(std::pair(WireID(iCLayer, iWire).getEWire(), eff));
    }
 
    unsigned short getEntries() const
    {
      return m_wires.size();
    }
 
    std::map<unsigned short, float> getWires() const
    {
      return m_wires;
    }
 
    bool isBadWire(const WireID& wid) const
    {
      std::map<unsigned short, float>::const_iterator it = m_wires.find(wid.getEWire());
      if (it != m_wires.end() && it->second == 0.) {
        return true;
      } else {
        return false;
      }
    }
 
    bool isDeadWire(const WireID& wid, double& eff) const
    {
      std::map<unsigned short, float>::const_iterator it = m_wires.find(wid.getEWire());
      if (it != m_wires.end() && 0. <= it->second && it->second < 1.) {
        eff = it->second;
        return true;
      } else {
        return false;
      }
    }
 
    bool isHotWire(const WireID& wid) const
    {
      std::map<unsigned short, float>::const_iterator it = m_wires.find(wid.getEWire());
      if (it != m_wires.end() && it->second > 1.) {
        return true;
      } else {
        return false;
      }
    }
 
    void dump() const
    {
      std::cout << " " << std::endl;
      std::cout << "Bad wire list" << std::endl;
      std::cout << "#entries= " << m_wires.size() << std::endl;
      std::cout << "in order of ICLayer and IWire and EWire" << std::endl;
 
      for (auto const ent : m_wires) {
        std::cout << WireID(ent.first).getICLayer() << " " << WireID(ent.first).getIWire() << " " << ent.second << std::endl;
      }
    }
 
    void outputToFile(std::string fileName) const
    {
      std::ofstream fout(fileName);
 
      if (fout.bad()) {
        B2ERROR("Specified output file could not be opened!");
      } else {
        for (auto const ent : m_wires) {
          fout << std::setw(2) << std::right << WireID(ent.first).getICLayer() << "  " << std::setw(3) << WireID(
                 ent.first).getIWire() << " " << ent.second << std::endl;
        }
        fout.close();
      }
    }
 
  private:
    std::map<unsigned short, float> m_wires; 
    ClassDef(CDCBadWires, 2); 
  };
 
} // end namespace Belle2