diff --git a/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.cc b/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.cc index a4b0cb9a75fb574ec281637d2b4a6a295ae900e9..7513237ccbca68256a62feadc531671d2e7ad471 100644 --- a/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.cc +++ b/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.cc @@ -333,9 +333,6 @@ complex<double> AbsPsi2STo2K2PiGamLh::chiToK0K2Amp(Psi2STo2K2PiGamData::Psi2STo2 - // result+=theMag*expiphi* - // (ampK0ToKpPi0*ampK2ToKmPi1+ampK0ToKpPi1*ampK2ToKmPi0+ampK0ToKmPi0+ampK2ToKpPi1+ampK0ToKmPi1*ampK2ToKpPi0); - result+=theMag*expiphi* (ampK0ToKpPi0*ampK2ToKmPi1+ampK0ToKpPi1*ampK2ToKmPi0+ampK0ToKmPi0*ampK2ToKpPi1+ampK0ToKmPi1*ampK2ToKpPi0); @@ -345,6 +342,49 @@ complex<double> AbsPsi2STo2K2PiGamLh::chiToK0K2Amp(Psi2STo2K2PiGamData::Psi2STo2 } +complex<double> AbsPsi2STo2K2PiGamLh::chiToK0K1Amp(Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& ChiToK0K1, double K0Mass, double K0Width, double K1Mass, double K1Width){ + complex<double> result(0.,0.); + + Vector4<double> KpPi0(theData->KpPi0_HeliChic0_4V.E(), theData->KpPi0_HeliChic0_4V.Px(), + theData->KpPi0_HeliChic0_4V.Py(), theData->KpPi0_HeliChic0_4V.Pz()); + + Vector4<double> KpPi1(theData->KpPi1_HeliChic0_4V.E(), theData->KpPi1_HeliChic0_4V.Px(), + theData->KpPi1_HeliChic0_4V.Py(), theData->KpPi1_HeliChic0_4V.Pz()); + + Vector4<double> KmPi0(theData->KmPi0_HeliChic0_4V.E(), theData->KmPi0_HeliChic0_4V.Px(), + theData->KmPi0_HeliChic0_4V.Py(), theData->KmPi0_HeliChic0_4V.Pz()); + + Vector4<double> KmPi1(theData->KmPi1_HeliChic0_4V.E(), theData->KmPi1_HeliChic0_4V.Px(), + theData->KmPi1_HeliChic0_4V.Py(), theData->KmPi1_HeliChic0_4V.Pz()); + + complex<double> ampK0ToKpPi0= BreitWigner(KpPi0, K0Mass, K0Width); + complex<double> ampK0ToKpPi1= BreitWigner(KpPi1, K0Mass, K0Width); + complex<double> ampK0ToKmPi0= BreitWigner(KmPi0, K0Mass, K0Width); + complex<double> ampK0ToKmPi1= BreitWigner(KmPi1, K0Mass, K0Width); + complex<double> ampK1ToKpPi0= sqrt(3.)*(conj(theData->DfKst1pToKpPi0[1][0][0])*BreitWignerBlattW(KpPi0, 0.493677, 0.1349766, K1Mass, K1Width, 1)); + complex<double> ampK1ToKpPi1= sqrt(3.)*(conj(theData->DfKst1pToKpPi1[1][0][0])*BreitWignerBlattW(KpPi1, 0.493677, 0.1349766, K1Mass, K1Width, 1)); + complex<double> ampK1ToKmPi0= sqrt(3.)*(conj(theData->DfKst1mToKmPi0[1][0][0])*BreitWignerBlattW(KmPi0, 0.493677, 0.1349766, K1Mass, K1Width, 1)); + complex<double> ampK1ToKmPi1= sqrt(3.)*(conj(theData->DfKst1mToKmPi1[1][0][0])*BreitWignerBlattW(KmPi1, 0.493677, 0.1349766, K1Mass, K1Width, 1)); + + std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >::iterator it; + for ( it=ChiToK0K1.begin(); it!=ChiToK0K1.end(); ++it){ + + boost::shared_ptr<const JPCLS> ChiToK0K1_State=it->first; + double theMag=it->second.first; + double thePhi=it->second.second; + complex<double> expiphi(cos(thePhi), sin(thePhi)); + + + result+=theMag*expiphi*sqrt(2.*ChiToK0K1_State->L+1.)* + (ampK0ToKpPi0*ampK1ToKmPi1+ampK0ToKpPi1*ampK1ToKmPi0+ampK0ToKmPi0*ampK1ToKpPi1+ampK0ToKmPi1*ampK1ToKpPi0); + + } + + return result; + +} + + complex<double> AbsPsi2STo2K2PiGamLh::chiTo2K_0_Amp(Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& ChiTo2K_0, double K_0_Mass0, double K_0_Width0, double K_0_Mass1, double K_0_Width1){ complex<double> result(0.,0.); diff --git a/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.hh b/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.hh index 7f8335790e29a49d04363be566d6ac9f2131b341..f8616ca177ec396a567852ac3aca0cc4573c0a23 100644 --- a/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.hh +++ b/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.hh @@ -77,7 +77,7 @@ protected: virtual complex<double> chiToK0K2Amp(Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& ChiToK0K2, double K0Mass, double K0Width, double K2Mass, double K2Width); -// virtual complex<double> chiToK1ToK1piAmp(Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& K1ToK1dPi, double K1Mass, double K1Width, double K1dMass, double K1dWidth); + virtual complex<double> chiToK0K1Amp(Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& ChiToK0K1, double K0Mass, double K0Width, double K1Mass, double K1Width); virtual complex<double> chiToK1ToK1piAmp(Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& K1Prod, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& K1Dec, double K1Mass, double K1Width, double K1dMass, double K1dWidth); diff --git a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc index c8d11733339ccdffd33d4cb8fcd04ff8fcbd1971..f2d236ba0114df1c617a1b2b669c4da16688a525 100644 --- a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc +++ b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc @@ -62,6 +62,7 @@ FitParams2K2PiGam::FitParams2K2PiGam() _jpclsMap[paramEnum2K2PiGam::K_1_2400Tof980K]=theStates.K1pTof0KStates(); _jpclsMap[paramEnum2K2PiGam::K892K_1_1680]=theStates.ChiTo2K892States(); _jpclsMap[paramEnum2K2PiGam::ChiToK1680K1680]=theStates.ChiTo2K892States(); + _jpclsMap[paramEnum2K2PiGam::ChiToK1680K_0_1430]=theStates.ChiToKst0Kst1States(); _jpclsMap[paramEnum2K2PiGam::K892K_1_2300]=theStates.ChiTo2K892States(); _jpclsMap[paramEnum2K2PiGam::sigmaf980]=theStates.ChiTof0f0States(); _jpclsMap[paramEnum2K2PiGam::ChiToSigmaf1370]=theStates.ChiTof0f0States(); @@ -161,6 +162,7 @@ std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collecti else if (index==paramEnum2K2PiGam::K_1_2400Tof980K) return params.K_1_2400Tof980K; else if (index==paramEnum2K2PiGam::K892K_1_1680) return params.ChiToK892K1680; else if (index==paramEnum2K2PiGam::ChiToK1680K1680) return params.ChiToK1680K1680; + else if (index==paramEnum2K2PiGam::ChiToK1680K_0_1430) return params.ChiToK1680K_0_1430; else if (index==paramEnum2K2PiGam::K892K_1_2300) return params.ChiToK892K2300; else if (index==paramEnum2K2PiGam::sigmaf980) return params.ChiToSigmaf980; else if (index==paramEnum2K2PiGam::ChiToSigmaf1370) return params.ChiToSigmaf1370; diff --git a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh index dca54f65954f60fa72c3e4cb06e308a842c60278..4562c121b1d419343dca6dab168a906bcb1a200b 100644 --- a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh +++ b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh @@ -76,6 +76,7 @@ struct param2K2PiGam std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToKappaK_0_1950; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK892K1680; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK1680K1680; + std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK1680K_0_1430; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK892K2300; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > K_0_1460ToK892Pi; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > K_0_1460ToK_0_1430Pi; @@ -150,7 +151,7 @@ struct paramEnum2K2PiGam{ ChiToK_1_2400K, K_1_2400Tof980K, f980_pif1710_k, f980_kf1710_pi, 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, K892K_1_2300, + K_0_2400KToKf980, ChiToK_0_2400ToKf_0_1710, K892K_1_1680, ChiToK1680K1680, ChiToK1680K_0_1430, K892K_1_2300, sigmaf980, ChiToSigmaf1370, sigmaf1710, sigmaf2200, 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, 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, @@ -177,7 +178,7 @@ struct paramEnum2K2PiGam{ "ChiToK_1_2400K", "K_1_2400Tof980K", "f980_pif1710_k", "f980_kf1710_pi", "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", "K892K_1_2300", + "K_0_2400KToKf980", "ChiToK_0_2400ToKf_0_1710", "K892K_1_1680", "ChiToK1680K1680", "ChiToK1680K_0_1430", "K892K_1_2300", "sigmaf980", "ChiToSigmaf1370", "sigmaf1710", "sigmaf2200", "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", "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/Hyp7Lh.cc b/Examples/Psi2STo2K2PiGam/Hyp7Lh.cc index bed233a934289332e9c0c78e833d361f61d61110..b8f993dfc4b986a645071598c5926745108c7881 100644 --- a/Examples/Psi2STo2K2PiGam/Hyp7Lh.cc +++ b/Examples/Psi2STo2K2PiGam/Hyp7Lh.cc @@ -12,6 +12,7 @@ Hyp7Lh::Hyp7Lh(boost::shared_ptr<const Psi2STo2K2PiGamEvtList> theEvtList, const ,_KappaHyp(true) ,_K1_1680Hyp(true) ,_K1_1680K1_1680Hyp7(true) + ,_K1_1680K0_1430Hyp7(true) ,_K1_2300Hyp(true) ,_nFitParams(0) { @@ -23,6 +24,7 @@ Hyp7Lh::Hyp7Lh( boost::shared_ptr<AbsPsi2STo2K2PiGamLh> theLhPtr, const std::map ,_KappaHyp(true) ,_K1_1680Hyp(true) ,_K1_1680K1_1680Hyp7(true) + ,_K1_1680K0_1430Hyp7(true) ,_K1_2300Hyp(true) ,_nFitParams(0) { @@ -49,7 +51,7 @@ complex<double> Hyp7Lh::chi0DecAmps(const param2K2PiGam& theParamVal, Psi2STo2K2 result+=chiTo2K_0_Amp(theData, ChiTo2Kappa, KappaMass, KappaWidth, KappaMass, KappaWidth); } - if(_K1_1680Hyp || _K1_1680K1_1680Hyp7){ + if(_K1_1680Hyp || _K1_1680K1_1680Hyp7 || _K1_1680K0_1430Hyp7){ double K_1_1680Mass=theParamVal.BwK_1_1680.first; double K_1_1680Width=theParamVal.BwK_1_1680.second; @@ -64,6 +66,14 @@ complex<double> Hyp7Lh::chi0DecAmps(const param2K2PiGam& theParamVal, Psi2STo2K2 result+=chiToK1K1Amp(theData, ChiToK1680K1680, K_1_1680Mass, K_1_1680Width, K_1_1680Mass, K_1_1680Width); } + if(_K1_1680K0_1430Hyp7){ + + double K_0_1430Mass=theParamVal.BwK_0_1430.first; + double K_0_1430Width=theParamVal.BwK_0_1430.second; + std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK1680K_0_1430=theParamVal.ChiToK1680K_0_1430; + result+=chiToK0K1Amp(theData, ChiToK1680K_0_1430, K_0_1430Mass, K_0_1430Width, K_1_1680Mass, K_1_1680Width); + } + } if(_K1_2300Hyp){ @@ -220,6 +230,16 @@ void Hyp7Lh::setUp(const std::map<const std::string, bool>& hypMap){ } else Alert << "hypothesis K1_1680K1_1680Hyp7 not set!!!" <<endmsg; + + iter= hypMap.find("K1_1680K0_1430Hyp7"); + if (iter !=hypMap.end()){ + _K1_1680K0_1430Hyp7= iter->second; + Info<< "hypothesis " << iter->first << "\t" << _K1_1680K0_1430Hyp7 <<endmsg; + _hypMap[iter->first]= iter->second; + } + else Alert << "hypothesis K1_1680K0_1430Hyp7 not set!!!" <<endmsg; + + iter= hypMap.find("K1_2300Hyp7"); if (iter !=hypMap.end()){ @@ -236,7 +256,7 @@ void Hyp7Lh::setUp(const std::map<const std::string, bool>& hypMap){ if(!_KappaK_0_1950Hyp6) _massVec.push_back(paramEnum2K2PiGam::Kappa); } - if(_K1_1680Hyp || _K1_1680K1_1680Hyp7){ + if(_K1_1680Hyp || _K1_1680K1_1680Hyp7 || _K1_1680K0_1430Hyp7){ _massVec.push_back(paramEnum2K2PiGam::K_1_1680); if(_K1_1680Hyp){ @@ -245,6 +265,11 @@ void Hyp7Lh::setUp(const std::map<const std::string, bool>& hypMap){ if(_K1_1680K1_1680Hyp7){ _ampVec.push_back(paramEnum2K2PiGam::ChiToK1680K1680); } + + if(_K1_1680K0_1430Hyp7){ + _ampVec.push_back(paramEnum2K2PiGam::ChiToK1680K_0_1430); + } + } if(_K1_2300Hyp){ diff --git a/Examples/Psi2STo2K2PiGam/Hyp7Lh.hh b/Examples/Psi2STo2K2PiGam/Hyp7Lh.hh index 6204a1401704e02587da83b53101281538628b2f..27ea4ebc2b26836c671f5f87d263a2b1ded11f69 100644 --- a/Examples/Psi2STo2K2PiGam/Hyp7Lh.hh +++ b/Examples/Psi2STo2K2PiGam/Hyp7Lh.hh @@ -53,6 +53,7 @@ protected: bool _KappaHyp; bool _K1_1680Hyp; bool _K1_1680K1_1680Hyp7; + bool _K1_1680K0_1430Hyp7; bool _K1_2300Hyp; virtual complex<double> chi0DecAmps(const param2K2PiGam& theParamVal, Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData); diff --git a/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc b/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc index 1b865065d9e84e3df88ffd803cdfdbeba5f07667..fe21b404dcc5b2a568088c53eb5c4e79317a2c6c 100644 --- a/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc +++ b/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc @@ -174,6 +174,7 @@ int main(int __argc,char *__argv[]){ hypMap["KappaK_0_1950Hyp6"]=true; hypMap["K1_1680Hyp7"]=true; hypMap["K1_1680K1_1680Hyp7"]=true; + hypMap["K1_1680K0_1430Hyp7"]=true; hypMap["K1_2300Hyp7"]=true; hypMap["KappaHyp7"]=true; hypMap["K_0_1460ToKstPiHyp8"]=true; diff --git a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.cc b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.cc index 6f90d3383d7d8129d2178ec19b6168368a5cb331..4cbc7561590cd2fd18019372ba66bb9d7fdb7d82 100644 --- a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.cc +++ b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.cc @@ -38,6 +38,9 @@ Psi2STo2K2PiGamStates::Psi2STo2K2PiGamStates() : //LS combinations for the Chi_c0 decay to K*0(1430) K*2(1430) fillJPCLS(_chic0JPC, _Kst0JPC, _Kst2JPC, _ChiToK0K2_JPCLS); + //LS combinations for the Chi_c0 decay to K*0 K*1 + fillJPCLS(_chic0JPC, _Kst0JPC, _Kst1JPC, _ChiToKst0Kst1JPCLS); + //LS combinations for the Chi_c0 decay to f0 f0 fillJPCLS(_chic0JPC, _f0JPC, _f0JPC,_ChiTof0f0JPCLS); @@ -168,6 +171,12 @@ void Psi2STo2K2PiGamStates::print(std::ostream& os) const os << "\n" << std::endl; } + os << "*** Chi_c0: LS combinations for the decay to K*0 K*1 *** "<< std::endl; + for ( itJPCLS=_ChiToKst0Kst1JPCLS.begin(); itJPCLS!=_ChiToKst0Kst1JPCLS.end(); ++itJPCLS){ + (*itJPCLS)->print(os); + os << "\n" << std::endl; + } + os << "*** Chi_c0: LS combinations for the decay to f0 f0 *** "<< std::endl; for ( itJPCLS=_ChiTof0f0JPCLS.begin(); itJPCLS!=_ChiTof0f0JPCLS.end(); ++itJPCLS){ (*itJPCLS)->print(os); diff --git a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.hh b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.hh index b85172d2ae17b33f4d83c0012b66c871ae0a309c..748e7e8b42c76b9135ce6738e9f590a489b5cd6b 100644 --- a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.hh +++ b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.hh @@ -35,6 +35,7 @@ public: std::vector< boost::shared_ptr<const JPCLS> > ChiTof0f2States() const {return _ChiTof0f2JPCLS;} std::vector< boost::shared_ptr<const JPCLS> > K1400ToKst1PiStates() const {return _K1400ToKst1PiJPCLS;} std::vector< boost::shared_ptr<const JPCLS> > ChiToK0K0States() const {return _ChiToK0K0JPCLS;} + std::vector< boost::shared_ptr<const JPCLS> > ChiToKst0Kst1States() const {return _ChiToKst0Kst1JPCLS;} std::vector< boost::shared_ptr<const JPCLS> > K1ToK0PiStates() const {return _K1ToK0PiJPCLS;} std::vector< boost::shared_ptr<const JPCLS> > ChiToPi_2PiStates() const {return _ChiToPi_2PiJPCLS;} std::vector< boost::shared_ptr<const JPCLS> > Pi_2Tof_2PiStates() const {return _Pi_2Tof_2PiJPCLS;} @@ -77,6 +78,7 @@ private: std::vector< boost::shared_ptr<const JPCLS> > _ChiTof0f0JPCLS; std::vector< boost::shared_ptr<const JPCLS> > _ChiTof0f2JPCLS; std::vector< boost::shared_ptr<const JPCLS> > _ChiToK0K0JPCLS; + std::vector< boost::shared_ptr<const JPCLS> > _ChiToKst0Kst1JPCLS; std::vector< boost::shared_ptr<const JPCLS> > _K1ToK0PiJPCLS; std::vector< boost::shared_ptr<const JPCLS> > _K1400ToKst1PiJPCLS; std::vector< boost::shared_ptr<const JPCLS> > _ChiToPi_2PiJPCLS;