GammaGammaLeptons.hh

// -*- C++ -*-
#ifndef RIVET_GammaGammaLeptons_HH
#define RIVET_GammaGammaLeptons_HH
 
#include "Rivet/Event.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Projections/Beam.hh"
#include "Rivet/Projections/HadronicFinalState.hh"
#include "Rivet/Projections/LeptonFinder.hh"
#include "Rivet/Projections/PromptFinalState.hh"
#include "Rivet/Projections/UndressBeamLeptons.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