RdGlobalFit.h

Go to the documentation of this file.

#ifndef _RdGlobalFit_h_
#define _RdGlobalFit_h_

#include <fwk/VModule.h>
#include <string>
#include <rdet/RDetector.h>
#include <revt/REvent.h>
#include <utl/CoordinateSystemPtr.h>
#include <evt/ShowerRRecData.h>
#include <vector>
#include <algorithm>
#include "math.h"

#include "TF2.h"

#include <revt/Station.h>
#include <revt/Channel.h>
#include <string>
#include <vector>
#include <utl/CoordinateSystem.h>

#include <utl/Point.h>
#include "FitModels.h"
#include "Minuit2/MnUserCovariance.h"
#include "Minuit2/MnUserParameters.h"

#include "Minuit2/Minuit2Minimizer.h"
#include "Minuit2/MnUserParameters.h"
#include "Minuit2/MnMinimize.h"
#include "Minuit2/MnContours.h"
#include "Minuit2/MnMigrad.h"
#include "Minuit2/MnContours.h"
#include "Minuit2/ContoursError.h"
#include "Minuit2/FunctionMinimum.h"
#include "Minuit2/MinimumParameters.h"
#include "Minuit2/MnPrint.h"
#include "Minuit2/MnScan.h"
#include "Minuit2/MnPlot.h"
#include "Minuit2/MnUserCovariance.h"
namespace revt {
class Channel;
}

namespace evt {
class Event;
}

namespace revt {
class Station;
}

namespace RdGlobalFit {

class RdGlobalFit : public fwk::VModule {

 public:
  // If you define a constructor, you MUST provide an implementation
  // for it, otherwise the REGISTER_MODULE macro will not register
  // your module. (Here, implementation is in MyModule.cc)
  RdGlobalFit();
  virtual ~RdGlobalFit();

  fwk::VModule::ResultFlag Init();
  fwk::VModule::ResultFlag Run(evt::Event& event);
  fwk::VModule::ResultFlag Finish();
  std::vector<std::pair<std::vector<std::pair<double, double> >, std::string> > scanParameters(
      RdGlobalFitMinimizationCriterion,
      ROOT::Minuit2::FunctionMinimum/*,  ROOT::Minuit2::MnStrategy*/, std::vector<std::string>);
  void printPlotsVec(std::vector<std::pair<std::vector<std::pair<double, double> >, std::string> >);
  utl::CoordinateSystemPtr getlocalCSPtr(evt::Event& event);
  utl::Point adaptSDCoreTofLocalCS(evt::Event& event);
  void setGlFitDataStartValues(evt::Event& event);
  void setGlFitDataCore(evt::Event& event);
  void setGlFitDataDirection(evt::Event& event);
  int calcNDF(ROOT::Minuit2::MnUserParameters upar, int nrSigSt);
  double XmaxFromRmax(evt::Event& event, double rmax, double zenithInLocalCS);
  utl::Point getCoreInCSPlane(utl::Point core, utl::Vector& axis,
                              utl::CoordinateSystemPtr targetCS);
 private:
  // The REGISTER_MODULE macro makes the framework aware of your module.
  // You need to put this macro at the very end of your module class definition.
  // Give your module a name (first argument) and specify its
  // class name (second argument).  In the ModuleSequenceExample.xml file,
  // you can see this module being refered to by the name we give it here.
  utl::CoordinateSystemPtr fLocalCS;
  enum InfoLevel {
    eNone = 0,
    eFinal = 1,
    eIntermediate = 2,
    eObscure = 3,
    eDebug = 4
  };

  int fInfoLevel;

  FitConfig fFitConfig;
  FitConfig fOrigFitConfig;
  RdGlobalFitData fGlobalFitData;
  LDFConstsTable fLDFConstsTable;
  calcLDFConsts fcalcLDFConsts;

REGISTER_MODULE("RdGlobalFit", RdGlobalFit);
};

}

#endif