RecoRunSummary.cc

Go to the documentation of this file.

#include "DQM/HcalTasks/interface/RecoRunSummary.h"

namespace hcaldqm
{
    RecoRunSummary::RecoRunSummary(std::string const& name, 
        std::string const& taskname, edm::ParameterSet const& ps) :
        DQClient(name, taskname, ps)
    {
        _thresh_unihf = ps.getUntrackedParameter<double>("thresh_unihf",
            0.2);
        _thresh_tcds = ps.getUntrackedParameter<double>("thresh_tcds",
            1.5);
    }

    /* virtual */ void RecoRunSummary::beginRun(edm::Run const& r,
        edm::EventSetup const& es)
    {
        DQClient::beginRun(r,es);
    }

    /*
     *
     */
    /* virtual */ void RecoRunSummary::endLuminosityBlock(DQMStore::IBooker& ib,
        DQMStore::IGetter& ig, edm::LuminosityBlock const& lb,
        edm::EventSetup const& es)
    {
        DQClient::endLuminosityBlock(ib, ig, lb, es);
    }

    /*
     *
     */
    /* virtual */ std::vector<flag::Flag> RecoRunSummary::endJob(
        DQMStore::IBooker& ib, DQMStore::IGetter& ig)
    {

        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;
    }
}