159 new_egammas->setBXRange(firstBX, lastBX);
160 new_rlxtaus->setBXRange(firstBX, lastBX);
161 new_isotaus->setBXRange(firstBX, lastBX);
162 new_jets->setBXRange(firstBX, lastBX);
163 new_preGtJets->setBXRange(firstBX, lastBX);
164 new_etsums->setBXRange(firstBX, lastBX);
165 new_hfsums->setBXRange(firstBX, lastBX);
166 new_hfcounts->setBXRange(firstBX, lastBX);
168 for(
int bx = firstBX; bx <= lastBX; ++bx)
171 itEGamma != old_egammas->
end(bx); ++itEGamma)
174 const double et = emScale->
et( itEGamma->hwPt() );
175 const double eta = getPhysicalEta(itEGamma->hwEta());
176 const double phi = getPhysicalPhi(itEGamma->hwPhi());
180 itEGamma->hwEta(), itEGamma->hwPhi(),
181 itEGamma->hwQual(), itEGamma->hwIso());
182 new_egammas->push_back(bx, *&eg);
188 itTau != old_rlxtaus->
end(bx); ++itTau)
198 const double et = jetScale->
et( itTau->hwPt() );
200 const double eta = getPhysicalEta(itTau->hwEta());
201 const double phi = getPhysicalPhi(itTau->hwPhi());
205 itTau->hwEta(), itTau->hwPhi(),
206 itTau->hwQual(), itTau->hwIso());
207 new_rlxtaus->push_back(bx, *&
tau);
212 itTau != old_isotaus->
end(bx); ++itTau)
222 const double et = jetScale->
et( itTau->hwPt() );
224 const double eta = getPhysicalEta(itTau->hwEta());
225 const double phi = getPhysicalPhi(itTau->hwPhi());
229 itTau->hwEta(), itTau->hwPhi(),
230 itTau->hwQual(), itTau->hwIso());
231 new_isotaus->push_back(bx, *&
tau);
236 itJet != old_jets->
end(bx); ++itJet)
246 const double et = jetScale->
et( itJet->hwPt() );
248 const bool forward = ((itJet->hwQual() & 0x2) != 0);
249 const double eta = getPhysicalEta(itJet->hwEta(), forward);
250 const double phi = getPhysicalPhi(itJet->hwPhi());
254 itJet->hwEta(), itJet->hwPhi(),
256 new_jets->push_back(bx, *&
jet);
261 itJet != old_preGtJets->
end(bx); ++itJet)
268 const double et = itJet->hwPt() * emScale->
linearLsb();
273 const bool forward = ((itJet->hwQual() & 0x2) != 0);
274 const double eta = getPhysicalEta(itJet->hwEta(), forward);
275 const double phi = getPhysicalPhi(itJet->hwPhi());
279 itJet->hwEta(), itJet->hwPhi(),
281 new_preGtJets->push_back(bx, *&
jet);
287 itEtSum != old_etsums->
end(bx); ++itEtSum)
289 double et = itEtSum->hwPt() * emScale->
linearLsb();
294 const double eta = getPhysicalEta(itEtSum->hwEta());
295 double phi = getPhysicalPhi(itEtSum->hwPhi());
297 et = htMissScale->
et( itEtSum->hwPt() );
299 phi=phi+(regionPhiWidth/2.);
305 EtSum eg(*&
p4, sumType, itEtSum->hwPt(),
306 itEtSum->hwEta(), itEtSum->hwPhi(),
308 new_etsums->push_back(bx, *&eg);
314 itCaloSpare != old_hfsums->
end(bx); ++itCaloSpare)
318 new_hfsums->push_back(bx, *itCaloSpare);
322 itCaloSpare != old_hfcounts->
end(bx); ++itCaloSpare)
326 new_hfcounts->push_back(bx, *itCaloSpare);
edm::EDGetToken JetToken_
const_iterator end(int bx) const
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
edm::EDGetToken HfSumsToken_
bool getByToken(EDGetToken token, Handle< PROD > &result) const
edm::EDGetToken HfCountsToken_
edm::EDGetToken EtSumToken_
double et(const uint16_t rank) const
convert from rank to physically meaningful quantity
double linearLsb() const
get LSB of linear input scale
PtEtaPhiMLorentzVectorD PtEtaPhiMLorentzVector
Lorentz vector with cartesian internal representation.
edm::EDGetToken EGammaToken_
edm::EDGetToken RlxTauToken_
edm::EDGetToken preGtJetToken_
et
define resolution functions of each parameter
edm::EDGetToken IsoTauToken_
static const unsigned N_PHI
const_iterator begin(int bx) const
std::vector< EGamma >::const_iterator const_iterator