StMCFilter/StExampleFilter.cxx

// @(#)root/eg:$Id: StExampleFilter.cxx,v 1.7 2010/07/29 19:19:41 jwebb Exp $
// Author: Victor Perev  17/03/2009

//______________________________________________________________________________
#include "stdlib.h"
#include "math.h"

#include "StExampleFilter.h"
#include "StGenParticle.h"


static StExampleFilter qwerty;




//______________________________________________________________________________
int StExampleFilter::RejectEG(const StGenParticleMaster &ptl) const
{
  ptl.Print("************** In RejectEG ************** ");
// Condition: number of tracks etaGate[0]<eta<= etaGate[1] must be bigger  etaGate[2]
const static double etaGate[3]={0.8,1.2, 3};
  const StGenParticle *tk=0;
  int n = ptl.Size();
  int ntk=0;
  for (int i=0;i<n;i++) {
    tk = ptl(i); if (!tk)       continue;
    if (tk->GetStatusCode()!=1) continue;
    if (tk->Eta() < etaGate[0]) continue;
    if (tk->Eta() > etaGate[1]) continue;
    ntk++;
  }

  if (ntk<etaGate[2]) return 1;
  return 0;
}
//______________________________________________________________________________
int StExampleFilter::RejectGT(const StGenParticleMaster &ptl) const
{
  ptl.Print("************** In RejectGT ************** ");
  return 0;
}
//______________________________________________________________________________
int StExampleFilter::RejectGE(const StGenParticleMaster &ptl) const
{
  ptl.Print("************** In RejectGE ************** ");
  return 0;
}

---