BFieldCutter.h

Go to the documentation of this file.

//*************************************************************************
//*                                                                       *
//*                     class BFieldCutter  GenParticle Selector          *
//*                                                                       *
//*  Select GenParticles above pt threshold to reach barrel calorimeter   *
//*  Threshold is determined by the magnetic field and GenParticle charge *
//*                                                                       *
//*************************************************************************  
#ifndef HEPMCHELPER_BFIELDCUTTER_H
#define HEPMCHELPER_BFIELDCUTTER_H

#ifndef HEPMCHELPER_IMCSELECTOR_H
#include "AtlfastUtils/HepMC_helper/IMCselector.h"
#endif

#ifndef ALTFAST_CHARGESERVICE_H
#include "AtlfastEvent/ChargeService.h"
#endif

#ifndef CLHEP_UNITS_SYSTEMOFUNITS_H
#include "CLHEP/Units/SystemOfUnits.h"
#define CLHEP_UNITS_SYSTEMOFUNITS_H
#endif

namespace HepMC_helper{

  class BFieldCutter: public IMCselector {
  public:
    
    BFieldCutter(bool fieldOn){
      try{
        m_chargeService = new Atlfast::ChargeService();
      }catch(std::string errMsg){
        throw errMsg;
      }catch(...){
        throw;
      }
      m_ptCut = (fieldOn)? 0.5*GeV:0.0*GeV;
    }
   
    IMCselector* create() const;
    bool operator()( const Particle* const )const;
    bool operator()( const Particle& p ) const;
    
  private:
    double m_ptCut;
    Atlfast::ChargeService* m_chargeService;    
  };
}//namespace
#endif

---