hcaldqm::RecoRunSummary Class Reference

#include <RecoRunSummary.h>

Inheritance diagram for hcaldqm::RecoRunSummary:

Public Member Functions
void	beginRun (edm::Run const &, edm::EventSetup const &) override
std::vector< flag::Flag >	endJob (DQMStore::IBooker &, DQMStore::IGetter &) override
void	endLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &) override
	RecoRunSummary (std::string const &, std::string const &, edm::ParameterSet const &)
	~RecoRunSummary () override
Public Member Functions inherited from hcaldqm::DQClient
virtual void	beginLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &lb, edm::EventSetup const &)
	DQClient (std::string const &, std::string const &, edm::ParameterSet const &)
	~DQClient () override
Public Member Functions inherited from hcaldqm::DQModule
	DQModule (edm::ParameterSet const &)
virtual	~DQModule ()(false)

Protected Types
enum	RecoFlag { fTCDS = 0, fUniSlotHF = 1, nRecoFlag = 2 }

Protected Attributes
double	_thresh_tcds
double	_thresh_unihf
Protected Attributes inherited from hcaldqm::DQClient
HcalElectronicsMap const *	_emap
int	_maxProcessedLS
std::string	_taskname
int	_totalLS
std::vector< uint32_t >	_vcdaqEids
std::vector< int >	_vCrates
std::vector< int >	_vFEDs
std::vector< uint32_t >	_vhashCrates
std::vector< uint32_t >	_vhashFEDs
ContainerXXX< uint32_t >	_xQuality
Protected Attributes inherited from hcaldqm::DQModule
int	_currentLS
int	_debug
int	_evsPerLS
int	_evsTotal
Logger	_logger
int	_maxLS
ModuleType	_mtype
std::string	_name
ProcessingType	_ptype
std::string	_runkeyName
int	_runkeyVal
std::string	_subsystem

Detailed Description

Definition at line 8 of file RecoRunSummary.h.

Member Enumeration Documentation

enum hcaldqm::RecoRunSummary::RecoFlag

protected

Enumerator
fTCDS
fUniSlotHF
nRecoFlag

Definition at line 23 of file RecoRunSummary.h.

{ fTCDS = 0, fUniSlotHF = 1, nRecoFlag = 2 };

Constructor & Destructor Documentation

hcaldqm::RecoRunSummary::RecoRunSummary	(	std::string const &	name,
		std::string const &	taskname,
		edm::ParameterSet const &	ps
	)

Definition at line 6 of file RecoRunSummary.cc.

References _thresh_tcds, _thresh_unihf, and edm::ParameterSet::getUntrackedParameter().

      : DQClient(name, taskname, ps) {
    _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf", 0.2);
    _thresh_tcds = ps.getUntrackedParameter<double>("thresh_tcds", 1.5);
  }

hcaldqm::RecoRunSummary::~RecoRunSummary ( )

inlineoverride

Definition at line 11 of file RecoRunSummary.h.

References beginRun(), endJob(), and endLuminosityBlock().

{}

Member Function Documentation

void hcaldqm::RecoRunSummary::beginRun ( edm::Run const &  r, edm::EventSetup const &  es  )

overridevirtual

Reimplemented from hcaldqm::DQClient.

Definition at line 12 of file RecoRunSummary.cc.

References hcaldqm::DQClient::beginRun().

Referenced by ~RecoRunSummary().

                                                                          {
    DQClient::beginRun(r, es);
  }

std::vector< flag::Flag > hcaldqm::RecoRunSummary::endJob ( DQMStore::IBooker &  ib, DQMStore::IGetter &  ig  )

overridevirtual

Reimplemented from hcaldqm::DQClient.

Definition at line 29 of file RecoRunSummary.cc.

