Build a collection of piZeros from objects in the input jet.
225 edm::LogPrint(
"RecoTauPiZeroStripPlugin3") <<
"<RecoTauPiZeroStripPlugin3::operator()>:";
241 for (CandPtrs::iterator
cand = candsVector.begin();
cand != candsVector.end(); ++
cand) {
244 <<
"PFGamma #" <<
idx <<
" (" <<
cand->id() <<
":" <<
cand->key() <<
"): Et = " << (*cand)->et()
245 <<
", eta = " << (*cand)->eta() <<
", phi = " << (*cand)->phi();
249 edm::LogPrint(
"RecoTauPiZeroStripPlugin3") <<
"--> assigning seedCandId = " << seedCands.size();
253 <<
"track: Pt = " << candTrack->
pt() <<
" eta = " << candTrack->
eta()
254 <<
", phi = " << candTrack->
phi() <<
", charge = " << candTrack->
charge();
261 <<
", dPt/Pt = " << (candTrack->
ptError() / candTrack->
pt()) <<
")";
264 seedCands.push_back(*
cand);
267 edm::LogPrint(
"RecoTauPiZeroStripPlugin3") <<
"--> assigning addCandId = " << addCands.size();
269 addCands.push_back(*
cand);
274 std::vector<bool> seedCandFlags(seedCands.size());
275 std::vector<bool> addCandFlags(addCands.size());
277 std::set<size_t> seedCandIdsCurrentStrip;
278 std::set<size_t> addCandIdsCurrentStrip;
281 while (idxSeed < seedCands.size()) {
283 edm::LogPrint(
"RecoTauPiZeroStripPlugin3") <<
"processing seed #" << idxSeed;
285 seedCandIdsCurrentStrip.clear();
286 addCandIdsCurrentStrip.clear();
289 strip->addDaughter(seedCands[idxSeed]);
290 seedCandIdsCurrentStrip.insert(idxSeed);
293 int stripBuildIteration = 0;
305 ++stripBuildIteration;
311 <<
"Building strip: Et = " <<
strip->et() <<
", eta = " <<
strip->eta() <<
", phi = " <<
strip->phi();
314 if (
strip->daughterPtr(0).isNonnull())
315 strip->setVertex(
strip->daughterPtr(0)->vertex());
319 markCandsInStrip(seedCandFlags, seedCandIdsCurrentStrip);
320 markCandsInStrip(addCandFlags, addCandIdsCurrentStrip);
324 <<
"Discarding strip: Et = " <<
strip->et() <<
", eta = " <<
strip->eta() <<
", phi = " <<
strip->phi();
328 while (idxSeed < seedCands.size() && seedCandFlags[idxSeed]) {
337 std::sort(
output.begin(),
output.end(), [&](
auto& arg1,
auto& arg2) {
return arg1->pt() > arg2->pt(); });
342 for (PiZeroVector::const_iterator firstIter =
output.begin(); firstIter != end_iter - 1; ++firstIter) {
343 for (PiZeroVector::const_iterator secondIter = firstIter + 1; secondIter != end_iter; ++secondIter) {
344 auto const&
first = *firstIter;
345 auto const&
second = *secondIter;
353 auto combinedStrips =
354 std::make_unique<RecoTauPiZero>(0,
364 for (
auto const&
gamma :
first->daughterPtrVector()) {
365 combinedStrips->addDaughter(
gamma);
367 for (
auto const&
gamma :
second->daughterPtrVector()) {
368 combinedStrips->addDaughter(
gamma);
371 if (combinedStrips->daughterPtr(0).isNonnull()) {
372 combinedStrips->setVertex(combinedStrips->daughterPtr(0)->vertex());
376 stripCombinations.push_back(
std::move(combinedStrips));
381 std::move(stripCombinations.begin(), stripCombinations.end(), std::back_inserter(
output));
387 double bendCorrEta = 0.;
388 double bendCorrPhi = 0.;
390 for (
auto const&
gamma :
strip->daughterPtrVector()) {
400 strip->setBendCorrEta(bendCorrEta);
401 strip->setBendCorrPhi(bendCorrPhi);
int numberOfValidPixelHits() const
std::vector< CandidatePtr > pfCandidates(const Jet &jet, int particleId, bool sort=true)
InputIterator takeNElements(const InputIterator &begin, const InputIterator &end, size_t N)
void set(reco::Candidate &c) const
set up a candidate
double ptError() const
error on Pt (set to 1000 TeV if charge==0 for safety)
bool isNonnull() const
Checks for non-null.
void addCandsToStrip(RecoTauPiZero &, CandPtrs &, const std::vector< bool > &, std::set< size_t > &, bool &) const
bool updateStripAfterEachDaughter_
reco::VertexRef associatedVertex(const Jet &jet) const
U second(std::pair< T, U > const &p)
double pt() const
track transverse momentum
int charge() const
track electric charge
double dz() const
dz parameter (= dsz/cos(lambda)). This is the track z0 w.r.t (0,0,0) only if the refPoint is close to...
std::vector< reco::CandidatePtr > CandPtrs
double phi() const
azimuthal angle of momentum vector
std::unique_ptr< const TFormula > etaAssociationDistance_
Log< level::Warning, true > LogPrint
std::unique_ptr< RecoTauQualityCuts > qcuts_
static const TrackGhostTrackState * getTrack(const BasicGhostTrackState *basic)
double eta() const
pseudorapidity of momentum vector
double minGammaEtStripSeed_
const HitPattern & hitPattern() const
Access the hit pattern, indicating in which Tracker layers the track has hits.
int maxStripBuildIterations_
double normalizedChi2() const
chi-squared divided by n.d.o.f. (or chi-squared * 1e6 if n.d.o.f. is zero)
std::vector< std::unique_ptr< RecoTauPiZero > > PiZeroVector
double combinatoricStripMassHypo_
math::XYZTLorentzVector LorentzVector
Lorentz vector.
std::unique_ptr< const TFormula > phiAssociationDistance_
RecoTauVertexAssociator vertexAssociator_
AddFourMomenta p4Builder_
std::vector< int > inputParticleIds_
double minGammaEtStripAdd_
math::XYZPoint Point
point in the space
double energySum(const DataFrame &df, int fs, int ls)
int numberOfValidTrackerHits() const
double dxy() const
dxy parameter. (This is the transverse impact parameter w.r.t. to (0,0,0) ONLY if refPoint is close t...