diff --git a/analyses/pluginALICE/ALICE_2012_I1126966.cc b/analyses/pluginALICE/ALICE_2012_I1126966.cc --- a/analyses/pluginALICE/ALICE_2012_I1126966.cc +++ b/analyses/pluginALICE/ALICE_2012_I1126966.cc @@ -1,151 +1,151 @@ //-*- C++ -*- #include "Rivet/Projections/CentralityProjection.hh" #include "Rivet/Projections/AliceCommon.hh" #include "Rivet/Tools/AliceCommon.hh" #include "Rivet/Tools/Cuts.hh" namespace Rivet { /// Pion, Kaon, and Proton Production in 0-5% /// central Pb--Pb Collisions at 2.76 TeV class ALICE_2012_I1126966 : public Analysis { public: /// Constructor DEFAULT_RIVET_ANALYSIS_CTOR(ALICE_2012_I1126966); /// Book histograms and initialise projections before the run void init() { // Particles of interest. declare(ALICE::PrimaryParticles(Cuts::absrap < 0.5),"CFS"); // The event trigger. declare(ALICE::V0AndTrigger(), "V0-AND"); // The centrality projection. declareCentrality(ALICE::V0MMultiplicity(), "ALICE_2015_PBPBCentrality", "V0M", "V0M"); // Invariant pT distributions. _histPtPi = bookHisto1D("d01-x01-y01"); //pi+ _histPtPibar = bookHisto1D("d01-x01-y02");// pi- _histPtKaon = bookHisto1D("d02-x01-y01"); //K+ _histPtKaonbar = bookHisto1D("d02-x01-y02"); //K- _histPtProton = bookHisto1D("d03-x01-y01"); //P+ _histPtProtonbar = bookHisto1D("d03-x01-y02"); //P- // Yield histograms. _histNpi = bookHisto1D("d04-x01-y01"); _histNpibar = bookHisto1D("d04-x01-y02"); _histNKaon = bookHisto1D("d04-x01-y03"); _histNKaonbar = bookHisto1D("d04-x01-y04"); _histNproton = bookHisto1D("d04-x01-y05"); _histNprotonbar =bookHisto1D("d04-x01-y06"); // Sum of weights of triggered events. sow = bookCounter("sow"); } /// Perform the per-event analysis void analyze(const Event& event) { // Event weight. const double weight = event.weight(); // Analysis only considers 0-5% central events if (apply(event,"V0M")() > 5.0) vetoEvent; // Event trigger. if (!apply(event, "V0-AND")() ) vetoEvent; sow->fill(weight); // ID particles counters for this event. int Npi=0; int Npibar=0; int NKaon=0; int NKaonbar=0; int Nproton=0; int Nprotonbar=0; for (const auto& p : apply(event,"CFS").particles()) { const double pWeight = weight / p.pT() / 2. / M_PI; switch (p.pid()) { case 211: // pi+ Npi++; _histPtPi->fill(p.pT()/GeV, pWeight); break; case -211: //pi- Npibar++; _histPtPibar->fill(p.pT()/GeV, pWeight); break; case 2212: // proton Nproton++; _histPtProton->fill(p.pT()/GeV, pWeight); break; case -2212: // p-bar Nprotonbar++; _histPtProtonbar->fill(p.pT()/GeV, pWeight); break; case 321: // K+ NKaon++; _histPtKaon->fill(p.pT()/GeV, pWeight); break; case -321: // K- NKaonbar++; _histPtKaonbar->fill(p.pT()/GeV, pWeight); break; } } // Particle loop ends. // Fill yield histograms. - _histNpi->fill(Npi, weight); - _histNpibar->fill(Npibar, weight); - _histNKaon->fill(NKaon, weight); - _histNKaonbar->fill(NKaonbar, weight); - _histNproton->fill(Nproton, weight); - _histNprotonbar->fill(Nprotonbar, weight); + _histNpi->fill(0.0, Npi, weight); + _histNpibar->fill(0.0, Npibar, weight); + _histNKaon->fill(0.0, NKaon, weight); + _histNKaonbar->fill(0.0, NKaonbar, weight); + _histNproton->fill(0.0, Nproton, weight); + _histNprotonbar->fill(0.0, Nprotonbar, weight); } void finalize() { const double s = 1./sow->sumW(); _histPtPi->scaleW(s); _histPtPibar->scaleW(s); _histPtKaon->scaleW(s); _histPtKaonbar->scaleW(s); _histPtProton->scaleW(s); _histPtProtonbar->scaleW(s); _histNpi->scaleW(s); _histNpibar->scaleW(s); _histNKaon->scaleW(s); _histNKaonbar->scaleW(s); _histNproton->scaleW(s); _histNprotonbar->scaleW(s); } private: // pT histograms Histo1DPtr _histPtPi; Histo1DPtr _histPtKaon; Histo1DPtr _histPtProton; Histo1DPtr _histPtPibar; Histo1DPtr _histPtKaonbar; Histo1DPtr _histPtProtonbar; Histo1DPtr _histNpi; Histo1DPtr _histNpibar; Histo1DPtr _histNKaon; Histo1DPtr _histNKaonbar; Histo1DPtr _histNproton; Histo1DPtr _histNprotonbar; CounterPtr sow; }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(ALICE_2012_I1126966); } diff --git a/analyses/pluginALICE/ALICE_2013_I1225979.plot b/analyses/pluginALICE/ALICE_2013_I1225979.plot --- a/analyses/pluginALICE/ALICE_2013_I1225979.plot +++ b/analyses/pluginALICE/ALICE_2013_I1225979.plot @@ -1,8 +1,32 @@ # BEGIN PLOT /ALICE_2013_I1225979/d01-x01-y01 -#Title=[Uncomment and insert title for histogram d01-x01-y01 here] -#XLabel=[Uncomment and insert x-axis label for histogram d01-x01-y01 here] -#YLabel=[Uncomment and insert y-axis label for histogram d01-x01-y01 here] -# + any additional plot settings you might like, see make-plots documentation +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 0-5 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 # END PLOT -# ... add more histograms as you need them ... +# BEGIN PLOT /ALICE_2013_I1225979/d01-x01-y02 +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 5-10 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 +# END PLOT + +# BEGIN PLOT /ALICE_2013_I1225979/d01-x01-y03 +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 10-20 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 +# END PLOT + +# BEGIN PLOT /ALICE_2013_I1225979/d01-x01-y04 +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 20-30 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 +# END PLOT + diff --git a/analyses/pluginALICE/ALICE_2014_I1243865.plot b/analyses/pluginALICE/ALICE_2014_I1243865.plot --- a/analyses/pluginALICE/ALICE_2014_I1243865.plot +++ b/analyses/pluginALICE/ALICE_2014_I1243865.plot @@ -1,134 +1,134 @@ #BEGIN PLOT /ALICE_2014_I1243865/d01-x01-y01 Title=$\Xi^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=0-10$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d01-x01-y02 Title=$\bar{\Xi}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=0-10$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d02-x01-y01 Title=$\Xi^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=10-20$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d02-x01-y02 Title=$\bar{\Xi}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=10-20$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d03-x01-y01 Title=$\Xi^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=20-40$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d03-x01-y02 Title=$\bar{\Xi}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=20-40$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d04-x01-y01 -Title=$\Xi^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=20-40$ \%. +Title=$\Xi^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=40-60$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d04-x01-y02 Title=$\bar{\Xi}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=40-60$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d05-x01-y01 Title=$\Xi^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=60-80$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d05-x01-y02 Title=$\bar{\Xi}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=60-80$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d06-x01-y01 Title=$\Omega^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=0-10$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d06-x01-y02 Title=$\bar{\Omega}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=0-10$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d07-x01-y01 Title=$\Omega^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=10-20$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d07-x01-y02 -Title=$\bar{\Omega}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=1-20$ \%. +Title=$\bar{\Omega}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=10-20$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d08-x01-y01 Title=$\Omega^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=20-40$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d08-x01-y02 Title=$\bar{\Omega}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=20-40$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d09-x01-y01 -Title=$\Omega^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=20-40$ \%. +Title=$\Omega^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=40-60$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d09-x01-y02 Title=$\bar{\Omega}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=40-60$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d10-x01-y01 Title=$\Omega^-$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=60-80$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}} \frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d10-x01-y02 Title=$\bar{\Omega}$ $p_\perp$ in $|\eta| < 0.5$ with centrality $=60-80$ \%. XLabel=$p_{\rm T}$ (GeV/$c$) YLabel=$\frac{1}{N_{\rm inel}}\frac{{\rm d}^2N}{{\rm d}p_{\rm T}{\rm d}\eta}$ ($c$/GeV) #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d14-x01-y01 Title=$\Xi/\pi$ yield ratios as function of $\langle N_{part} \rangle$. XLabel=$\langle N_{part} \rangle$ YLabel=$N(\Xi) / N(\pi)$ LogX=0 LogY=0 #END PLOT #BEGIN PLOT /ALICE_2014_I1243865/d14-x01-y02 Title=$\Omega/\pi$ yield ratios as function of $\langle N_{part} \rangle$. XLabel=$\langle N_{part} \rangle$ YLabel=$N(\Omega) / N(\pi)$ LogX=0 LogY=0 #END PLOT diff --git a/analyses/pluginALICE/ALICE_2016_I1394676.plot b/analyses/pluginALICE/ALICE_2016_I1394676.plot --- a/analyses/pluginALICE/ALICE_2016_I1394676.plot +++ b/analyses/pluginALICE/ALICE_2016_I1394676.plot @@ -1,8 +1,47 @@ # BEGIN PLOT /ALICE_2016_I1394676/d01-x01-y01 -#Title=[Uncomment and insert title for histogram d01-x01-y01 here] -#XLabel=[Uncomment and insert x-axis label for histogram d01-x01-y01 here] -#YLabel=[Uncomment and insert y-axis label for histogram d01-x01-y01 here] -# + any additional plot settings you might like, see make-plots documentation +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 30-40 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 # END PLOT -# ... add more histograms as you need them ... +# BEGIN PLOT /ALICE_2016_I1394676/d01-x01-y02 +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 40-50 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 +# END PLOT + +# BEGIN PLOT /ALICE_2016_I1394676/d01-x01-y03 +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 50-60 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 +# END PLOT + +# BEGIN PLOT /ALICE_2016_I1394676/d01-x01-y04 +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 60-70 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 +# END PLOT + +# BEGIN PLOT /ALICE_2016_I1394676/d01-x01-y05 +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 70-80 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 +# END PLOT + +# BEGIN PLOT /ALICE_2016_I1394676/d01-x01-y06 +Title=d$N_\text{ch}$/d$\eta$, $\sqrt{s}=2.76$ TeV/nn, Centrality: 80-90 \% +XLabel=$\eta$ +YLabel=$\mathrm{d}N_\text{ch} / \mathrm{d}\eta$ +LogY=0 +LogX=0 +# END PLOT