/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_14/src/GeneratorInterface/HiGenCommon/interface/HadronDecayGenEvtSelector.h

Go to the documentation of this file.

#ifndef _HI_HadronDecayGenEvtSelector_h__
#define _HI_HadronDecayGenEvtSelector_h__

#include <vector>
#include "GeneratorInterface/HiGenCommon/interface/BaseHiGenEvtSelector.h"

class HadronDecayGenEvtSelector : public BaseHiGenEvtSelector 
{
 public:
  HadronDecayGenEvtSelector(const edm::ParameterSet& pset);
  virtual ~HadronDecayGenEvtSelector(){;}
  
  bool filter(HepMC::GenEvent *);
  bool selectParticle(HepMC::GenParticle* par, int status, int pdg /*Absolute*/, double etaMax, double etaMin, double pMin, double ptMax, double ptMin){
    return (par->status() == status &&
            abs(par->pdg_id()) == pdg &&
            par->momentum().eta() < etaMax && par->momentum().eta() > etaMin &&
            par->momentum().rho() > pMin &&
            par->momentum().perp() < ptMax && par->momentum().perp() > ptMin);
  }

 private:

  std::vector<int>    hadronId_;
  std::vector<int>    hadronStatus_;
  std::vector<double> hadronEtaMax_;
  std::vector<double> hadronEtaMin_;
  std::vector<double> hadronPMin_;
  std::vector<double> hadronPtMax_;
  std::vector<double> hadronPtMin_;

  int    decayId_;
  int    decayStatus_;
  double decayEtaMax_;
  double decayEtaMin_;
  double decayPMin_;
  double decayPtMax_;
  double decayPtMin_;
  int    decayNtrig_;

};

#endif