From a5c8d4265d530f2d8e2fad947eff6caff61520ed Mon Sep 17 00:00:00 2001
From: Bertram Kopf <bertram@pc14.(none)>
Date: Fri, 2 Mar 2012 00:29:38 +0100
Subject: [PATCH] Psi2STo2K2PiGam: added sigma f2_1270
---
Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc | 6 ++--
Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh | 5 ++--
Examples/Psi2STo2K2PiGam/Hyp4Lh.cc | 30 +++++++++++++++----
Examples/Psi2STo2K2PiGam/Hyp4Lh.hh | 1 +
.../Mpsi2STo2K2PiGamTestApp.cc | 1 +
5 files changed, 34 insertions(+), 9 deletions(-)
diff --git a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc
index bd68d73e..f06117ce 100644
--- a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc
+++ b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc
@@ -91,8 +91,9 @@ void FitParams2K2PiGam::filljpclsMap(){
_jpclsMap[paramEnum2K2PiGam::ChiToSigmaf1370]=_statesPtr->ChiTof0f0States();
_jpclsMap[paramEnum2K2PiGam::sigmaf1710]=_statesPtr->ChiTof0f0States();
_jpclsMap[paramEnum2K2PiGam::sigmaf2200]=_statesPtr->ChiTof0f0States();
- _jpclsMap[paramEnum2K2PiGam::ChiToSigmaf_2_1525]=_statesPtr->ChiTof0f0States();
- _jpclsMap[paramEnum2K2PiGam::ChiToSigmaf_2_1950]=_statesPtr->ChiTof0f0States();
+ _jpclsMap[paramEnum2K2PiGam::ChiToSigmaf_2_1270]=_statesPtr->ChiTof0f2States();
+ _jpclsMap[paramEnum2K2PiGam::ChiToSigmaf_2_1525]=_statesPtr->ChiTof0f2States();
+ _jpclsMap[paramEnum2K2PiGam::ChiToSigmaf_2_1950]=_statesPtr->ChiTof0f2States();
_jpclsMap[paramEnum2K2PiGam::K_0_1460ToK892Pi]=_statesPtr->ChiToK0K0States();
_jpclsMap[paramEnum2K2PiGam::K_0_1460ToK_0_1430Pi]=_statesPtr->ChiToK0K0States();
_jpclsMap[paramEnum2K2PiGam::K_0_1830ToK892Pi]=_statesPtr->ChiToK0K0States();
@@ -203,6 +204,7 @@ std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collecti
else if (index==paramEnum2K2PiGam::ChiToSigmaf1370) return params.ChiToSigmaf1370;
else if (index==paramEnum2K2PiGam::sigmaf1710) return params.ChiToSigmaf1710;
else if (index==paramEnum2K2PiGam::sigmaf2200) return params.ChiToSigmaf2200;
+ else if (index==paramEnum2K2PiGam::ChiToSigmaf_2_1270) return params.ChiToSigmaf_2_1270;
else if (index==paramEnum2K2PiGam::ChiToSigmaf_2_1525) return params.ChiToSigmaf_2_1525;
else if (index==paramEnum2K2PiGam::ChiToSigmaf_2_1950) return params.ChiToSigmaf_2_1950;
else if (index==paramEnum2K2PiGam::K_0_1460ToK892Pi) return params.K_0_1460ToK892Pi;
diff --git a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh
index 5c003104..88763e69 100644
--- a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh
+++ b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh
@@ -59,6 +59,7 @@ struct param2K2PiGam
std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToSigmaf1370;
std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToSigmaf1710;
std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToSigmaf2200;
+ std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToSigmaf_2_1270;
std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToSigmaf_2_1525;
std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToSigmaf_2_1950;
std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > f980_pif1370_k;
@@ -173,7 +174,7 @@ struct paramEnum2K2PiGam{
f980_pif1710_k, f980_kf1710_pi, ChiTof1710f1710, f980f980, f980f2200, ChiTof980f_2_2300, ChiTof_2_2300sigma,
f980_pif1370_k, f980_kf1370_pi, f980_pif1500_k, f980_kf1500_pi, f1710_pif1370_k, f1710_kf1370_pi,
K_0_2400KToKf980, ChiToK_0_2400ToKf_0_1710, K892K_1_1680, ChiToK1680K1680, ChiToK1680K_0_1430, K892K_1_2300,
- sigmaf980, ChiToSigmaf1370, sigmaf1710, sigmaf2200, ChiToSigmaf_2_1525, ChiToSigmaf_2_1950,
+ sigmaf980, ChiToSigmaf1370, sigmaf1710, sigmaf2200, ChiToSigmaf_2_1270, ChiToSigmaf_2_1525, ChiToSigmaf_2_1950,
K_0_1460ToK892Pi, K_0_1460ToK_0_1430Pi, K_0_1830ToK892Pi, K_0_1830ToK_0_1430Pi,
f980_pif_2_1270_k, f980_kf_2_1270_pi, f980_pif_2_1430_k, f980_kf_2_1430_pi, f980_pif_2_1525_k, f980_kf_2_1525_pi, f980_pif_2_1950_k, f980_kf_2_1950_pi, f980_pif_2_2010_k, f980_kf_2_2010_pi,
f1500_pif_2_1525_k, f1500_kf_2_1525_pi, f1710_pif_2_1430_k, f1710_kf_2_1430_pi, f1710_pif_2_1950_k, f1710_kf_2_1950_pi,
@@ -203,7 +204,7 @@ struct paramEnum2K2PiGam{
"f980_pif1710_k", "f980_kf1710_pi", "ChiTof1710f1710", "f980f980", "f980f2200", "ChiTof980f_2_2300","ChiTof_2_2300sigma",
"f980_pif1370_k", "f980_kf1370_pi", "f980_pif1500_k", "f980_kf1500_pi", "f1710_pif1370_k", "f1710_kf1370_pi",
"K_0_2400KToKf980", "ChiToK_0_2400ToKf_0_1710", "K892K_1_1680", "ChiToK1680K1680", "ChiToK1680K_0_1430", "K892K_1_2300",
- "sigmaf980", "ChiToSigmaf1370", "sigmaf1710", "sigmaf2200", "ChiToSigmaf_2_1525", "ChiToSigmaf_2_1950",
+ "sigmaf980", "ChiToSigmaf1370", "sigmaf1710", "sigmaf2200", "ChiToSigmaf_2_1270", "ChiToSigmaf_2_1525", "ChiToSigmaf_2_1950",
"K_0_1460ToK892Pi", "K_0_1460ToK_0_1430Pi","K_0_1830ToK892Pi", "K_0_1830ToK_0_1430Pi",
"f980_pif_2_1270_k", "f980_kf_2_1270_pi","f980_pif_2_1430_k", "f980_kf_2_1430_pi","f980_pif_2_1525_k", "f980_kf_2_1525_pi", "f980_pif_2_1950_k", "f980_kf_2_1950_pi", "f980_pif_2_2010_k", "f980_kf_2_2010_pi",
"f1500_pif_2_1525_k", "f1500_kf_2_1525_pi", "f1710_pif_2_1430_k", "f1710_kf_2_1430_pi", "f1710_pif_2_1950_k", "f1710_kf_2_1950_pi",
diff --git a/Examples/Psi2STo2K2PiGam/Hyp4Lh.cc b/Examples/Psi2STo2K2PiGam/Hyp4Lh.cc
index 9dfd205e..dbe466cb 100644
--- a/Examples/Psi2STo2K2PiGam/Hyp4Lh.cc
+++ b/Examples/Psi2STo2K2PiGam/Hyp4Lh.cc
@@ -20,6 +20,7 @@ Hyp4Lh::Hyp4Lh(boost::shared_ptr<const Psi2STo2K2PiGamEvtList> theEvtList, const
,_f1500f_2_1525Hyp4(true)
,_f1710f_2_1430Hyp4(true)
,_f1710f_2_1950Hyp4(true)
+ ,_sigmaf_2_1270Hyp4(true)
,_sigmaf_2_1525Hyp4(true)
,_sigmaf_2_1950Hyp4(true)
,_doHyp4(true)
@@ -44,6 +45,7 @@ Hyp4Lh::Hyp4Lh( boost::shared_ptr<AbsPsi2STo2K2PiGamLh> theLhPtr, const std::map
,_f1500f_2_1525Hyp4(true)
,_f1710f_2_1430Hyp4(true)
,_f1710f_2_1950Hyp4(true)
+ ,_sigmaf_2_1270Hyp4(true)
,_sigmaf_2_1525Hyp4(true)
,_sigmaf_2_1950Hyp4(true)
,_doHyp4(true)
@@ -224,6 +226,15 @@ complex<double> Hyp4Lh::chi0DecAmps(const param2K2PiGam& theParamVal, Psi2STo2K2
currentResult+=chiTof2_pif0_kAmp(theData, ChiTof1710_kf2_pi, f_2_1950Mass, f_2_1950Width, f1710Mass, f1710Width);
}
+ if(_sigmaf_2_1270Hyp4){
+ std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToSigmaf_2_1270=theParamVal.ChiToSigmaf_2_1270;
+ double f_2_1270Mass=theParamVal.Bwf_2_1270.first;
+ double f_2_1270Width=theParamVal.Bwf_2_1270.second;
+ double sigmaMass=theParamVal.BwSigma.first;
+ double sigmaWidth=theParamVal.BwSigma.second;
+
+ currentResult+=chiTof0_pif2_kAmp(theData, ChiToSigmaf_2_1270, sigmaMass, sigmaWidth, f_2_1270Mass, f_2_1270Width);
+ }
if (_sigmaf_2_1525Hyp4){
std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToSigmaf_2_1525=theParamVal.ChiToSigmaf_2_1525;
@@ -413,6 +424,9 @@ void Hyp4Lh::setUp(const std::map<const std::string, bool>& hypMap){
theKey="sigmaf1370Hyp4";
setHyps( hypMap, _sigmaf1370Hyp4, theKey);
+ theKey="sigmaf_2_1270Hyp4";
+ setHyps( hypMap, _sigmaf_2_1270Hyp4, theKey);
+
theKey="sigmaf_2_1525Hyp4";
setHyps( hypMap, _sigmaf_2_1525Hyp4, theKey);
@@ -420,7 +434,7 @@ void Hyp4Lh::setUp(const std::map<const std::string, bool>& hypMap){
setHyps( hypMap, _sigmaf_2_1950Hyp4, theKey);
- if(!_f980f1370Hyp4 && !_f1710f1370Hyp4 && !_f980f_2_1430Hyp4 && !_f1710f_2_1430Hyp4 && !_f980f_2_1950Hyp4 && !_f980f_2_2010Hyp4 && !_f1710f_2_1950Hyp4 && !_sigmaf1370Hyp4 && !_sigmaf_2_1525Hyp4 && !_sigmaf_2_1950Hyp4) _doHyp4=false;
+ if(!_f980f1370Hyp4 && !_f1710f1370Hyp4 && !_f980f_2_1430Hyp4 && !_f1710f_2_1430Hyp4 && !_f980f_2_1950Hyp4 && !_f980f_2_2010Hyp4 && !_f1710f_2_1950Hyp4 && !_sigmaf1370Hyp4 && !_sigmaf_2_1270Hyp4 && !_sigmaf_2_1525Hyp4 && !_sigmaf_2_1950Hyp4) _doHyp4=false;
if (!_doHyp4) return;
@@ -438,10 +452,16 @@ void Hyp4Lh::setUp(const std::map<const std::string, bool>& hypMap){
_ampVec.push_back(paramEnum2K2PiGam::f1710_pif1370_k);
_ampVec.push_back(paramEnum2K2PiGam::f1710_kf1370_pi);
}
- if(_f980f_2_1270Hyp4){
- _ampVec.push_back(paramEnum2K2PiGam::f980_pif_2_1270_k);
- _ampVec.push_back(paramEnum2K2PiGam::f980_kf_2_1270_pi);
+ if(_f980f_2_1270Hyp4 || _sigmaf_2_1270Hyp4){
_massVec.push_back(paramEnum2K2PiGam::f_2_1270);
+
+ if(_f980f_2_1270Hyp4){
+ _ampVec.push_back(paramEnum2K2PiGam::f980_pif_2_1270_k);
+ _ampVec.push_back(paramEnum2K2PiGam::f980_kf_2_1270_pi);
+ }
+ if(_sigmaf_2_1270Hyp4){
+ _ampVec.push_back(paramEnum2K2PiGam::ChiToSigmaf_2_1270);
+ }
}
if (_f980f_2_1430Hyp4){
_ampVec.push_back(paramEnum2K2PiGam::f980_pif_2_1430_k);
@@ -489,7 +509,7 @@ void Hyp4Lh::setUp(const std::map<const std::string, bool>& hypMap){
if (_f980f_2_1525Hyp4 || _f1500f_2_1525Hyp4 || _sigmaf_2_1525Hyp4) _massVec.push_back(paramEnum2K2PiGam::f_2_1525);
// fill all other resonances
- if (_sigmaf1370Hyp4 || _sigmaf_2_1525Hyp4 || _sigmaf_2_1950Hyp4){
+ if (_sigmaf1370Hyp4 || _sigmaf_2_1270Hyp4 || _sigmaf_2_1525Hyp4 || _sigmaf_2_1950Hyp4){
_massVecRemain.push_back(paramEnum2K2PiGam::sigma);
}
diff --git a/Examples/Psi2STo2K2PiGam/Hyp4Lh.hh b/Examples/Psi2STo2K2PiGam/Hyp4Lh.hh
index 4643ca2e..8836d02f 100644
--- a/Examples/Psi2STo2K2PiGam/Hyp4Lh.hh
+++ b/Examples/Psi2STo2K2PiGam/Hyp4Lh.hh
@@ -62,6 +62,7 @@ protected:
bool _f1500f_2_1525Hyp4;
bool _f1710f_2_1430Hyp4;
bool _f1710f_2_1950Hyp4;
+ bool _sigmaf_2_1270Hyp4;
bool _sigmaf_2_1525Hyp4;
bool _sigmaf_2_1950Hyp4;
bool _doHyp4;
diff --git a/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc b/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc
index b0035afa..9b53baf0 100644
--- a/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc
+++ b/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc
@@ -173,6 +173,7 @@ int main(int __argc,char *__argv[]){
hypMap["f1710f_2_1430Hyp4"]=true;
hypMap["f1710f_2_1950Hyp4"]=true;
hypMap["sigmaf1370Hyp4"]=true;
+ hypMap["sigmaf_2_1270Hyp4"]=true;
hypMap["sigmaf_2_1525Hyp4"]=true;
hypMap["sigmaf_2_1950Hyp4"]=true;
hypMap["doHyp5"]=true;
--
GitLab