ALICE_2010_S8624100.cc

Go to the documentation of this file.

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/RivetAIDA.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/ChargedFinalState.hh"

namespace Rivet {


  class ALICE_2010_S8624100 : public Analysis {
  public:

    /// @name Constructors etc.
    //@{

    /// Constructor
    ALICE_2010_S8624100()
      : Analysis("ALICE_2010_S8624100")
    {
      /// @todo Set whether your finalize method needs the generator cross section
      setNeedsCrossSection(false);
      setBeams(PROTON, PROTON);
    }

    //@}


  public:

    /// @name Analysis methods
    //@{

    /// Book histograms and initialise projections before the run
    void init() {

      ChargedFinalState cfs05(-0.5, 0.5);
      ChargedFinalState cfs10(-1.0, 1.0);
      ChargedFinalState cfs13(-1.3, 1.3);
      addProjection(cfs05, "CFS05");
      addProjection(cfs10, "CFS10");
      addProjection(cfs13, "CFS13");

      if (fuzzyEquals(sqrtS()/GeV, 900, 1E-3)) {
        _h_dN_dNch_05    = bookHistogram1D(11, 1, 1);
        _h_dN_dNch_10    = bookHistogram1D(12, 1, 1);
        _h_dN_dNch_13    = bookHistogram1D(13, 1, 1);
      } else if (fuzzyEquals(sqrtS()/GeV, 2360, 1E-3)) {
        _h_dN_dNch_05    = bookHistogram1D(17, 1, 1);
        _h_dN_dNch_10    = bookHistogram1D(18, 1, 1);
        _h_dN_dNch_13    = bookHistogram1D(19, 1, 1);
      }

    }


    /// Perform the per-event analysis
    void analyze(const Event& event) {
      const double weight = event.weight();

      const ChargedFinalState& charged_05 = applyProjection<ChargedFinalState>(event, "CFS05");
      const ChargedFinalState& charged_10 = applyProjection<ChargedFinalState>(event, "CFS10");
      const ChargedFinalState& charged_13 = applyProjection<ChargedFinalState>(event, "CFS13");
      
      _h_dN_dNch_05->fill(charged_05.size(), weight);
      _h_dN_dNch_10->fill(charged_10.size(), weight);
      _h_dN_dNch_13->fill(charged_13.size(), weight);
    }


    /// Normalise histograms etc., after the run
    void finalize() {

      normalize(_h_dN_dNch_05);
      normalize(_h_dN_dNch_10);
      normalize(_h_dN_dNch_13);

    }

    //@}


  private:

    /// @name Histograms
    //@{

    AIDA::IHistogram1D *_h_dN_dNch_05;
    AIDA::IHistogram1D *_h_dN_dNch_10;
    AIDA::IHistogram1D *_h_dN_dNch_13;
    //@}


  };



  // This global object acts as a hook for the plugin system
  AnalysisBuilder<ALICE_2010_S8624100> plugin_ALICE_2010_S8624100;


}