DTHVCheckByAbsoluteValues.cc

Go to the documentation of this file.

/*
 *  See header file for a description of this class.
 *
 *  $Date: 2010/04/15 20:38:42 $
 *  $Revision: 1.4 $
 *  \author Paolo Ronchese INFN Padova
 *
 */
 
//-----------------------
// This Class' Header --
//-----------------------
#include "CondTools/DT/plugins/DTHVCheckByAbsoluteValues.h"
 
//-------------------------------
// Collaborating Class Headers --
//-------------------------------
#include "DataFormats/MuonDetId/interface/DTWireId.h"
#include "FWCore/PluginManager/interface/ModuleDef.h"
#include "FWCore/ServiceRegistry/interface/ServiceMaker.h"
 
//---------------
// C++ Headers --
//---------------
#include <iostream>
 
namespace cond {
  namespace service {
 
    //-------------------
    // Initializations --
    //-------------------
 
    //----------------
    // Constructors --
    //----------------
    DTHVCheckByAbsoluteValues::DTHVCheckByAbsoluteValues(const edm::ParameterSet& iConfig, edm::ActivityRegistry& iAR) {
      if (instance == nullptr) {
        std::cout << "create DTHVCheckByAbsoluteValues" << std::endl;
        minHV = new float[4];
        maxHV = new float[4];
        //    minHV[0] = 3500.0;
        //    minHV[1] = 3500.0;
        //    minHV[2] = 1400.0;
        //    minHV[3] =  800.0;
        minHV[0] = 3500.0;
        minHV[1] = 3500.0;
        minHV[2] = 1700.0;
        minHV[3] = 1100.0;
        maxHV[0] = 4000.0;
        maxHV[1] = 4000.0;
        maxHV[2] = 2200.0;
        maxHV[3] = 1600.0;
        maxCurrent = 30.0;
        instance = this;
      }
    }
 
    //--------------
    // Destructor --
    //--------------
    DTHVCheckByAbsoluteValues::~DTHVCheckByAbsoluteValues() {}
 
    //--------------
    // Operations --
    //--------------
    DTHVAbstractCheck::flag DTHVCheckByAbsoluteValues::checkCurrentStatus(
        int rawId,
        int type,
        float valueA,
        float valueC,
        float valueS,
        const std::map<int, timedMeasurement>& snapshotValues,
        const std::map<int, int>& aliasMap,
        const std::map<int, int>& layerMap) {
      // find all values for this channel
      //  ind dpid = 0;
      //  std::map<int,int>::const_iterator lpartIter;
      //  std::map<int,int>::const_iterator lpartIend = layerMap.end();
      //  if ( ( layerIter = layerMap.find( chp0 ) ) != layerIend )
      //      dpid = layerIter.second;
      //  std::map<int,timedMeasurement>::const_iterator snapIter;
      //  std::map<int,timedMeasurement>::const_iterator snapIend =
      //                                                 snapshotValues.end();
      //  float val1 = -999999.0;
      //  float val2 = -999999.0;
      //  int chan = dpId * 10;
      //  if ( ( snapIter = snapshotValues.find( chan + 1 ) ) != snapIend )
      //      val1 = snapIter->second.second;
      //  if ( ( snapIter = snapshotValues.find( chan + 2 ) ) != snapIend )
      //      val2 = snapIter->second.second;
 
      // find dp identifier for all channels in this layer
      //  DTLayerId lay = chlId.layerId();
      //  int chp0 = DTWireId( lay, 10 ).rawId();
      //  int chp1 = DTWireId( lay, 11 ).rawId();
      //  int chp2 = DTWireId( lay, 12 ).rawId();
      //  int chp3 = DTWireId( lay, 13 ).rawId();
      //  ind dpi0 = 0;
      //  ind dpi1 = 0;
      //  ind dpi2 = 0;
      //  ind dpi3 = 0;
      //  std::map<int,int>::const_iterator layerIter;
      //  std::map<int,int>::const_iterator layerIend = layerMap.end();
      //  if ( ( layerIter = layerMap.find( chp0 ) ) != layerIend )
      //      dpi0 = layerIter.second;
      //  if ( ( layerIter = layerMap.find( chp1 ) ) != layerIend )
      //      dpi1 = layerIter.second;
      //  if ( ( layerIter = layerMap.find( chp2 ) ) != layerIend )
      //      dpi2 = layerIter.second;
      //  if ( ( layerIter = layerMap.find( chp3 ) ) != layerIend )
      //      dpi3 = layerIter.second;
 
      DTWireId chlId(rawId);
      int part = chlId.wire() - 10;
      DTHVAbstractCheck::flag flag;
      flag.a = flag.c = flag.s = 0;
      if (type == 1) {
        if (valueA < minHV[part])
          flag.a = 1;
        if (valueA > maxHV[part])
          flag.a = 2;
        if (valueS < minHV[2])
          flag.s = 1;
        if (valueS > maxHV[2])
          flag.s = 2;
        if (valueC < minHV[3])
          flag.c = 1;
        if (valueC > maxHV[3])
          flag.c = 2;
      }
      if (type == 2) {
        float voltA = 0.0;
        float voltS = 0.0;
        float voltC = 0.0;
        DTLayerId lay = chlId.layerId();
        int l_p = chlId.wire();
        DTWireId chA(lay, l_p);
        DTWireId chS(lay, 12);
        DTWireId chC(lay, 13);
        std::map<int, int>::const_iterator layerIter;
        std::map<int, int>::const_iterator layerIend = layerMap.end();
        std::map<int, timedMeasurement>::const_iterator snapIter;
        std::map<int, timedMeasurement>::const_iterator snapIend = snapshotValues.end();
        int chan;
        if ((layerIter = layerMap.find(chA.rawId())) != layerIend) {
          chan = (layerIter->second * 10) + l_p;
          if ((snapIter = snapshotValues.find(chan)) != snapIend) {
            voltA = snapIter->second.second;
          }
        }
        if ((layerIter = layerMap.find(chS.rawId())) != layerIend) {
          chan = (layerIter->second * 10) + 2;
          if ((snapIter = snapshotValues.find(chan)) != snapIend) {
            voltS = snapIter->second.second;
          }
        }
        if ((layerIter = layerMap.find(chC.rawId())) != layerIend) {
          chan = (layerIter->second * 10) + 3;
          if ((snapIter = snapshotValues.find(chan)) != snapIend) {
            voltC = snapIter->second.second;
          }
        }
        if ((valueA > maxCurrent) && (voltA >= minHV[part]))
          flag.a = 4;
        if ((valueS > maxCurrent) && (voltS >= minHV[2]))
          flag.s = 4;
        if ((valueC > maxCurrent) && (voltC >= minHV[3]))
          flag.c = 4;
      }
      return flag;
    }
 
    DEFINE_FWK_SERVICE(DTHVCheckByAbsoluteValues);
  }  // namespace service
}  // namespace cond