Hyp4Lh.cc

#include <getopt.h>
#include <fstream>
#include <string>

#include "Examples/Psi2STo2K2PiGam/Hyp4Lh.hh"
#include "Examples/Psi2STo2K2PiGam/Psi2STo2K2PiGamEvtList.hh"
#include "ErrLogger/ErrLogger.hh"

Hyp4Lh::Hyp4Lh(boost::shared_ptr<const Psi2STo2K2PiGamEvtList> theEvtList, const std::map<const std::string, bool>& hypMap, bool cacheAmps ) :
  Hyp3Lh(theEvtList, hypMap, cacheAmps)
  ,_f980f1370Hyp4(true)
  ,_f980f1500Hyp4(true)
  ,_sigmaf1370Hyp4(true)
  ,_f1710f1370Hyp4(true)
  ,_f980f_2_1270Hyp4(true)
  ,_f980f_2_1430Hyp4(true)
  ,_f980f_2_1525Hyp4(true)
  ,_f980f_2_1950Hyp4(true)
  ,_f980f_2_2010Hyp4(true)
  ,_f1500f_2_1525Hyp4(true)
  ,_f1710f_2_1430Hyp4(true)
  ,_f1710f_2_1950Hyp4(true)
  ,_doHyp4(true)
  ,_nFitParams(0)
  ,_f980FlatteRemain(false)
  ,_evtCounter(0)
  ,_equalParameter(false)
{
  setUp(hypMap);
}

Hyp4Lh::Hyp4Lh( boost::shared_ptr<AbsPsi2STo2K2PiGamLh> theLhPtr, const std::map<const std::string, bool>& hypMap, bool cacheAmps  ) :
  Hyp3Lh(theLhPtr->getEventList(), hypMap, cacheAmps)
  ,_f980f1370Hyp4(true)
  ,_f980f1500Hyp4(true)
  ,_sigmaf1370Hyp4(true)
  ,_f1710f1370Hyp4(true)
  ,_f980f_2_1270Hyp4(true)
  ,_f980f_2_1430Hyp4(true)
  ,_f980f_2_1525Hyp4(true)
  ,_f980f_2_1950Hyp4(true)
  ,_f980f_2_2010Hyp4(true)
  ,_f1500f_2_1525Hyp4(true)
  ,_f1710f_2_1430Hyp4(true)
  ,_f1710f_2_1950Hyp4(true)
  ,_doHyp4(true)
  ,_nFitParams(0)
  ,_f980FlatteRemain(false)
  ,_evtCounter(0)
  ,_equalParameter(false)
{
  setUp(hypMap);
}

Hyp4Lh::~Hyp4Lh()
{;
}

bool  Hyp4Lh::equalChic0DecParams(){
  bool result=false; 
  bool equalRemainHyps=Hyp3Lh::equalChic0DecParams();
  if(!_doHyp4) return equalRemainHyps;

  _equalParameter=equalParams();
  DebugMsg << "equal parameter: "<< _equalParameter << endmsg;

  if(_equalParameter && equalRemainHyps) result=true;
  return result;
}