References hcaldqm::DQClient::_emap, hcaldqm::DQModule::_name, hcaldqm::DQModule::_ptype, hcaldqm::DQModule::_subsystem, hcaldqm::DQClient::_taskname, _thresh_tcds, _thresh_unihf, hcaldqm::DQClient::_vhashCrates, a, HcalElectronicsMap::allPrecisionId(), b, hcaldqm::ContainerXXX< STDTYPE >::begin(), hcaldqm::ContainerXXX< STDTYPE >::book(), hcaldqm::ContainerSingle2D::book(), hcaldqm::Container1D::book(), EnergyCorrector::c, runTauDisplay::eid, hcaldqm::ContainerXXX< STDTYPE >::end(), hcaldqm::flag::fBAD, hcaldqm::hashfunctions::fCrate, hcaldqm::hashfunctions::fCrateSlot, hcaldqm::electronicsmap::fD2EHashMap, hcaldqm::hashfunctions::fdepth, hcaldqm::flag::fGOOD, hcaldqm::hashfunctions::fHBHEPartition, hcaldqm::constants::FIBER_uTCA_MIN1, hcaldqm::constants::FIBERCH_MIN, hcaldqm::quantity::fieta, hcaldqm::quantity::fiphi, nanoDQM_cff::Flag, hcaldqm::quantity::fN, hcaldqm::fOffline, hcaldqm::filter::fPreserver, hcaldqm::quantity::fState, fTCDS, hcaldqm::quantity::fTiming_ns, fUniSlotHF, hcaldqm::ContainerXXX< STDTYPE >::get(), hcaldqm::Container2D::getBinContent(), hcaldqm::Container1D::getMean(), HcalBarrel, HcalEndcap, HcalForward, hcaldqm::ContainerXXX< STDTYPE >::initialize(), hcaldqm::filter::HashFilter::initialize(), hcaldqm::Container2D::initialize(), hcaldqm::ContainerSingle2D::initialize(), hcaldqm::electronicsmap::ElectronicsMap::initialize(), hcaldqm::Container1D::initialize(), hcaldqm::Container1D::load(), hcaldqm::electronicsmap::ElectronicsMap::lookup(), HcalElectronicsMap::lookup(), nRecoFlag, hcaldqm::flag::Flag::reset(), hcaldqm::ContainerSingle2D::setBinContent(), hcaldqm::constants::SLOT_uTCA_MIN, HcalDetId::subdet(), globals_cff::x1, and globals_cff::x2.

