@@ -2825,8 +2825,8 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
28252825 // polarization parameters
28262826 bool useHE = fgUsedVars[kCosThetaHE ] || fgUsedVars[kPhiHE ]; // helicity frame
28272827 bool useCS = fgUsedVars[kCosThetaCS ] || fgUsedVars[kPhiCS ]; // Collins-Soper frame
2828- bool usePP = fgUsedVars[kCosThetaPP ]; // production plane frame
2829- bool useRM = fgUsedVars[kCosThetaRM ]; // Random frame
2828+ bool usePP = fgUsedVars[kCosThetaPP ]; // production plane frame
2829+ bool useRM = fgUsedVars[kCosThetaRM ]; // Random frame
28302830
28312831 if (useHE || useCS || usePP || useRM) {
28322832 // TO DO: get the correct values from CCDB
@@ -2847,7 +2847,8 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
28472847 ROOT::Math::XYZVectorF zaxis_HE{(v12.Vect ()).Unit ()};
28482848 ROOT::Math::XYZVectorF yaxis_HE{(Beam1_CM.Cross (Beam2_CM)).Unit ()};
28492849 ROOT::Math::XYZVectorF xaxis_HE{(yaxis_HE.Cross (zaxis_HE)).Unit ()};
2850- if (fgUsedVars[kCosThetaHE ]) values[kCosThetaHE ] = zaxis_HE.Dot (v_CM);
2850+ if (fgUsedVars[kCosThetaHE ])
2851+ values[kCosThetaHE ] = zaxis_HE.Dot (v_CM);
28512852 if (fgUsedVars[kPhiHE ]) {
28522853 values[kPhiHE ] = TMath::ATan2 (yaxis_HE.Dot (v_CM), xaxis_HE.Dot (v_CM));
28532854 if (values[kPhiHE ] < 0 ) {
@@ -2877,7 +2878,8 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
28772878 ROOT::Math::XYZVectorF zaxis_CS{(Beam1_CM - Beam2_CM).Unit ()};
28782879 ROOT::Math::XYZVectorF yaxis_CS{(Beam1_CM.Cross (Beam2_CM)).Unit ()};
28792880 ROOT::Math::XYZVectorF xaxis_CS{(yaxis_CS.Cross (zaxis_CS)).Unit ()};
2880- if (fgUsedVars[kCosThetaCS ]) values[kCosThetaCS ] = zaxis_CS.Dot (v_CM);
2881+ if (fgUsedVars[kCosThetaCS ])
2882+ values[kCosThetaCS ] = zaxis_CS.Dot (v_CM);
28812883 if (fgUsedVars[kPhiCS ]) {
28822884 values[kPhiCS ] = TMath::ATan2 (yaxis_CS.Dot (v_CM), xaxis_CS.Dot (v_CM));
28832885 if (values[kPhiCS ] < 0 ) {
@@ -2939,7 +2941,8 @@ void VarManager::FillPair(T1 const& t1, T2 const& t2, float* values)
29392941 double randomCostheta = gRandom ->Uniform (-1 ., 1 .);
29402942 double randomPhi = gRandom ->Uniform (0 ., 2 . * TMath::Pi ());
29412943 ROOT::Math::XYZVectorF zaxis_RM (randomCostheta, std::sqrt (1 - randomCostheta * randomCostheta) * std::cos (randomPhi), std::sqrt (1 - randomCostheta * randomCostheta) * std::sin (randomPhi));
2942- if (fgUsedVars[kCosThetaRM ]) values[kCosThetaRM ] = zaxis_RM.Dot (v_CM);
2944+ if (fgUsedVars[kCosThetaRM ])
2945+ values[kCosThetaRM ] = zaxis_RM.Dot (v_CM);
29432946 }
29442947 }
29452948
@@ -3334,11 +3337,11 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
33343337 values[kMCPhi ] = v12.Phi ();
33353338 values[kMCY ] = -v12.Rapidity ();
33363339
3337- // polarization parameters
3340+ // polarization parameters
33383341 bool useHE = fgUsedVars[kMCCosThetaHE ] || fgUsedVars[kMCPhiHE ]; // helicity frame
33393342 bool useCS = fgUsedVars[kMCCosThetaCS ] || fgUsedVars[kMCPhiCS ]; // Collins-Soper frame
3340- bool usePP = fgUsedVars[kMCCosThetaPP ]; // production plane frame
3341- bool useRM = fgUsedVars[kMCCosThetaRM ]; // Random frame
3343+ bool usePP = fgUsedVars[kMCCosThetaPP ]; // production plane frame
3344+ bool useRM = fgUsedVars[kMCCosThetaRM ]; // Random frame
33423345
33433346 if (useHE || useCS || usePP || useRM) {
33443347 // TO DO: get the correct values from CCDB
@@ -3359,7 +3362,8 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
33593362 ROOT::Math::XYZVectorF zaxis_HE{(v12.Vect ()).Unit ()};
33603363 ROOT::Math::XYZVectorF yaxis_HE{(Beam1_CM.Cross (Beam2_CM)).Unit ()};
33613364 ROOT::Math::XYZVectorF xaxis_HE{(yaxis_HE.Cross (zaxis_HE)).Unit ()};
3362- if (fgUsedVars[kMCCosThetaHE ]) values[kMCCosThetaHE ] = zaxis_HE.Dot (v_CM);
3365+ if (fgUsedVars[kMCCosThetaHE ])
3366+ values[kMCCosThetaHE ] = zaxis_HE.Dot (v_CM);
33633367 if (fgUsedVars[kMCPhiHE ]) {
33643368 values[kMCPhiHE ] = TMath::ATan2 (yaxis_HE.Dot (v_CM), xaxis_HE.Dot (v_CM));
33653369 if (values[kMCPhiHE ] < 0 ) {
@@ -3389,7 +3393,8 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
33893393 ROOT::Math::XYZVectorF zaxis_CS{(Beam1_CM - Beam2_CM).Unit ()};
33903394 ROOT::Math::XYZVectorF yaxis_CS{(Beam1_CM.Cross (Beam2_CM)).Unit ()};
33913395 ROOT::Math::XYZVectorF xaxis_CS{(yaxis_CS.Cross (zaxis_CS)).Unit ()};
3392- if (fgUsedVars[kMCCosThetaCS ]) values[kMCCosThetaCS ] = zaxis_CS.Dot (v_CM);
3396+ if (fgUsedVars[kMCCosThetaCS ])
3397+ values[kMCCosThetaCS ] = zaxis_CS.Dot (v_CM);
33933398 if (fgUsedVars[kMCPhiCS ]) {
33943399 values[kMCPhiCS ] = TMath::ATan2 (yaxis_CS.Dot (v_CM), xaxis_CS.Dot (v_CM));
33953400 if (values[kMCPhiCS ] < 0 ) {
@@ -3451,7 +3456,8 @@ void VarManager::FillPairMC(T1 const& t1, T2 const& t2, float* values)
34513456 double randomCostheta = gRandom ->Uniform (-1 ., 1 .);
34523457 double randomPhi = gRandom ->Uniform (0 ., 2 . * TMath::Pi ());
34533458 ROOT::Math::XYZVectorF zaxis_RM (randomCostheta, std::sqrt (1 - randomCostheta * randomCostheta) * std::cos (randomPhi), std::sqrt (1 - randomCostheta * randomCostheta) * std::sin (randomPhi));
3454- if (fgUsedVars[kMCCosThetaRM ]) values[kMCCosThetaRM ] = zaxis_RM.Dot (v_CM);
3459+ if (fgUsedVars[kMCCosThetaRM ])
3460+ values[kMCCosThetaRM ] = zaxis_RM.Dot (v_CM);
34553461 }
34563462 }
34573463}
0 commit comments