AtlasCommon.hh

// -*- C++ -*-
#ifndef RIVET_ATLAS_COMMON_HH
#define RIVET_ATLAS_COMMON_HH
 
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/SingleValueProjection.hh"
#include "Rivet/Projections/TriggerProjection.hh"
 
namespace Rivet {
 
  namespace ATLAS {
 

    class SumET_PB_Centrality : public SingleValueProjection {
    public:

      SumET_PB_Centrality() {
        declare(FinalState(Cuts::eta < -3.2 && Cuts::eta > -4.9 && Cuts::pT > 0.1 * GeV),
                "SumET_PB_Centrality");
      }

      RIVET_DEFAULT_PROJ_CLONE(SumET_PB_Centrality);

      using Projection::operator=;
 
    protected:

      void project(const Event& e) {
        clear();
        const FinalState& fsfwd = apply<FinalState>(e, "SumET_PB_Centrality");
        double estimate = 0.0;
        for (const Particle& p : fsfwd.particles()) {
          estimate += p.Et();
        }
        setValue(estimate);
      }

      CmpState compare(const Projection& p) const {
        return mkNamedPCmp(p, "SumET_PB_Centrality");
      }
    };
 

    class SumET_PBPB_Centrality : public SingleValueProjection {
    public:

      SumET_PBPB_Centrality() {
        declare(FinalState(Cuts::abseta > 3.2 && Cuts::abseta < 4.9 && Cuts::pT > 0.1 * GeV),
                "SumET_PBPB_Centrality");
      }

      RIVET_DEFAULT_PROJ_CLONE(SumET_PBPB_Centrality);

      using Projection::operator=;
 
    protected:

      void project(const Event& e) {
        clear();
        const FinalState& fsfwd = apply<FinalState>(e, "SumET_PBPB_Centrality");
        double estimate = 0.0;
        for (const Particle& p : fsfwd.particles()) {
          estimate += p.Et();
        }
        setValue(estimate);
      }

      CmpState compare(const Projection& p) const {
        return mkNamedPCmp(p, "SumET_PBPB_Centrality");
      }
    };
 

    class MinBiasTrigger : public TriggerProjection {
    public:

      MinBiasTrigger() {
        declare(ChargedFinalState(Cuts::eta > 2.09 && Cuts::eta < 3.84 && Cuts::pT > 0.1 * GeV), "MBB");
        declare(ChargedFinalState(Cuts::eta < -2.09 && Cuts::eta > -3.84 && Cuts::pT > 0.1 * GeV), "MBF");
      }

      RIVET_DEFAULT_PROJ_CLONE(MinBiasTrigger);

      using Projection::operator=;
 
    protected:

      void project(const Event& event) {
        pass();
        if (apply<FinalState>(event, "MBF").particles().empty()
            || apply<FinalState>(event, "MBB").particles().empty())
          fail();
      }

      CmpState compare(const Projection& p) const {
        return mkNamedPCmp(p, "MBF") || mkNamedPCmp(p, "MBB");
      }
    };
 
 
  }
}
 
#endif