ZDCTask Class Reference

#include <ZDCTask.h>

Inheritance diagram for ZDCTask:

Public Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	ZDCTask (edm::ParameterSet const &)
	~ZDCTask () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Public Attributes
	cut
	debug
	name
	ped
	runkeyName
	runkeyVal
	tagQIE10
	zdcQIE10Task
	zdcTask

Protected Member Functions
void	analyze (edm::Event const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes
MonitorElement *	_cADC
std::map< std::string, MonitorElement * >	_cADC_EChannel
MonitorElement *	_cADC_vs_TS
std::map< std::string, MonitorElement * >	_cADC_vs_TS_EChannel
MonitorElement *	_cShape
std::map< std::string, MonitorElement * >	_cShape_EChannel
double	_cut
int	_ped
edm::InputTag	_tagQIE10
edm::EDGetTokenT< ZDCDigiCollection >	_tokQIE10
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T... >
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Detailed Description

Definition at line 22 of file ZDCTask.h.

Constructor & Destructor Documentation

ZDCTask()

ZDCTask::ZDCTask

(

edm::ParameterSet const &

ps

)

Definition at line 7 of file ZDCTask.cc.

References edm::ParameterSet::getUntrackedParameter(), and ProducerED_cfi::InputTag.

                                          {
  //    tags
  _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagQIE10", edm::InputTag("hcalDigis"));
  _tokQIE10 = consumes<ZDCDigiCollection>(_tagQIE10);

  //    cuts
  _cut = ps.getUntrackedParameter<double>("cut", 50.0);
  _ped = ps.getUntrackedParameter<int>("ped", 4);
}

~ZDCTask()

ZDCTask::~ZDCTask ( )

inlineoverride

Definition at line 25 of file ZDCTask.h.

{}

Member Function Documentation

analyze()

void ZDCTask::analyze ( edm::Event const &  e, edm::EventSetup const &    )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 181 of file ZDCTask.cc.

References gpuClustering::adc, MillePedeFileConverter_cfg::e, runTauDisplay::eid, mps_fire::i, dqmiolumiharvest::j, and edm::SortedCollection< T, SORT >::size().

                                                               {
  edm::Handle<ZDCDigiCollection> cqie10;
  if (!e.getByToken(_tokQIE10, cqie10))
    edm::LogError("Collection ZDCDigiCollection isn't available" + _tagQIE10.label() + " " + _tagQIE10.instance());

  for (uint32_t i = 0; i < cqie10->size(); i++) {
    ZDCDataFrame frame = static_cast<ZDCDataFrame>((*cqie10)[i]);
    HcalElectronicsId eid = frame.elecId();

    char histoname[300];
    sprintf(histoname, "%d_%d_%d_%d", eid.fiberChanId(), eid.fiberIndex(), eid.spigot(), eid.dccid());

    //  compute the signal, ped subracted
    //double q = hcaldqm::utilities::sumQ_v10<ZDCDataFrame>(frame, constants::adc2fC[_ped], 0, frame.size()-1);

    //  iterate thru all TS and fill
    for (int j = 0; j < frame.size(); j++) {
      _cShape_EChannel[histoname]->Fill(j, frame[j].nominal_fC());
      _cShape->Fill(j, frame[j].nominal_fC());

      _cADC_EChannel[histoname]->Fill(frame[j].adc());
      _cADC->Fill(frame[j].adc());

      _cADC_vs_TS_EChannel[histoname]->Fill(j, frame[j].adc());
      _cADC_vs_TS->Fill(j, frame[j].adc());
    }
  }
}

bookHistograms()

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 17 of file ZDCTask.cc.

References hcaldqm::constants::BIT_OFFSET, runTauDisplay::eid, hcaldqm::quantity::fADC_128, hcaldqm::quantity::ffC_10000, hcaldqm::quantity::fTiming_TS, hcaldqm::quantity::fYAxis, cuy::ib, hcaldqm::quantity::ValueQuantity::max(), hcaldqm::quantity::ValueQuantity::min(), and hcaldqm::quantity::ValueQuantity::nbins().

                                                                                            {
  //############################## hardcode manually the zdc mapping #############################
  //############################# this follows from https://github.com/cms-sw/cmssw/blob/CMSSW_8_0_X/EventFilter/CastorRawToDigi/src/ZdcUnpacker.cc#L118
  //##############################################################################################
  //PZDC
  std::map<HcalElectronicsId, DetId> myEMap;
  HcalElectronicsId eid = HcalElectronicsId(0, 1, 0, 3);
  eid.setHTR(18, 8, 1);
  myEMap[eid] = DetId(0x54000051);  //PZDC EM1

  eid = HcalElectronicsId(1, 1, 0, 3);
  eid.setHTR(18, 8, 1);
  myEMap[eid] = DetId(0x54000052);  //PZDC EM2

  eid = HcalElectronicsId(2, 1, 0, 3);
  eid.setHTR(18, 8, 1);
  myEMap[eid] = DetId(0x54000053);  //PZDC EM3

  eid = HcalElectronicsId(0, 2, 0, 3);
  eid.setHTR(18, 8, 1);
  myEMap[eid] = DetId(0x54000061);  //PZDC HAD1

  eid = HcalElectronicsId(1, 2, 0, 3);
  eid.setHTR(18, 8, 1);
  myEMap[eid] = DetId(0x54000054);  //PZDC EM4

  eid = HcalElectronicsId(2, 2, 0, 3);
  eid.setHTR(18, 8, 1);
  myEMap[eid] = DetId(0x54000055);  //PZDC EM5

  eid = HcalElectronicsId(0, 3, 0, 3);
  eid.setHTR(18, 8, 1);
  myEMap[eid] = DetId(0x54000062);  //PZDC HAD2

  eid = HcalElectronicsId(1, 3, 0, 3);
  eid.setHTR(18, 8, 1);
  myEMap[eid] = DetId(0x54000063);  //PZDC HAD3

  eid = HcalElectronicsId(2, 3, 0, 3);
  eid.setHTR(18, 8, 1);
  myEMap[eid] = DetId(0x54000064);  //PZDC HAD4

  //NZDC
  eid = HcalElectronicsId(0, 1, 1, 3);
  eid.setHTR(18, 8, 0);
  myEMap[eid] = DetId(0x54000011);  //NZDC EM1

  eid = HcalElectronicsId(1, 1, 1, 3);
  eid.setHTR(18, 8, 0);
  myEMap[eid] = DetId(0x54000012);  //NZDC EM2

  eid = HcalElectronicsId(2, 1, 1, 3);
  eid.setHTR(18, 8, 0);
  myEMap[eid] = DetId(0x54000013);  //NZDC EM3

  eid = HcalElectronicsId(0, 2, 1, 3);
  eid.setHTR(18, 8, 0);
  myEMap[eid] = DetId(0x54000015);  //NZDC EM5

  eid = HcalElectronicsId(1, 2, 1, 3);
  eid.setHTR(18, 8, 0);
  myEMap[eid] = DetId(0x54000021);  //NZDC HAD1

  eid = HcalElectronicsId(2, 2, 1, 3);
  eid.setHTR(18, 8, 0);
  myEMap[eid] = DetId(0x54000014);  //NZDC EM4

  eid = HcalElectronicsId(0, 3, 1, 3);
  eid.setHTR(18, 8, 0);
  myEMap[eid] = DetId(0x54000022);  //NZDC HAD2

  eid = HcalElectronicsId(1, 3, 1, 3);
  eid.setHTR(18, 8, 0);
  myEMap[eid] = DetId(0x54000023);  //NZDC HAD3

  eid = HcalElectronicsId(2, 3, 1, 3);
  eid.setHTR(18, 8, 0);
  myEMap[eid] = DetId(0x54000024);  //NZDC HAD4
  //##################################### end hardcoding ###################################

  ib.cd();

  //quantities for axis
  hcaldqm::quantity::ValueQuantity xAxisShape(hcaldqm::quantity::fTiming_TS);
  hcaldqm::quantity::ValueQuantity yAxisShape(hcaldqm::quantity::ffC_10000);

  hcaldqm::quantity::ValueQuantity xAxisADC(hcaldqm::quantity::fADC_128);

  //book histos per channel
  for (std::map<HcalElectronicsId, DetId>::const_iterator itr = myEMap.begin(); itr != myEMap.end(); ++itr) {
    char histoname[300];

    sprintf(histoname,
            "%d_%d_%d_%d",
            itr->first.fiberChanId(),
            itr->first.fiberIndex(),
            itr->first.spigot(),
            itr->first.dccid());

    ib.setCurrentFolder("Hcal/ZDCTask/Shape_perChannel");
    _cShape_EChannel[histoname] = ib.bookProfile(histoname,
                                                 histoname,
                                                 xAxisShape.nbins(),
                                                 xAxisShape.min(),
                                                 xAxisShape.max(),
                                                 yAxisShape.nbins(),
                                                 yAxisShape.min(),
                                                 yAxisShape.max());
    _cShape_EChannel[histoname]->setAxisTitle("Timing", 1);
    _cShape_EChannel[histoname]->setAxisTitle("fC QIE8", 2);

    ib.setCurrentFolder("Hcal/ZDCTask/ADC_perChannel");
    _cADC_EChannel[histoname] = ib.book1DD(histoname, histoname, xAxisADC.nbins(), xAxisADC.min(), xAxisADC.max());
    _cADC_EChannel[histoname]->getTH1()->SetBit(
        BIT(hcaldqm::constants::BIT_OFFSET + hcaldqm::quantity::AxisType::fYAxis));
    _cADC_EChannel[histoname]->setAxisTitle("ADC QIE8", 1);

    ib.setCurrentFolder("Hcal/ZDCTask/ADC_vs_TS_perChannel");
    _cADC_vs_TS_EChannel[histoname] = ib.book2D(histoname,
                                                histoname,
                                                xAxisShape.nbins(),
                                                xAxisShape.min(),
                                                xAxisShape.max(),
                                                xAxisADC.nbins(),
                                                xAxisADC.min(),
                                                xAxisADC.max());
    _cADC_vs_TS_EChannel[histoname]->getTH1()->SetBit(
        BIT(hcaldqm::constants::BIT_OFFSET + hcaldqm::quantity::AxisType::fYAxis));
    _cADC_vs_TS_EChannel[histoname]->setAxisTitle("Timing", 1);
    _cADC_vs_TS_EChannel[histoname]->setAxisTitle("ADC QIE8", 2);
  }

  //book global histos
  ib.setCurrentFolder("Hcal/ZDCTask");

  _cShape = ib.bookProfile("Shape",
                           "Shape",
                           xAxisShape.nbins(),
                           xAxisShape.min(),
                           xAxisShape.max(),
                           yAxisShape.nbins(),
                           yAxisShape.min(),
                           yAxisShape.max());
  _cShape->setAxisTitle("Timing", 1);
  _cShape->setAxisTitle("fC QIE8", 2);

  _cADC = ib.book1DD("ADC", "ADC", xAxisADC.nbins(), xAxisADC.min(), xAxisADC.max());
  _cADC->getTH1()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET + hcaldqm::quantity::AxisType::fYAxis));
  _cADC->setAxisTitle("ADC QIE8", 1);

  _cADC_vs_TS = ib.book2D("ADC_vs_TS",
                          "ADC_vs_TS",
                          xAxisShape.nbins(),
                          xAxisShape.min(),
                          xAxisShape.max(),
                          xAxisADC.nbins(),
                          xAxisADC.min(),
                          xAxisADC.max());
  _cADC_vs_TS->getTH1()->SetBit(BIT(hcaldqm::constants::BIT_OFFSET + hcaldqm::quantity::AxisType::fYAxis));
  _cADC_vs_TS->setAxisTitle("Timing", 1);
  _cADC_vs_TS->setAxisTitle("ADC QIE8", 2);
}

