#include <QIE10Task.h>

Inheritance diagram for QIE10Task:

Public Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

	QIE10Task (edm::ParameterSet const &)

	~QIE10Task () override

Public Member Functions inherited from hcaldqm::DQTask
void	analyze (edm::Event const &, edm::EventSetup const &) override

void	beginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (edm::Run const &, edm::EventSetup const &) override

	DQTask (edm::ParameterSet const &)

void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override

	~DQTask () override

Public Member Functions inherited from one::DQMEDAnalyzer< one::DQMLuminosityBlockElements >
	DQMEDAnalyzer ()=default

	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete

	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete

	~DQMEDAnalyzer () override=default

Public Member Functions inherited from hcaldqm::DQModule
	DQModule (edm::ParameterSet const &)

virtual	~DQModule ()(false)

Protected Member Functions
void	_process (edm::Event const &, edm::EventSetup const &) override

void	_resetMonitors (hcaldqm::UpdateFreq) override

Protected Member Functions inherited from hcaldqm::DQTask
virtual int	_getCalibType (edm::Event const &)

virtual bool	_isApplicable (edm::Event const &)

Protected Attributes
hcaldqm::ContainerSingle1D	_cADC

hcaldqm::Container1D	_cADC_EChannel [10][36]

hcaldqm::ContainerSingle1D	_cLETDC

hcaldqm::Container1D	_cLETDC_EChannel [10][36]

hcaldqm::Container1D	_cLETDCTime_EChannel [36]

hcaldqm::ContainerSingle2D	_cLETDCTimevsADC

hcaldqm::ContainerSingle2D	_cLETDCvsADC

hcaldqm::Container2D	_cLETDCvsADC_EChannel [10][36]

hcaldqm::Container2D	_cLETDCvsTS_EChannel [36]

hcaldqm::Container2D	_cOccupancy_Crate

hcaldqm::Container2D	_cOccupancy_CrateSlot

hcaldqm::Container2D	_cOccupancy_depth

hcaldqm::ContainerSingleProf1D	_cShapeCut

hcaldqm::ContainerProf1D	_cShapeCut_EChannel [36]

double	_cut

hcaldqm::electronicsmap::ElectronicsMap	_ehashmap

hcaldqm::filter::HashFilter	_filter_slot [36]

int	_ped

edm::InputTag	_tagHF

edm::InputTag	_tagQIE10

edm::EDGetTokenT< HFDigiCollection >	_tokHF

edm::EDGetTokenT< QIE10DigiCollection >	_tokQIE10

Protected Attributes inherited from hcaldqm::DQTask
ContainerI	_cEvsPerLS

ContainerI	_cEvsTotal

ContainerS	_cProcessingTypeName

ContainerS	_cRunKeyName

ContainerI	_cRunKeyVal

edm::ESHandle< HcalDbService >	_dbService

HcalElectronicsMap const *	_emap

int	_procLSs

edm::InputTag	_tagRaw

edm::EDGetTokenT< FEDRawDataCollection >	_tokRaw

std::vector< uint32_t >	_vcdaqEids

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 22 of file QIE10Task.h.

Constructor & Destructor Documentation

QIE10Task::QIE10Task ( edm::ParameterSet const & ps )

Definition at line 8 of file QIE10Task.cc.

References _cut, _ped, _tagHF, _tagQIE10, _tokHF, _tokQIE10, and edm::ParameterSet::getUntrackedParameter().

                                               : DQTask(ps) {
   //    tags
   _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagQIE10", edm::InputTag("hcalDigis"));
   _tagHF = ps.getUntrackedParameter<edm::InputTag>("tagHF", edm::InputTag("hcalDigis"));
   _tokQIE10 = consumes<QIE10DigiCollection>(_tagQIE10);
   _tokHF = consumes<HFDigiCollection>(_tagHF);
 
   //    cuts
   _cut = ps.getUntrackedParameter<double>("cut", 50.0);
   _ped = ps.getUntrackedParameter<int>("ped", 4);
 }