complex<double> Hyp4Lh::chi0DecAmps(const param2K2PiGam& theParamVal, Psi2STo2K2PiGamData::Psi2STo2K2PiGamEvtData* theData){

  complex<double> result=Hyp3Lh::chi0DecAmps(theParamVal, theData);

  if(!_doHyp4) return result;

  if(_equalParameter){
    result+=_currentResultHyp4[_evtCounter];
    _evtCounter++;
    return result;
  }


  complex<double> currentResult(0.,0.);
  
  double f1370Mass=theParamVal.Bwf1370.first;
  double f1370Width=theParamVal.Bwf1370.second;

  double f_2_1430Mass=theParamVal.Bwf_2_1430.first;
  double f_2_1430Width=theParamVal.Bwf_2_1430.second;

  double f_2_1525Mass=theParamVal.Bwf_2_1525.first;
  double f_2_1525Width=theParamVal.Bwf_2_1525.second;

  double f_2_1950Mass=theParamVal.Bwf_2_1950.first;
  double f_2_1950Width=theParamVal.Bwf_2_1950.second;

  double f_2_2010Mass=theParamVal.Bwf_2_2010.first;
  double f_2_2010Width=theParamVal.Bwf_2_2010.second;

  if (_f980f1370Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_pif1370_k=theParamVal.f980_pif1370_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_kf1370_pi=theParamVal.f980_kf1370_pi;
    double f980_Mass=theParamVal.Flatf980;
    double f980_gPiPi=theParamVal.Flatf980gPiPi;
    double f980_gKK=theParamVal.Flatf980gKK;
    
    //Chi_c0 decay to f0(980) f0(1370) -> (pi0 pi0) (K K) 
    currentResult+=chiTof980_pif0_kAmp(theData, ChiTof980_pif1370_k, f980_Mass, f980_gPiPi,  f980_gKK, f1370Mass, f1370Width);
    currentResult+=chiTof980_kf0_piAmp(theData, ChiTof980_kf1370_pi, f980_Mass, f980_gKK,  f980_gPiPi, f1370Mass, f1370Width);
  }

  double f1500Mass=theParamVal.Bwf1500.first;
  double f1500Width=theParamVal.Bwf1500.second;

  if (_f980f1500Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_pif1500_k=theParamVal.f980_pif1500_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_kf1500_pi=theParamVal.f980_kf1500_pi;
    double f980_Mass=theParamVal.Flatf980;
    double f980_gPiPi=theParamVal.Flatf980gPiPi;
    double f980_gKK=theParamVal.Flatf980gKK;
    
    //Chi_c0 decay to f0(980) f0(1500) -> (pi0 pi0) (K K) 
    currentResult+=chiTof980_pif0_kAmp(theData, ChiTof980_pif1500_k, f980_Mass, f980_gPiPi,  f980_gKK, f1500Mass, f1500Width);
    //Chi_c0 decay to f0(980) f0(1500) -> (K K) (pi0 pi0) 
    currentResult+=chiTof980_kf0_piAmp(theData, ChiTof980_kf1500_pi, f980_Mass, f980_gKK,  f980_gPiPi, f1500Mass, f1500Width);
  }

  if (_f1710f1370Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof1710_pif1370_k=theParamVal.f1710_pif1370_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof1710_kf1370_pi=theParamVal.f1710_kf1370_pi;
    double f1710Mass=theParamVal.Bwf1710.first;
    double f1710Width=theParamVal.Bwf1710.second;

    currentResult+=chiTof0_pif0_kAmp(theData, ChiTof1710_pif1370_k, f1710Mass, f1710Width,  f1370Mass, f1370Width);
    currentResult+=chiTof0_pif0_kAmp(theData, ChiTof1710_kf1370_pi, f1370Mass, f1370Width,  f1710Mass, f1710Width);
  }

  if(_f980f_2_1270Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > f980_pif_2_1270_k=theParamVal.f980_pif_2_1270_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > f980_kf_2_1270_pi=theParamVal.f980_kf_2_1270_pi;

    double f980_Mass=theParamVal.Flatf980;
    double f980_gPiPi=theParamVal.Flatf980gPiPi;
    double f980_gKK=theParamVal.Flatf980gKK;
    double f_2_1270Mass=theParamVal.Bwf_2_1270.first;
    double f_2_1270Width=theParamVal.Bwf_2_1270.second;
    currentResult+=chiTof980_pif2_kAmp(theData, f980_pif_2_1270_k, f980_Mass, f980_gPiPi, f980_gKK, f_2_1270Mass, f_2_1270Width);
    currentResult+=chiTof980_kf2_piAmp(theData, f980_kf_2_1270_pi, f980_Mass, f980_gKK, f980_gPiPi, f_2_1270Mass, f_2_1270Width);
  }

  if (_f980f_2_1430Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_pif2_k=theParamVal.f980_pif_2_1430_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_kf2_pi=theParamVal.f980_kf_2_1430_pi;

    double f980_Mass=theParamVal.Flatf980;
    double f980_gPiPi=theParamVal.Flatf980gPiPi;
    double f980_gKK=theParamVal.Flatf980gKK;
    currentResult+=chiTof980_pif2_kAmp(theData, ChiTof980_pif2_k, f980_Mass, f980_gPiPi, f980_gKK, f_2_1430Mass, f_2_1430Width);
    currentResult+=chiTof980_kf2_piAmp(theData, ChiTof980_kf2_pi, f980_Mass, f980_gKK, f980_gPiPi, f_2_1430Mass, f_2_1430Width);
  }

  if (_f980f_2_1525Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_pif2_k=theParamVal.f980_pif_2_1525_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_kf2_pi=theParamVal.f980_kf_2_1525_pi;

    double f980_Mass=theParamVal.Flatf980;
    double f980_gPiPi=theParamVal.Flatf980gPiPi;
    double f980_gKK=theParamVal.Flatf980gKK;
    currentResult+=chiTof980_pif2_kAmp(theData, ChiTof980_pif2_k, f980_Mass, f980_gPiPi, f980_gKK, f_2_1525Mass, f_2_1525Width);
    currentResult+=chiTof980_kf2_piAmp(theData, ChiTof980_kf2_pi, f980_Mass, f980_gKK, f980_gPiPi, f_2_1525Mass, f_2_1525Width);
  }

  if (_f980f_2_1950Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_pif2_k=theParamVal.f980_pif_2_1950_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_kf2_pi=theParamVal.f980_kf_2_1950_pi;

    double f980_Mass=theParamVal.Flatf980;
    double f980_gPiPi=theParamVal.Flatf980gPiPi;
    double f980_gKK=theParamVal.Flatf980gKK;
    currentResult+=chiTof980_pif2_kAmp(theData, ChiTof980_pif2_k, f980_Mass, f980_gPiPi, f980_gKK, f_2_1950Mass, f_2_1950Width);
    currentResult+=chiTof980_kf2_piAmp(theData, ChiTof980_kf2_pi, f980_Mass, f980_gKK, f980_gPiPi, f_2_1950Mass, f_2_1950Width);
  }

  if (_f980f_2_2010Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_pif2_k=theParamVal.f980_pif_2_2010_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof980_kf2_pi=theParamVal.f980_kf_2_2010_pi;

    double f980_Mass=theParamVal.Flatf980;
    double f980_gPiPi=theParamVal.Flatf980gPiPi;
    double f980_gKK=theParamVal.Flatf980gKK;
    currentResult+=chiTof980_pif2_kAmp(theData, ChiTof980_pif2_k, f980_Mass, f980_gPiPi, f980_gKK, f_2_2010Mass, f_2_2010Width);
    currentResult+=chiTof980_kf2_piAmp(theData, ChiTof980_kf2_pi, f980_Mass, f980_gKK, f980_gPiPi, f_2_2010Mass, f_2_2010Width);
  }

  if (_f1500f_2_1525Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof1500_pif2_k=theParamVal.f1500_pif_2_1525_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof1500_kf2_pi=theParamVal.f1500_kf_2_1525_pi;

  currentResult+=chiTof0_pif2_kAmp(theData, ChiTof1500_pif2_k, f1500Mass, f1500Width, f_2_1525Mass, f_2_1525Width);
  currentResult+=chiTof2_pif0_kAmp(theData, ChiTof1500_kf2_pi, f_2_1525Mass, f_2_1525Width, f1500Mass, f1500Width);
  }

  if (_f1710f_2_1430Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof1710_pif2_k=theParamVal.f1710_pif_2_1430_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof1710_kf2_pi=theParamVal.f1710_kf_2_1430_pi;

  double f1710Mass=theParamVal.Bwf1710.first;
  double f1710Width=theParamVal.Bwf1710.second;

  currentResult+=chiTof0_pif2_kAmp(theData, ChiTof1710_pif2_k, f1710Mass, f1710Width, f_2_1430Mass, f_2_1430Width);
  currentResult+=chiTof2_pif0_kAmp(theData, ChiTof1710_kf2_pi, f_2_1430Mass, f_2_1430Width, f1710Mass, f1710Width);
  }

  if (_f1710f_2_1950Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof1710_pif2_k=theParamVal.f1710_pif_2_1950_k;
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiTof1710_kf2_pi=theParamVal.f1710_kf_2_1950_pi;

  double f1710Mass=theParamVal.Bwf1710.first;
  double f1710Width=theParamVal.Bwf1710.second;

  currentResult+=chiTof0_pif2_kAmp(theData, ChiTof1710_pif2_k, f1710Mass, f1710Width, f_2_1950Mass, f_2_1950Width);
  currentResult+=chiTof2_pif0_kAmp(theData, ChiTof1710_kf2_pi, f_2_1950Mass, f_2_1950Width, f1710Mass, f1710Width);
  }

  if (_sigmaf1370Hyp4){
    std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > ChiToSigmaf1370=theParamVal.ChiToSigmaf1370;
    
    double sigmaMass=theParamVal.BwSigma.first;
    double sigmaWidth=theParamVal.BwSigma.second;
    
    currentResult+=chiTof0_pif0_kAmp(theData, ChiToSigmaf1370, sigmaMass, sigmaWidth,  f1370Mass, f1370Width);
    
  }

  _currentResultHyp4[_evtCounter]=currentResult; 
  _evtCounter++;
  result+=currentResult;
  return result;
}