Member Data Documentation

_cADC

MonitorElement* ZDCTask::_cADC

protected

Definition at line 47 of file ZDCTask.h.

_cADC_EChannel

std::map<std::string, MonitorElement*> ZDCTask::_cADC_EChannel

protected

Definition at line 42 of file ZDCTask.h.

_cADC_vs_TS

MonitorElement* ZDCTask::_cADC_vs_TS

protected

Definition at line 48 of file ZDCTask.h.

_cADC_vs_TS_EChannel

std::map<std::string, MonitorElement*> ZDCTask::_cADC_vs_TS_EChannel

protected

Definition at line 43 of file ZDCTask.h.

_cShape

MonitorElement* ZDCTask::_cShape

protected

Definition at line 46 of file ZDCTask.h.

_cShape_EChannel

std::map<std::string, MonitorElement*> ZDCTask::_cShape_EChannel

protected

Definition at line 41 of file ZDCTask.h.

_cut

double ZDCTask::_cut

protected

Definition at line 37 of file ZDCTask.h.

_ped

int ZDCTask::_ped

protected

Definition at line 38 of file ZDCTask.h.

_tagQIE10

edm::InputTag ZDCTask::_tagQIE10

protected

Definition at line 33 of file ZDCTask.h.

_tokQIE10

edm::EDGetTokenT<ZDCDigiCollection> ZDCTask::_tokQIE10

