#ifdef __CINT__

#else

#include "TROOT.h"
#include "TTree.h"

Int_t simple_check(Float_t *, Int_t);
Int_t calo_check(Float_t *, Int_t);
Int_t fiducial_check(Float_t * , Int_t);

#endif

#define NO_CUT 0
#define SIMPLE_METHOD 1
#define CALO_METHOD 2

Int_t fiducial_check(Float_t * data, Int_t method)
{
  switch (method)
    {
    case SIMPLE_METHOD:
      return simple_check(data);
      break;
    case CALO_METHOD:
      return calo_check(data);
      break;
    default:
      return 1;
      break;
    }

}

Int_t simple_check(Float_t * data)
{
  // Data array is:
  // 0: thx
  // 1: thy

  if ( ( fabs(data[0]) < 0.170 ) 
       && ( fabs(data[1]) < 0.140 ) &&
       ( fabs(data[1]) > 0.040 ) )
    return 1;
  else
    return 0;
}

Int_t calo_check(Float_t * data)
{
  // Data array is:
  // 0: calX
  // 1: calY
  // 2: frontClampX
  // 3: frontClampY
  // 4: rearClampX
  // 5: rearClampY
  // 6: rearClampYOther

  if (  fabs(data[0]) <= 175.0  
	&& fabs(data[1]) >= 30.0  
	&& fabs(data[1]) <= 108.00  
	&& fabs(data[2]) < 31.0
	&& fabs(data[3]) > 7.0
	&& fabs(data[4]) <= 100.0
	&& fabs(data[5]) <= 54.0
	&& fabs(data[6]) <= 54.0 )
    return 1;
  else
    return 0;
}