Page MenuHomeHEPForge
Files F8310394

STAR_2006_S6500200.cc
No OneTemporary
Actions

STAR_2006_S6500200.cc
View Options

// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Projections/ChargedFinalState.hh"
#include "Rivet/Projections/IdentifiedFinalState.hh"
#include "Rivet/RivetAIDA.hh"
namespace Rivet {
/// @brief STAR identified hadron spectra in pp at 200 GeV
class STAR_2006_S6500200 : public Analysis {
public:
/// Constructor
STAR_2006_S6500200()
: Analysis("STAR_2006_S6500200"),
_sumWeightSelected(0.0)
{ }
/// Book projections and histograms
void init() {
ChargedFinalState bbc1(-5.0,-3.3, 0.0*GeV); // beam-beam-counter trigger
ChargedFinalState bbc2( 3.3, 5.0, 0.0*GeV); // beam-beam-counter trigger
addProjection(bbc1, "BBC1");
addProjection(bbc2, "BBC2");
IdentifiedFinalState pionfs(-2.5, 2.5, 0.3*GeV);
IdentifiedFinalState protonfs(-2.5, 2.5, 0.4*GeV);
pionfs.acceptIdPair(PIPLUS);
protonfs.acceptIdPair(PROTON);
addProjection(pionfs, "PIONFS");
addProjection(protonfs, "PROTONFS");
_h_pT_piplus = bookHistogram1D(1, 1, 1);
_h_pT_piminus = bookHistogram1D(1, 2, 1);
_h_pT_proton = bookHistogram1D(1, 3, 1);
_h_pT_antiproton = bookHistogram1D(1, 4, 1);
}
/// Do the analysis
void analyze(const Event& event) {
const ChargedFinalState& bbc1 = applyProjection<ChargedFinalState>(event, "BBC1");
const ChargedFinalState& bbc2 = applyProjection<ChargedFinalState>(event, "BBC2");
if (bbc1.size()<1 || bbc2.size()<1) {
MSG_DEBUG("Failed beam-beam-counter trigger");
vetoEvent;
}
const double weight = event.weight();
const IdentifiedFinalState& pionfs = applyProjection<IdentifiedFinalState>(event, "PIONFS");
foreach (const Particle& p, pionfs.particles()) {
if (fabs(p.momentum().rapidity()) < 0.5) {
const double pT = p.momentum().pT() / GeV;
if (p.pdgId()>0) {
_h_pT_piplus->fill(pT, weight/pT);
}
else {
_h_pT_piminus->fill(pT, weight/pT);
}
}
}
const IdentifiedFinalState& protonfs = applyProjection<IdentifiedFinalState>(event, "PROTONFS");
foreach (const Particle& p, protonfs.particles()) {
if (fabs(p.momentum().rapidity()) < 0.5) {
const double pT = p.momentum().pT() / GeV;
if (p.pdgId()>0) {
_h_pT_proton->fill(pT, weight/pT);
}
else {
_h_pT_antiproton->fill(pT, weight/pT);
}
}
}
_sumWeightSelected += event.weight();
}
/// Finalize
void finalize() {
AIDA::IHistogramFactory& hf = histogramFactory();
const string dir = histoDir();
hf.divide(dir + "/d02-x01-y01", *_h_pT_piminus, *_h_pT_piplus);
hf.divide(dir + "/d02-x02-y01", *_h_pT_antiproton, *_h_pT_proton);
hf.divide(dir + "/d02-x03-y01", *_h_pT_proton, *_h_pT_piplus);
hf.divide(dir + "/d02-x04-y01", *_h_pT_antiproton, *_h_pT_piminus);
scale(_h_pT_piplus, 1./(2*M_PI*_sumWeightSelected));
scale(_h_pT_piminus, 1./(2*M_PI*_sumWeightSelected));
scale(_h_pT_proton, 1./(2*M_PI*_sumWeightSelected));
scale(_h_pT_antiproton, 1./(2*M_PI*_sumWeightSelected));
MSG_DEBUG("sumOfWeights() = " << sumOfWeights());
MSG_DEBUG("_sumWeightSelected = " << _sumWeightSelected);
}
private:
double _sumWeightSelected;
AIDA::IHistogram1D * _h_pT_piplus;
AIDA::IHistogram1D * _h_pT_piminus;
AIDA::IHistogram1D * _h_pT_proton;
AIDA::IHistogram1D * _h_pT_antiproton;
};
// The hook for the plugin system
DECLARE_RIVET_PLUGIN(STAR_2006_S6500200);
}

File Metadata

Mime Type
text/x-c
Expires
Sat, Dec 21, 6:21 PM (7 h, 57 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
4023778
Default Alt Text
STAR_2006_S6500200.cc (3 KB)

Event Timeline