Skip to content

PWGHF: make uniform event selection #5261

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
fgrosa merged 20 commits into from
Mar 20, 2024
Merged

PWGHF: make uniform event selection #5261

Show file tree
Hide file tree
Changes from 3 commits
Commits
Show all changes
20 commits
Select commit Hold shift + click to select a range
ea65985
PWGHF: make uniform event selection
fgrosa Mar 20, 2024
044c2d7
Please consider the following formatting changes
alibuild Mar 20, 2024
af0d988
Merge pull request #69 from alibuild/alibot-cleanup-5261
fgrosa Mar 20, 2024
0ff35f5
Implement Vit's comments and restore possibility to select Run2 trigg…
fgrosa Mar 20, 2024
f0dede0
Please consider the following formatting changes
alibuild Mar 20, 2024
94e1d8d
Merge pull request #70 from alibuild/alibot-cleanup-5261
fgrosa Mar 20, 2024
fe5be41
Implement Vit's comments
fgrosa Mar 20, 2024
548ac82
Please consider the following formatting changes
alibuild Mar 20, 2024
2f6a327
Merge pull request #71 from alibuild/alibot-cleanup-5261
fgrosa Mar 20, 2024
799cf0c
Update PWGHF/Utils/utilsEvSelHf.h
fgrosa Mar 20, 2024
a8a39aa
Fix copy-paste mistake
fgrosa Mar 20, 2024
ae9d07f
Trivial commit
fgrosa Mar 20, 2024
b7f0602
Change type in data model
fgrosa Mar 20, 2024
b2a7c71
Improve histo title
fgrosa Mar 20, 2024
3f4d802
Document template parameters
fgrosa Mar 20, 2024
c7247cf
Update PWGHF/Utils/utilsEvSelHf.h
fgrosa Mar 20, 2024
8a23f02
Update PWGHF/Utils/utilsEvSelHf.h
fgrosa Mar 20, 2024
86976ca
Please consider the following formatting changes
alibuild Mar 20, 2024
64566f7
Merge pull request #72 from alibuild/alibot-cleanup-5261
fgrosa Mar 20, 2024
2ce2b5c
Update PWGHF/Utils/utilsEvSelHf.h
fgrosa Mar 20, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
130 changes: 72 additions & 58 deletions PWGHF/TableProducer/candidateCreator2Prong.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -31,6 +31,7 @@
#include "CommonConstants/PhysicsConstants.h"
#include "DCAFitter/DCAFitterN.h"
#include "Framework/AnalysisTask.h"
#include "Framework/HistogramRegistry.h"
#include "Framework/runDataProcessing.h"
#include "ReconstructionDataFormats/DCA.h"

