VMinMethodFunctor.h

Go to the documentation of this file.

#ifndef _MdLDFFinderAG_VMinMethodFunctor_h_
#define _MdLDFFinderAG_VMinMethodFunctor_h_

#include <mevt/MEvent.h>
#include <mdet/MDetector.h>
#include <utl/AxialVector.h>

#include <boost/tuple/tuple.hpp>
#include <vector>
#include <ostream>

#include "VLDFFunctor.h"

#include <Minuit2/MnMinimize.h>
#include <Minuit2/MnHesse.h>
#include <Minuit2/FCNBase.h>
#include <Minuit2/FunctionMinimum.h>
#include <Minuit2/MnUserParameters.h>
#include <Minuit2/MnUserCovariance.h>
#include <Minuit2/MnPrint.h>

namespace MdLDFFinderAG {

   class VLDFFunctor;

//   typedef std::vector< boost::tuple<double, unsigned int , std::vector< unsigned int > > > ModuleList; // area projected in the shower plane, number of strips, strips with signals in each time window
  typedef std::vector<mevt::Counter*> CounterListBase;

  class CounterList: public CounterListBase {
    public:
      friend std::ostream& operator<<(std::ostream& os, const CounterList& dt);
  };

  class VMinMethodFunctor: public ROOT::Minuit2::FCNBase {
    
  public:

  enum FitParameter {
    eShowerSize = 0,
    eBeta = 1,
    eCoreX = 2,
    eCoreY = 3
  };

  virtual double operator()( const std::vector<double>& par ) const = 0;
  virtual double Up() const = 0;
  ROOT::Minuit2::FunctionMinimum Minimize( std::vector<double>& par, std::vector<double>& epar, const std::vector<bool>& fpar ) const;  
  virtual ~VMinMethodFunctor();

  void SetZcore(double z) {zcore = z;}
  double GetZcore() const { return zcore; }

  void SetAxis(utl::Vector axis);
  utl::Vector GetAxis() const { return rAxis; }

  double DistanceInShowerPlane(const utl::Vector v) const {
      return utl::cross( rAxis, v ).GetMag() ;
  }

  void SetCandidateCounterList(CounterList cl) { fCandidateCounters=cl; }
  void SetSilentCounterList(CounterList cl) { fSilentCounters=cl; } 
      
  protected:

   double zcore;  // z coordinate of the core position (fixed parameter for the reconstruction)
   utl::Vector rAxis;
   double cosTheta;

   const VLDFFunctor* fLDFunction;

   bool fUseSilent;
   unsigned int fSilentLimit;

   CounterList fCandidateCounters;
   CounterList fSilentCounters;

   utl::CoordinateSystemPtr siteCS;

   double CalculateSilentLikelihood(const std::vector<double>& par) const;

  };

}
#endif // _MdLDFFinderAG_VMinMethodFunctor_h_