CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HLTGetDigi.cc
Go to the documentation of this file.
1 
12 #include "HLTrigger/HLTanalyzers/interface/HLTGetDigi.h"
13 
14 #include "DataFormats/Common/interface/Handle.h"
15 #include "FWCore/Framework/interface/EventSetup.h"
16 #include "FWCore/Framework/interface/ESHandle.h"
17 
18 // system include files
19 #include <memory>
20 #include <vector>
21 #include <map>
22 
23 #include "FWCore/MessageLogger/interface/MessageLogger.h"
24 #include "DataFormats/Common/interface/DetSetVector.h"
25 #include "DataFormats/SiPixelDigi/interface/PixelDigi.h"
26 #include "FWCore/Utilities/interface/InputTag.h"
27 
28 #include "DataFormats/SiStripDigi/interface/SiStripDigi.h"
29 #include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
30 #include "DataFormats/EcalDigi/interface/EBDataFrame.h"
31 #include "DataFormats/EcalDigi/interface/EEDataFrame.h"
32 #include "DataFormats/EcalDigi/interface/ESDataFrame.h"
33 #include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
34 #include "DataFormats/HcalDigi/interface/HBHEDataFrame.h"
35 #include "DataFormats/HcalDigi/interface/HODataFrame.h"
36 #include "DataFormats/HcalDigi/interface/HFDataFrame.h"
37 #include "DataFormats/MuonData/interface/MuonDigiCollection.h"
38 #include "DataFormats/CSCDigi/interface/CSCWireDigi.h"
39 #include "DataFormats/CSCDigi/interface/CSCWireDigiCollection.h"
40 #include "DataFormats/CSCDigi/interface/CSCStripDigi.h"
41 #include "DataFormats/CSCDigi/interface/CSCStripDigiCollection.h"
42 #include "DataFormats/DTDigi/interface/DTDigi.h"
43 #include "DataFormats/DTDigi/interface/DTDigiCollection.h"
44 #include "DataFormats/RPCDigi/interface/RPCDigi.h"
45 #include "DataFormats/RPCDigi/interface/RPCDigiCollection.h"
46 #include "DataFormats/L1GlobalCaloTrigger/interface/L1GctEtSums.h"
47 #include "DataFormats/L1GlobalCaloTrigger/interface/L1GctCollections.h"
48 #include "DataFormats/L1GlobalCaloTrigger/interface/L1GctJetCounts.h"
49 #include "DataFormats/L1CaloTrigger/interface/L1CaloCollections.h"
50 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
51 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerEvmReadoutRecord.h"
52 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapRecord.h"
53 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMap.h"
54 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTReadoutCollection.h"
55 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTCand.h"
56 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTExtendedCand.h"
57 #include "CondFormats/L1TObjects/interface/L1GtParameters.h"
58 #include "CondFormats/DataRecord/interface/L1GtParametersRcd.h"
59 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTCand.h"
60 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTExtendedCand.h"
61 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTReadoutCollection.h"
62 
63 using namespace edm;
64 using namespace std;
65 
66 //
67 // constructors and destructor
68 //
69 HLTGetDigi::HLTGetDigi(const edm::ParameterSet& ps)
70 {
71  EBdigiCollection_ = ps.getParameter<edm::InputTag>("EBdigiCollection");
72  EEdigiCollection_ = ps.getParameter<edm::InputTag>("EEdigiCollection");
73  ESdigiCollection_ = ps.getParameter<edm::InputTag>("ESdigiCollection");
74  HBHEdigiCollection_ = ps.getParameter<edm::InputTag>("HBHEdigiCollection");
75  HOdigiCollection_ = ps.getParameter<edm::InputTag>("HOdigiCollection");
76  HFdigiCollection_ = ps.getParameter<edm::InputTag>("HFdigiCollection");
77  PXLdigiCollection_ = ps.getParameter<edm::InputTag>("SiPixeldigiCollection");
78  SSTdigiCollection_ = ps.getParameter<edm::InputTag>("SiStripdigiCollection");
79  CSCStripdigiCollection_ = ps.getParameter<edm::InputTag>("CSCStripdigiCollection");
80  CSCWiredigiCollection_ = ps.getParameter<edm::InputTag>("CSCWiredigiCollection");
81  DTdigiCollection_ = ps.getParameter<edm::InputTag>("DTdigiCollection");
82  RPCdigiCollection_ = ps.getParameter<edm::InputTag>("RPCdigiCollection");
83 
84  GctCaloEmLabel_ = ps.getParameter<edm::InputTag>("L1CaloEmCollection");
85  GctCaloRegionLabel_ = ps.getParameter<edm::InputTag>("L1CaloRegionCollection");
86 
87  GctIsoEmLabel_ = ps.getParameter<edm::InputTag>("GctIsoEmCollection");
88  GctNonIsoEmLabel_ = ps.getParameter<edm::InputTag>("GctNonIsoEmCollection");
89 
90  GctCenJetLabel_ = ps.getParameter<edm::InputTag>("GctCenJetCollection");
91  GctForJetLabel_ = ps.getParameter<edm::InputTag>("GctForJetCollection");
92  GctTauJetLabel_ = ps.getParameter<edm::InputTag>("GctTauJetCollection");
93  GctJetCountsLabel_ = ps.getParameter<edm::InputTag>("GctJetCounts");
94 
95  GctEtHadLabel_ = ps.getParameter<edm::InputTag>("GctEtHadCollection");
96  GctEtMissLabel_ = ps.getParameter<edm::InputTag>("GctEtMissCollection");
97  GctEtTotalLabel_ = ps.getParameter<edm::InputTag>("GctEtTotalCollection");
98 
99  GtEvmRRLabel_ = ps.getParameter<edm::InputTag>("GtEvmReadoutRecord");
100  GtObjectMapLabel_ = ps.getParameter<edm::InputTag>("GtObjectMapRecord");
101  GtRRLabel_ = ps.getParameter<edm::InputTag>("GtReadoutRecord");
102 
103  GmtCandsLabel_ = ps.getParameter<edm::InputTag>("GmtCands");
104  GmtReadoutCollection_ = ps.getParameter<edm::InputTag>("GmtReadoutCollection");
105 
106  //--- Define which digis we want ---//
107  getEcalDigis_ = ps.getUntrackedParameter<bool>("getEcal",true) ;
108  getEcalESDigis_ = ps.getUntrackedParameter<bool>("getEcalES",true) ;
109  getHcalDigis_ = ps.getUntrackedParameter<bool>("getHcal",true) ;
110  getPixelDigis_ = ps.getUntrackedParameter<bool>("getPixels",true) ;
111  getStripDigis_ = ps.getUntrackedParameter<bool>("getStrips",true) ;
112  getCSCDigis_ = ps.getUntrackedParameter<bool>("getCSC",true) ;
113  getDTDigis_ = ps.getUntrackedParameter<bool>("getDT",true) ;
114  getRPCDigis_ = ps.getUntrackedParameter<bool>("getRPC",true) ;
115  getGctEmDigis_ = ps.getUntrackedParameter<bool>("getGctEm",true) ;
116  getGctJetDigis_ = ps.getUntrackedParameter<bool>("getGctJet",true) ;
117  getGctJetCounts_ = ps.getUntrackedParameter<bool>("getGctJetCounts",true) ;
118  getGctEtDigis_ = ps.getUntrackedParameter<bool>("getGctEt",true) ;
119  getL1Calo_ = ps.getUntrackedParameter<bool>("getL1Calo",true) ;
120  getGtEvmRR_ = ps.getUntrackedParameter<bool>("getGtEvmRR",true) ;
121  getGtObjectMap_ = ps.getUntrackedParameter<bool>("getGtObjectMap",true) ;
122  getGtRR_ = ps.getUntrackedParameter<bool>("getGtReadoutRecord",true) ;
123  getGmtCands_ = ps.getUntrackedParameter<bool>("getGmtCands",true) ;
124  getGmtRC_ = ps.getUntrackedParameter<bool>("getGmtReadout",true) ;
125 
126 }
127 
128 HLTGetDigi::~HLTGetDigi()
129 { }
130 
131 //
132 // member functions
133 //
134 
135 // ------------ method called to produce the data ------------
136 void
137 HLTGetDigi::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
138 {
139  using namespace edm;
140 
141  //--- L1 GCT and GT Digis ---//
142  edm::Handle<L1GctEtHad> GctEtHad ;
143  edm::Handle<L1GctEtMiss> GctEtMiss ;
144  edm::Handle<L1GctEtTotal> GctEtTotal ;
145 
146  const L1GctEtHad* etHad = 0 ;
147  const L1GctEtMiss* etMiss = 0 ;
148  const L1GctEtTotal* etTotal = 0 ;
149 
150  if (getGctEtDigis_) {
151  iEvent.getByLabel(GctEtHadLabel_,GctEtHad) ;
152  iEvent.getByLabel(GctEtMissLabel_,GctEtMiss) ;
153  iEvent.getByLabel(GctEtTotalLabel_,GctEtTotal) ;
154  etHad = GctEtHad.product() ;
155  etMiss = GctEtMiss.product() ;
156  etTotal = GctEtTotal.product() ;
157 
158  LogDebug("DigiInfo") << "Value of L1GctEtHad::et(): " << etHad->et() ;
159  LogDebug("DigiInfo") << "Value of L1GctEtMiss::et(): " << etMiss->et() << ", phi(): " << etMiss->phi() ;
160  LogDebug("DigiInfo") << "Value of L1GctEtTotal::et(): " << etTotal->et() ;
161  }
162 
163  edm::Handle<L1GctEmCandCollection> GctIsoEM ;
164  edm::Handle<L1GctEmCandCollection> GctNonIsoEM ;
165 
166  const L1GctEmCandCollection* isoEMdigis = 0 ;
167  const L1GctEmCandCollection* nonIsoEMdigis = 0 ;
168  if (getGctEmDigis_) {
169  iEvent.getByLabel(GctIsoEmLabel_,GctIsoEM) ;
170  isoEMdigis = GctIsoEM.product() ;
171  iEvent.getByLabel(GctNonIsoEmLabel_,GctNonIsoEM) ;
172  nonIsoEMdigis = GctNonIsoEM.product() ;
173  LogDebug("DigiInfo") << "total # Gct Iso EM digis: " << isoEMdigis->size() ;
174  LogDebug("DigiInfo") << "total # Gct non-Iso EM digis: " << nonIsoEMdigis->size() ;
175  }
176 
177  edm::Handle<L1GctJetCandCollection> GctCenJets ;
178  edm::Handle<L1GctJetCandCollection> GctForJets ;
179  edm::Handle<L1GctJetCandCollection> GctTauJets ;
180  edm::Handle<L1GctJetCounts> GctJetCounts ;
181 
182  const L1GctJetCandCollection* cenJetDigis = 0 ;
183  const L1GctJetCandCollection* forJetDigis = 0 ;
184  const L1GctJetCandCollection* tauJetDigis = 0 ;
185  std::auto_ptr<L1GctJetCounts> newCounts( new L1GctJetCounts() ) ;
186  L1GctJetCounts* counts = newCounts.get() ;
187 
188  if (getGctJetDigis_) {
189  iEvent.getByLabel(GctCenJetLabel_,GctCenJets) ;
190  cenJetDigis = GctCenJets.product() ;
191  iEvent.getByLabel(GctForJetLabel_,GctForJets) ;
192  forJetDigis = GctForJets.product() ;
193  iEvent.getByLabel(GctTauJetLabel_,GctTauJets) ;
194  tauJetDigis = GctTauJets.product() ;
195  LogDebug("DigiInfo") << "total # Gct central Jet digis: " << cenJetDigis->size() ;
196  LogDebug("DigiInfo") << "total # Gct forward Jet digis: " << forJetDigis->size() ;
197  LogDebug("DigiInfo") << "total # Gct tau Jet digis: " << tauJetDigis->size() ;
198  }
199 
200  if (getGctJetCounts_) {
201  iEvent.getByLabel(GctJetCountsLabel_,GctJetCounts) ;
202  *counts = *GctJetCounts.product() ;
203  }
204 
205  edm::Handle<L1CaloEmCollection> GctCaloEm ;
206  edm::Handle<L1CaloRegionCollection> GctCaloRegion ;
207 
208  const L1CaloEmCollection* caloEm = 0 ;
209  const L1CaloRegionCollection* caloRegion = 0 ;
210 
211  if (getL1Calo_) {
212  iEvent.getByLabel(GctCaloEmLabel_,GctCaloEm) ;
213  iEvent.getByLabel(GctCaloRegionLabel_,GctCaloRegion) ;
214 
215  caloEm = GctCaloEm.product() ;
216  caloRegion = GctCaloRegion.product() ;
217 
218  LogDebug("DigiInfo") << "Calo EM size: " << caloEm->size() ;
219  LogDebug("DigiInfo") << "Calo region size: " << caloRegion->size() ;
220  }
221 
222  edm::Handle<L1GlobalTriggerEvmReadoutRecord> gtEvmRR ;
223  edm::Handle<L1GlobalTriggerObjectMapRecord> gtMap ;
224  edm::Handle<L1GlobalTriggerReadoutRecord> gtRR ;
225 
226  edm::ESHandle< L1GtParameters > l1GtPar ;
227  iSetup.get< L1GtParametersRcd >().get( l1GtPar ) ;
228  int nBx = l1GtPar->gtTotalBxInEvent() ;
229 
230  std::auto_ptr<L1GlobalTriggerEvmReadoutRecord> newGtEvm( new L1GlobalTriggerEvmReadoutRecord(nBx) ) ;
231  std::auto_ptr<L1GlobalTriggerObjectMapRecord> newGtMap( new L1GlobalTriggerObjectMapRecord() ) ;
232  std::auto_ptr<L1GlobalTriggerReadoutRecord> newGtRR( new L1GlobalTriggerReadoutRecord(nBx) ) ;
233  L1GlobalTriggerEvmReadoutRecord* evm = newGtEvm.get() ;
234  L1GlobalTriggerObjectMapRecord* map = newGtMap.get() ;
235  L1GlobalTriggerReadoutRecord* rr = newGtRR.get() ;
236 
237  if (getGtEvmRR_) {
238  iEvent.getByLabel(GtEvmRRLabel_, gtEvmRR) ;
239  *evm = *gtEvmRR.product() ;
240  }
241  if (getGtObjectMap_) {
242  iEvent.getByLabel(GtObjectMapLabel_, gtMap) ;
243  *map = *gtMap.product() ;
244  }
245  if (getGtRR_) {
246  iEvent.getByLabel(GtRRLabel_, gtRR) ;
247  *rr = *gtRR.product() ;
248  }
249 
250  edm::Handle< std::vector<L1MuGMTCand> > GmtCands ;
251  edm::Handle<L1MuGMTReadoutCollection> GmtMuCollection ;
252  std::auto_ptr<std::vector<L1MuGMTCand> > cands( new std::vector<L1MuGMTCand> ) ;
253  std::auto_ptr<L1MuGMTReadoutCollection> muCollection(new L1MuGMTReadoutCollection(nBx)) ;
254 
255  if (getGmtCands_) {
256  iEvent.getByLabel(GmtCandsLabel_,GmtCands) ;
257  *cands = *GmtCands.product() ;
258  }
259  if (getGmtRC_) {
260  iEvent.getByLabel(GmtReadoutCollection_,GmtMuCollection) ;
261  *muCollection = *GmtMuCollection.product() ;
262  std::vector<L1MuGMTExtendedCand> muons = muCollection->getRecord().getGMTCands() ;
263  LogDebug("DigiInfo") << "GMT muons present: " << muons.size() ;
264  }
265 
266  edm::Handle< DetSetVector<PixelDigi> > input;
267  auto_ptr<DetSetVector<PixelDigi> > NewPixelDigi(new DetSetVector<PixelDigi> );
268  DetSetVector<PixelDigi>* tt = NewPixelDigi.get();
269  if (getPixelDigis_) {
270  iEvent.getByLabel(PXLdigiCollection_, input);
271  *tt = *input.product();
272  }
273 
274  edm::Handle< edm::DetSetVector<SiStripDigi> > input2;
275  auto_ptr<DetSetVector<SiStripDigi> > NewSiDigi(new DetSetVector<SiStripDigi> );
276  DetSetVector<SiStripDigi>* uu = NewSiDigi.get();
277  if (getStripDigis_) {
278  iEvent.getByLabel(SSTdigiCollection_, input2);
279  *uu = *input2.product();
280  }
281 
282  Handle<EBDigiCollection> EcalDigiEB;
283  Handle<EEDigiCollection> EcalDigiEE;
284  Handle<ESDigiCollection> EcalDigiES;
285  const EBDigiCollection* EBdigis = 0;
286  const EEDigiCollection* EEdigis = 0;
287  const ESDigiCollection* ESdigis = 0;
288 
289  if (getEcalDigis_) {
290  iEvent.getByLabel( EBdigiCollection_, EcalDigiEB );
291  EBdigis = EcalDigiEB.product();
292  LogDebug("DigiInfo") << "total # EBdigis: " << EBdigis->size() ;
293 
294  iEvent.getByLabel( EEdigiCollection_, EcalDigiEE );
295  EEdigis = EcalDigiEE.product();
296  LogDebug("DigiInfo") << "total # EEdigis: " << EEdigis->size() ;
297  }
298 
299  if (getEcalESDigis_) {
300  iEvent.getByLabel( ESdigiCollection_, EcalDigiES );
301  ESdigis = EcalDigiES.product();
302  LogDebug("DigiInfo") << "total # ESdigis: " << ESdigis->size() ;
303  }
304 
305  Handle<HBHEDigiCollection> HcalDigiHBHE ;
306  Handle<HODigiCollection> HcalDigiHO ;
307  Handle<HFDigiCollection> HcalDigiHF ;
308  const HBHEDigiCollection* HBHEdigis = 0 ;
309  const HODigiCollection* HOdigis = 0 ;
310  const HFDigiCollection* HFdigis = 0 ;
311 
312  if (getHcalDigis_) {
313  iEvent.getByLabel( HBHEdigiCollection_, HcalDigiHBHE );
314  HBHEdigis = HcalDigiHBHE.product();
315  LogDebug("DigiInfo") << "total # HBHEdigis: " << HBHEdigis->size() ;
316 
317  iEvent.getByLabel( HOdigiCollection_, HcalDigiHO );
318  HOdigis = HcalDigiHO.product();
319  LogDebug("DigiInfo") << "total # HOdigis: " << HOdigis->size() ;
320 
321  iEvent.getByLabel( HFdigiCollection_, HcalDigiHF );
322  HFdigis = HcalDigiHF.product();
323  LogDebug("DigiInfo") << "total # HFdigis: " << HFdigis->size() ;
324  }
325 
326  Handle<CSCStripDigiCollection> CSCDigiStrip ;
327  Handle<CSCWireDigiCollection> CSCDigiWire ;
328 
329  if (getCSCDigis_) {
330  iEvent.getByLabel( CSCStripdigiCollection_, CSCDigiStrip );
331  iEvent.getByLabel( CSCWiredigiCollection_, CSCDigiWire );
332 
333  int numDigis = 0 ;
334  for (CSCStripDigiCollection::DigiRangeIterator iter=CSCDigiStrip->begin();
335  iter!=CSCDigiStrip->end(); iter++) {
336  for ( vector<CSCStripDigi>::const_iterator digiIter = (*iter).second.first;
337  digiIter != (*iter).second.second; digiIter++) numDigis++ ;
338  }
339  LogDebug("DigiInfo") << "total # CSCstripdigis: " << numDigis ;
340  numDigis = 0 ;
341  for (CSCWireDigiCollection::DigiRangeIterator iter=CSCDigiWire->begin();
342  iter!=CSCDigiWire->end(); iter++) {
343  for ( vector<CSCWireDigi>::const_iterator digiIter = (*iter).second.first;
344  digiIter != (*iter).second.second; digiIter++) numDigis++ ;
345  }
346  LogDebug("DigiInfo") << "total # CSCwiredigis: " << numDigis ;
347  }
348 
349  Handle<DTDigiCollection> DTDigiHandle ;
350 
351  if (getDTDigis_) {
352  iEvent.getByLabel( DTdigiCollection_, DTDigiHandle );
353 
354  int numDigis = 0 ;
355  for (DTDigiCollection::DigiRangeIterator iter=DTDigiHandle->begin();
356  iter!=DTDigiHandle->end(); iter++) {
357  for ( vector<DTDigi>::const_iterator digiIter = (*iter).second.first;
358  digiIter != (*iter).second.second; digiIter++) numDigis++ ;
359  }
360  LogDebug("DigiInfo") << "total # DTdigis: " << numDigis ;
361  }
362 
363  Handle<RPCDigiCollection> RPCDigiHandle ;
364 
365  if (getRPCDigis_) {
366  iEvent.getByLabel( RPCdigiCollection_, RPCDigiHandle );
367 
368  int numDigis = 0 ;
369  for (RPCDigiCollection::DigiRangeIterator iter=RPCDigiHandle->begin();
370  iter!=RPCDigiHandle->end(); iter++) {
371  for ( vector<RPCDigi>::const_iterator digiIter = (*iter).second.first;
372  digiIter != (*iter).second.second; digiIter++) numDigis++ ;
373  }
374  LogDebug("DigiInfo") << "total # RPCdigis: " << numDigis ;
375  }
376 
377  LogDebug("DigiInfo") << "***--------------- End of Event -----------------***" ;
378 
379 }
LogDebug
#define LogDebug(id)
Definition: MessageLogger.h:501
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
L1CaloEmCollection
std::vector< L1CaloEmCand > L1CaloEmCollection
Definition: L1CaloCollections.h:10
LaserDQM_cfg.input
tuple input
Definition: LaserDQM_cfg.py:38
L1GctEtMiss
Persistable copy of missing Et measured at Level-1.
Definition: L1GctEtMiss.h:18
input2
#define input2
Definition: AMPTWrapper.h:149
Association.map
dictionary map
Definition: Association.py:205
L1GctJetCandCollection
std::vector< L1GctJetCand > L1GctJetCandCollection
Definition: L1GctCollections.h:31
iEvent
int iEvent
Definition: GenABIO.cc:243
L1GctEtTotal
Persistable copy of total Et measured at Level-1.
Definition: L1GctEtTotal.h:18
L1GctEtHad
Persistable copy of total Ht measured at Level-1.
Definition: L1GctEtHad.h:18
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:361
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
HLTGetDigi::analyze
void analyze(const edm::Event &, const edm::EventSetup &)
Definition: HLTGetDigi.cc:137
patZpeak.muons
tuple muons
Definition: patZpeak.py:38
HLTGetDigi::HLTGetDigi
HLTGetDigi(const edm::ParameterSet &)
Definition: HLTGetDigi.cc:69
L1CaloRegionCollection
std::vector< L1CaloRegion > L1CaloRegionCollection
Definition: L1CaloCollections.h:11
L1GctEmCandCollection
std::vector< L1GctEmCand > L1GctEmCandCollection
Definition: L1GctCollections.h:30
Generated for CMSSW Reference Manual by doxygen 1.8.5