139 int firstBX = old_egammas->getFirstBX();
140 int lastBX = old_egammas->getLastBX();
142 new_egammas->setBXRange(firstBX, lastBX);
143 new_rlxtaus->setBXRange(firstBX, lastBX);
144 new_isotaus->setBXRange(firstBX, lastBX);
145 new_jets->setBXRange(firstBX, lastBX);
146 new_preGtJets->setBXRange(firstBX, lastBX);
147 new_etsums->setBXRange(firstBX, lastBX);
148 new_hfsums->setBXRange(firstBX, lastBX);
149 new_hfcounts->setBXRange(firstBX, lastBX);
151 for (
int bx = firstBX;
bx <= lastBX; ++
bx) {
155 const double et = emScale->et(itEGamma->hwPt());
156 const double eta = getPhysicalEta(itEGamma->hwEta());
157 const double phi = getPhysicalPhi(itEGamma->hwPhi());
160 EGamma eg(*&p4, itEGamma->hwPt(), itEGamma->hwEta(), itEGamma->hwPhi(), itEGamma->hwQual(), itEGamma->hwIso());
161 new_egammas->push_back(
bx, *&eg);
173 const double et = jetScale->et(itTau->hwPt());
175 const double eta = getPhysicalEta(itTau->hwEta());
176 const double phi = getPhysicalPhi(itTau->hwPhi());
179 Tau tau(*&p4, itTau->hwPt(), itTau->hwEta(), itTau->hwPhi(), itTau->hwQual(), itTau->hwIso());
180 new_rlxtaus->push_back(
bx, *&
tau);
192 const double et = jetScale->et(itTau->hwPt());
194 const double eta = getPhysicalEta(itTau->hwEta());
195 const double phi = getPhysicalPhi(itTau->hwPhi());
198 Tau tau(*&p4, itTau->hwPt(), itTau->hwEta(), itTau->hwPhi(), itTau->hwQual(), itTau->hwIso());
199 new_isotaus->push_back(
bx, *&
tau);
211 const double et = jetScale->et(itJet->hwPt());
213 const bool forward = ((itJet->hwQual() & 0x2) != 0);
214 const double eta = getPhysicalEta(itJet->hwEta(), forward);
215 const double phi = getPhysicalPhi(itJet->hwPhi());
218 Jet jet(*&p4, itJet->hwPt(), itJet->hwEta(), itJet->hwPhi(), itJet->hwQual());
219 new_jets->push_back(
bx, *&
jet);
228 const double et = itJet->hwPt() * emScale->linearLsb();
233 const bool forward = ((itJet->hwQual() & 0x2) != 0);
234 const double eta = getPhysicalEta(itJet->hwEta(), forward);
235 const double phi = getPhysicalPhi(itJet->hwPhi());
238 Jet jet(*&p4, itJet->hwPt(), itJet->hwEta(), itJet->hwPhi(), itJet->hwQual());
239 new_preGtJets->push_back(
bx, *&
jet);
243 double et = itEtSum->hwPt() * emScale->linearLsb();
248 const double eta = getPhysicalEta(itEtSum->hwEta());
249 double phi = getPhysicalPhi(itEtSum->hwPhi());
251 et = htMissScale->et(itEtSum->hwPt());
253 phi = phi + (regionPhiWidth / 2.);
258 EtSum eg(*&p4, sumType, itEtSum->hwPt(), itEtSum->hwEta(), itEtSum->hwPhi(), itEtSum->hwQual());
259 new_etsums->push_back(
bx, *&eg);
266 new_hfsums->push_back(
bx, *itCaloSpare);
270 itCaloSpare != old_hfcounts->end(
bx);
274 new_hfcounts->push_back(
bx, *itCaloSpare);
edm::EDGetToken JetToken_
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
edm::EDGetToken HfSumsToken_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::ESGetToken< L1CaloEtScale, L1HtMissScaleRcd > htMissScaleToken_
edm::EDGetToken HfCountsToken_
edm::EDGetToken EtSumToken_
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
std::vector< T >::const_iterator const_iterator
edm::EDGetToken EGammaToken_
edm::EDGetToken RlxTauToken_
edm::EDGetToken preGtJetToken_
edm::EDGetToken IsoTauToken_
static const unsigned N_PHI
ESHandle< T > getHandle(const ESGetToken< T, R > &iToken) const
edm::ESGetToken< L1CaloEtScale, L1EmEtScaleRcd > emScaleToken_
edm::ESGetToken< L1CaloEtScale, L1JetEtScaleRcd > jetScaleToken_