void Hyp4Lh::setMnUsrParams(MnUserParameters& upar, param2K2PiGam& startVal, param2K2PiGam& errVal){

  Hyp3Lh::setMnUsrParams(upar, startVal, errVal);

  if (!_doHyp4) return;

  std::vector<unsigned int>::const_iterator itAmps;
  for ( itAmps=_ampVec.begin(); itAmps!=_ampVec.end(); ++itAmps){

    _fitParams2K2PiGam.setMnUsrParamsDec(upar, startVal, errVal, (*itAmps));
  }

  std::vector<unsigned int>::const_iterator itMasses;
  for ( itMasses=_massVec.begin(); itMasses!=_massVec.end(); ++itMasses){
    _fitParams2K2PiGam.setMnUsrParamsMass(upar, startVal, errVal, (*itMasses) );  
  }

}



int Hyp4Lh::setFitParamVal(param2K2PiGam& theParamVal, const std::vector<double>& par){

  if (par.size() != nFitParams() ) {
    Alert << "size of parameters wrong!!! par.size()=" << par.size() << 
      "\t it should be" << nFitParams() << endmsg;
    exit(1);
  }  

  int counter=Hyp3Lh::setFitParamVal(theParamVal, par);

  if (!_doHyp4) return counter;

  std::vector<unsigned int>::const_iterator itAmps;
  for ( itAmps=_ampVec.begin(); itAmps!=_ampVec.end(); ++itAmps){
    counter=_fitParams2K2PiGam.setFitParamValDec(theParamVal, par, counter, (*itAmps));
  } 

  std::vector<unsigned int>::const_iterator itMasses;
  for ( itMasses=_massVec.begin(); itMasses!=_massVec.end(); ++itMasses){
    counter=_fitParams2K2PiGam.setFitParamValMass(theParamVal, par, counter, (*itMasses) ); 
  }

  return counter;
}