Expand Down Expand Up @@ -90,21 +91,8 @@ struct HfCandidateCreator2Prong {
double massKPi{0.};
double bz{0.};

OutputObj<TH1F> hMass2{TH1F("hMass2", "2-prong candidates;inv. mass (#pi K) (GeV/#it{c}^{2});entries", 500, 0., 5.)};
OutputObj<TH1F> hCovPVXX{TH1F("hCovPVXX", "2-prong candidates;XX element of cov. matrix of prim. vtx. position (cm^{2});entries", 100, 0., 1.e-4)};
OutputObj<TH1F> hCovSVXX{TH1F("hCovSVXX", "2-prong candidates;XX element of cov. matrix of sec. vtx. position (cm^{2});entries", 100, 0., 0.2)};
OutputObj<TH1F> hCovPVYY{TH1F("hCovPVYY", "2-prong candidates;YY element of cov. matrix of prim. vtx. position (cm^{2});entries", 100, 0., 1.e-4)};
OutputObj<TH1F> hCovSVYY{TH1F("hCovSVYY", "2-prong candidates;YY element of cov. matrix of sec. vtx. position (cm^{2});entries", 100, 0., 0.2)};
OutputObj<TH1F> hCovPVXZ{TH1F("hCovPVXZ", "2-prong candidates;XZ element of cov. matrix of prim. vtx. position (cm^{2});entries", 100, -1.e-4, 1.e-4)};
OutputObj<TH1F> hCovSVXZ{TH1F("hCovSVXZ", "2-prong candidates;XZ element of cov. matrix of sec. vtx. position (cm^{2});entries", 100, -1.e-4, 0.2)};
OutputObj<TH1F> hCovPVZZ{TH1F("hCovPVZZ", "2-prong candidates;ZZ element of cov. matrix of prim. vtx. position (cm^{2});entries", 100, 0., 1.e-4)};
OutputObj<TH1F> hCovSVZZ{TH1F("hCovSVZZ", "2-prong candidates;ZZ element of cov. matrix of sec. vtx. position (cm^{2});entries", 100, 0., 0.2)};
OutputObj<TH2F> hDcaXYProngs{TH2F("hDcaXYProngs", "DCAxy of 2-prong candidates;#it{p}_{T} (GeV/#it{c};#it{d}_{xy}) (#mum);entries", 100, 0., 20., 200, -500., 500.)};
OutputObj<TH2F> hDcaZProngs{TH2F("hDcaZProngs", "DCAz of 2-prong candidates;#it{p}_{T} (GeV/#it{c};#it{d}_{z}) (#mum);entries", 100, 0., 20., 200, -500., 500.)};
OutputObj<TH1F> hVertexerType{TH1F("hVertexerType", "Use KF or DCAFitterN;Vertexer type;entries", 2, -0.5, 1.5)}; // See o2::aod::hf_cand::VertexerType
OutputObj<TH1D> hCollisions{TH1D("hCollisions", "HF event counter", ValuesEvSel::NEvSel, -0.5f, static_cast<float>(ValuesEvSel::NEvSel) - 0.5f)};
OutputObj<TH1D> hPosZBeforeEvSel{TH1D("hPosZBeforeEvSel", "PV position Z before ev. selection;posZ (cm);entries", 400, -20, 20)};
OutputObj<TH1D> hPosZAfterEvSel{TH1D("hPosZAfterEvSel", "PV position Z after ev. selection;posZ (cm);entries", 400, -20, 20)};
std::shared_ptr<TH1> hCollisions, hPosZBeforeEvSel, hPosZAfterEvSel, hPosXAfterEvSel, hPosYAfterEvSel, hNumPvContributorsAfterSel;
HistogramRegistry registry{"registry"};

void init(InitContext const&)
{
Expand Down Expand Up @@ -135,11 +123,31 @@ struct HfCandidateCreator2Prong {
}
}

// histograms
registry.add("hMass2", "2-prong candidates;inv. mass (#pi K) (GeV/#it{c}^{2});entries", {HistType::kTH1F, {{500, 1.6, 2.1}}});
registry.add("hCovPVXX", "2-prong candidates;XX element of cov. matrix of prim. vtx. position (cm^{2});entries", {HistType::kTH1F, {{100, 0., 1.e-4}}});
registry.add("hCovSVXX", "2-prong candidates;XX element of cov. matrix of sec. vtx. position (cm^{2});entries", {HistType::kTH1F, {{100, 0., 0.2}}});
registry.add("hCovPVYY", "2-prong candidates;YY element of cov. matrix of prim. vtx. position (cm^{2});entries", {HistType::kTH1F, {{100, 0., 1.e-4}}});
registry.add("hCovSVYY", "2-prong candidates;YY element of cov. matrix of sec. vtx. position (cm^{2});entries", {HistType::kTH1F, {{100, 0., 0.2}}});
registry.add("hCovPVXZ", "2-prong candidates;XZ element of cov. matrix of prim. vtx. position (cm^{2});entries", {HistType::kTH1F, {{100, -1.e-4, 1.e-4}}});
registry.add("hCovSVXZ", "2-prong candidates;XZ element of cov. matrix of sec. vtx. position (cm^{2});entries", {HistType::kTH1F, {{100, -1.e-4, 0.2}}});
registry.add("hCovPVZZ", "2-prong candidates;ZZ element of cov. matrix of prim. vtx. position (cm^{2});entries", {HistType::kTH1F, {{100, 0., 1.e-4}}});
registry.add("hCovSVZZ", "2-prong candidates;ZZ element of cov. matrix of sec. vtx. position (cm^{2});entries", {HistType::kTH1F, {{100, 0., 0.2}}});
registry.add("hDcaXYProngs", "DCAxy of 2-prong candidate daughters;#it{p}_{T} (GeV/#it{c};#it{d}_{xy}) (#mum);entries", {HistType::kTH2F, {{100, 0., 20.}, {200, -500., 500.}}});
registry.add("hDcaZProngs", "DCAz of 2-prong candidate daughters;#it{p}_{T} (GeV/#it{c};#it{d}_{z}) (#mum);entries", {HistType::kTH2F, {{100, 0., 20.}, {200, -500., 500.}}});
registry.add("hVertexerType", "Use KF or DCAFitterN;Vertexer type;entries", {HistType::kTH1D, {{2, -0.5, 1.5}}}); // See o2::aod::hf_cand::VertexerType
hCollisions = registry.add<TH1>("hCollisions", "HF event counter;;entries", {HistType::kTH1D, {{ValuesEvSel::NEvSel, -0.5f, static_cast<float>(ValuesEvSel::NEvSel) - 0.5f}}});
hPosZBeforeEvSel = registry.add<TH1>("hPosZBeforeEvSel", "all events;#it{z}_{prim. vtx.} (cm);entries", {HistType::kTH1D, {{400, -20., 20.}}});
hPosZAfterEvSel = registry.add<TH1>("hPosZAfterEvSel", "selected events;#it{z}_{prim. vtx.} (cm);entries", {HistType::kTH1D, {{400, -20., 20.}}});
hPosXAfterEvSel = registry.add<TH1>("hPosXAfterEvSel", "selected events;#it{x}_{prim. vtx.} (cm);entries", {HistType::kTH1D, {{200, -0.5, 0.5}}});
hPosYAfterEvSel = registry.add<TH1>("hPosYAfterEvSel", "selected events;#it{y}_{prim. vtx.} (cm);entries", {HistType::kTH1D, {{200, -0.5, 0.5}}});
hNumPvContributorsAfterSel = registry.add<TH1>("hNumPvContributorsAfterSel", "selected events;#it{y}_{prim. vtx.} (cm);entries", {HistType::kTH1D, {{500, -0.5, 499.5}}});

massPi = MassPiPlus;
massK = MassKPlus;

if (std::accumulate(doprocessDF.begin(), doprocessDF.end(), 0) == 1) {
hVertexerType->Fill(aod::hf_cand::VertexerType::DCAFitter);
registry.fill(HIST("hVertexerType"), aod::hf_cand::VertexerType::DCAFitter);
// Configure DCAFitterN
// df.setBz(bz);
df.setPropagateToPCA(propagateToPCA);
Expand All @@ -151,7 +159,7 @@ struct HfCandidateCreator2Prong {
df.setWeightedFinalPCA(useWeightedFinalPCA);
}
if (std::accumulate(doprocessKF.begin(), doprocessKF.end(), 0) == 1) {
hVertexerType->Fill(aod::hf_cand::VertexerType::KfParticle);
registry.fill(HIST("hVertexerType"), aod::hf_cand::VertexerType::KfParticle);
}

ccdb->setURL(ccdbUrl);
Expand All @@ -161,7 +169,7 @@ struct HfCandidateCreator2Prong {
runNumber = 0;

/// collision monitoring
setLabelHistoEvSel(hCollisions.object);
setLabelHistoEvSel(hCollisions);
}

template <bool doPvRefit, o2::aod::hf_collision_centrality::CentralityEstimator centEstimator, typename Coll, typename CandType, typename TTracks>
Expand All @@ -175,7 +183,8 @@ struct HfCandidateCreator2Prong {

/// reject candidates not satisfying the event selections
auto collision = rowTrackIndexProng2.template collision_as<Coll>();
const auto rejectionMask = getHfCollisionRejectionMask<centEstimator>(collision, centralityMin, centralityMax, useSel8Trigger, maxPvPosZ, useTimeFrameBorderCut);
float centrality{-1.f};
const auto rejectionMask = getHfCollisionRejectionMask<true, centEstimator>(collision, centrality, centralityMin, centralityMax, useSel8Trigger, useTimeFrameBorderCut, -maxPvPosZ, maxPvPosZ, 0, -1.f);
if (rejectionMask != 0) {
/// at least one event selection not satisfied --> reject the candidate
continue;
Expand Down Expand Up @@ -207,10 +216,10 @@ struct HfCandidateCreator2Prong {
const auto& secondaryVertex = df.getPCACandidate();
auto chi2PCA = df.getChi2AtPCACandidate();
auto covMatrixPCA = df.calcPCACovMatrixFlat();
hCovSVXX->Fill(covMatrixPCA[0]); // FIXME: Calculation of errorDecayLength(XY) gives wrong values without this line.
hCovSVYY->Fill(covMatrixPCA[2]);
hCovSVXZ->Fill(covMatrixPCA[3]);
hCovSVZZ->Fill(covMatrixPCA[5]);
registry.fill(HIST("hCovSVXX"), covMatrixPCA[0]); // FIXME: Calculation of errorDecayLength(XY) gives wrong values without this line.
registry.fill(HIST("hCovSVYY"), covMatrixPCA[2]);
registry.fill(HIST("hCovSVXZ"), covMatrixPCA[3]);
registry.fill(HIST("hCovSVZZ"), covMatrixPCA[5]);
auto trackParVar0 = df.getTrack(0);
auto trackParVar1 = df.getTrack(1);

Expand Down Expand Up @@ -240,18 +249,18 @@ struct HfCandidateCreator2Prong {
primaryVertex.setSigmaZ2(rowTrackIndexProng2.pvRefitSigmaZ2());
covMatrixPV = primaryVertex.getCov();
}
hCovPVXX->Fill(covMatrixPV[0]);
hCovPVYY->Fill(covMatrixPV[2]);
hCovPVXZ->Fill(covMatrixPV[3]);
hCovPVZZ->Fill(covMatrixPV[5]);
registry.fill(HIST("hCovPVXX"), covMatrixPV[0]);
registry.fill(HIST("hCovPVYY"), covMatrixPV[2]);
registry.fill(HIST("hCovPVXZ"), covMatrixPV[3]);
registry.fill(HIST("hCovPVZZ"), covMatrixPV[5]);
o2::dataformats::DCA impactParameter0;
o2::dataformats::DCA impactParameter1;
trackParVar0.propagateToDCA(primaryVertex, bz, &impactParameter0);
trackParVar1.propagateToDCA(primaryVertex, bz, &impactParameter1);
hDcaXYProngs->Fill(track0.pt(), impactParameter0.getY() * toMicrometers);
hDcaXYProngs->Fill(track1.pt(), impactParameter1.getY() * toMicrometers);
hDcaZProngs->Fill(track0.pt(), impactParameter0.getZ() * toMicrometers);
hDcaZProngs->Fill(track1.pt(), impactParameter1.getZ() * toMicrometers);
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0.getY() * toMicrometers);
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1.getY() * toMicrometers);
registry.fill(HIST("hDcaZProngs"), track0.pt(), impactParameter0.getZ() * toMicrometers);
registry.fill(HIST("hDcaZProngs"), track1.pt(), impactParameter1.getZ() * toMicrometers);

// get uncertainty of the decay length
double phi, theta;
Expand Down Expand Up @@ -287,8 +296,8 @@ struct HfCandidateCreator2Prong {
auto arrayMomenta = std::array{pvec0, pvec1};
massPiK = RecoDecay::m(arrayMomenta, std::array{massPi, massK});
massKPi = RecoDecay::m(arrayMomenta, std::array{massK, massPi});
hMass2->Fill(massPiK);
hMass2->Fill(massKPi);
registry.fill(HIST("hMass2"), massPiK);
registry.fill(HIST("hMass2"), massKPi);
}
}
}
Expand All @@ -304,7 +313,8 @@ struct HfCandidateCreator2Prong {

/// reject candidates in collisions not satisfying the event selections
auto collision = rowTrackIndexProng2.template collision_as<Coll>();
const auto rejectionMask = getHfCollisionRejectionMask<centEstimator>(collision, centralityMin, centralityMax, useSel8Trigger, maxPvPosZ, useTimeFrameBorderCut);
float centrality{-1.f};
const auto rejectionMask = getHfCollisionRejectionMask<true, centEstimator>(collision, centrality, centralityMin, centralityMax, useSel8Trigger, useTimeFrameBorderCut, -maxPvPosZ, maxPvPosZ, 0, -1.f);
if (rejectionMask != 0) {
/// at least one event selection not satisfied --> reject the candidate
continue;
Expand Down Expand Up @@ -340,10 +350,10 @@ struct HfCandidateCreator2Prong {
}
kfpVertex.GetCovarianceMatrix(covMatrixPV);
KFParticle KFPV(kfpVertex);
hCovPVXX->Fill(covMatrixPV[0]);
hCovPVYY->Fill(covMatrixPV[2]);
hCovPVXZ->Fill(covMatrixPV[3]);
hCovPVZZ->Fill(covMatrixPV[5]);
registry.fill(HIST("hCovPVXX"), covMatrixPV[0]);
registry.fill(HIST("hCovPVYY"), covMatrixPV[2]);
registry.fill(HIST("hCovPVXZ"), covMatrixPV[3]);
registry.fill(HIST("hCovPVZZ"), covMatrixPV[5]);

KFPTrack kfpTrack0 = createKFPTrackFromTrack(track0);
KFPTrack kfpTrack1 = createKFPTrackFromTrack(track1);
Expand All @@ -355,18 +365,18 @@ struct HfCandidateCreator2Prong {

float impactParameter0XY = 0., errImpactParameter0XY = 0., impactParameter1XY = 0., errImpactParameter1XY = 0.;
if (!kfPosPion.GetDistanceFromVertexXY(KFPV, impactParameter0XY, errImpactParameter0XY)) {
hDcaXYProngs->Fill(track0.pt(), impactParameter0XY * toMicrometers);
hDcaZProngs->Fill(track0.pt(), std::sqrt(kfPosPion.GetDistanceFromVertex(KFPV) * kfPosPion.GetDistanceFromVertex(KFPV) - impactParameter0XY * impactParameter0XY) * toMicrometers);
registry.fill(HIST("hDcaXYProngs"), track0.pt(), impactParameter0XY * toMicrometers);
registry.fill(HIST("hDcaZProngs"), track0.pt(), std::sqrt(kfPosPion.GetDistanceFromVertex(KFPV) * kfPosPion.GetDistanceFromVertex(KFPV) - impactParameter0XY * impactParameter0XY) * toMicrometers);
} else {
hDcaXYProngs->Fill(track0.pt(), -999.);
hDcaZProngs->Fill(track0.pt(), -999.);
registry.fill(HIST("hDcaXYProngs"), track0.pt(), -999.f);
registry.fill(HIST("hDcaZProngs"), track0.pt(), -999.f);
}
if (!kfNegPion.GetDistanceFromVertexXY(KFPV, impactParameter1XY, errImpactParameter1XY)) {
hDcaXYProngs->Fill(track1.pt(), impactParameter1XY * toMicrometers);
hDcaZProngs->Fill(track1.pt(), std::sqrt(kfNegPion.GetDistanceFromVertex(KFPV) * kfNegPion.GetDistanceFromVertex(KFPV) - impactParameter1XY * impactParameter1XY) * toMicrometers);
registry.fill(HIST("hDcaXYProngs"), track1.pt(), impactParameter1XY * toMicrometers);
registry.fill(HIST("hDcaZProngs"), track1.pt(), std::sqrt(kfPosPion.GetDistanceFromVertex(KFPV) * kfPosPion.GetDistanceFromVertex(KFPV) - impactParameter1XY * impactParameter1XY) * toMicrometers);
} else {
hDcaXYProngs->Fill(track1.pt(), -999.);
hDcaZProngs->Fill(track1.pt(), -999.);
registry.fill(HIST("hDcaXYProngs"), track1.pt(), -999.f);
registry.fill(HIST("hDcaZProngs"), track1.pt(), -999.f);
}

KFParticle kfCandD0;
Expand All @@ -381,10 +391,11 @@ struct HfCandidateCreator2Prong {
auto massD0 = kfCandD0.GetMass();
auto massD0bar = kfCandD0bar.GetMass();

hCovSVXX->Fill(kfCandD0.Covariance(0, 0));
hCovSVYY->Fill(kfCandD0.Covariance(1, 1));
hCovSVXZ->Fill(kfCandD0.Covariance(2, 0));
hCovSVZZ->Fill(kfCandD0.Covariance(2, 2));
registry.fill(HIST("hCovSVXX"), kfCandD0.Covariance(0, 0));
registry.fill(HIST("hCovSVYY"), kfCandD0.Covariance(1, 1));
registry.fill(HIST("hCovSVXZ"), kfCandD0.Covariance(2, 0));
registry.fill(HIST("hCovSVZZ"), kfCandD0.Covariance(2, 2));

auto covMatrixSV = kfCandD0.CovarianceMatrix();

double phi, theta;
Expand Down Expand Up @@ -426,8 +437,8 @@ struct HfCandidateCreator2Prong {

// fill histograms
if (fillHistograms) {
hMass2->Fill(massD0);
hMass2->Fill(massD0bar);
registry.fill(HIST("hMass2"), massD0);
registry.fill(HIST("hMass2"), massD0bar);
}
}
}
Expand Down Expand Up @@ -583,10 +594,11 @@ struct HfCandidateCreator2Prong {
for (const auto& collision : collisions) {

/// bitmask with event. selection info
const auto rejectionMask = getHfCollisionRejectionMask<CentralityEstimator::None>(collision, centralityMin, centralityMax, useSel8Trigger, maxPvPosZ, useTimeFrameBorderCut);
float centrality{-1.f};
const auto rejectionMask = getHfCollisionRejectionMask<true, CentralityEstimator::None>(collision, centrality, centralityMin, centralityMax, useSel8Trigger, useTimeFrameBorderCut, -maxPvPosZ, maxPvPosZ, 0, -1.f);

/// monitor the satisfied event selections
monitorCollision(collision, rejectionMask, hCollisions.object, hPosZBeforeEvSel.object, hPosZAfterEvSel.object);
monitorCollision(collision, rejectionMask, hCollisions, hPosZBeforeEvSel, hPosZAfterEvSel, hPosXAfterEvSel, hPosYAfterEvSel, hNumPvContributorsAfterSel);

} /// end loop over collisions
}
Expand All @@ -599,10 +611,11 @@ struct HfCandidateCreator2Prong {
for (const auto& collision : collisions) {

/// bitmask with event. selection info
const auto rejectionMask = getHfCollisionRejectionMask<CentralityEstimator::FT0C>(collision, centralityMin, centralityMax, useSel8Trigger, maxPvPosZ, useTimeFrameBorderCut);
float centrality{-1.f};
const auto rejectionMask = getHfCollisionRejectionMask<true, CentralityEstimator::FT0C>(collision, centrality, centralityMin, centralityMax, useSel8Trigger, useTimeFrameBorderCut, -maxPvPosZ, maxPvPosZ, 0, -1.f);

/// monitor the satisfied event selections
monitorCollision(collision, rejectionMask, hCollisions.object, hPosZBeforeEvSel.object, hPosZAfterEvSel.object);
monitorCollision(collision, rejectionMask, hCollisions, hPosZBeforeEvSel, hPosZAfterEvSel, hPosXAfterEvSel, hPosYAfterEvSel, hNumPvContributorsAfterSel);

} /// end loop over collisions
}
Expand All @@ -615,10 +628,11 @@ struct HfCandidateCreator2Prong {
for (const auto& collision : collisions) {

/// bitmask with event. selection info
const auto rejectionMask = getHfCollisionRejectionMask<CentralityEstimator::FT0M>(collision, centralityMin, centralityMax, useSel8Trigger, maxPvPosZ, useTimeFrameBorderCut);
float centrality{-1.f};
const auto rejectionMask = getHfCollisionRejectionMask<true, CentralityEstimator::FT0M>(collision, centrality, centralityMin, centralityMax, useSel8Trigger, useTimeFrameBorderCut, -maxPvPosZ, maxPvPosZ, 0, -1.f);

/// monitor the satisfied event selections
monitorCollision(collision, rejectionMask, hCollisions.object, hPosZBeforeEvSel.object, hPosZAfterEvSel.object);
monitorCollision(collision, rejectionMask, hCollisions, hPosZBeforeEvSel, hPosZAfterEvSel, hPosXAfterEvSel, hPosYAfterEvSel, hNumPvContributorsAfterSel);

} /// end loop over collisions
}
Expand Down
Loading