//************************************************************************//
//                                                                        //
//  Copyright 2013 Bertram Kopf (bertram@ep1.rub.de)                      //
//                 Julian Pychy (julian@ep1.rub.de)                       //
//                 - Ruhr-Universität Bochum                              //
//                                                                        //
//  This file is part of Pawian.                                          //
//                                                                        //
//  Pawian is free software: you can redistribute it and/or modify        //
//  it under the terms of the GNU General Public License as published by  //
//  the Free Software Foundation, either version 3 of the License, or     //
//  (at your option) any later version.                                   //
//                                                                        //
//  Pawian is distributed in the hope that it will be useful,             //
//  but WITHOUT ANY WARRANTY; without even the implied warranty of        //
//  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the         //
//  GNU General Public License for more details.                          //
//                                                                        //
//  You should have received a copy of the GNU General Public License     //
//  along with Pawian.  If not, see <http://www.gnu.org/licenses/>.       //
//                                                                        //
//************************************************************************//

// WaveContribution class definition file. -*- C++ -*-
// Copyright 2013 Julian Pychy

#pragma once

#include <memory>
//#include "FitParams/FitParColBase.hh"
#include "FitParams/AbsPawianParameters.hh"

struct calcContributionData {
   calcContributionData(std::string contribName, std::vector<std::string>& contribZeroAmpVec) :
     _contribName(contribName)
    ,_contribZeroAmpVec(contribZeroAmpVec)
  {}

  std::string _contribName;
  std::vector<std::string> _contribZeroAmpVec;

  virtual bool operator==(const calcContributionData& compare) const {
    bool result=false;
    if ( _contribName==compare._contribName) result=true;
    return result;
  }

 virtual bool operator<(const calcContributionData& compare) const {
   bool result=false;
   if(_contribName < compare._contribName) result=true;
    return result;
  }
};

class AbsLh;
//class fitParCol;
class EvtData;
class PwaCovMatrix;


class WaveContribution{
  
  public:
    // WaveContribution(std::shared_ptr<AbsLh> theLh, fitParCol& theFitParams);
    // WaveContribution(std::shared_ptr<AbsLh> theLh, fitParCol& theFitParams,
    //  		     std::shared_ptr<PwaCovMatrix> thePwaCovMatrix);

  WaveContribution(std::shared_ptr<AbsLh> theLh, std::shared_ptr<AbsPawianParameters> theFitParams);
  WaveContribution(std::shared_ptr<AbsLh> theLh, std::shared_ptr<AbsPawianParameters> theFitParams,
		     std::shared_ptr<PwaCovMatrix> thePwaCovMatrix);

    std::string GetContributionName(unsigned int index);
    std::shared_ptr<AbsPawianParameters> GetParametersForContribution(unsigned int index);
    unsigned int NoOfContributions();
    std::pair<double,double> CalcContribution();
    std::vector<std::pair<std::string,std::pair<double,double>>> CalcSingleContributions();
    double CalcError(double result, std::shared_ptr<AbsPawianParameters> currentParameters);

   private:
    bool _calcError;
    std::shared_ptr<AbsLh> _theLh;
    std::shared_ptr<PwaCovMatrix> _thePwaCovMatrix;
    // FitParColBase _theFitParColBase;
    // fitParCol* _theFitParamsOriginal;
    std::vector<EvtData*> _MCDataList;

  std::shared_ptr<AbsPawianParameters> _fitParamsOriginal;
 
    // double CalcContribution(fitParCol& theFitParams);
  double CalcContribution(std::shared_ptr<AbsPawianParameters> theParams);
};