Skip to content

Commit 6f4b2e5

Browse files
committed
Fill vertexing chi2PCA for triplets
1 parent 4c43f6f commit 6f4b2e5

File tree

1 file changed

+5
-0
lines changed

1 file changed

+5
-0
lines changed

PWGDQ/Core/VarManager.h

Lines changed: 5 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -3803,6 +3803,7 @@ void VarManager::FillTripletVertexing(C const& collision, T const& t1, T const&
38033803
values[VarManager::kVertexingProcCode] = procCode;
38043804
if (procCode == 0) {
38053805
// TODO: set the other variables to appropriate values and return
3806+
values[kVertexingChi2PCA] = -999.;
38063807
values[kVertexingLxy] = -999.;
38073808
values[kVertexingLxyz] = -999.;
38083809
values[kVertexingLz] = -999.;
@@ -3837,6 +3838,10 @@ void VarManager::FillTripletVertexing(C const& collision, T const& t1, T const&
38373838
o2::dataformats::VertexBase primaryVertex = {std::move(vtxXYZ), std::move(vtxCov)};
38383839
auto covMatrixPV = primaryVertex.getCov();
38393840

3841+
auto chi2PCA = fgFitterThreeProngBarrel.getChi2AtPCACandidate();
3842+
if (fgUsedVars[kVertexingChi2PCA])
3843+
values[VarManager::kVertexingChi2PCA] = chi2PCA;
3844+
38403845
double phi = std::atan2(secondaryVertex[1] - collision.posY(), secondaryVertex[0] - collision.posX());
38413846
double theta = std::atan2(secondaryVertex[2] - collision.posZ(),
38423847
std::sqrt((secondaryVertex[0] - collision.posX()) * (secondaryVertex[0] - collision.posX()) +

0 commit comments

Comments
 (0)