@@ -695,6 +695,8 @@ class VarManager : public TObject
695695 kWV2SP ,
696696 kWV2EP ,
697697 kU2Q2 ,
698+ kU2Q2TPCPOS ,
699+ kU2Q2TPCNEG ,
698700 kU3Q3 ,
699701 kQ42XA ,
700702 kQ42YA ,
@@ -4888,6 +4890,8 @@ void VarManager::FillPairVn(T1 const& t1, T2 const& t2, float* values)
48884890 // Compute the scalar product UQ using Q-vector from A, for second and third harmonic
48894891 // Dilepton vn could be accessible after dividing this product with the R factor
48904892 values[kU2Q2 ] = values[kQ2X0A ] * TMath::Cos (2 * v12.Phi ()) + values[kQ2Y0A ] * TMath::Sin (2 * v12.Phi ());
4893+ values[kU2Q2TPCPOS ] = values[kQ2X0APOS ] * TMath::Cos (2 * v12.Phi ()) + values[kQ2Y0APOS ] * TMath::Sin (2 * v12.Phi ());
4894+ values[kU2Q2TPCNEG ] = values[kQ2X0ANEG ] * TMath::Cos (2 * v12.Phi ()) + values[kQ2Y0ANEG ] * TMath::Sin (2 * v12.Phi ());
48914895 values[kU3Q3 ] = values[kQ3X0A ] * TMath::Cos (3 * v12.Phi ()) + values[kQ3Y0A ] * TMath::Sin (3 * v12.Phi ());
48924896 values[kR2SP_AB ] = (values[kQ2X0A ] * values[kQ2X0B ] + values[kQ2Y0A ] * values[kQ2Y0B ]);
48934897 values[kR2SP_AC ] = (values[kQ2X0A ] * values[kQ2X0C ] + values[kQ2Y0A ] * values[kQ2Y0C ]);
@@ -4924,6 +4928,15 @@ void VarManager::FillPairVn(T1 const& t1, T2 const& t2, float* values)
49244928 values[kR2SP_AC ] = -999 .;
49254929 values[kR2SP_BC ] = -999 .;
49264930 }
4931+
4932+ if (std::isnan (VarManager::fgValues[VarManager::kU2Q2TPCPOS ]) == true ) {
4933+ values[kU2Q2TPCPOS ] = -999 .;
4934+ }
4935+
4936+ if (std::isnan (VarManager::fgValues[VarManager::kU2Q2TPCNEG ]) == true ) {
4937+ values[kU2Q2TPCNEG ] = -999 .;
4938+ }
4939+
49274940 if (std::isnan (VarManager::fgValues[VarManager::kU3Q3 ]) == true ) {
49284941 values[kU3Q3 ] = -999 .;
49294942 values[kR3SP ] = -999 .;
@@ -4934,6 +4947,15 @@ void VarManager::FillPairVn(T1 const& t1, T2 const& t2, float* values)
49344947 values[kR2EP_AC ] = -999 .;
49354948 values[kR2EP_BC ] = -999 .;
49364949 }
4950+
4951+ if (std::isnan (VarManager::fgValues[VarManager::kCos2DeltaPhiTPCPOS ]) == true ) {
4952+ values[kCos2DeltaPhiTPCPOS ] = -999 .;
4953+ }
4954+
4955+ if (std::isnan (VarManager::fgValues[VarManager::kCos2DeltaPhiTPCNEG ]) == true ) {
4956+ values[kCos2DeltaPhiTPCNEG ] = -999 .;
4957+ }
4958+
49374959 if (std::isnan (VarManager::fgValues[VarManager::kCos3DeltaPhi ]) == true ) {
49384960 values[kCos3DeltaPhi ] = -999 .;
49394961 values[kR3EP ] = -999 .;
0 commit comments