HitDistributor.h

#ifndef HITDISTRIBUTOR_H
#define HITDISTRIBUTOR_H 1

//MARLINTPC
#include "TPCVoxel.h"

//C++
#include <list>
#include <stdexcept>

//LCIO
#include <EVENT/SimTrackerHit.h>

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


namespace marlintpc{

class GlobalPadIndex : public gear::GlobalPadIndex
{
public:
       GlobalPadIndex(int padIndex, int moduleID)
       :gear::GlobalPadIndex(padIndex, moduleID)
       {}

       bool operator == (const  GlobalPadIndex & cmp) const
       {
          return (_moduleID==cmp._moduleID && _padIndex==cmp._padIndex);
       }
   
       inline void setPadIndex(int padindex){_padIndex = padindex;}
       inline void setModuleID(int moduleID){_moduleID = moduleID;}
};


class HitDistributor{

 public:

  HitDistributor(const gear::TPCParameters*, double zBinLength, double inductionFactor=0);

  ~HitDistributor();

  std::list<TPCVoxel> distributeHit(SimTrackerHit* hit, double sigmaTrans, double sigmaLong = 0);

 protected:
 
  std::list<GlobalPadIndex> getPads( double sigmaTrans );

  void getPadsOnTheLeft( std::list<GlobalPadIndex> & padList, GlobalPadIndex centerIndex, double sigmaTrans);
  void getPadsOnTheRight( std::list<GlobalPadIndex> & padList, GlobalPadIndex centerIndex, double sigmaTrans );

  bool checkDistance( GlobalPadIndex globalPadIndex, double sigmaTrans );

  std::list<TPCVoxel> distributeHitOnPadsPolar( std::list<GlobalPadIndex> padIndexList, SimTrackerHit* hit, double sigmaTrans, double sigmaLong );
  std::list<TPCVoxel> distributeHitOnPadsRect(  std::list<GlobalPadIndex> padIndexList, SimTrackerHit* hit, double sigmaTrans, double sigmaLong );

  double CalculateGausIntegral(double sigma, double start, double end, double center);
  void ApplyInduction(double hitTime, double sigmaLong,double charge,GlobalPadIndex gpadIndex, std::list<TPCVoxel> &voxelList);
  void FillZBins(double hitTime, double sigmaLong,double charge,GlobalPadIndex padIndex, std::list<TPCVoxel> &voxelList );

  const gear::TPCParameters* _tpcParameters;

  double _zBinLength;

  double _sigmaTrans,_sigmaLong;

  double _c0Hit, _c1Hit;

  int _padLayoutType;
  int _tpcCoordinateType;
  double _inductionFactor;
};

} //end of namespace marlintpc



#endif // HITDISTRIBUTOR_H