hcaldqm::RecoRunSummary Class Reference

#include <RecoRunSummary.h>

Inheritance diagram for hcaldqm::RecoRunSummary:

Public Member Functions
virtual void	beginRun (edm::Run const &, edm::EventSetup const &)
virtual std::vector< flag::Flag >	endJob (DQMStore::IBooker &, DQMStore::IGetter &)
virtual void	endLuminosityBlock (DQMStore::IBooker &, DQMStore::IGetter &, edm::LuminosityBlock const &, edm::EventSetup const &)
	RecoRunSummary (std::string const &, std::string const &, edm::ParameterSet const &)
virtual	~RecoRunSummary ()
Public Member Functions inherited from hcaldqm::DQClient
	DQClient (std::string const &, std::string const &, edm::ParameterSet const &)
virtual	~DQClient ()
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 >	_vFEDs
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 9 of file RecoRunSummary.h.

Member Enumeration Documentation

enum hcaldqm::RecoRunSummary::RecoFlag

protected

Enumerator
fTCDS
fUniSlotHF
nRecoFlag

Definition at line 26 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 5 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);
    }

virtual hcaldqm::RecoRunSummary::~RecoRunSummary ( )

inlinevirtual

Definition at line 14 of file RecoRunSummary.h.

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

{}

Member Function Documentation

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

virtual

Reimplemented from hcaldqm::DQClient.

Definition at line 15 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  )

virtual

Reimplemented from hcaldqm::DQClient.

Definition at line 34 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::_vcdaqEids, hcaldqm::DQClient::_vFEDs, hcaldqm::DQClient::_vhashFEDs, a, HcalElectronicsMap::allPrecisionId(), b, hcaldqm::ContainerXXX< STDTYPE >::begin(), hcaldqm::ContainerXXX< STDTYPE >::book(), hcaldqm::ContainerSingle2D::book(), hcaldqm::Container1D::book(), EnergyCorrector::c, hcaldqm::constants::CRATE_uTCA_MIN, hcaldqm::constants::CRATE_VME_MIN, runTauDisplay::eid, hcaldqm::ContainerXXX< STDTYPE >::end(), hcaldqm::flag::fBAD, hcaldqm::electronicsmap::fD2EHashMap, hcaldqm::hashfunctions::fdepth, hcaldqm::hashfunctions::fElectronics, hcaldqm::hashfunctions::fFED, hcaldqm::hashfunctions::fFEDSlot, hcaldqm::filter::fFilter, hcaldqm::flag::fGOOD, hcaldqm::hashfunctions::fHBHEPartition, hcaldqm::constants::FIBER_uTCA_MIN1, hcaldqm::constants::FIBER_VME_MIN, hcaldqm::constants::FIBERCH_MIN, hcaldqm::quantity::fieta, spr::find(), hcaldqm::quantity::fiphi, hcaldqm::quantity::fN, hcaldqm::flag::fNCDAQ, 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, HcalForward, hcaldqm::Container2D::initialize(), hcaldqm::ContainerXXX< STDTYPE >::initialize(), hcaldqm::ContainerSingle2D::initialize(), hcaldqm::filter::HashFilter::initialize(), hcaldqm::electronicsmap::ElectronicsMap::initialize(), hcaldqm::Container1D::initialize(), hcaldqm::utilities::isFEDHBHE(), hcaldqm::utilities::isFEDHF(), hcaldqm::Container1D::load(), hcaldqm::electronicsmap::ElectronicsMap::lookup(), nRecoFlag, hcaldqm::flag::Flag::reset(), hcaldqm::ContainerSingle2D::setBinContent(), hcaldqm::constants::SLOT_uTCA_MIN, and hcaldqm::constants::SPIGOT_MIN.

Referenced by ~RecoRunSummary().

    {

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

        // FILTERS, some useful vectors, hash maps
        std::vector<uint32_t> vhashFEDHF;
        vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN,
            FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
        vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN,
            FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
        vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN,
            FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
        vhashFEDHF.push_back(HcalElectronicsId(22, SLOT_uTCA_MIN+6,
            FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
        vhashFEDHF.push_back(HcalElectronicsId(29, SLOT_uTCA_MIN+6,
            FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
        vhashFEDHF.push_back(HcalElectronicsId(32, SLOT_uTCA_MIN+6,
            FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
        filter::HashFilter filter_FEDHF;
        filter_FEDHF.initialize(filter::fPreserver, hashfunctions::fFED,
            vhashFEDHF);    // preserve only HF FEDs
        electronicsmap::ElectronicsMap ehashmap;
        ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
        bool tcdsshift = false;
        filter::HashFilter filter_VME;
        filter::HashFilter filter_uTCA;
        std::vector<uint32_t> vhashVME,vhashuTCA;
        vhashVME.push_back(HcalElectronicsId(constants::FIBERCH_MIN,
            constants::FIBER_VME_MIN, SPIGOT_MIN, CRATE_VME_MIN).rawId());
        vhashuTCA.push_back(HcalElectronicsId(CRATE_uTCA_MIN, SLOT_uTCA_MIN,
            FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
        filter_uTCA.initialize(filter::fFilter, hashfunctions::fElectronics,
            vhashuTCA);
        filter_VME.initialize(filter::fFilter, hashfunctions::fElectronics,
            vhashVME);
        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::fFED);
        xUniHF.initialize(hashfunctions::fFEDSlot);
        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::FEDQuantity(_vFEDs),
            new quantity::FlagQuantity(vflags),
            new quantity::ValueQuantity(quantity::fState),0);

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


        //  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 ifed=0;
        for (std::vector<uint32_t>::const_iterator it=_vhashFEDs.begin();
            it!=_vhashFEDs.end(); ++it)
        {
            flag::Flag fSum("RECO");
            HcalElectronicsId eid(*it);

            std::vector<uint32_t>::const_iterator cit=std::find(
                _vcdaqEids.begin(), _vcdaqEids.end(),*it);
            if (cit==_vcdaqEids.end())
            {
                //  not registered @cDAQ
                sumflags.push_back(flag::Flag("RECO", flag::fNCDAQ));
                ifed++;
                continue;
            }

            //  registered @cDAQ
            if (utilities::isFEDHBHE(eid))
            {
                if (tcdsshift)
                    vflags[fTCDS]._state = flag::fBAD;
                else
                    vflags[fTCDS]._state = flag::fGOOD;
            }
            if (utilities::isFEDHF(eid))
            {
                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);
            ifed++;
        }

        return sumflags;
    }

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

virtual

Reimplemented from hcaldqm::DQClient.

Definition at line 24 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 24 of file RecoRunSummary.h.

Referenced by endJob(), and RecoRunSummary().

double hcaldqm::RecoRunSummary::_thresh_unihf

protected

Definition at line 24 of file RecoRunSummary.h.

Referenced by endJob(), and RecoRunSummary().