unsigned int  Hyp4Lh::nFitParams(){
  unsigned int remainingFitParams=Hyp3Lh::nFitParams();
  return _nFitParams+remainingFitParams;
}

void Hyp4Lh::print(std::ostream& os) const{
  os << "Hyp4Lh::print\n";
}

void Hyp4Lh::printCurrentFitResult(param2K2PiGam& theParamVal){
  Hyp3Lh::printCurrentFitResult(theParamVal);

  if (!_doHyp4) return;

  std::vector<unsigned int>::const_iterator itAmps;
  for ( itAmps=_ampVec.begin(); itAmps!=_ampVec.end(); ++itAmps){
    std::vector< boost::shared_ptr<const JPCLS> > JPCLSs=_fitParams2K2PiGam.jpclsVec(*itAmps);

    std::vector< boost::shared_ptr<const JPCLS> >::const_iterator itJPCLS;    
    
    for ( itJPCLS=JPCLSs.begin(); itJPCLS!=JPCLSs.end(); ++itJPCLS){
      DebugMsg<< (*itJPCLS)->name()<< paramEnum2K2PiGam::name(*itAmps) << endmsg;
      std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > currentMap=_fitParams2K2PiGam.ampMap(theParamVal, *itAmps);
      std::pair<double, double> tmpParam=currentMap[(*itJPCLS)];
      DebugMsg <<"\t mag:" << tmpParam.first <<"\t phi:" << tmpParam.second  << endmsg;
    }  
  }
  
  std::vector<unsigned int>::const_iterator itMasses;
  for ( itMasses=_massVec.begin(); itMasses!=_massVec.end(); ++itMasses){
    DebugMsg<< paramEnum2K2PiGam::name(*itMasses) << endmsg;
    std::pair<double, double> tmpParam=_fitParams2K2PiGam.massPair(theParamVal, *itMasses);
    DebugMsg <<"\t mag:" << tmpParam.first <<"\t phi:" << tmpParam.second  << endmsg;
  }

}