protected

Definition at line 34 of file ZDCTask.h.

cut

ZDCTask.cut

Definition at line 17 of file ZDCTask.py.

Referenced by postprocessor.PostProcessor.run().

debug

ZDCTask.debug

Definition at line 9 of file ZDCTask.py.

Referenced by util.rrapi.RRApi.dprint(), rrapi.RRApi.dprint(), pkg.AbstractPkg.generate(), util.rrapi.RRApi.get(), rrapi.RRApi.get(), pkg.AbstractPkg.get_kwds(), runTauIdMVA.TauIDEmbedder.loadMVA_WPs_run2_2017(), runTauIdMVA.TauIDEmbedder.runTauID(), and pkg.AbstractPkg.write().

name

ZDCTask.name

Definition at line 8 of file ZDCTask.py.

Referenced by ElectronMVAID.ElectronMVAID.__call__(), FWLite.ElectronMVAID.__call__(), dirstructure.Directory.__create_pie_image(), DisplayManager.DisplayManager.__del__(), dqm_interfaces.DirID.__eq__(), dirstructure.Directory.__get_full_path(), dirstructure.Comparison.__get_img_name(), dirstructure.Comparison.__make_image(), core.autovars.NTupleVariable.__repr__(), core.autovars.NTupleObjectType.__repr__(), core.autovars.NTupleObject.__repr__(), core.autovars.NTupleCollection.__repr__(), dirstructure.Directory.__repr__(), dqm_interfaces.DirID.__repr__(), dirstructure.Comparison.__repr__(), config.Service.__setattr__(), config.CFG.__str__(), counter.Counter.__str__(), average.Average.__str__(), FWLite.WorkingPoints._reformat_cut_definitions(), core.autovars.NTupleObjectType.addSubObjects(), core.autovars.NTupleObjectType.addVariables(), core.autovars.NTupleObjectType.allVars(), dataset.CMSDataset.buildListOfFiles(), dataset.LocalDataset.buildListOfFiles(), dataset.CMSDataset.buildListOfFilesDBS(), dirstructure.Directory.calcStats(), validation.Sample.digest(), python.rootplot.utilities.Hist.divide(), python.rootplot.utilities.Hist.divide_wilson(), DisplayManager.DisplayManager.Draw(), TreeCrawler.Package.dump(), core.autovars.NTupleVariable.fillBranch(), core.autovars.NTupleObject.fillBranches(), core.autovars.NTupleCollection.fillBranchesScalar(), core.autovars.NTupleCollection.fillBranchesVector(), core.autovars.NTupleCollection.get_cpp_declaration(), core.autovars.NTupleCollection.get_cpp_wrapper_class(), core.autovars.NTupleCollection.get_py_wrapper_class(), utils.StatisticalTest.get_status(), production_tasks.Task.getname(), dataset.CMSDataset.getPrimaryDatasetEntries(), dataset.PrivateDataset.getPrimaryDatasetEntries(), VIDSelectorBase.VIDSelectorBase.initialize(), personalPlayback.Applet.log(), core.autovars.NTupleVariable.makeBranch(), core.autovars.NTupleObject.makeBranches(), core.autovars.NTupleCollection.makeBranchesScalar(), core.autovars.NTupleCollection.makeBranchesVector(), dirstructure.Directory.print_report(), dataset.BaseDataset.printInfo(), dataset.Dataset.printInfo(), production_tasks.MonitorJobs.run(), python.rootplot.utilities.Hist.TGraph(), python.rootplot.utilities.Hist.TH1F(), counter.Counter.write(), and average.Average.write().

ped

ZDCTask.ped

Definition at line 18 of file ZDCTask.py.

runkeyName

ZDCTask.runkeyName

Definition at line 11 of file ZDCTask.py.

runkeyVal

ZDCTask.runkeyVal

Definition at line 10 of file ZDCTask.py.

tagQIE10

ZDCTask.tagQIE10

Definition at line 14 of file ZDCTask.py.

zdcQIE10Task

ZDCTask.zdcQIE10Task

Definition at line 21 of file ZDCTask.py.

zdcTask

ZDCTask.zdcTask

Definition at line 4 of file ZDCTask.py.