PhotoelectricClusterTools.h

/*
 * PhotoelectricClusterTools.h
 *
 *  Created on: Apr 24, 2014
 *      Author: zenker
 */

#ifndef PHOTOELECTRICCLUSTERTOOLS_H_
#define PHOTOELECTRICCLUSTERTOOLS_H_

#include <set>
#include <vector>

// CLHEP
#include <CLHEP/Vector/TwoVector.h>

// LCIO
#include <lcio.h>
#include <EVENT/TrackerPulse.h>
#include <IMPL/TrackerHitImpl.h>

// GEAR
#include <gear/TPCParameters.h>
#include <gearimpl/GlobalPadIndex.h>

/* NOTE this uses the C++ version of Minuit - not the Fortran version */
#include <Minuit2/FCNBase.h>

namespace marlintpc{

CLHEP::Hep2Vector calculateSeed(const std::set<lcio::TrackerPulse*> &pulses, const gear::TPCParameters& tpc_parameters, IMPL::TrackerHitImpl *seed_photodot = 0);

bool doLikelihoodFit(const lcio::TrackerHit &seed, lcio::TrackerHitImpl* fit_result, double sigma, bool fitSigma = false);

class PhotodotLikelihood : public ROOT::Minuit2::FCNBase
{

  public:

    PhotodotLikelihood(const lcio::TrackerHit &photodot_data);

    virtual double operator()(const std::vector<double>& x) const;

    virtual double Up() const;

  private:

    const lcio::TrackerHit &_photodot_data;
};
}
#endif /* PHOTOELECTRICCLUSTERTOOLS_H_ */