StMyPointMaker.h

#ifndef __StMyPointMaker_h__
#define __StMyPointMaker_h__

#include "StEEmcGenericPointMaker.h"

class StMyPointMaker : public StEEmcGenericPointMaker
{
public:
  
  StMyPointMaker(const Char_t *name="EEmcPointMaker", const StEEmcA2EMaker *a2e=NULL, const StEEmcGenericClusterMaker *cl=NULL );
  virtual ~StMyPointMaker(){ /* nada */ };

  virtual Int_t Init();
  virtual Int_t Make();
  virtual void  Clear(Option_t *opts="");

  void setSplit(Bool_t s = true){ mAllowSplitting=s; }
  void setSplitMinimumET( Float_t et ){ mSplitMinimumET=et; }
  
  void setSmdMinFraction( Float_t f ){ mSmdMinFrac=f; }

protected:

  Bool_t mAllowSplitting;
  Float_t mSplitMinimumET;

  Float_t mSmdMinFrac; 
  Int_t mMaxClusterId; // maximal cluster ID assigned by a cluster maker. our new clusters will begin from this.

  Bool_t split( const StEEmcSmdCluster &in1, const StEEmcSmdCluster &in2, StEEmcSmdCluster &out1, StEEmcSmdCluster &out2, Float_t &chi2 );

  Bool_t AssociateClusters( const StEEmcSmdClusterVec_t &c1, StEEmcSmdClusterVec_t &c2 );

  Bool_t SplitClusters( StEEmcSmdClusterVec_t &c1, const StEEmcSmdClusterVec_t &c2 );

  Float_t energyChi2( const StEEmcSmdCluster &c1, const StEEmcSmdCluster &c2 ) const;
  Float_t energyChi2( const StEEmcSmdCluster &c1, const StEEmcSmdCluster &c2, const StEEmcSmdCluster &c3 ) const;

  ClassDef(StMyPointMaker,1);

};

#endif