QIE10Task::~QIE10Task ( )

inlineoverride

Definition at line 25 of file QIE10Task.h.

References _process(), _resetMonitors(), bookHistograms(), and endLuminosityBlock().

25 {}

Member Function Documentation

void QIE10Task::_process	(	edm::Event const &	e,
		edm::EventSetup const &
	)

overrideprotectedvirtual

Implements hcaldqm::DQTask.

Definition at line 162 of file QIE10Task.cc.

References _cADC, _cADC_EChannel, _cLETDC, _cLETDC_EChannel, _cLETDCTime_EChannel, _cLETDCTimevsADC, _cLETDCvsADC, _cLETDCvsADC_EChannel, _cLETDCvsTS_EChannel, _cOccupancy_Crate, _cOccupancy_CrateSlot, _cOccupancy_depth, _cShapeCut, _cShapeCut_EChannel, _cut, hcaldqm::DQTask::_dbService, _ehashmap, hcaldqm::DQModule::_logger, _tagHF, _tokHF, _tokQIE10, ecalMGPA::adc(), HcalElectronicsId::crateId(), QIE10DataFrame::detid(), Logger::dqmthrow(), runTauDisplay::eid, hcaldqm::ContainerSingle1D::fill(), hcaldqm::Container2D::fill(), hcaldqm::ContainerSingle2D::fill(), hcaldqm::Container1D::fill(), edm::Event::getByToken(), HcalForward, mps_fire::i, edm::InputTag::instance(), edm::InputTag::label(), hcaldqm::electronicsmap::ElectronicsMap::lookup(), lumiQueryAPI::q, QIE10DataFrame::samples(), edm::DataFrameContainer::size(), HcalElectronicsId::slot(), HcalDetId::subdet(), hcaldqm::utilities::sumQ(), and ntuplemaker::time.

Referenced by ~QIE10Task().

                                                                   {
   edm::Handle<QIE10DigiCollection> cqie10;
   edm::Handle<HFDigiCollection> chf;
   if (!e.getByToken(_tokQIE10, cqie10))
     return;
   if (!e.getByToken(_tokHF, chf))
     _logger.dqmthrow("Collection HFDigiCollection isn't available" + _tagHF.label() + " " + _tagHF.instance());
 
   std::map<uint32_t, QIE10DataFrame> mqie10;
   for (uint32_t i = 0; i < cqie10->size(); i++) {
     QIE10DataFrame frame = static_cast<QIE10DataFrame>((*cqie10)[i]);
     HcalDetId did = frame.detid();
     HcalElectronicsId eid = HcalElectronicsId(_ehashmap.lookup(did));
     if (did.subdet() != HcalForward) {
       continue;
     }
 
     // Compute index for EChannel plots
     int fakecrate = -1;
     if (eid.crateId() == 22) {
       fakecrate = 0;
     } else if (eid.crateId() == 29) {
       fakecrate = 1;
     } else if (eid.crateId() == 32) {
       fakecrate = 2;
     } else {
       // Unknown crate, skip digi
       continue;
     }
     int index = fakecrate * 12 + eid.slot() - 1;
 
     //  compute the signal, ped subracted
     CaloSamples digi_fC = hcaldqm::utilities::loadADC2fCDB<QIE10DataFrame>(_dbService, did, frame);
     double sumQ = hcaldqm::utilities::sumQDB<QIE10DataFrame>(_dbService, digi_fC, did, frame, 0, frame.samples() - 1);
     // double sumQ = hcaldqm::utilities::sumQ_v10<QIE10DataFrame>(frame, constants::adc2fC[_ped], 0, frame.samples()-1);
 
     _cOccupancy_Crate.fill(eid);
     _cOccupancy_CrateSlot.fill(eid);
     _cOccupancy_depth.fill(did);
 
     //  iterate thru all TS and fill
     for (int j = 0; j < frame.samples(); j++) {
       //    shapes are after the cut
       if (sumQ > _cut) {
         double q = hcaldqm::utilities::adc2fCDBMinusPedestal<QIE10DataFrame>(_dbService, digi_fC, did, frame, j);
         _cShapeCut_EChannel[index].fill(eid, j, q);
         _cShapeCut.fill(j, q);
       }
       //    w/o a cut
       _cLETDCvsADC_EChannel[j][index].fill(eid, frame[j].adc(), frame[j].le_tdc());
       _cLETDCvsADC.fill(frame[j].adc(), frame[j].le_tdc());
       _cLETDC_EChannel[j][index].fill(eid, frame[j].le_tdc());
       _cLETDC.fill(frame[j].le_tdc());
       _cADC_EChannel[j][index].fill(eid, frame[j].adc());
       _cADC.fill(frame[j].adc());
       _cLETDCvsTS_EChannel[index].fill(eid, j, frame[j].le_tdc());
 
       // TDC conversion to time
       if (frame[j].le_tdc() < 50) {
         // Each TDC count is 0.5 ns.
         // tdc == 62 or 63 means value was below or above threshold for whole time slice.
         double time = j * 25. + (frame[j].le_tdc() / 2.);
         _cLETDCTime_EChannel[index].fill(eid, time);
         _cLETDCTimevsADC.fill(frame[j].adc(), time);
       }
     }
 
     mqie10[did.rawId()] = frame;
   }
 }

