108 auto pOut = std::make_unique<reco::JPTJetCollection>();
109 auto pOut1 = std::make_unique<reco::CaloJetCollection>();
111 double scaleJPT = 1.;
115 for (
auto const& oldjet : addjets_h) {
124 auto const& mycalo = addjets_h[icalo];
125 std::vector<edm::Ptr<reco::Track> > tracksinjet =
jet.tracks();
128 for (
auto const& itrack : tracksinjet) {
129 for (
auto const& ixtrp : iExtrapolations) {
130 if (ixtrp.positions().empty())
133 double mydphi =
deltaPhi(ixtrp.track()->phi(), itrack->phi());
134 if (fabs(ixtrp.track()->pt() - itrack->pt()) > 0.001 || fabs(ixtrp.track()->eta() - itrack->eta()) > 0.001 ||
138 if (itrack.id() != ixtrp.track().id() || itrack.key() != ixtrp.track().key())
157 mJPTalgo->
correction(corrected, mycalo, iEvent, iSetup, tracksinvert, tracksincalo, p4, pions, muons, elecs);
172 reco::JPTJet fJet(p4,
jet.primaryVertex()->position(), jptspe, mycalo.getJetConstituents());
173 pOut->push_back(fJet);
174 pOut1->push_back(mycalo);
179 for (
auto const& oldjet : jets_h) {
183 double factorZSP = 1.;
199 scaleJPT =
mJPTalgo->
correction(corrected, oldjet, iEvent, iSetup, pions, muons, elecs, ok);
201 corrected.
py() * scaleJPT,
202 corrected.
pz() * scaleJPT,
203 corrected.
energy() * scaleJPT);
205 scaleJPT =
mJPTalgo->
correction(corrected, oldjet, iEvent, iSetup, p4, pions, muons, elecs, ok);
235 double deEta2Tr = 0.;
236 double dePhi2Tr = 0.;
240 double denominator_tracks = 0.;
245 double deR =
deltaR((*it)->eta(), (*it)->phi(), p4.eta(), p4.phi());
246 double deEta = (*it)->eta() - p4.eta();
247 double dePhi =
deltaPhi((*it)->phi(), p4.phi());
248 if ((**it).ptError() / (**it).pt() < 0.1) {
249 deR2Tr = deR2Tr + deR * deR * (*it)->pt();
250 deEta2Tr = deEta2Tr + deEta * deEta * (*it)->pt();
251 dePhi2Tr = dePhi2Tr + dePhi * dePhi * (*it)->pt();
252 denominator_tracks = denominator_tracks + (*it)->pt();
253 Zch = Zch + (*it)->pt();
255 Pout2 = Pout2 + (**it).p() * (**it).p() - (Zch * p4.P()) * (Zch * p4.P());
262 double deR =
deltaR((*it)->eta(), (*it)->phi(), p4.eta(), p4.phi());
263 double deEta = (*it)->eta() - p4.eta();
264 double dePhi =
deltaPhi((*it)->phi(), p4.phi());
265 if ((**it).ptError() / (**it).pt() < 0.1) {
266 deR2Tr = deR2Tr + deR * deR * (*it)->pt();
267 deEta2Tr = deEta2Tr + deEta * deEta * (*it)->pt();
268 dePhi2Tr = dePhi2Tr + dePhi * dePhi * (*it)->pt();
269 denominator_tracks = denominator_tracks + (*it)->pt();
270 Zch = Zch + (*it)->pt();
272 Pout2 = Pout2 + (**it).p() * (**it).p() - (Zch * p4.P()) * (Zch * p4.P());
278 double deR =
deltaR((*it)->eta(), (*it)->phi(), p4.eta(), p4.phi());
279 double deEta = (*it)->eta() - p4.eta();
280 double dePhi =
deltaPhi((*it)->phi(), p4.phi());
281 if ((**it).ptError() / (**it).pt() < 0.1) {
282 deR2Tr = deR2Tr + deR * deR * (*it)->pt();
283 deEta2Tr = deEta2Tr + deEta * deEta * (*it)->pt();
284 dePhi2Tr = dePhi2Tr + dePhi * dePhi * (*it)->pt();
285 denominator_tracks = denominator_tracks + (*it)->pt();
286 Zch = Zch + (*it)->pt();
288 Pout2 = Pout2 + (**it).p() * (**it).p() - (Zch * p4.P()) * (Zch * p4.P());
295 Zch = Zch + (*it)->pt();
300 Zch = Zch + (*it)->pt();
305 Zch = Zch + (*it)->pt();
312 Pout =
sqrt(fabs(Pout2)) / ntracks;
314 if (denominator_tracks != 0) {
315 deR2Tr = deR2Tr / denominator_tracks;
316 deEta2Tr = deEta2Tr / denominator_tracks;
317 dePhi2Tr = dePhi2Tr / denominator_tracks;
323 specific.
Pout = Pout;
333 if (pvCollection.
isValid() && !pvCollection->empty())
334 vertex_ = pvCollection->begin()->position();
341 pOut->push_back(fJet);
343 std::sort(pOut->begin(), pOut->end(),
sort_by_pt);
reco::TrackRefVector muonsInVertexOutCalo
float mSumEnergyOfChargedWithoutEff
std::remove_cv< typename std::remove_reference< argument_type >::type >::type key_type
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
double pz() const final
z coordinate of momentum vector
Jets made from CaloTowers.
reco::TrackRefVector muonsInVertexInCalo
virtual void scaleEnergy(double fScale)
scale energy of the jet
bool getByToken(EDGetToken token, Handle< PROD > &result) const
float mChargedHadronEnergy
float mSumPtOfChargedWithEff
reco::TrackRefVector inVertexInCalo_
reco::TrackRefVector muonsOutVertexInCalo
reco::TrackRefVector inVertexOutOfCalo_
virtual double correction(const reco::Jet &, const edm::Event &, const edm::EventSetup &) const
apply correction using Event information
double getResponseOfChargedWithoutEff()
float mResponseOfChargedWithEff
edm::EDGetTokenT< edm::View< reco::CaloJet > > input_jets_token_
reco::TrackRefVector elecsOutVertexInCalo
virtual Constituents getJetConstituents() const
list of constituents
reco::TrackRefVector pionsInVertexOutCalo
const_iterator end() const
Termination of iteration.
const_iterator begin() const
Initialize an iterator over the RefVector.
XYZTLorentzVectorD XYZTLorentzVector
Lorentz vector with cylindrical internal representation using pseudorapidity.
double px() const final
x coordinate of momentum vector
ZSPJPTJetCorrector * mZSPalgo
edm::EDGetTokenT< edm::View< reco::CaloJet > > input_addjets_token_
Jets made from CaloJets corrected for ZSP and tracks.
double correction(const reco::Jet &, const reco::Jet &, const edm::Event &, const edm::EventSetup &, P4 &, jpt::MatchedTracks &pions, jpt::MatchedTracks &muons, jpt::MatchedTracks &elecs, bool &validMatches)
Vectorial correction method (corrected 4-momentum passed by reference)
double getResponseOfChargedWithEff()
Get responses/sumPT/SumEnergy with and without Efficiency correction.
reco::TrackRefVector pionsOutVertexInCalo
float mSumEnergyOfChargedWithEff
math::XYZPoint Point
point in the space
bool sort_by_pt(const reco::JPTJet &a, const reco::JPTJet &b)
reco::TrackRefVector outOfVertexInCalo_
bool get(ProductID const &oid, Handle< PROD > &result) const
reco::TrackRefVector elecsInVertexInCalo
math::XYZPoint Point
point in the space
edm::EDGetTokenT< std::vector< reco::TrackExtrapolation > > mExtrapolations_
double py() const final
y coordinate of momentum vector
double getSumPtWithoutEff()
RefProd< PROD > getRefBeforePut()
edm::RefToBase< reco::Jet > theCaloJetRef
float mResponseOfChargedWithoutEff
reco::TrackRefVector pionsInVertexInCalo
constexpr auto deltaR2(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
float mSumPtOfChargedWithoutEff
reco::TrackRefVector elecsInVertexOutCalo
Particles matched to tracks that are in/in, in/out, out/in at Vertex and CaloFace.
JetPlusTrackCorrector * mJPTalgo
void push_back(value_type const &ref)
Add a Ref<C, T> to the RefVector.
double getSumEnergyWithoutEff()
edm::EDGetTokenT< reco::VertexCollection > input_vertex_token_
double getSumEnergyWithEff()
*vegas h *****************************************************used in the default bin number in original ***version of VEGAS is ***a higher bin number might help to derive a more precise ***grade subtle point
edm::EDGetTokenT< edm::View< reco::TrackJet > > input_trackjets_token_
std::vector< CaloJet > CaloJetCollection
collection of CaloJet objects
math::XYZPoint Point
point in the space
double energy() const final
energy