Referenced by ~RecoRunSummary().

                                                                                       {
    if (_ptype != fOffline)
      return std::vector<flag::Flag>();

    // FILTERS, some useful vectors, hash maps
    std::vector<uint32_t> vhashCrateHF;
    vhashCrateHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashCrateHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    vhashCrateHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
    filter::HashFilter filter_CrateHF;
    filter_CrateHF.initialize(filter::fPreserver, hashfunctions::fCrate,
                              vhashCrateHF);  // preserve only HF crates
    electronicsmap::ElectronicsMap ehashmap;
    ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
    bool tcdsshift = false;
    std::vector<flag::Flag> vflags;
    vflags.resize(nRecoFlag);
    vflags[fUniSlotHF] = flag::Flag("UniSlotHF");
    vflags[fTCDS] = flag::Flag("TCDS");

    //  INITIALIZE
    Container2D cOccupancy_depth, cOccupancyCut_depth;
    ContainerSingle2D cSummary;
    Container1D cTimingCut_HBHEPartition;
    ContainerXXX<double> xUniHF, xUni;
    xUni.initialize(hashfunctions::fCrate);
    xUniHF.initialize(hashfunctions::fCrateSlot);
    cOccupancy_depth.initialize(_taskname,
                                "Occupancy",
                                hashfunctions::fdepth,
                                new quantity::DetectorQuantity(quantity::fieta),
                                new quantity::DetectorQuantity(quantity::fiphi),
                                new quantity::ValueQuantity(quantity::fN),
                                0);
    cOccupancyCut_depth.initialize(_taskname,
                                   "OccupancyCut",
                                   hashfunctions::fdepth,
                                   new quantity::DetectorQuantity(quantity::fieta),
                                   new quantity::DetectorQuantity(quantity::fiphi),
                                   new quantity::ValueQuantity(quantity::fN),
                                   0);
    cTimingCut_HBHEPartition.initialize(_taskname,
                                        "TimingCut",
                                        hashfunctions::fHBHEPartition,
                                        new quantity::ValueQuantity(quantity::fTiming_ns),
                                        new quantity::ValueQuantity(quantity::fN),
                                        0);

    cSummary.initialize(_name,
                        "Summary",
                        new quantity::CrateQuantity(_emap),
                        new quantity::FlagQuantity(vflags),
                        new quantity::ValueQuantity(quantity::fState),
                        0);

    //  BOOK
    xUniHF.book(_emap, filter_CrateHF);

    //  LOAD
    cOccupancy_depth.load(ig, _emap, _subsystem);
    cOccupancyCut_depth.load(ig, _emap, _subsystem);
    cTimingCut_HBHEPartition.book(ib, _emap, _subsystem);
    cSummary.book(ib, _subsystem);

    //  iterate over all channels
    std::vector<HcalGenericDetId> gids = _emap->allPrecisionId();
    for (std::vector<HcalGenericDetId>::const_iterator it = gids.begin(); it != gids.end(); ++it) {
      if (!it->isHcalDetId())
        continue;

      HcalDetId did(it->rawId());
      HcalElectronicsId eid = HcalElectronicsId(ehashmap.lookup(did));

      if (did.subdet() == HcalForward)
        xUniHF.get(eid) += cOccupancyCut_depth.getBinContent(did);
    }

    //  iphi/slot HF non uniformity
    for (doubleCompactMap::const_iterator it = xUniHF.begin(); it != xUniHF.end(); ++it) {
      uint32_t hash1 = it->first;
      HcalElectronicsId eid1(hash1);
      double x1 = it->second;
      for (doubleCompactMap::const_iterator jt = xUniHF.begin(); jt != xUniHF.end(); ++jt) {
        if (jt == it)
          continue;

        double x2 = jt->second;
        if (x2 == 0)
          continue;
        if (x1 / x2 < _thresh_unihf)
          xUni.get(eid1)++;
      }
    }

    //  TCDS shift
    double a = cTimingCut_HBHEPartition.getMean(HcalDetId(HcalBarrel, 1, 5, 1));
    double b = cTimingCut_HBHEPartition.getMean(HcalDetId(HcalBarrel, 1, 30, 1));
    double c = cTimingCut_HBHEPartition.getMean(HcalDetId(HcalBarrel, 1, 55, 1));
    double dab = fabs(a - b);
    double dac = fabs(a - c);
    double dbc = fabs(b - c);
    if (dab >= _thresh_tcds || dac >= _thresh_tcds || dbc >= _thresh_tcds)
      tcdsshift = true;

    //  summary flags
    std::vector<flag::Flag> sumflags;
    int icrate = 0;
    for (std::vector<uint32_t>::const_iterator it = _vhashCrates.begin(); it != _vhashCrates.end(); ++it) {
      flag::Flag fSum("RECO");
      HcalElectronicsId eid(*it);
      HcalDetId did = HcalDetId(_emap->lookup(eid));

      //    registered @cDAQ
      if (did.subdet() == HcalBarrel || did.subdet() == HcalEndcap) {
        if (tcdsshift)
          vflags[fTCDS]._state = flag::fBAD;
        else
          vflags[fTCDS]._state = flag::fGOOD;
      }
      if (did.subdet() == HcalForward) {
        if (xUni.get(eid) > 0)
          vflags[fUniSlotHF]._state = flag::fBAD;
        else
          vflags[fUniSlotHF]._state = flag::fGOOD;
      }

      //    combine
      int iflag = 0;
      for (std::vector<flag::Flag>::iterator ft = vflags.begin(); ft != vflags.end(); ++ft) {
        cSummary.setBinContent(eid, iflag, ft->_state);
        fSum += (*ft);
        iflag++;
        ft->reset();
      }
      sumflags.push_back(fSum);
      icrate++;
    }

    return sumflags;
  }

void hcaldqm::RecoRunSummary::endLuminosityBlock ( DQMStore::IBooker &  ib, DQMStore::IGetter &  ig, edm::LuminosityBlock const &  lb, edm::EventSetup const &  es  )

overridevirtual

Reimplemented from hcaldqm::DQClient.

Definition at line 19 of file RecoRunSummary.cc.

References hcaldqm::DQClient::endLuminosityBlock().

Referenced by ~RecoRunSummary().

                                                                                 {
    DQClient::endLuminosityBlock(ib, ig, lb, es);
  }

Member Data Documentation

double hcaldqm::RecoRunSummary::_thresh_tcds

protected

Definition at line 21 of file RecoRunSummary.h.

Referenced by endJob(), and RecoRunSummary().

double hcaldqm::RecoRunSummary::_thresh_unihf

protected

Definition at line 21 of file RecoRunSummary.h.

Referenced by endJob(), and RecoRunSummary().