void Hyp4Lh::dumpCurrentResult(std::ostream& os, param2K2PiGam& theParamVal, std::string& suffix){
 
  if ( suffix.compare("Val") != 0 && suffix.compare("Err") !=0 ){
    Warning << "suffix " << suffix << " not supported!!! Use Val or Err" << endmsg;
    return;
  }
  Hyp3Lh::dumpCurrentResult(os, theParamVal, suffix);

  if (!_doHyp4) return;

  std::vector<unsigned int>::const_iterator itAmps;
  for ( itAmps=_ampVec.begin(); itAmps!=_ampVec.end(); ++itAmps){
    std::vector< boost::shared_ptr<const JPCLS> > JPCLSs=_fitParams2K2PiGam.jpclsVec(*itAmps);

    std::vector< boost::shared_ptr<const JPCLS> >::const_iterator itJPCLS;    
    
    for ( itJPCLS=JPCLSs.begin(); itJPCLS!=JPCLSs.end(); ++itJPCLS){
      std::map< boost::shared_ptr<const JPCLS>, pair<double, double>, pawian::Collection::SharedPtrLess > currentMap=_fitParams2K2PiGam.ampMap(theParamVal, *itAmps);
      std::pair<double, double> tmpParam=currentMap[(*itJPCLS)];

      std::string tmpStringDec=(*itJPCLS)->name()+paramEnum2K2PiGam::name(*itAmps)+suffix;
      os << tmpStringDec << "\t" << tmpParam.first  << "\t" << tmpParam.second << std::endl;
    }  
  }

  std::vector<unsigned int>::const_iterator itMasses;
  for ( itMasses=_massVec.begin(); itMasses!=_massVec.end(); ++itMasses){
    std::string tmpStringMass=paramEnum2K2PiGam::name(*itMasses)+suffix;

    std::pair<double, double> tmpParam=_fitParams2K2PiGam.massPair(theParamVal, *itMasses);
    os << tmpStringMass << "\t" << tmpParam.first  << "\t" << tmpParam.second << std::endl;
  }

}

