#ifndef _Hyp8Lh_H #define _Hyp8Lh_H #include <iostream> #include <fstream> #include <string> #include <vector> #include <complex> #include <cassert> #include <boost/shared_ptr.hpp> #include "TROOT.h" // #include <TSystem.h> #include "qft++/topincludes/relativistic-quantum-mechanics.hh" #include "Examples/Psi2STo2K2PiGam/Hyp7Lh.hh" #include "Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamData.hh" #include "PwaUtils/DataUtils.hh" #include "Minuit2/MnUserParameters.h" class Hyp8Lh : public Hyp7Lh{ public: // create/copy/destroy: ///Constructor Hyp8Lh(boost::shared_ptr<const Psi2STo2K2PiGamEvtList>, const std::map<const std::string, bool>& hypMap, boost::shared_ptr<Psi2STo2K2PiGamStates> theStatesPtr, bool cacheAmps=false); Hyp8Lh(boost::shared_ptr<AbsPsi2STo2K2PiGamLh>, const std::map<const std::string, bool>& hypMap, boost::shared_ptr<Psi2STo2K2PiGamStates> theStatesPtr, bool cacheAmps=false); /** Destructor */ virtual ~Hyp8Lh(); virtual AbsPsi2STo2K2PiGamLh* clone_() { return new Hyp8Lh(_Psi2STo2K2PiGamEvtListPtr, _hypMap, _fitParams2K2PiGam.states(), _cacheAmps); } // Getters: virtual void setMnUsrParams(MnUserParameters& upar, param2K2PiGam& startVal, param2K2PiGam& errVal); virtual int setFitParamVal(param2K2PiGam& theParamVal, const std::vector<double>& par); virtual unsigned int nFitParams(); virtual void print(std::ostream& os) const; virtual void printCurrentFitResult(param2K2PiGam& theParamVal); virtual void dumpCurrentResult(std::ostream& os, param2K2PiGam& theParamVal, std::string& suffix); protected: bool _K_0_1460ToKstPiHyp8; bool _K_0_1460ToSigmaKHyp8; bool _K_0_1830ToKstPiHyp8; bool _K_1_1650Hyp8; bool _doHyp8; virtual complex<double> chi0DecAmps(const param2K2PiGam& theParamVal, Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData); virtual bool equalChic0DecParams(); virtual void copyCurrentVals(Hyp8Lh* theLh); std::map<unsigned int, complex<double> > _currentResultHyp8; private: unsigned int _nFitParams; std::vector<unsigned int> _ampVec; std::vector<unsigned int> _massVec; std::vector<unsigned int> _massVecRemain; void setUp(const std::map<const std::string, bool>& hypMap); bool _equalParameter; bool equalParams(); }; #endif