DemoMCFilter.h

Go to the documentation of this file.

/***********************************************************************
  $Id: DemoMCFilter.h,v 1.2 2002/02/17 22:50:21 messier Exp $

  A job module to select certain events based on the MC information

  messier@huhepl.harvard.edu
***********************************************************************/
#ifndef DEMOMCFILTER
#define DEMOMCFILTER

#include "JobControl/JobCModule.h"
class REROOT_NeuKin;

class DemoMCFilter : public JobCModule 
{
public:
  DemoMCFilter();

  const Registry& DefaultConfig() const;

  JobCResult      Ana(const MomNavigator* mom);
  void            Config(const Registry& r);
  void            Report();

  // These are obsolete and will eventually go away
  void HandleCommand(JobCommand* cmd);
  void Help();


  // Categories of neutrino interactions modes - selected by flavor
  enum ENeuInt {
    kCCnue      = 0x0001,
    kCCnuebar   = 0x0002,
    kCCnumu     = 0x0004,
    kCCnumubar  = 0x0008,
    kCCnutau    = 0x0010,
    kCCnutaubar = 0x0020,
    kNC         = 0x0040,
    kAllNuInt   = 0xFFFF
  };
  typedef int NuInt_t;

  // Categories of neutrino interaction modes - selected by cross-section type
  enum Emode {
    kElastic  = 0x0001,
    kResonant = 0x0002,
    kDIS      = 0x0003,
    kAllMode  = 0xFFFF
  };
  typedef int Mode_t;
  
  void SetNuInt(NuInt_t i);
  void SetNuMode(Mode_t i);

  void SetEnuMin(double e);
  void SetEnuMax(double e);
  
  void SetYmin(double y);
  void SetYmax(double y);

  void SetPrint(bool doPrint=true);
  
private:
  NuInt_t GetNuInt(const REROOT_NeuKin* nukin);
  Mode_t  GetMode(const REROOT_NeuKin* nukin);
  double  GetEnu(const REROOT_NeuKin* nukin);
  double  GetY(const REROOT_NeuKin* nukin);
  void    Print(const REROOT_NeuKin* nukin);

private:
  double fEnuMin;       // Lower limit on neutrino energies
  double fEnuMax;       // Upper limit on neutrino energies
  double fYmin;         // Minimum y value to allow
  double fYmax;         // Maximum y value to allow
  int    fNuIntSelect;  // Which neutrino interactions to select
  int    fModeSelect;   // Which interaction modes to select
  bool   fDoPrint;      // Print stuff out for each event?
};

#endif // DEMOMCFILTER

---