RHICCommon.hh

// -*- C++ -*-
#ifndef RIVET_RHIC_COMMON_HH
#define RIVET_RHIC_COMMON_HH
 
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/SingleValueProjection.hh"
#include "Rivet/Projections/TriggerProjection.hh"

namespace Rivet {
 
 
  class STAR_BES_Centrality : public SingleValueProjection {
  public:
 
    STAR_BES_Centrality() {
      declare(ChargedFinalState(Cuts::abseta < 0.5 &&
                                Cuts::absrap < 0.1 && Cuts::pT > 0.2 * GeV),
              "STAR_BES_Centrality");
    }
 
    // Destructor
    virtual ~STAR_BES_Centrality() {}

    RIVET_DEFAULT_PROJ_CLONE(STAR_BES_Centrality);

    using Projection::operator =;
 
  protected:
 
    void project(const Event& e) {
      clear();
      double estimate =
        apply<FinalState>(e, "STAR_BES_Centrality").particles().size();
      setValue(estimate);
    }

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

  class BRAHMSCentrality : public SingleValueProjection {
  public:
 
    // Constructor
    BRAHMSCentrality() : SingleValueProjection() {
      // Using here the BRAHMS reaction centrality from eg. 1602.01183, which
      // might not be correct.
      declare(ChargedFinalState(Cuts::pT > 0.1*GeV && Cuts::abseta < 2.2),
              "ChargedFinalState");
    }
 
    // Destructor
    virtual ~BRAHMSCentrality() {}
 
    // Clone on the heap.
    RIVET_DEFAULT_PROJ_CLONE(BRAHMSCentrality);

    using Projection::operator =;
 
  protected:
 
    // Do the projection. Count the number of charged particles in
    // the specified range.
    virtual void project(const Event& e) {
      clear();
      setValue(apply<ChargedFinalState>(e, "ChargedFinalState").particles().size());
    }
 
    // Compare to another projection.
    virtual CmpState compare(const Projection& p) const {
      return mkNamedPCmp(p, "ChargedFinalState");
    }
 
  };
 
 
}
 
#endif