Implements edm::EDAnalyzer.
Definition at line 137 of file HLTGetDigi.cc.
References CSCStripdigiCollection_, CSCWiredigiCollection_, DTdigiCollection_, EBdigiCollection_, EEdigiCollection_, ESdigiCollection_, L1GctEtTotal::et(), L1GctEtHad::et(), L1GctEtMiss::et(), etHad, etMiss, GctCaloEmLabel_, GctCaloRegionLabel_, GctCenJetLabel_, GctEtHadLabel_, GctEtMissLabel_, GctEtTotalLabel_, GctForJetLabel_, GctIsoEmLabel_, GctJetCountsLabel_, GctNonIsoEmLabel_, GctTauJetLabel_, edm::EventSetup::get(), edm::Event::getByLabel(), getCSCDigis_, getDTDigis_, getEcalDigis_, getEcalESDigis_, getGctEmDigis_, getGctEtDigis_, getGctJetCounts_, getGctJetDigis_, getGmtCands_, getGmtRC_, getGtEvmRR_, getGtObjectMap_, getGtRR_, getHcalDigis_, getL1Calo_, getPixelDigis_, getRPCDigis_, getStripDigis_, GmtCandsLabel_, GmtReadoutCollection_, GtEvmRRLabel_, GtObjectMapLabel_, GtRRLabel_, HBHEdigiCollection_, HFdigiCollection_, HOdigiCollection_, iggi_31X_cfg::input, iter, LogDebug, muons_cfi::muons, L1GctEtMiss::phi(), edm::Handle< T >::product(), PXLdigiCollection_, RPCdigiCollection_, edm::SortedCollection< T, SORT >::size(), edm::DataFrameContainer::size(), SSTdigiCollection_, and tt.
00138 {
00139 using namespace edm;
00140
00141
00142 edm::Handle<L1GctEtHad> GctEtHad ;
00143 edm::Handle<L1GctEtMiss> GctEtMiss ;
00144 edm::Handle<L1GctEtTotal> GctEtTotal ;
00145
00146 const L1GctEtHad* etHad = 0 ;
00147 const L1GctEtMiss* etMiss = 0 ;
00148 const L1GctEtTotal* etTotal = 0 ;
00149
00150 if (getGctEtDigis_) {
00151 iEvent.getByLabel(GctEtHadLabel_,GctEtHad) ;
00152 iEvent.getByLabel(GctEtMissLabel_,GctEtMiss) ;
00153 iEvent.getByLabel(GctEtTotalLabel_,GctEtTotal) ;
00154 etHad = GctEtHad.product() ;
00155 etMiss = GctEtMiss.product() ;
00156 etTotal = GctEtTotal.product() ;
00157
00158 LogDebug("DigiInfo") << "Value of L1GctEtHad::et(): " << etHad->et() ;
00159 LogDebug("DigiInfo") << "Value of L1GctEtMiss::et(): " << etMiss->et() << ", phi(): " << etMiss->phi() ;
00160 LogDebug("DigiInfo") << "Value of L1GctEtTotal::et(): " << etTotal->et() ;
00161 }
00162
00163 edm::Handle<L1GctEmCandCollection> GctIsoEM ;
00164 edm::Handle<L1GctEmCandCollection> GctNonIsoEM ;
00165
00166 const L1GctEmCandCollection* isoEMdigis = 0 ;
00167 const L1GctEmCandCollection* nonIsoEMdigis = 0 ;
00168 if (getGctEmDigis_) {
00169 iEvent.getByLabel(GctIsoEmLabel_,GctIsoEM) ;
00170 isoEMdigis = GctIsoEM.product() ;
00171 iEvent.getByLabel(GctNonIsoEmLabel_,GctNonIsoEM) ;
00172 nonIsoEMdigis = GctNonIsoEM.product() ;
00173 LogDebug("DigiInfo") << "total # Gct Iso EM digis: " << isoEMdigis->size() ;
00174 LogDebug("DigiInfo") << "total # Gct non-Iso EM digis: " << nonIsoEMdigis->size() ;
00175 }
00176
00177 edm::Handle<L1GctJetCandCollection> GctCenJets ;
00178 edm::Handle<L1GctJetCandCollection> GctForJets ;
00179 edm::Handle<L1GctJetCandCollection> GctTauJets ;
00180 edm::Handle<L1GctJetCounts> GctJetCounts ;
00181
00182 const L1GctJetCandCollection* cenJetDigis = 0 ;
00183 const L1GctJetCandCollection* forJetDigis = 0 ;
00184 const L1GctJetCandCollection* tauJetDigis = 0 ;
00185 std::auto_ptr<L1GctJetCounts> newCounts( new L1GctJetCounts() ) ;
00186 L1GctJetCounts* counts = newCounts.get() ;
00187
00188 if (getGctJetDigis_) {
00189 iEvent.getByLabel(GctCenJetLabel_,GctCenJets) ;
00190 cenJetDigis = GctCenJets.product() ;
00191 iEvent.getByLabel(GctForJetLabel_,GctForJets) ;
00192 forJetDigis = GctForJets.product() ;
00193 iEvent.getByLabel(GctTauJetLabel_,GctTauJets) ;
00194 tauJetDigis = GctTauJets.product() ;
00195 LogDebug("DigiInfo") << "total # Gct central Jet digis: " << cenJetDigis->size() ;
00196 LogDebug("DigiInfo") << "total # Gct forward Jet digis: " << forJetDigis->size() ;
00197 LogDebug("DigiInfo") << "total # Gct tau Jet digis: " << tauJetDigis->size() ;
00198 }
00199
00200 if (getGctJetCounts_) {
00201 iEvent.getByLabel(GctJetCountsLabel_,GctJetCounts) ;
00202 *counts = *GctJetCounts.product() ;
00203 }
00204
00205 edm::Handle<L1CaloEmCollection> GctCaloEm ;
00206 edm::Handle<L1CaloRegionCollection> GctCaloRegion ;
00207
00208 const L1CaloEmCollection* caloEm = 0 ;
00209 const L1CaloRegionCollection* caloRegion = 0 ;
00210
00211 if (getL1Calo_) {
00212 iEvent.getByLabel(GctCaloEmLabel_,GctCaloEm) ;
00213 iEvent.getByLabel(GctCaloRegionLabel_,GctCaloRegion) ;
00214
00215 caloEm = GctCaloEm.product() ;
00216 caloRegion = GctCaloRegion.product() ;
00217
00218 LogDebug("DigiInfo") << "Calo EM size: " << caloEm->size() ;
00219 LogDebug("DigiInfo") << "Calo region size: " << caloRegion->size() ;
00220 }
00221
00222 edm::Handle<L1GlobalTriggerEvmReadoutRecord> gtEvmRR ;
00223 edm::Handle<L1GlobalTriggerObjectMapRecord> gtMap ;
00224 edm::Handle<L1GlobalTriggerReadoutRecord> gtRR ;
00225
00226 edm::ESHandle< L1GtParameters > l1GtPar ;
00227 iSetup.get< L1GtParametersRcd >().get( l1GtPar ) ;
00228 int nBx = l1GtPar->gtTotalBxInEvent() ;
00229
00230 std::auto_ptr<L1GlobalTriggerEvmReadoutRecord> newGtEvm( new L1GlobalTriggerEvmReadoutRecord(nBx) ) ;
00231 std::auto_ptr<L1GlobalTriggerObjectMapRecord> newGtMap( new L1GlobalTriggerObjectMapRecord() ) ;
00232 std::auto_ptr<L1GlobalTriggerReadoutRecord> newGtRR( new L1GlobalTriggerReadoutRecord(nBx) ) ;
00233 L1GlobalTriggerEvmReadoutRecord* evm = newGtEvm.get() ;
00234 L1GlobalTriggerObjectMapRecord* map = newGtMap.get() ;
00235 L1GlobalTriggerReadoutRecord* rr = newGtRR.get() ;
00236
00237 if (getGtEvmRR_) {
00238 iEvent.getByLabel(GtEvmRRLabel_, gtEvmRR) ;
00239 *evm = *gtEvmRR.product() ;
00240 }
00241 if (getGtObjectMap_) {
00242 iEvent.getByLabel(GtObjectMapLabel_, gtMap) ;
00243 *map = *gtMap.product() ;
00244 }
00245 if (getGtRR_) {
00246 iEvent.getByLabel(GtRRLabel_, gtRR) ;
00247 *rr = *gtRR.product() ;
00248 }
00249
00250 edm::Handle< std::vector<L1MuGMTCand> > GmtCands ;
00251 edm::Handle<L1MuGMTReadoutCollection> GmtMuCollection ;
00252 std::auto_ptr<std::vector<L1MuGMTCand> > cands( new std::vector<L1MuGMTCand> ) ;
00253 std::auto_ptr<L1MuGMTReadoutCollection> muCollection(new L1MuGMTReadoutCollection(nBx)) ;
00254
00255 if (getGmtCands_) {
00256 iEvent.getByLabel(GmtCandsLabel_,GmtCands) ;
00257 *cands = *GmtCands.product() ;
00258 }
00259 if (getGmtRC_) {
00260 iEvent.getByLabel(GmtReadoutCollection_,GmtMuCollection) ;
00261 *muCollection = *GmtMuCollection.product() ;
00262 std::vector<L1MuGMTExtendedCand> muons = muCollection->getRecord().getGMTCands() ;
00263 LogDebug("DigiInfo") << "GMT muons present: " << muons.size() ;
00264 }
00265
00266 edm::Handle< DetSetVector<PixelDigi> > input;
00267 auto_ptr<DetSetVector<PixelDigi> > NewPixelDigi(new DetSetVector<PixelDigi> );
00268 DetSetVector<PixelDigi>* tt = NewPixelDigi.get();
00269 if (getPixelDigis_) {
00270 iEvent.getByLabel(PXLdigiCollection_, input);
00271 *tt = *input.product();
00272 }
00273
00274 edm::Handle< edm::DetSetVector<SiStripDigi> > input2;
00275 auto_ptr<DetSetVector<SiStripDigi> > NewSiDigi(new DetSetVector<SiStripDigi> );
00276 DetSetVector<SiStripDigi>* uu = NewSiDigi.get();
00277 if (getStripDigis_) {
00278 iEvent.getByLabel(SSTdigiCollection_, input2);
00279 *uu = *input2.product();
00280 }
00281
00282 Handle<EBDigiCollection> EcalDigiEB;
00283 Handle<EEDigiCollection> EcalDigiEE;
00284 Handle<ESDigiCollection> EcalDigiES;
00285 const EBDigiCollection* EBdigis = 0;
00286 const EEDigiCollection* EEdigis = 0;
00287 const ESDigiCollection* ESdigis = 0;
00288
00289 if (getEcalDigis_) {
00290 iEvent.getByLabel( EBdigiCollection_, EcalDigiEB );
00291 EBdigis = EcalDigiEB.product();
00292 LogDebug("DigiInfo") << "total # EBdigis: " << EBdigis->size() ;
00293
00294 iEvent.getByLabel( EEdigiCollection_, EcalDigiEE );
00295 EEdigis = EcalDigiEE.product();
00296 LogDebug("DigiInfo") << "total # EEdigis: " << EEdigis->size() ;
00297 }
00298
00299 if (getEcalESDigis_) {
00300 iEvent.getByLabel( ESdigiCollection_, EcalDigiES );
00301 ESdigis = EcalDigiES.product();
00302 LogDebug("DigiInfo") << "total # ESdigis: " << ESdigis->size() ;
00303 }
00304
00305 Handle<HBHEDigiCollection> HcalDigiHBHE ;
00306 Handle<HODigiCollection> HcalDigiHO ;
00307 Handle<HFDigiCollection> HcalDigiHF ;
00308 const HBHEDigiCollection* HBHEdigis = 0 ;
00309 const HODigiCollection* HOdigis = 0 ;
00310 const HFDigiCollection* HFdigis = 0 ;
00311
00312 if (getHcalDigis_) {
00313 iEvent.getByLabel( HBHEdigiCollection_, HcalDigiHBHE );
00314 HBHEdigis = HcalDigiHBHE.product();
00315 LogDebug("DigiInfo") << "total # HBHEdigis: " << HBHEdigis->size() ;
00316
00317 iEvent.getByLabel( HOdigiCollection_, HcalDigiHO );
00318 HOdigis = HcalDigiHO.product();
00319 LogDebug("DigiInfo") << "total # HOdigis: " << HOdigis->size() ;
00320
00321 iEvent.getByLabel( HFdigiCollection_, HcalDigiHF );
00322 HFdigis = HcalDigiHF.product();
00323 LogDebug("DigiInfo") << "total # HFdigis: " << HFdigis->size() ;
00324 }
00325
00326 Handle<CSCStripDigiCollection> CSCDigiStrip ;
00327 Handle<CSCWireDigiCollection> CSCDigiWire ;
00328 const CSCStripDigiCollection* CSCstripdigis ;
00329 const CSCWireDigiCollection* CSCwiredigis ;
00330
00331 if (getCSCDigis_) {
00332 iEvent.getByLabel( CSCStripdigiCollection_, CSCDigiStrip );
00333 CSCstripdigis = CSCDigiStrip.product();
00334 iEvent.getByLabel( CSCWiredigiCollection_, CSCDigiWire );
00335 CSCwiredigis = CSCDigiWire.product();
00336
00337 int numDigis = 0 ;
00338 for (CSCStripDigiCollection::DigiRangeIterator iter=CSCDigiStrip->begin();
00339 iter!=CSCDigiStrip->end(); iter++) {
00340 for ( vector<CSCStripDigi>::const_iterator digiIter = (*iter).second.first;
00341 digiIter != (*iter).second.second; digiIter++) numDigis++ ;
00342 }
00343 LogDebug("DigiInfo") << "total # CSCstripdigis: " << numDigis ;
00344 numDigis = 0 ;
00345 for (CSCWireDigiCollection::DigiRangeIterator iter=CSCDigiWire->begin();
00346 iter!=CSCDigiWire->end(); iter++) {
00347 for ( vector<CSCWireDigi>::const_iterator digiIter = (*iter).second.first;
00348 digiIter != (*iter).second.second; digiIter++) numDigis++ ;
00349 }
00350 LogDebug("DigiInfo") << "total # CSCwiredigis: " << numDigis ;
00351 }
00352
00353 Handle<DTDigiCollection> DTDigiHandle ;
00354 const DTDigiCollection* DTdigis ;
00355
00356 if (getDTDigis_) {
00357 iEvent.getByLabel( DTdigiCollection_, DTDigiHandle );
00358 DTdigis = DTDigiHandle.product();
00359
00360 int numDigis = 0 ;
00361 for (DTDigiCollection::DigiRangeIterator iter=DTDigiHandle->begin();
00362 iter!=DTDigiHandle->end(); iter++) {
00363 for ( vector<DTDigi>::const_iterator digiIter = (*iter).second.first;
00364 digiIter != (*iter).second.second; digiIter++) numDigis++ ;
00365 }
00366 LogDebug("DigiInfo") << "total # DTdigis: " << numDigis ;
00367 }
00368
00369 Handle<RPCDigiCollection> RPCDigiHandle ;
00370 const RPCDigiCollection* RPCdigis ;
00371
00372 if (getRPCDigis_) {
00373 iEvent.getByLabel( RPCdigiCollection_, RPCDigiHandle );
00374 RPCdigis = RPCDigiHandle.product();
00375
00376 int numDigis = 0 ;
00377 for (RPCDigiCollection::DigiRangeIterator iter=RPCDigiHandle->begin();
00378 iter!=RPCDigiHandle->end(); iter++) {
00379 for ( vector<RPCDigi>::const_iterator digiIter = (*iter).second.first;
00380 digiIter != (*iter).second.second; digiIter++) numDigis++ ;
00381 }
00382 LogDebug("DigiInfo") << "total # RPCdigis: " << numDigis ;
00383 }
00384
00385 LogDebug("DigiInfo") << "***--------------- End of Event -----------------***" ;
00386
00387 }
1.5.4