PadResponseFunctionGEM.h

#ifndef PRFGEM_H
#define PRFGEM_H 1

//LCIO
#include "lcio.h"
#include <EVENT/TrackerHit.h>
#include <EVENT/Track.h>

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

//ROOT
#include "TGraph.h"
#include "TF1.h"


namespace marlintpc
{
class PadResponseFunctionGEM
{
 public:
  PadResponseFunctionGEM(double defocussing, double diffusion, const gear::TPCParameters& gearparameter);
  bool  EvaluateHit(lcio::TrackerHit *hit, lcio::Track *track=NULL);
  double* GetPRFPosition();
  double GetPRFCharge(){return _charge;}
  double GetPRFChargeError(){return _chargeErr;}
  bool GetPRFSucceeded(){return _fitOK;}
  void SetPRCOnly(int value){_PRConly=value;}


 protected:
  void CalculateFixedWidth(lcio::TrackerHit *hit, lcio::Track *track=NULL);
  void PrepareHitData(lcio::TrackerHit *hit);
  void FitPRF(lcio::TrackerHit *hit, bool fixCharge=false);
  void CalculatePRC(lcio::TrackerHit *hit);
  void KeepCentreOfGravityInfo(lcio::TrackerHit *hit);


  double _defocussing;
  double _diffusion;
  double _fixedWidth;
  double _position;
  double _positionRow;
  double _prfpos[2];
  double _positionRowLocal;
  double _charge;
  double _chargeErr;
  bool _status;
  bool _fitOK;
  TGraph *_graph;
  double _maxCharge;
  double _maxPos;
  double _min,_max;
  const gear::TPCParameters& _tpcParameters;
  TF1 *_PRFShift;
  TF1 *_PRF;
  int _PRConly;

};

}// namespace marlintpc
#endif