Page MenuHomeHEPForge
Files F19189532

D91.1759065340.diff
No OneTemporary
Actions

Size
3 KB
Referenced Files
None
Subscribers
None

D91.1759065340.diff
View Options

diff --git a/History.md b/History.md
--- a/History.md
+++ b/History.md
@@ -11,6 +11,10 @@
===
## R02-02-0X
+3 Feb 2023 John Back
+* D91: Check for non-zero momentum for EvtSLBaryonAmp parent spin density matrix.
+ Print out integrals of JSON test histograms.
+
16 Dec 2022 John Back
* D89: Added probabilities for B_c -> V pi+ and V pi+ pi0 EvtBcVHad modes.
diff --git a/src/EvtGenModels/EvtSLBaryonAmp.cpp b/src/EvtGenModels/EvtSLBaryonAmp.cpp
--- a/src/EvtGenModels/EvtSLBaryonAmp.cpp
+++ b/src/EvtGenModels/EvtSLBaryonAmp.cpp
@@ -427,7 +427,7 @@
EvtVector4R vector4P = parent->getP4Lab();
double pmag = vector4P.d3mag();
- double cosTheta = vector4P.get( 3 ) / pmag;
+ double cosTheta = pmag > 0.0 ? vector4P.get( 3 ) / pmag : 1.0;
double theta = acos( cosTheta );
double phi = atan2( vector4P.get( 2 ), vector4P.get( 1 ) );
@@ -461,7 +461,7 @@
FormFactors->getdiracff( par_num, bar_num, q2, baryonmass, &f1, &f2,
&f3, &g1, &g2, &g3 );
- const double form_fact[6] = {f1, f2, f3, g1, g2, g3};
+ const double form_fact[6] = { f1, f2, f3, g1, g2, g3 };
EvtVector4C b11, b12, b21, b22, l1, l2;
@@ -584,7 +584,7 @@
FormFactors->getraritaff( par_num, bar_num, q2, baryonmass, &f1, &f2,
&f3, &f4, &g1, &g2, &g3, &g4 );
- const double form_fact[8] = {f1, f2, f3, f4, g1, g2, g3, g4};
+ const double form_fact[8] = { f1, f2, f3, f4, g1, g2, g3, g4 };
EvtId l_num = parent->getDaug( 1 )->getId();
diff --git a/test/testDecayModel.cc b/test/testDecayModel.cc
--- a/test/testDecayModel.cc
+++ b/test/testDecayModel.cc
@@ -1064,9 +1064,12 @@
double chiSq{ 0.0 };
int nDof{ 0 };
int iGood{ 0 };
- double pValue{ refHist->Chi2TestX( hist, chiSq, nDof, iGood, "WW" ) };
+ const double pValue{
+ refHist->Chi2TestX( hist, chiSq, nDof, iGood, "WW" ) };
+ const double integral{ refHist->Integral() };
std::cout << "Histogram " << histName << " chiSq/nDof = " << chiSq
- << "/" << nDof << ", pValue = " << pValue << std::endl;
+ << "/" << nDof << ", pValue = " << pValue
+ << ", integral = " << integral << std::endl;
} else {
std::cerr << "Could not find reference histogram " << histName
@@ -1084,9 +1087,12 @@
double chiSq{ 0.0 };
int nDof{ 0 };
int iGood{ 0 };
- double pValue{ refHist->Chi2TestX( hist, chiSq, nDof, iGood, "WW" ) };
+ const double pValue{
+ refHist->Chi2TestX( hist, chiSq, nDof, iGood, "WW" ) };
+ const double integral{ refHist->Integral() };
std::cout << "Histogram " << histName << " chiSq/nDof = " << chiSq
- << "/" << nDof << ", pValue = " << pValue << std::endl;
+ << "/" << nDof << ", pValue = " << pValue
+ << ", integral = " << integral << std::endl;
} else {
std::cerr << "Could not find reference histogram " << histName

File Metadata

Mime Type
text/plain
Expires
Sun, Sep 28, 2:15 PM (2 s ago)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
6531084
Default Alt Text
D91.1759065340.diff (3 KB)

Event Timeline

HEPForge.org · .