#ifndef __USERFUNC__
#define __USERFUNC__

#include "hparticlecand.h"
#include "hparticlecandsim.h"
#include "hparticledef.h"

#include "hades.h"
#include "hruntimedb.h"
#include "hgeomtransform.h"
#include "hgeomvector.h"
#include "hmdcsizescells.h"
#include "hrichgeometrypar.h"
#include "hrichpadtab.h"
#include "hrichpad.h"
#include "hrichpadcorner.h"
#include "hrichframe.h"
#include "hrichframecorner.h"


#include "TString.h"
#include "TRandom.h"
#include "TStyle.h"
#include "TVector3.h"
#include "TMath.h"
#include "TEveManager.h"
#include "TEveQuadSet.h"
#include "TEveTrans.h"
#include "TEveFrameBox.h"
#include "TEveRGBAPalette.h"
#include "TGLCameraOverlay.h"
#include "TGLViewer.h"

using namespace std;

//---------------------------------------------------------
//             USER FUNC
// This macro provides some helper functions to be
// used in nextEven() function
//---------------------------------------------------------


static Bool_t selectLeptonsBeta(HParticleCand* pcand)
{
    //  selection function for lepton candidates.

    Bool_t test = kFALSE;
    if( pcand->isFlagAND(5,
			 Particle::kIsAcceptedHitRICH,
			 Particle::kIsAcceptedHitInnerMDC,
			 Particle::kIsAcceptedHitOuterMDC,
			 Particle::kIsAcceptedHitMETA,
			 Particle::kIsAcceptedRK
			)
       &&
       pcand->getInnerSegmentChi2() > 0
       &&
       pcand->getChi2()             < 1000      // RK
      ) test = kTRUE;

    if(test)
    {
      //if (fabs(pcand->getDeltaTheta()) < 4. && fabs(pcand->getDeltaPhi()) < 4. ){
      if (pcand->getRichMatchingQuality() < 1.5){
          if(pcand->getBeta() > 0.9) {
	      if(pcand->getMomentum() < 1600) {
		  if(pcand->getMdcdEdx() < 5 ) {
		      if(pcand->getMetaMatchQuality() < 4 ) {
			  return kTRUE;
		      } else return kFALSE;
		  } else return kFALSE;
	      } else return kFALSE;
	  } else return kFALSE;
      }else return kFALSE;
    }
    return kFALSE;
}



#endif // __USERFUNC__