CMS 3D CMS Logo

FCDTask.cc
Go to the documentation of this file.
1 
3 
4 bool operator==(const FCDTask::FCDChannel& lhs, const FCDTask::FCDChannel& rhs) {
5  return ((lhs.crate == rhs.crate) && (lhs.slot == rhs.slot) && (lhs.fiber == rhs.fiber) &&
6  (lhs.fiberChannel == rhs.fiberChannel));
7 }
8 
10  // tags
11  _tagQIE10 = ps.getUntrackedParameter<edm::InputTag>("tagQIE10", edm::InputTag("hcalDigis", "ZDC"));
12  _tokQIE10 = consumes<QIE10DigiCollection>(_tagQIE10);
13 
14  // channels
15  edm::ParameterSet channelPSet = ps.getParameter<edm::ParameterSet>("fcdChannels");
16  std::vector<int32_t> crates = channelPSet.getUntrackedParameter<std::vector<int32_t> >("crate");
17  std::vector<int32_t> slots = channelPSet.getUntrackedParameter<std::vector<int32_t> >("slot");
18  std::vector<int32_t> fibers = channelPSet.getUntrackedParameter<std::vector<int32_t> >("fiber");
19  std::vector<int32_t> fiberChannels = channelPSet.getUntrackedParameter<std::vector<int32_t> >("fiber_channel");
20  for (unsigned int i = 0; i < crates.size(); ++i) {
21  _channels.push_back({crates[i], slots[i], fibers[i], fiberChannels[i]});
22  }
23 }
24 
25 /* virtual */ void FCDTask::bookHistograms(DQMStore::IBooker& ib, edm::Run const& r, edm::EventSetup const& es) {
27  es.get<HcalDbRecord>().get(dbService);
28  _emap = dbService->getHcalMapping();
30 
31  ib.cd();
32 
33  //book histos per channel
34  std::string histoname;
35  std::vector<HcalGenericDetId> gids = _emap->allPrecisionId();
36  for (auto& it_gid : gids) {
37  if (it_gid.genericSubdet() != HcalGenericDetId::HcalGenZDC) {
38  continue;
39  }
41  for (auto& it_channel : _channels) {
42  if ((eid.crateId() == it_channel.crate) && (eid.slot() == it_channel.slot) &&
43  (eid.fiberIndex() == it_channel.fiber) && (eid.fiberChanId() == it_channel.fiberChannel)) {
44  _fcd_eids.push_back(eid);
45  }
46  }
47  }
48  for (auto& it_eid : _fcd_eids) {
49  // EM Pos
50  histoname = std::to_string(it_eid.crateId()) + "-" + std::to_string(it_eid.slot()) + "-" +
51  std::to_string(it_eid.fiberIndex()) + "-" + std::to_string(it_eid.fiberChanId());
52  ib.setCurrentFolder("Hcal/FCDTask/ADC");
53  _cADC[it_eid] = ib.book1D(histoname.c_str(), histoname.c_str(), 256, 0, 256);
54  _cADC[it_eid]->setAxisTitle("ADC", 1);
55  _cADC[it_eid]->setAxisTitle("N", 2);
56 
57  ib.setCurrentFolder("Hcal/FCDTask/ADC_vs_TS"),
58  _cADC_vs_TS[it_eid] = ib.book2D(histoname.c_str(), histoname.c_str(), 10, 0, 10, 64, 0, 256);
59  _cADC_vs_TS[it_eid]->setAxisTitle("TS", 1);
60  _cADC_vs_TS[it_eid]->setAxisTitle("ADC", 2);
61 
62  ib.setCurrentFolder("Hcal/FCDTask/TDCTime");
63  _cTDCTime[it_eid] = ib.book1D(histoname.c_str(), histoname.c_str(), 500, 0., 250.);
64  _cTDCTime[it_eid]->setAxisTitle("TDC time [ns]", 1);
65 
66  ib.setCurrentFolder("Hcal/FCDTask/TDC");
67  _cTDC[it_eid] = ib.book1D(histoname.c_str(), histoname.c_str(), 64, -0.5, 63.5);
68  _cTDC[it_eid]->setAxisTitle("TDC", 1);
69  }
70 }
71 
72 /* virtual */ void FCDTask::analyze(edm::Event const& e, edm::EventSetup const&) {
74  if (!e.getByToken(_tokQIE10, digis))
75  edm::LogError("Collection QIE10DigiCollection for ZDC isn't available" + _tagQIE10.label() + " " +
77 
78  for (auto it = digis->begin(); it != digis->end(); it++) {
79  const QIE10DataFrame digi = static_cast<const QIE10DataFrame>(*it);
80  HcalGenericDetId const& gdid = digi.detid();
82  if (std::find(_fcd_eids.begin(), _fcd_eids.end(), eid) == _fcd_eids.end()) {
83  continue;
84  }
85 
86  for (int i = 0; i < digi.samples(); i++) {
87  // iter over all samples
88  _cADC[eid]->Fill(digi[i].adc());
89  _cADC_vs_TS[eid]->Fill(i, digi[i].adc());
90  _cTDC[eid]->Fill(digi[i].le_tdc());
91  if (digi[i].le_tdc() <= 50.) {
92  double tdctime = 25. * i + 0.5 * digi[i].le_tdc();
93  _cTDCTime[eid]->Fill(tdctime);
94  }
95  }
96  }
97 }
98 
FCDTask.h
FCDTask::_cADC_vs_TS
std::map< HcalElectronicsId, MonitorElement * > _cADC_vs_TS
Definition: FCDTask.h:37
mps_fire.i
i
Definition: mps_fire.py:428
hcaldqm::electronicsmap::ElectronicsMap::lookup
uint32_t lookup(DetId const &)
Definition: ElectronicsMap.cc:122
FCDTask::FCDChannel::fiber
int fiber
Definition: FCDTask.h:18
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
HcalGenericDetId
Definition: HcalGenericDetId.h:15
QIE10DataFrame::samples
constexpr int samples() const
total number of samples in the digi
Definition: QIE10DataFrame.h:86
edm::Run
Definition: Run.h:45
FCDTask::_channels
std::vector< FCDChannel > _channels
Definition: FCDTask.h:42
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89287
hcaldqm::electronicsmap::ElectronicsMap::initialize
void initialize(HcalElectronicsMap const *, ElectronicsMapType etype=fHcalElectronicsMap)
Definition: ElectronicsMap.cc:6
FCDTask::FCDChannel::slot
int slot
Definition: FCDTask.h:17
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
spr::find
void find(edm::Handle< EcalRecHitCollection > &hits, DetId thisDet, std::vector< EcalRecHitCollection::const_iterator > &hit, bool debug=false)
Definition: FindCaloHit.cc:19
edm::Handle
Definition: AssociativeIterator.h:50
FCDTask::_tagQIE10
edm::InputTag _tagQIE10
Definition: FCDTask.h:32
HcalDbService::getHcalMapping
const HcalElectronicsMap * getHcalMapping() const
Definition: HcalDbService.cc:343
FCDTask::FCDChannel
Definition: FCDTask.h:15
FCDTask::FCDChannel::crate
int crate
Definition: FCDTask.h:16
ecalLiteDTU::adc
constexpr int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalLiteDTUSample.h:12
HcalElectronicsMap::allPrecisionId
std::vector< HcalGenericDetId > allPrecisionId() const
Definition: HcalElectronicsMap.cc:139
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
HcalGenericDetId::HcalGenZDC
Definition: HcalGenericDetId.h:24
edm::EventSetup::get
T get() const
Definition: EventSetup.h:80
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
QIE10DataFrame::detid
constexpr DetId detid() const
Get the detector id.
Definition: QIE10DataFrame.h:76
operator==
bool operator==(const FCDTask::FCDChannel &lhs, const FCDTask::FCDChannel &rhs)
Definition: FCDTask.cc:4
FCDTask::_cTDC
std::map< HcalElectronicsId, MonitorElement * > _cTDC
Definition: FCDTask.h:38
edm::ESHandle< HcalDbService >
HcalElectronicsId
Readout chain identification for Hcal.
Definition: HcalElectronicsId.h:32
FCDTask::analyze
void analyze(edm::Event const &, edm::EventSetup const &) override
Definition: FCDTask.cc:72
FCDTask::_tokQIE10
edm::EDGetTokenT< QIE10DigiCollection > _tokQIE10
Definition: FCDTask.h:33
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
FCDTask::_cADC
std::map< HcalElectronicsId, MonitorElement * > _cADC
Definition: FCDTask.h:36
edm::ParameterSet
Definition: ParameterSet.h:47
FCDTask::_emap
HcalElectronicsMap const * _emap
Definition: FCDTask.h:43
runTauDisplay.eid
eid
Definition: runTauDisplay.py:298
edm::DataFrameContainer::begin
const_iterator begin() const
The iterator returned can not safely be used across threads.
Definition: DataFrameContainer.h:149
cuy.ib
ib
Definition: cuy.py:662
edm::EventSetup
Definition: EventSetup.h:57
get
#define get
FCDTask::_ehashmap
hcaldqm::electronicsmap::ElectronicsMap _ehashmap
Definition: FCDTask.h:44
FCDTask::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: FCDTask.cc:25
alignCSCRings.r
r
Definition: alignCSCRings.py:93
FCDTask::FCDTask
FCDTask(edm::ParameterSet const &)
Definition: FCDTask.cc:9
FCDTask
Definition: FCDTask.py:1
hcaldqm::electronicsmap::fD2EHashMap
Definition: ElectronicsMap.h:25
QIE10DataFrame
Definition: QIE10DataFrame.h:11
FCDTask::FCDChannel::fiberChannel
int fiberChannel
Definition: FCDTask.h:19
FCDTask::_cTDCTime
std::map< HcalElectronicsId, MonitorElement * > _cTDCTime
Definition: FCDTask.h:39
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
dqm::implementation::IBooker
Definition: DQMStore.h:43
HcalDbRecord
Definition: HcalDbRecord.h:30
edm::Event
Definition: Event.h:73
edm::Log
Definition: MessageLogger.h:70
FCDTask::_fcd_eids
std::vector< HcalElectronicsId > _fcd_eids
Definition: FCDTask.h:41
edm::DataFrameContainer::end
const_iterator end() const
Definition: DataFrameContainer.h:152
edm::InputTag
Definition: InputTag.h:15
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37