L1TNtuples/src/L1AnalysisRCT.cc

0001 #include "L1Trigger/L1TNtuples/interface/L1AnalysisRCT.h"
0002
0003 // need of maxRCTREG ??
0004
0005 L1Analysis::L1AnalysisRCT::L1AnalysisRCT() {}
0006
0007 L1Analysis::L1AnalysisRCT::L1AnalysisRCT(int maxRCTREG) {
0008   rct_.maxRCTREG_ = maxRCTREG;
0009   rct_.Reset();
0010 }
0011
0012 L1Analysis::L1AnalysisRCT::~L1AnalysisRCT() {}
0013
0014 void L1Analysis::L1AnalysisRCT::SetHdRCT(const edm::Handle<L1CaloRegionCollection> rgn) {
0015   // Regions
0016   rct_.RegSize = rgn->size();
0017   for (L1CaloRegionCollection::const_iterator ireg = rgn->begin(); ireg != rgn->end(); ireg++) {
0018     // local eta phi
0019     rct_.RegEta.push_back(ireg->rctEta());
0020     rct_.RegPhi.push_back(ireg->rctPhi());
0021     // global eta phi
0022     rct_.RegGEta.push_back(ireg->gctEta());
0023     rct_.RegGPhi.push_back(ireg->gctPhi());
0024
0025     rct_.RegRnk.push_back(ireg->et());
0026     rct_.RegVeto.push_back(ireg->tauVeto());
0027     rct_.RegBx.push_back(ireg->bx());
0028     rct_.RegOverFlow.push_back(ireg->overFlow());
0029     rct_.RegMip.push_back(ireg->mip());
0030     rct_.RegFGrain.push_back(ireg->fineGrain());
0031   }
0032 }
0033
0034 void L1Analysis::L1AnalysisRCT::SetEmRCT(const edm::Handle<L1CaloEmCollection> em) {
0035   // Isolated and non-isolated EM
0036   rct_.EmSize = em->size();
0037   for (L1CaloEmCollection::const_iterator emit = em->begin(); emit != em->end(); emit++) {
0038     rct_.IsIsoEm.push_back(emit->isolated());
0039     rct_.EmEta.push_back(emit->regionId().ieta());
0040     rct_.EmPhi.push_back(emit->regionId().iphi());
0041     rct_.EmRnk.push_back(emit->rank());
0042     rct_.EmBx.push_back(emit->bx());
0043   }
0044 }