void Hyp4Lh::setUp(const std::map<const std::string, bool>& hypMap){

  std::map<const std::string, bool>::const_iterator iter= hypMap.find("f980f1370Hyp4");

  if (iter !=hypMap.end()){
    _f980f1370Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f980f1370Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f980f1370Hyp4 not set!!!" <<endmsg;
    exit(0);
  } 

  iter= hypMap.find("f980f1500Hyp4");

  if (iter !=hypMap.end()){
    _f980f1500Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f980f1500Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f980f1500Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

  iter= hypMap.find("f1710f1370Hyp4");

  if (iter !=hypMap.end()){
    _f1710f1370Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f1710f1370Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f1710f1370Hyp4 not set!!!" <<endmsg;
    exit(0);
  }


  iter= hypMap.find("f980f_2_1270Hyp4");

  if (iter !=hypMap.end()){
    _f980f_2_1270Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f980f_2_1270Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f980f_2_1270Hyp4 not set!!!" <<endmsg;
    exit(0);
  }


  iter= hypMap.find("f980f_2_1430Hyp4");

  if (iter !=hypMap.end()){
    _f980f_2_1430Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f980f_2_1430Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f980f_2_1430Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

  iter= hypMap.find("f980f_2_1525Hyp4");
  if (iter !=hypMap.end()){
    _f980f_2_1525Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f980f_2_1525Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f980f_2_1525Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

  iter= hypMap.find("f980f_2_1950Hyp4");
  if (iter !=hypMap.end()){
    _f980f_2_1950Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f980f_2_1950Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f980f_2_1950Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

  iter= hypMap.find("f980f_2_2010Hyp4");
  if (iter !=hypMap.end()){
    _f980f_2_2010Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f980f_2_2010Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f980f_2_2010Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

  iter= hypMap.find("f980f1500Hyp4");
  if (iter !=hypMap.end()){
    _f980f1500Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f980f1500Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f980f1500Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

  iter= hypMap.find("f1500f_2_1525Hyp4");

  if (iter !=hypMap.end()){
    _f1500f_2_1525Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f1500f_2_1525Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f1500f_2_1525Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

  iter= hypMap.find("f1710f_2_1430Hyp4");

  if (iter !=hypMap.end()){
    _f1710f_2_1430Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f1710f_2_1430Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f1710f_2_1430Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

  iter= hypMap.find("f1710f_2_1950Hyp4");

  if (iter !=hypMap.end()){
    _f1710f_2_1950Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _f1710f_2_1950Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis f1710f_2_1950Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

 iter= hypMap.find("sigmaf1370Hyp4");

  if (iter !=hypMap.end()){
    _sigmaf1370Hyp4= iter->second;
    Info<< "hypothesis " << iter->first << "\t" << _sigmaf1370Hyp4 <<endmsg;
    _hypMap[iter->first]= iter->second;
  }
  else{
    Alert << "hypothesis sigmaf1370Hyp4 not set!!!" <<endmsg;
    exit(0);
  }

  if(!_f980f1370Hyp4 && !_f1710f1370Hyp4 && !_f980f_2_1430Hyp4 && !_f1710f_2_1430Hyp4 && !_f980f_2_1950Hyp4 && !_f980f_2_2010Hyp4 && !_f1710f_2_1950Hyp4 && !_sigmaf1370Hyp4) _doHyp4=false; 

  if (!_doHyp4) return;

  if (_f980f1370Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::f980_pif1370_k);
    _ampVec.push_back(paramEnum2K2PiGam::f980_kf1370_pi);
  }

  if (_f980f1500Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::f980_pif1500_k);
    _ampVec.push_back(paramEnum2K2PiGam::f980_kf1500_pi);
  }

  if (_f1710f1370Hyp4){
    _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);
    _massVec.push_back(paramEnum2K2PiGam::f_2_1270);
  }
  if (_f980f_2_1430Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::f980_pif_2_1430_k);
    _ampVec.push_back(paramEnum2K2PiGam::f980_kf_2_1430_pi);
  }
  if (_f980f_2_1525Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::f980_pif_2_1525_k);
    _ampVec.push_back(paramEnum2K2PiGam::f980_kf_2_1525_pi);
  }
  if (_f980f_2_1950Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::f980_pif_2_1950_k);
    _ampVec.push_back(paramEnum2K2PiGam::f980_kf_2_1950_pi);
  }
  if (_f980f_2_2010Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::f980_pif_2_2010_k);
    _ampVec.push_back(paramEnum2K2PiGam::f980_kf_2_2010_pi);
  }
  if (_f1500f_2_1525Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::f1500_pif_2_1525_k);
    _ampVec.push_back(paramEnum2K2PiGam::f1500_kf_2_1525_pi);
  }
  if (_f1710f_2_1430Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::f1710_pif_2_1430_k);
    _ampVec.push_back(paramEnum2K2PiGam::f1710_kf_2_1430_pi);
  }
  if (_f1710f_2_1950Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::f1710_pif_2_1950_k);
    _ampVec.push_back(paramEnum2K2PiGam::f1710_kf_2_1950_pi);
  }
  if(_sigmaf1370Hyp4){
    _ampVec.push_back(paramEnum2K2PiGam::ChiToSigmaf1370);
  }
  if (_f980f1370Hyp4 || _f1710f1370Hyp4 || _sigmaf1370Hyp4) _massVec.push_back(paramEnum2K2PiGam::f1370);
  if (_f980f1500Hyp4 || _f1500f_2_1525Hyp4) _massVec.push_back(paramEnum2K2PiGam::f1500);
  if (_f980f_2_1430Hyp4 || _f1710f_2_1430Hyp4) _massVec.push_back(paramEnum2K2PiGam::f_2_1430);
  if (_f980f_2_1950Hyp4 || _f1710f_2_1950Hyp4) _massVec.push_back(paramEnum2K2PiGam::f_2_1950);
  if (_f980f_2_2010Hyp4) _massVec.push_back(paramEnum2K2PiGam::f_2_2010);
  if (_f980f_2_1525Hyp4 || _f1500f_2_1525Hyp4) _massVec.push_back(paramEnum2K2PiGam::f_2_1525);

  // fill all other resonances
  if (_sigmaf1370Hyp4){
    _massVecRemain.push_back(paramEnum2K2PiGam::sigma);
  }

  if (_f1710f1370Hyp4 || _f1710f_2_1430Hyp4 || _f1710f_2_1950Hyp4){
      _massVecRemain.push_back(paramEnum2K2PiGam::f1710);
  }

  if(_f980f1370Hyp4 || _f980f1500Hyp4 || _f980f_2_1270Hyp4 || _f980f_2_1430Hyp4 || _f980f_2_1525Hyp4 || _f980f_2_1950Hyp4 || _f980f_2_2010Hyp4){
    _f980FlatteRemain=true;
  }


  std::vector<unsigned int>::iterator ampIt;
  for (ampIt=_ampVec.begin(); ampIt!=_ampVec.end(); ++ampIt){
    std::vector< boost::shared_ptr<const JPCLS> > JPCLSs=_fitParams2K2PiGam.jpclsVec(*ampIt);
    _nFitParams+=2*JPCLSs.size();
  }

  std::vector<unsigned int>::iterator massIt; 
  for (massIt=_massVec.begin(); massIt!=_massVec.end(); ++massIt){
    _nFitParams+=2;
  }
}



void Hyp4Lh::copyCurrentVals(Hyp4Lh* theLh){
  Hyp3Lh::copyCurrentVals(theLh);
  std::map<unsigned int, complex<double> > newResult; 
  std::map<unsigned int, complex<double> >::iterator it;
  for (it= _currentResultHyp4.begin(); it!= _currentResultHyp4.end(); ++it){
    newResult[it->first]=it->second;
  }
  theLh->_currentResultHyp4=newResult;
  
}



bool Hyp4Lh::equalParams(){

  bool result=true;
  std::vector< boost::shared_ptr<const JPCLS> >::const_iterator itJPCLS;

  if (!compAmpParms( _ampVec )) return false;
  if (!compMassParms(_massVec)) return false;
  if (!compMassParms(_massVecRemain)) return false;
  if(_f980FlatteRemain){
    if (!compFlatteParms()) return false;
  }
  return result;

}