diff --git a/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.cc b/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.cc index 7513237ccbca68256a62feadc531671d2e7ad471..2db61c3a665c40914c377401035a1aa2bee48fe7 100644 --- a/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.cc +++ b/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.cc @@ -385,6 +385,63 @@ complex<double> AbsPsi2STo2K2PiGamLh::chiToK0K1Amp(Psi2STo2K2PiGamData::Psi2STo2 } + +complex<double> AbsPsi2STo2K2PiGamLh::chiToK1K2Amp(Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& ChiToK1K2, double K1Mass, double K1Width, double K2Mass, double K2Width){ + 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()); + + + std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >::iterator it; + for ( it=ChiToK1K2.begin(); it!=ChiToK1K2.end(); ++it){ + + boost::shared_ptr<const JPCLS> ChiToK1K2_State=it->first; + double theMag=it->second.first; + double thePhi=it->second.second; + complex<double> expiphi(cos(thePhi), sin(thePhi)); + + complex<double> tmpResult=theMag*expiphi*sqrt(2.*ChiToK1K2_State->L+1.)*sqrt(3.)*sqrt(5.); + for (Spin lamK1=-1; lamK1<=1; ++lamK1){ + Spin lamK2=lamK1; + Spin lamKK=lamK2-lamK1; + + complex<double> ampK1ToKpPi0= conj(theData->DfKst1pToKpPi0[1][lamK1][0])*BreitWignerBlattW(KpPi0, 0.493677, 0.1349766, K1Mass, K1Width, 1); + complex<double> ampK1ToKpPi1= conj(theData->DfKst1pToKpPi1[1][lamK1][0])*BreitWignerBlattW(KpPi1, 0.493677, 0.1349766, K1Mass, K1Width, 1); + complex<double> ampK1ToKmPi0= conj(theData->DfKst1mToKmPi0[1][lamK1][0])*BreitWignerBlattW(KmPi0, 0.493677, 0.1349766, K1Mass, K1Width, 1); + complex<double> ampK1ToKmPi1= conj(theData->DfKst1mToKmPi1[1][lamK1][0])*BreitWignerBlattW(KmPi1, 0.493677, 0.1349766, K1Mass, K1Width, 1); + + complex<double> ampK2ToKpPi0= conj(theData->DfKst2pToKpPi0[2][lamK2][0])*BreitWignerBlattW(KpPi0, 0.493677, 0.1349766, K2Mass, K2Width, 2); + complex<double> ampK2ToKpPi1= conj(theData->DfKst2pToKpPi1[2][lamK2][0])*BreitWignerBlattW(KpPi1, 0.493677, 0.1349766, K2Mass, K2Width, 2); + complex<double> ampK2ToKmPi0= conj(theData->DfKst2mToKmPi0[2][lamK2][0])*BreitWignerBlattW(KmPi0, 0.493677, 0.1349766, K2Mass, K2Width, 2); + complex<double> ampK2ToKmPi1= conj(theData->DfKst2mToKmPi1[2][lamK2][0])*BreitWignerBlattW(KmPi1, 0.493677, 0.1349766, K2Mass, K2Width, 2); + + tmpResult*=Clebsch(ChiToK1K2_State->L, 0, ChiToK1K2_State->S, lamKK, ChiToK1K2_State->J, lamKK)*Clebsch(2,lamK2, 1, -lamK1, ChiToK1K2_State->S, lamKK)* + ( + ampK1ToKpPi0*ampK2ToKmPi1 + +ampK1ToKpPi1*ampK2ToKmPi0 + +ampK1ToKmPi0*ampK2ToKpPi1 + +ampK1ToKmPi1*ampK2ToKpPi0 ); + + } + result+=tmpResult; + } + + 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 f8616ca177ec396a567852ac3aca0cc4573c0a23..d16312d8ef1f3c20e14f3016bec12d22771a88ab 100644 --- a/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.hh +++ b/Examples/Psi2STo2K2PiGam/AbsPsi2STo2K2PiGamLh.hh @@ -79,6 +79,8 @@ protected: 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> chiToK1K2Amp(Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& ChiToK1K2, double K1Mass, double K1Width, double K2Mass, double K2Width); + 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); virtual complex<double> chiToK0KT0KpipiKAmp(Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData, std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess >& K0ToK0Pi, double K0pMass, double K0pWidth, double K0dMass, double K0dWidth); diff --git a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc index 8ded6267cc00380c8d1c6c838a6691000f08f472..802654e111a45cebd12c7281febb36903ff6fc4c 100644 --- a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc +++ b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.cc @@ -21,6 +21,7 @@ FitParams2K2PiGam::FitParams2K2PiGam() _jpclsMap[paramEnum2K2PiGam::ChiToK_0_1430_K892]=theStates.ChiToKst0Kst1States(); _jpclsMap[paramEnum2K2PiGam::K_1_1410K_1_1410]=theStates.ChiTo2K892States(); _jpclsMap[paramEnum2K2PiGam::ChiToK_1_1410_K892]=theStates.ChiTo2K892States(); + _jpclsMap[paramEnum2K2PiGam::ChiToK_2_1430_K892]=theStates.ChiToKst1Kst2States(); _jpclsMap[paramEnum2K2PiGam::ChiToK_1_1400K]=theStates.ChiToK1400ToK892piStates(); _jpclsMap[paramEnum2K2PiGam::K_1_1400ToK892Pi]=theStates.K1400ToKst1PiStates(); _jpclsMap[paramEnum2K2PiGam::ChiToK_1_1270_K]=theStates.ChiToK1400ToK892piStates(); @@ -124,6 +125,7 @@ std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collecti else if (index==paramEnum2K2PiGam::ChiToK_0_1430_K892) return params.ChiToK_0_1430_K892; else if (index==paramEnum2K2PiGam::K_1_1410K_1_1410) return params.ChiToK_1_1410_K_1_1410; else if (index==paramEnum2K2PiGam::ChiToK_1_1410_K892) return params.ChiToK_1_1410_K892; + else if (index==paramEnum2K2PiGam::ChiToK_2_1430_K892) return params.ChiToK_2_1430_K892; else if (index==paramEnum2K2PiGam::ChiToK_1_1400K) return params.ChiToK_1_1400K; else if (index==paramEnum2K2PiGam::K_1_1400ToK892Pi) return params.K1400ToK892Pi; else if (index==paramEnum2K2PiGam::ChiToK_1_1270_K) return params.ChiToK_1_1270_K; @@ -356,11 +358,11 @@ void FitParams2K2PiGam::setMnUsrParamsDec(MnUserParameters& upar, param2K2PiGam& double magErr=errPair.first; double phiErr=errPair.second; - double magMin=magVal-magErr; + double magMin=magVal-2.*magErr; if (magMin<0.) magMin=0.; - upar.Add(magStr, magVal, magErr, magMin, magVal+magErr); - upar.Add(phiStr, phiVal, phiErr, -3*M_PI, 3*M_PI); + upar.Add(magStr, magVal, magErr, magMin, magVal+2.*magErr); + upar.Add(phiStr, phiVal, phiErr, -3.*M_PI, 3.*M_PI); counter++; } diff --git a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh index a040a54172636d1851451a55fa091a51b402ebb4..52ebf66f5ff2647790182246bdb6a788b9589beb 100644 --- a/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh +++ b/Examples/Psi2STo2K2PiGam/FitParams2K2PiGam.hh @@ -34,6 +34,7 @@ struct param2K2PiGam std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK_0_1430_K892; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK_1_1410_K_1_1410; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK_1_1410_K892; + std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK_2_1430_K892; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK_1_1400K; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > K1400ToK892Pi; std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK_1_1270_K; @@ -149,7 +150,7 @@ struct param2K2PiGam struct paramEnum2K2PiGam{ enum {ChiGam=0, K892K892, KappaKappa, K_0_1430K_0_1430, ChiToK_0_1430_K892, - K_2_1430K_2_1430, K_0_1430K_2_1430, K_1_1410K_1_1410, ChiToK_1_1410_K892, + K_2_1430K_2_1430, K_0_1430K_2_1430, K_1_1410K_1_1410, ChiToK_1_1410_K892, ChiToK_2_1430_K892, ChiToK_1_1400K, K_1_1400ToK892Pi, ChiToK_1_1270_K, K_1_1270ToK892Pi, K_1_1270ToK_0_1430Pi, ChiToK_1_1650K, K_1_1650ToK892Pi, K_1_1650ToK_0_1430Pi, K_0_1430K_0_1950, KappaK_0_1950, ChiToKappaK_0_1430, ChiToK_1_2400K, K_1_2400Tof980K, @@ -177,7 +178,7 @@ struct paramEnum2K2PiGam{ { static std::string fitName[paramEnum2K2PiGam::nPhaseSpace] ={"ChiGam", "K892K892", "KappaKappa", "K_0_1430K_0_1430", "ChiToK_0_1430_K892", - "K_2_1430K_2_1430", "K_0_1430K_2_1430", "K_1_1410K_1_1410","ChiToK_1_1410_K892", + "K_2_1430K_2_1430", "K_0_1430K_2_1430", "K_1_1410K_1_1410","ChiToK_1_1410_K892","ChiToK_2_1430_K892", "ChiToK_1_1400K", "K_1_1400ToK892Pi", "ChiToK_1_1270_K", "K_1_1270ToK892Pi", "K_1_1270ToK_0_1430Pi", "ChiToK_1_1650K", "K_1_1650ToK892Pi", "K_1_1650ToK_0_1430Pi","K_0_1430K_0_1950", "KappaK_0_1950","ChiToKappaK_0_1430", "ChiToK_1_2400K", "K_1_2400Tof980K", diff --git a/Examples/Psi2STo2K2PiGam/Hyp1Lh.cc b/Examples/Psi2STo2K2PiGam/Hyp1Lh.cc index 59fbc5f00d018a8404d4b3dd1ef7f489b4c87c25..33f2e8fd9ab5f28a7e346b1ff66037b908efcf2c 100644 --- a/Examples/Psi2STo2K2PiGam/Hyp1Lh.cc +++ b/Examples/Psi2STo2K2PiGam/Hyp1Lh.cc @@ -13,6 +13,7 @@ Hyp1Lh::Hyp1Lh(boost::shared_ptr<const Psi2STo2K2PiGamEvtList> theEvtList, const ,_K2_1430_K2_1430Hyp(false) ,_K0_1430_K2_1430Hyp(false) ,_K0_1430_K892Hyp1(true) + ,_K2_1430_K892Hyp1(true) ,_K1_1410_K1_1410Hyp(false) ,_K1_1410_K892Hyp1(true) ,_f1710_f1710Hyp1(true) @@ -29,6 +30,7 @@ Hyp1Lh::Hyp1Lh( boost::shared_ptr<AbsPsi2STo2K2PiGamLh> theLhPtr, const std::map ,_K2_1430_K2_1430Hyp(false) ,_K0_1430_K2_1430Hyp(false) ,_K0_1430_K892Hyp1(true) + ,_K2_1430_K892Hyp1(true) ,_K1_1410_K1_1410Hyp(false) ,_K1_1410_K892Hyp1(true) ,_f1710_f1710Hyp1(true) @@ -89,6 +91,12 @@ complex<double> Hyp1Lh::chi0DecAmps(const param2K2PiGam& theParamVal, Psi2STo2K2 std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK_0_1430_K892=theParamVal.ChiToK_0_1430_K892; result+=chiToK0K1Amp(theData, ChiToK_0_1430_K892, K_0_1430Mass, K_0_1430Width, K892Mass, K892Width); } + + if(_K2_1430_K892Hyp1){ + std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToK_2_1430_K892=theParamVal.ChiToK_2_1430_K892; + result+=chiToK1K2Amp(theData, ChiToK_2_1430_K892, K892Mass, K892Width, K_2_1430Mass, K_2_1430Width); + } + if (_K1_1410_K1_1410Hyp || _K1_1410_K892Hyp1){ double K_1_1410Mass=theParamVal.BwK_1_1410.first; double K_1_1410Width=theParamVal.BwK_1_1410.second; @@ -338,6 +346,13 @@ void Hyp1Lh::setUp(const std::map<const std::string, bool>& hypMap){ } else Alert << "hypothesis K0_1430_K892Hyp1 not set!!!" <<endmsg; + iter= hypMap.find("K2_1430_K892Hyp1"); + if (iter !=hypMap.end()){ + _K2_1430_K892Hyp1= iter->second; + _hypMap[iter->first]= iter->second; + Info<< "hypothesis " << iter->first << "\t" << _K2_1430_K892Hyp1 <<endmsg; + } + else Alert << "hypothesis K2_1430_K892Hyp1 not set!!!" <<endmsg; iter= hypMap.find("K1_1410_K1_1410Hyp"); if (iter !=hypMap.end()){ @@ -383,6 +398,7 @@ void Hyp1Lh::setUp(const std::map<const std::string, bool>& hypMap){ if(_K0_1430_K0_1430Hyp) _ampVec.push_back(paramEnum2K2PiGam::K_0_1430K_0_1430); if(_K0_1430_K2_1430Hyp) _ampVec.push_back(paramEnum2K2PiGam::K_0_1430K_2_1430); if(_K0_1430_K892Hyp1) _ampVec.push_back(paramEnum2K2PiGam::ChiToK_0_1430_K892); + if(_K2_1430_K892Hyp1) _ampVec.push_back(paramEnum2K2PiGam::ChiToK_2_1430_K892); if(_K1_1410_K1_1410Hyp || _K1_1410_K892Hyp1){ _massVec.push_back(paramEnum2K2PiGam::K_1_1410); @@ -402,7 +418,7 @@ void Hyp1Lh::setUp(const std::map<const std::string, bool>& hypMap){ _massVec.push_back(paramEnum2K2PiGam::f1710); if(_K1_1270Hyp) _massVec.push_back(paramEnum2K2PiGam::K_1_1270); if(_K1_1400Hyp) _massVec.push_back(paramEnum2K2PiGam::K_1_1400); - if(_K2_1430_K2_1430Hyp || _K0_1430_K2_1430Hyp) _massVec.push_back(paramEnum2K2PiGam::K_2_1430); + if(_K2_1430_K2_1430Hyp || _K0_1430_K2_1430Hyp || _K2_1430_K892Hyp1) _massVec.push_back(paramEnum2K2PiGam::K_2_1430); if(_K0_1430_K0_1430Hyp || _K0_1430_K2_1430Hyp || _K1_1270Hyp || _K0_1430_K892Hyp1) _massVec.push_back(paramEnum2K2PiGam::K_0_1430); diff --git a/Examples/Psi2STo2K2PiGam/Hyp1Lh.hh b/Examples/Psi2STo2K2PiGam/Hyp1Lh.hh index d131592af72d78c9fefb39dbd654330193d6f512..616c2a4f8b8e44a00d76c7fc3095b11254387621 100644 --- a/Examples/Psi2STo2K2PiGam/Hyp1Lh.hh +++ b/Examples/Psi2STo2K2PiGam/Hyp1Lh.hh @@ -59,6 +59,7 @@ protected: bool _K2_1430_K2_1430Hyp; bool _K0_1430_K2_1430Hyp; bool _K0_1430_K892Hyp1; + bool _K2_1430_K892Hyp1; bool _K1_1410_K1_1410Hyp; bool _K1_1410_K892Hyp1; bool _f1710_f1710Hyp1; diff --git a/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc b/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc index ba037679f4b09b50ce8be560df792a57003827ab..be4b10171932dbed522782d7b9e7133cb231de1e 100644 --- a/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc +++ b/Examples/Psi2STo2K2PiGam/Mpsi2STo2K2PiGamTestApp.cc @@ -148,6 +148,7 @@ int main(int __argc,char *__argv[]){ hypMap["K2_1430_K2_1430Hyp"]=true; hypMap["K0_1430_K2_1430Hyp"]=true; hypMap["K0_1430_K892Hyp1"]=true; + hypMap["K2_1430_K892Hyp1"]=true; hypMap["K1_1410_K1_1410Hyp"]=true; hypMap["K1_1410_K892Hyp1"]=true; hypMap["f1710_f1710Hyp1"]=true; diff --git a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.cc b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.cc index 4cbc7561590cd2fd18019372ba66bb9d7fdb7d82..6bcb21c6ce25055f2fee240f5e79bf67dffe8561 100644 --- a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.cc +++ b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.cc @@ -41,6 +41,9 @@ Psi2STo2K2PiGamStates::Psi2STo2K2PiGamStates() : //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 K*1 K*2 + fillJPCLS(_chic0JPC, _Kst1JPC, _Kst2JPC, _ChiToKst1Kst2JPCLS); + //LS combinations for the Chi_c0 decay to f0 f0 fillJPCLS(_chic0JPC, _f0JPC, _f0JPC,_ChiTof0f0JPCLS); @@ -177,6 +180,12 @@ void Psi2STo2K2PiGamStates::print(std::ostream& os) const os << "\n" << std::endl; } + os << "*** Chi_c0: LS combinations for the decay to K*1 K*2 *** "<< std::endl; + for ( itJPCLS=_ChiToKst1Kst2JPCLS.begin(); itJPCLS!=_ChiToKst1Kst2JPCLS.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 748e7e8b42c76b9135ce6738e9f590a489b5cd6b..ea9090c48ea11f70e355f117ba332b5e7fc071d5 100644 --- a/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.hh +++ b/Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamStates.hh @@ -36,6 +36,7 @@ public: 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> > ChiToKst1Kst2States() const {return _ChiToKst1Kst2JPCLS;} 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;} @@ -79,6 +80,7 @@ private: 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> > _ChiToKst1Kst2JPCLS; 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;