GenMuonRadiationAlgorithm.h

Go to the documentation of this file.

#ifndef TauAnalysis_MCEmbeddingTools_GenMuonRadiationAlgorithm_h
#define TauAnalysis_MCEmbeddingTools_GenMuonRadiationAlgorithm_h

#include "DataFormats/Candidate/interface/Candidate.h"
#include "FWCore/Utilities/interface/StreamID.h"
#include "GeneratorInterface/Pythia6Interface/interface/Pythia6Service.h"
#include "GeneratorInterface/PhotosInterface/interface/PhotosInterfaceBase.h"
#include "GeneratorInterface/PhotosInterface/interface/PhotosFactory.h"

#include<string>

class myPythia6ServiceWithCallback;

class GenMuonRadiationAlgorithm
{
 public:
  explicit GenMuonRadiationAlgorithm(const edm::ParameterSet&);
  ~GenMuonRadiationAlgorithm();

  reco::Candidate::LorentzVector compFSR(const edm::StreamID& streamID, const reco::Candidate::LorentzVector&, int, const reco::Candidate::LorentzVector&, int&);

 private:
  double beamEnergy_;

  enum { kPYTHIA, kPHOTOS };
  int mode_;

  gen::PhotosInterfaceBase* photos_;
  static bool photos_isInitialized_;

  myPythia6ServiceWithCallback* pythia_;
  static bool pythia_isInitialized_;

  int verbosity_;
};

#endif