void QIE10Task::_resetMonitors ( hcaldqm::UpdateFreq )

overrideprotectedvirtual

Reimplemented from hcaldqm::DQTask.

Definition at line 233 of file QIE10Task.cc.

References DEFINE_FWK_MODULE.

Referenced by ~QIE10Task().

233 {}

void QIE10Task::bookHistograms	(	DQMStore::IBooker &	ib,
		edm::Run const &	r,
		edm::EventSetup const &	es
	)

override

Definition at line 19 of file QIE10Task.cc.

Referenced by ~QIE10Task().

                                                                                               {
   if (_ptype == fLocal)
     if (r.runAuxiliary().run() == 1)
       return;
 
   DQTask::bookHistograms(ib, r, es);
 
   //    GET WHAT YOU NEED
   edm::ESHandle<HcalDbService> dbs;
   es.get<HcalDbRecord>().get(dbs);
   _emap = dbs->getHcalMapping();
 
   unsigned int nTS = _ptype == fLocal ? 10 : 6;
 
   // create a slot filter and initialize what you need
   unsigned int itr = 0;
   for (auto& crate : constants::crateListHF) {
     for (unsigned int slot = SLOT_uTCA_MIN; slot <= SLOT_uTCA_MAX; ++slot) {
       std::vector<uint32_t> vhashSlot;
       vhashSlot.push_back(HcalElectronicsId(crate, slot, FIBER_uTCA_MIN1, FIBERCH_MIN, false).rawId());
       _filter_slot[itr].initialize(filter::fPreserver, hashfunctions::fCrateSlot, vhashSlot);
 
       _cShapeCut_EChannel[itr].initialize(_name,
                                           "ShapeCut",
                                           hcaldqm::hashfunctions::fEChannel,
                                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
                                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000));
       _cLETDCTime_EChannel[itr].initialize(_name,
                                            "LETDCTime",
                                            hcaldqm::hashfunctions::fEChannel,
                                            new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
                                            new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
       _cLETDCvsTS_EChannel[itr].initialize(_name,
                                            "TDCvsTS",
                                            hcaldqm::hashfunctions::fEChannel,
                                            new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
                                            new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
                                            new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
 
       for (unsigned int j = 0; j < nTS; j++) {
         _cLETDCvsADC_EChannel[j][itr].initialize(_name,
                                                  "LETDCvsADC",
                                                  hcaldqm::hashfunctions::fEChannel,
                                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
                                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
                                                  new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
         _cADC_EChannel[j][itr].initialize(_name,
                                           "ADC",
                                           hcaldqm::hashfunctions::fEChannel,
                                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
                                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
         _cLETDC_EChannel[j][itr].initialize(_name,
                                             "LETDC",
                                             hcaldqm::hashfunctions::fEChannel,
                                             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
                                             new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
       }
       ++itr;
     }
   }
 
   _cShapeCut.initialize(_name,
                         "ShapeCut",
                         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTiming_TS),
                         new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10fC_400000));
   _cLETDCvsADC.initialize(_name,
                           "LETDCvsADC",
                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
                           new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
   _cLETDCTimevsADC.initialize(_name,
                               "LETDCTimevsADC",
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fTime_ns_250),
                               new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
   _cLETDC.initialize(_name,
                      "LETDC",
                      new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10TDC_64),
                      new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
   _cADC.initialize(_name,
                    "ADC",
                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fQIE10ADC_256),
                    new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN, true));
 
   //OCCUPANCY PER CRATE/SLOT
   _cOccupancy_Crate.initialize(_name,
                                "Occupancy",
                                hashfunctions::fCrate,
                                new quantity::ElectronicsQuantity(quantity::fSlotuTCA),
                                new quantity::ElectronicsQuantity(quantity::fFiberuTCAFiberCh),
                                new quantity::ValueQuantity(quantity::fN));
   _cOccupancy_CrateSlot.initialize(_name,
                                    "Occupancy",
                                    hashfunctions::fCrateSlot,
                                    new quantity::ElectronicsQuantity(quantity::fFiberuTCA),
                                    new quantity::ElectronicsQuantity(quantity::fFiberCh),
                                    new quantity::ValueQuantity(quantity::fN));
 
   // OCCUPANCY IN DETECTOR COORDINATES
   _cOccupancy_depth.initialize(_name,
                                "Occupancy",
                                hcaldqm::hashfunctions::fdepth,
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fieta),
                                new hcaldqm::quantity::DetectorQuantity(hcaldqm::quantity::fiphi),
                                new hcaldqm::quantity::ValueQuantity(hcaldqm::quantity::fN));
 
   itr = 0;
   for (auto& crate : constants::crateListHF) {
     for (unsigned int slot = SLOT_uTCA_MIN; slot <= SLOT_uTCA_MAX; ++slot) {
       char aux[100];
       sprintf(aux, "/Crate%d_Slot%d", crate, slot);
       _cShapeCut_EChannel[itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux);
       _cLETDCTime_EChannel[itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux);
       _cLETDCvsTS_EChannel[itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux);
       for (unsigned int i = 0; i < nTS; i++) {
         char aux2[100];
         sprintf(aux2, "/Crate%d_Slot%d/TS%d", crate, slot, i);
         _cLETDCvsADC_EChannel[i][itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux);
         _cLETDC_EChannel[i][itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux);
         _cADC_EChannel[i][itr].book(ib, _emap, _filter_slot[itr], _subsystem, aux);
       }
       ++itr;
     }
   }
 
   _cShapeCut.book(ib, _subsystem);
   _cLETDCvsADC.book(ib, _subsystem);
   _cLETDCTimevsADC.book(ib, _subsystem);
   _cLETDC.book(ib, _subsystem);
   _cADC.book(ib, _subsystem);
 
   _cOccupancy_Crate.book(ib, _emap, _subsystem);
   _cOccupancy_CrateSlot.book(ib, _emap, _subsystem);
   _cOccupancy_depth.book(ib, _emap, _subsystem);
 
   _ehashmap.initialize(_emap, electronicsmap::fD2EHashMap);
 }

void QIE10Task::endLuminosityBlock	(	edm::LuminosityBlock const &	lb,
		edm::EventSetup const &	es
	)

override

Definition at line 157 of file QIE10Task.cc.

Referenced by ~QIE10Task().

                                                                                           {
   //    finish
   DQTask::endLuminosityBlock(lb, es);
 }