GammaGammaLeptons.hh

// -*- C++ -*-
#ifndef RIVET_GammaGammaLeptons_HH
#define RIVET_GammaGammaLeptons_HH
 
#include "Rivet/Projections/Beam.hh"
#include "Rivet/Projections/PromptFinalState.hh"
#include "Rivet/Projections/HadronicFinalState.hh"
#include "Rivet/Projections/LeptonFinder.hh"
#include "Rivet/Projections/UndressBeamLeptons.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Event.hh"
 
namespace Rivet {
 
 
  class GammaGammaLeptons : public Projection {
  public:
 
 
    GammaGammaLeptons(const Cut& cuts=Cuts::OPEN,
              LeptonReco lreco=LeptonReco::ALL, ObjOrdering lsort=ObjOrdering::ENERGY,
              double beamundresstheta=0.0, double isolDR=0.0, double dressDR=0.0)
      : _isolDR(isolDR), _lsort(lsort)
    {
      setName("GammaGammaLeptons");
      declare(HadronicFinalState(), "IFS");
 
      // Beam undressing
      if (beamundresstheta > 0.0) {
        declare(UndressBeamLeptons(beamundresstheta), "Beam");
      } else {
        declare(Beam(), "Beam");
      }
 
      // Lepton reco mode
      switch (lreco) {
      case LeptonReco::ALL:
        declare(FinalState(cuts), "LFS");
        break;
      case LeptonReco::ALL_DRESSED:
        declare(LeptonFinder(FinalState(), dressDR, cuts), "LFS");
        break;
      case LeptonReco::PROMPT_BARE:
        declare(PromptFinalState(cuts), "LFS");
        break;
      case LeptonReco::PROMPT_DRESSED:
        declare(LeptonFinder(PromptFinalState(), dressDR, cuts), "LFS");
        break;
      }
    }
 
 
    GammaGammaLeptons(Cut& cuts, LeptonReco lreco=LeptonReco::ALL,
              double beamundresstheta=0.0, double isolDR=0.0, double dressDR=0.0)
      : GammaGammaLeptons(cuts, lreco, ObjOrdering::ENERGY, beamundresstheta, isolDR, dressDR)
    {  }
 
    GammaGammaLeptons(LeptonReco lreco, ObjOrdering lsort=ObjOrdering::ENERGY,
          double beamundresstheta=0.0, double isolDR=0.0, double dressDR=0.0)
      : GammaGammaLeptons(Cuts::OPEN, lreco, lsort, beamundresstheta, isolDR, dressDR)
    {  }
 
    GammaGammaLeptons(LeptonReco lreco,
          double beamundresstheta=0.0, double isolDR=0.0, double dressDR=0.0)
      : GammaGammaLeptons(Cuts::OPEN, lreco, ObjOrdering::ENERGY, beamundresstheta, isolDR, dressDR)
    {  }
 
 
    GammaGammaLeptons(const FinalState& leptoncandidates,
              const Beam& beamproj=Beam(),
              const FinalState& isolationfs=FinalState(),
              double isolDR = 0.0, ObjOrdering lsort=ObjOrdering::ENERGY)
      : _isolDR(isolDR), _lsort(lsort)
    {
      declare(leptoncandidates, "LFS");
      declare(isolationfs, "IFS");
      declare(beamproj, "Beam");
    }
 
 
    RIVET_DEFAULT_PROJ_CLONE(GammaGammaLeptons);
 
 
 
    using Projection::operator =;
 
 
  protected:
 
    virtual void project(const Event& e);
 
    virtual CmpState compare(const Projection& p) const;
 
 
  public:
 
    const ParticlePair & in() const { return _incoming; }
 
    const ParticlePair & out() const { return _outgoing; }
 
 
  protected:
 
    ParticlePair _incoming;
 
    ParticlePair _outgoing;
 
    double _isolDR;
 
    ObjOrdering _lsort;
 
  };
 
 
}
 
#endif