diff --git a/analyses/pluginALICE/ALICE_2018_I1684320.cc b/analyses/pluginALICE/ALICE_2018_I1684320.cc deleted file mode 100644 --- a/analyses/pluginALICE/ALICE_2018_I1684320.cc +++ /dev/null @@ -1,226 +0,0 @@ -// -*- C++ -*- -#include "Rivet/Analysis.hh" -#include "Rivet/Projections/CentralityProjection.hh" -#include "Rivet/Projections/AliceCommon.hh" -#include "Rivet/Tools/AliceCommon.hh" -#include "Rivet/Projections/PrimaryParticles.hh" -namespace Rivet { - - /// Primary particles + resonances projection. - class PrimaryResonances : public Rivet::PrimaryParticles { - public: - PrimaryResonances(const Cut& c=Cuts::open()) - : Rivet::PrimaryParticles({},c) {} - - virtual int compare(const Projection& p) const { - return UNDEFINED; - } - - virtual std::unique_ptr clone() const { - return std::unique_ptr(new PrimaryResonances(*this)); - } - protected: - bool isPrimaryPID(const HepMC::GenParticle* p) const - { - int pdg = PID::abspid(p->pdg_id()); - if (pdg > 1000000000) return true; - - switch (pdg) { - case Rivet::PID::MUON: - case Rivet::PID::ELECTRON: - case Rivet::PID::GAMMA: - case Rivet::PID::PIPLUS: - case Rivet::PID::K0S: - case Rivet::PID::K0L: - case Rivet::PID::KPLUS: - case 313: // K*0 (892) - case 333: // phi - case Rivet::PID::PROTON: - case Rivet::PID::NEUTRON: - case Rivet::PID::LAMBDA: - case Rivet::PID::SIGMAMINUS: - case Rivet::PID::SIGMAPLUS: - case Rivet::PID::XIMINUS: - case Rivet::PID::XI0: - case Rivet::PID::OMEGAMINUS: - case Rivet::PID::NU_E: - case Rivet::PID::NU_MU: - case Rivet::PID::NU_TAU: - return true; - } - return false; - } - }; - - /// @brief Strangeness enhancement and in pp 7 TeV, no spectra. - class ALICE_2018_I1684320 : public Analysis { - public: - - /// Constructor - DEFAULT_RIVET_ANALYSIS_CTOR(ALICE_2018_I1684320); - - int profileIndex(vector cBins, double c) { - int index = 100; - if (c > 0 && c <= cBins[0]) return cBins.size() - 1; - for (size_t i = 0; i < cBins.size() - 1; ++i) { - if (c > cBins[i] && c <= cBins[i + 1]) { - index = i; - break; - } - } - // Catch low fluctuation. - return max(0, int(cBins.size() - index - 2)); - } - - /// Book histograms and initialise projections before the run - void init() { - // Centrality projection. - declareCentrality(ALICE::V0MMultiplicity(), - "ALICE_2015_PPCentrality","V0M","V0M"); - // Central primary particles - declare(PrimaryResonances(Cuts::abseta < 0.5),"PP"); - declare(PrimaryResonances(Cuts::absrap < 0.5),"PPy"); - centralityBins = {1.,5.,10.,15.,20., 30., 40., 50., 70., 100.}; - centralityBinsPhi = {1.,5.,10.,20., 30., 40., 50., 70., 100.}; - centralityBinsOmega = {5.,15.,30.,50.,100.}; - // Book histograms - - pipT = bookProfile1D(1,1,1); - kpmpT = bookProfile1D(2,1,1); - ppT = bookProfile1D(3,1,1); - kzeropT = bookProfile1D(4,1,1); - lambdapT = bookProfile1D(5,1,1); - kstarpT = bookProfile1D(6,1,1); - phipT = bookProfile1D(7,1,1); - xipT = bookProfile1D(8,1,1); - omegapT = bookProfile1D(9,1,1); - - piYield = bookProfile1D(1,2,1); - kpmYield = bookProfile1D(2,2,1); - kstarYield = bookProfile1D(3,2,1); - phiYield = bookProfile1D(4,2,1); - piRebinnedK = shared_ptr(kstarYield->newclone()); - piRebinnedK->setTitle("piRebinnedK"); - piRebinnedK->setPath("/" + name() + "/piRebinnedK"); - addAnalysisObject(piRebinnedK); - - piRebinnedP = shared_ptr(phiYield->newclone()); - piRebinnedP->setTitle("piRebinnedP"); - piRebinnedP->setPath("/" + name() + "/piRebinnedP"); - addAnalysisObject(piRebinnedP); - } - - - /// Perform the per-event analysis - void analyze(const Event& event) { - if (apply(event,"PP").particles().size() < 1) vetoEvent; - const PrimaryResonances& prim = apply(event,"PPy"); - const double weight = event.weight(); - const CentralityProjection& cent = apply(event,"V0M"); - double c = cent(); - // Fill the pt histograms and count yields. - int npi = 0, nkpm = 0, nkstar = 0, nphi = 0; - // Profile centrality index. - int index = profileIndex(centralityBins,c); - int indexPhi = profileIndex(centralityBinsPhi, c); - int indexOmega = profileIndex(centralityBinsOmega, c); - for (auto p : prim.particles()) { - const double pT = p.pT(); - const int pid = abs(p.pid()); - if (pid == 211) { - ++npi; - pipT->fillBin(index, pT, weight); - } - else if (pid == 321) { - ++nkpm; - kpmpT->fillBin(index, pT, weight); - } - else if (pid == 2212) { - ppT->fillBin(index, pT, weight); - } - else if (pid == 310) { - kzeropT->fillBin(index, pT, weight); - } - else if (pid == 3122) { - lambdapT->fillBin(index, pT, weight); - } - else if (pid == 313) { - ++nkstar; - kstarpT->fillBin(indexPhi, pT, weight); - } - else if (pid == 333) { - ++nphi; - phipT->fillBin(indexPhi, pT, weight); - } - else if (pid == 3312) { - xipT->fillBin(index, pT, weight); - } - else if (pid == 3334) { - omegapT->fillBin(indexOmega, pT, weight); - } - } - // Fill the profiles of yields. - piYield->fillBin(index, double(npi)/2., weight); - kpmYield->fillBin(index, double(nkpm)/2., weight); - kstarYield->fillBin(indexPhi, double(nkstar)/2., weight); - phiYield->fillBin(indexPhi, 2.*double(nphi), weight); - piRebinnedK->fillBin(indexPhi,double(npi)/2.,weight); - piRebinnedP->fillBin(indexPhi,double(npi)/2.,weight); - - } - - - /// Normalise histograms etc., after the run - void finalize() { - // Make the ratios - kpmpi = bookScatter2D(1, 1, 2, true); - kstarpi = bookScatter2D(2, 1, 2, true); - phipi = bookScatter2D(3, 1, 2, true); - - divide(kpmYield, piYield, kpmpi); - divide(kstarYield, piRebinnedK, kstarpi); - divide(phiYield, piRebinnedP, phipi); - } - - //@} - - - /// @name Histograms - //@{ - // Histograms ordered in centrality classes - vector centralityBins; - vector centralityBinsPhi; - vector centralityBinsOmega; - - // Average pT - Profile1DPtr pipT; - Profile1DPtr kpmpT; - Profile1DPtr ppT; - Profile1DPtr kzeropT; - Profile1DPtr lambdapT; - Profile1DPtr kstarpT; - Profile1DPtr phipT; - Profile1DPtr xipT; - Profile1DPtr omegapT; - - // Total yields - Profile1DPtr piYield; - Profile1DPtr kpmYield; - Profile1DPtr kstarYield; - Profile1DPtr phiYield; - Profile1DPtr piRebinnedK; - Profile1DPtr piRebinnedP; - - // Ratios - Scatter2DPtr kpmpi; - Scatter2DPtr kstarpi; - Scatter2DPtr phipi; - //@} - }; - - - // The hook for the plugin system - DECLARE_RIVET_PLUGIN(ALICE_2018_I1684320); - - -} diff --git a/analyses/pluginALICE/ALICE_2018_I1684320.info b/analyses/pluginALICE/ALICE_2018_I1684320.info deleted file mode 100644 --- a/analyses/pluginALICE/ALICE_2018_I1684320.info +++ /dev/null @@ -1,45 +0,0 @@ -Name: ALICE_2018_I1684320 -Year: 2018 -Summary: Mean pT of identified hadrons and phi/phi ratio -- DO NOT PUBLISH! -Experiment: ALICE -Collider: LHC -InspireID: 1684320 -Status: UNVALIDATED -Authors: - - Christian Bierlich -#References: -#- '' -#- '' -#- '' -RunInfo: -NeedCrossSection: no -#Beams: -#Energies: -#Luminosity_fb: -Beams: [p+, p+] -Energies: [7000] -Options: - - cent=REF,GEN,IMP,USR -Reentrant: True -Description: - ' 50\;\GeV$.>' -Keywords: [] -BibKey: Acharya:2018orn -BibTeX: '@article{Acharya:2018orn, - author = "Acharya, Shreyasi and others", - title = "{Multiplicity dependence of light-flavor hadron - production in pp collisions at $\sqrt{s}$ = 7 TeV}", - collaboration = "ALICE", - journal = "Submitted to: Phys. Rev.", - year = "2018", - eprint = "1807.11321", - archivePrefix = "arXiv", - primaryClass = "nucl-ex", - reportNumber = "CERN-EP-2018-209", - SLACcitation = "%%CITATION = ARXIV:1807.11321;%%" -}' -ToDo: - - Implement the analysis, test it, remove this ToDo, and mark as VALIDATED :-) - diff --git a/analyses/pluginALICE/ALICE_2018_I1684320.plot b/analyses/pluginALICE/ALICE_2018_I1684320.plot deleted file mode 100644 --- a/analyses/pluginALICE/ALICE_2018_I1684320.plot +++ /dev/null @@ -1,89 +0,0 @@ -# BEGIN PLOT /ALICE_2018_I1684320/d01-x01-y01 -Title=Average $p_\perp(N_{ch})$ for $\pi^+ + \pi^-$ -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$\langle p_\perp \rangle$ -LogY=0 -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d02-x01-y01 -Title=Average $p_\perp(N_{ch})$ for $K^+ + K^-$ -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$\langle p_\perp \rangle$ -LogY=0 -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d03-x01-y01 -Title=Average $p_\perp(N_{ch})$ for $p + \bar{p}$ -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$\langle p_\perp \rangle$ -LogY=0 -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d04-x01-y01 -Title=Average $p_\perp(N_{ch})$ for $K^0_s$ -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$\langle p_\perp \rangle$ -LogY=0 -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d05-x01-y01 -Title=Average $p_\perp(N_{ch})$ for $\Lambda + \bar{\Lambda}$ -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$\langle p_\perp \rangle$ -LogY=0 -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d06-x01-y01 -Title=Average $p_\perp(N_{ch})$ for $K^*(892) + \bar{K}^*(892)$ -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$\langle p_\perp \rangle$ -LogY=0 -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d07-x01-y01 -Title=Average $p_\perp(N_{ch})$ for $\phi$ -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$\langle p_\perp \rangle$ -LogY=0 -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d08-x01-y01 -Title=Average $p_\perp(N_{ch})$ for $\Xi^- + \bar{\Xi}^+$ -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$\langle p_\perp \rangle$ -LogY=0 -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d09-x01-y01 -Title=Average $p_\perp(N_{ch})$ for $\Omega^- + \bar{\Omega}^+$ -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$\langle p_\perp \rangle$ -LogY=0 -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d01-x02-y01 -Title=$\pi^+ \pi^-$ integrated yield -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$dN/dy$ -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d02-x02-y01 -Title=$K^+ K^-$ integrated yield -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$dN/dy$ -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d03-x02-y01 -Title=$K^*(892) + \bar{K}^*(892)$ integrated yield -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$dN/dy$ -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d04-x02-y01 -Title=$\phi$ integrated yield -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$dN/dy$ -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d01-x01-y02 -Title=$(K^+ + K^-)/(\pi^+ + \pi^-)$ ratio -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$N(K^+ + K^-)/N(\pi^+ + \pi^-)$ -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d02-x01-y02 -Title=$(K^*(892) + \bar{K}^*(892))/(\pi^+ + \pi^-)$ ratio -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$N(K^*(892) + \bar{K}^*(892))/N(\pi^+ + \pi^-)$ -# END PLOT -# BEGIN PLOT /ALICE_2018_I1684320/d03-x01-y02 -Title=$\phi/(\pi^+ + \pi^-)$ ratio -XLabel=$\langle \frac{dN}{d\eta} \rangle_{|\eta| = 0}$ -YLabel=$N(\phi)/N(\pi^+ + \pi^-)$ -# END PLOT diff --git a/analyses/pluginALICE/ALICE_2018_I1684320.yoda b/analyses/pluginALICE/ALICE_2018_I1684320.yoda deleted file mode 100644 --- a/analyses/pluginALICE/ALICE_2018_I1684320.yoda +++ /dev/null @@ -1,313 +0,0 @@ - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d01-x01-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d01-x01-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 3.904e-01 0.09e-01 0.09e-01 -3.896809e+00 1.169045e-01 1.169045e-01 4.237e-01 0.08e-01 0.08e-01 -5.397826e+00 1.619345e-01 1.619345e-01 4.427e-01 0.08e-01 0.08e-01 -6.724329e+00 2.017300e-01 2.017300e-01 4.557e-01 0.09e-01 0.09e-01 -8.447464e+00 2.534240e-01 2.534240e-01 4.693e-01 0.09e-01 0.09e-01 -1.007952e+01 3.023855e-01 3.023855e-01 4.806e-01 0.10e-01 0.10e-01 -1.150694e+01 3.452080e-01 3.452080e-01 4.898e-01 0.10e-01 0.10e-01 -1.345653e+01 4.036955e-01 4.036955e-01 5.008e-01 0.10e-01 0.10e-01 -1.651347e+01 4.954045e-01 4.954045e-01 5.158e-01 0.11e-01 0.11e-01 -2.129430e+01 6.388290e-01 6.388290e-01 5.359e-01 0.12e-01 0.12e-01 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d02-x01-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d02-x01-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 5.885e-01 0.16e-01 0.16e-01 -3.896809e+00 1.169045e-01 1.169045e-01 6.755e-01 0.15e-01 0.15e-01 -5.397826e+00 1.619345e-01 1.619345e-01 7.231e-01 0.13e-01 0.13e-01 -6.724329e+00 2.017300e-01 2.017300e-01 7.538e-01 0.13e-01 0.13e-01 -8.447464e+00 2.534240e-01 2.534240e-01 7.867e-01 0.13e-01 0.13e-01 -1.007952e+01 3.023855e-01 3.023855e-01 8.113e-01 0.13e-01 0.13e-01 -1.150694e+01 3.452080e-01 3.452080e-01 8.317e-01 0.13e-01 0.13e-01 -1.345653e+01 4.036955e-01 4.036955e-01 8.558e-01 0.14e-01 0.14e-01 -1.651347e+01 4.954045e-01 4.954045e-01 8.892e-01 0.13e-01 0.13e-01 -2.129430e+01 6.388290e-01 6.388290e-01 9.332e-01 0.14e-01 0.14e-01 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d03-x01-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d03-x01-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 6.647e-01 0.22e-01 0.22e-01 -3.896809e+00 1.169045e-01 1.169045e-01 7.716e-01 0.22e-01 0.22e-01 -5.397826e+00 1.619345e-01 1.619345e-01 8.285e-01 0.19e-01 0.19e-01 -6.724329e+00 2.017300e-01 2.017300e-01 8.725e-01 0.20e-01 0.20e-01 -8.447464e+00 2.534240e-01 2.534240e-01 9.170e-01 0.20e-01 0.20e-01 -1.007952e+01 3.023855e-01 3.023855e-01 9.539e-01 0.20e-01 0.20e-01 -1.150694e+01 3.452080e-01 3.452080e-01 9.815e-01 0.20e-01 0.20e-01 -1.345653e+01 4.036955e-01 4.036955e-01 1.020e+00 0.21e-01 0.21e-01 -1.651347e+01 4.954045e-01 4.954045e-01 1.071e+00 0.21e-01 0.21e-01 -2.129430e+01 6.388290e-01 6.388290e-01 1.135e+00 0.22e-01 0.22e-01 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d04-x01-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d04-x01-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 6.113e-01 0.07e-01 0.07e-01 -3.896809e+00 1.169045e-01 1.169045e-01 6.869e-01 0.08e-01 0.08e-01 -5.397826e+00 1.619345e-01 1.619345e-01 7.299e-01 0.09e-01 0.09e-01 -6.724329e+00 2.017300e-01 2.017300e-01 7.563e-01 0.09e-01 0.09e-01 -8.447464e+00 2.534240e-01 2.534240e-01 7.900e-01 0.10e-01 0.10e-01 -1.007952e+01 3.023855e-01 3.023855e-01 8.155e-01 0.10e-01 0.10e-01 -1.150694e+01 3.452080e-01 3.452080e-01 8.350e-01 0.10e-01 0.10e-01 -1.345653e+01 4.036955e-01 4.036955e-01 8.612e-01 0.11e-01 0.11e-01 -1.651347e+01 4.954045e-01 4.954045e-01 8.915e-01 0.11e-01 0.11e-01 -2.129430e+01 6.388290e-01 6.388290e-01 9.367e-01 0.12e-01 0.12e-01 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d05-x01-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d05-x01-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 7.669e-01 0.32e-01 0.32e-01 -3.896809e+00 1.169045e-01 1.169045e-01 8.708e-01 0.26e-01 0.26e-01 -5.397826e+00 1.619345e-01 1.619345e-01 9.309e-01 0.22e-01 0.22e-01 -6.724329e+00 2.017300e-01 2.017300e-01 9.759e-01 0.24e-01 0.24e-01 -8.447464e+00 2.534240e-01 2.534240e-01 1.027e+00 0.03e+00 0.03e+00 -1.007952e+01 3.023855e-01 3.023855e-01 1.069e+00 0.03e+00 0.03e+00 -1.150694e+01 3.452080e-01 3.452080e-01 1.103e+00 0.03e+00 0.03e+00 -1.345653e+01 4.036955e-01 4.036955e-01 1.144e+00 0.03e+00 0.03e+00 -1.651347e+01 4.954045e-01 4.954045e-01 1.201e+00 0.03e+00 0.03e+00 -2.129430e+01 6.388290e-01 6.388290e-01 1.277e+00 0.03e+00 0.03e+00 -END YODA_SCATTER2D_V2 - - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d06-x01-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d06-x01-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 7.185e-01 0.245e-01 0.245e-01 -3.896809e+00 1.169045e-01 1.169045e-01 8.656e-01 0.283e-01 0.283e-01 -5.397826e+00 1.619345e-01 1.619345e-01 9.557e-01 0.292e-01 0.292e-01 -6.724329e+00 2.017300e-01 2.017300e-01 1.014e+00 0.031e+00 0.031e+00 -8.447464e+00 2.534240e-01 2.534240e-01 1.064e+00 0.033e+00 0.033e+00 -1.079323e+01 3.237970e-01 3.237970e-01 1.150e+00 0.035e+00 0.035e+00 -1.345653e+01 4.036955e-01 4.036955e-01 1.211e+00 0.037e+00 0.037e+00 -1.651347e+01 4.954045e-01 4.954045e-01 1.295e+00 0.035e+00 0.035e+00 -2.129430e+01 6.388290e-01 6.388290e-01 1.307e+00 0.043e+00 0.043e+00 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d07-x01-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d07-x01-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 8.225e-01 0.374e-01 0.374e-01 -3.896809e+00 1.169045e-01 1.169045e-01 9.910e-01 0.268e-01 0.268e-01 -5.397826e+00 1.619345e-01 1.619345e-01 1.081e+00 0.031e+00 0.031e+00 -6.724329e+00 2.017300e-01 2.017300e-01 1.127e+00 0.033e+00 0.033e+00 -8.447464e+00 2.534240e-01 2.534240e-01 1.201e+00 0.034e+00 0.034e+00 -1.079323e+01 3.237970e-01 3.237970e-01 1.248e+00 0.030e+00 0.030e+00 -1.345653e+01 4.036955e-01 4.036955e-01 1.297e+00 0.036e+00 0.036e+00 -1.651347e+01 4.954045e-01 4.954045e-01 1.360e+00 0.034e+00 0.034e+00 -2.129430e+01 6.388290e-01 6.388290e-01 1.440e+00 0.037e+00 0.037e+00 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d08-x01-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d08-x01-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 8.975e-01 0.634e-01 0.634e-01 -3.896809e+00 1.169045e-01 1.169045e-01 1.012e+00 0.054e+00 0.054e+00 -5.397826e+00 1.619345e-01 1.619345e-01 1.092e+00 0.048e+00 0.048e+00 -6.724329e+00 2.017300e-01 2.017300e-01 1.141e+00 0.048e+00 0.048e+00 -8.447464e+00 2.534240e-01 2.534240e-01 1.188e+00 0.044e+00 0.044e+00 -1.007952e+01 3.023855e-01 3.023855e-01 1.232e+00 0.042e+00 0.042e+00 -1.150694e+01 3.452080e-01 3.452080e-01 1.268e+00 0.037e+00 0.037e+00 -1.345653e+01 4.036955e-01 4.036955e-01 1.293e+00 0.037e+00 0.037e+00 -1.651347e+01 4.954045e-01 4.954045e-01 1.382e+00 0.035e+00 0.035e+00 -2.129430e+01 6.388290e-01 6.388290e-01 1.463e+00 0.036e+00 0.036e+00 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d09-x01-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d09-x01-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.890000e+00 1.400000e-01 1.400000e-01 1.125e+00 0.16e+00 0.16e+00 -6.060000e+00 1.900000e-01 1.900000e-01 1.283e+00 0.11e+00 0.11e+00 -8.990000e+00 2.700000e-01 2.700000e-01 1.465e+00 0.09e+00 0.09e+00 -1.248000e+01 3.700000e-01 3.700000e-01 1.561e+00 0.10e+00 0.10e+00 -1.747000e+01 5.200000e-01 5.200000e-01 1.645e+00 0.08e+00 0.08e+00 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d01-x02-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d01-x02-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 12.100e-01 0.86e-01 0.09e-01 -3.896809e+00 1.169045e-01 1.169045e-01 19.887e-01 1.03e-01 0.08e-01 -5.397826e+00 1.619345e-01 1.619345e-01 2.699e+00 0.13e+00 0.13e+00 -6.724329e+00 2.017300e-01 2.017300e-01 3.326e+00 0.15e+00 0.15e+00 -8.447464e+00 2.534240e-01 2.534240e-01 4.138e+00 0.19e+00 0.19e+00 -1.007952e+01 3.023855e-01 3.023855e-01 4.892e+00 0.23e+00 0.23e+00 -1.150694e+01 3.452080e-01 3.452080e-01 5.554e+00 0.26e+00 0.26e+00 -1.345653e+01 4.036955e-01 4.036955e-01 6.459e+00 0.30e+00 0.30e+00 -1.651347e+01 4.954045e-01 4.954045e-01 7.878e+00 0.38e+00 0.38e+00 -2.129430e+01 6.388290e-01 6.388290e-01 10.035e+00 0.52e+00 0.52e+00 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d02-x02-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d02-x02-y01 -Title: -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 1.402e-01 0.12e-01 0.12e-01 -3.896809e+00 1.169045e-01 1.169045e-01 2.429e-01 0.16e-01 0.16e-01 -5.397826e+00 1.619345e-01 1.619345e-01 3.383e-01 0.21e-01 0.21e-01 -6.724329e+00 2.017300e-01 2.017300e-01 4.231e-01 0.26e-01 0.26e-01 -8.447464e+00 2.534240e-01 2.534240e-01 5.337e-01 0.33e-01 0.33e-01 -1.007952e+01 3.023855e-01 3.023855e-01 6.388e-01 0.40e-01 0.40e-01 -1.150694e+01 3.452080e-01 3.452080e-01 7.308e-01 0.45e-01 0.45e-01 -1.345653e+01 4.036955e-01 4.036955e-01 8.590e-01 0.54e-01 0.54e-01 -1.651347e+01 4.954045e-01 4.954045e-01 1.062e+00 0.07e+00 0.07e+00 -2.129430e+01 6.388290e-01 6.388290e-01 1.374e+00 0.10e+00 0.10e+00 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d03-x02-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d03-x02-y01 -Title: Table 96 -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 5.105e-02 0.83e-02 0.83e-02 -3.896809e+00 1.169045e-01 1.169045e-01 8.088e-02 1.14e-02 1.14e-02 -5.397826e+00 1.619345e-01 1.619345e-01 1.081e-01 0.14e-01 0.14e-01 -6.724329e+00 2.017300e-01 2.017300e-01 1.313e-01 0.18e-01 0.18e-01 -8.447464e+00 2.534240e-01 2.534240e-01 1.568e-01 0.21e-01 0.21e-01 -1.079323e+01 3.237970e-01 3.237970e-01 1.950e-01 0.28e-01 0.28e-01 -1.345653e+01 4.036955e-01 4.036955e-01 2.348e-01 0.35e-01 0.35e-01 -1.651347e+01 4.954045e-01 4.954045e-01 2.730e-01 0.40e-01 0.40e-01 -2.129430e+01 6.388290e-01 6.388290e-01 3.478e-01 0.58e-01 0.58e-01 -END YODA_SCATTER2D_V2 - -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d04-x02-y01 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d04-x02-y01 -Title: Table 96 -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261000e+00 1.170000e-01 1.170000e-01 1.729e-02 0.25e-02 0.25e-02 -3.897000e+00 1.430000e-01 1.430000e-01 2.905e-02 0.50e-02 0.50e-02 -5.398000e+00 1.730000e-01 1.730000e-01 3.988e-02 0.48e-02 0.48e-02 -6.724000e+00 2.070000e-01 2.070000e-01 5.178e-02 0.61e-02 0.61e-02 -8.447000e+00 2.540000e-01 2.540000e-01 6.536e-02 0.75e-02 0.75e-02 -1.079300e+01 3.240000e-01 3.240000e-01 8.711e-02 0.98e-02 0.98e-02 -1.345700e+01 4.040000e-01 4.040000e-01 1.098e-01 0.13e-01 0.13e-01 -1.651300e+01 4.950000e-01 4.950000e-01 1.311e-01 0.15e-01 0.15e-01 -2.129400e+01 6.390000e-01 6.390000e-01 1.697e-01 0.22e-01 0.22e-01 -END YODA_SCATTER2D_V2 - -# kpm / pi -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d01-x01-y02 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d01-x01-y02 -Title: Table 95 -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 1.158913e-01 1.311075e-02 1.311075e-02 -3.896809e+00 1.169045e-01 1.169045e-01 1.221590e-01 1.112001e-02 1.112001e-02 -5.397826e+00 1.619345e-01 1.619345e-01 1.253697e-01 1.053767e-02 1.053767e-02 -6.724329e+00 2.017300e-01 2.017300e-01 1.272244e-01 1.067623e-02 1.067623e-02 -8.447464e+00 2.534240e-01 2.534240e-01 1.289670e-01 1.084226e-02 1.084226e-02 -1.007952e+01 3.023855e-01 3.023855e-01 1.305718e-01 1.104409e-02 1.104409e-02 -1.150694e+01 3.452080e-01 3.452080e-01 1.315863e-01 1.116727e-02 1.116727e-02 -1.345653e+01 4.036955e-01 4.036955e-01 1.329932e-01 1.140330e-02 1.140330e-02 -1.651347e+01 4.954045e-01 4.954045e-01 1.347539e-01 1.160346e-02 1.160346e-02 -2.129430e+01 6.388290e-01 6.388290e-01 1.369054e-01 1.268176e-02 1.268176e-02 -END YODA_SCATTER2D_V2 - -# kstar / pi -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d02-x01-y02 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d02-x01-y02 -Title: Table 96 -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261159e+00 6.783450e-02 6.783450e-02 4.218826e-02 7.642046e-03 7.642046e-03 -3.896809e+00 1.169045e-01 1.169045e-01 4.066866e-02 6.487447e-03 6.487447e-03 -5.397826e+00 1.619345e-01 1.619345e-01 4.005038e-02 5.916134e-03 5.916134e-03 -6.724329e+00 2.017300e-01 2.017300e-01 3.946859e-02 5.947855e-03 5.947855e-03 -8.447464e+00 2.534240e-01 2.534240e-01 3.788867e-02 5.852208e-03 5.852208e-03 -1.079323e+01 3.237970e-01 3.237970e-01 3.745305e-02 5.485571e-03 5.485571e-03 -1.345653e+01 4.036955e-01 4.036955e-01 3.635241e-02 5.814801e-03 5.814801e-03 -1.651347e+01 4.954045e-01 4.954045e-01 3.464710e-02 5.214919e-03 5.214919e-03 -2.129430e+01 6.388290e-01 6.388290e-01 3.466082e-02 5.382984e-03 5.382984e-03 -END YODA_SCATTER2D_V2 - -# phi / pi -BEGIN YODA_SCATTER2D_V2 /REF/ALICE_2018_I1684320/d03-x01-y02 -Variations: [""] -IsRef: 1 -Path: /REF/ALICE_2018_I1684320/d03-x01-y02 -Title: Table 97 -Type: Scatter2D ---- -# xval xerr- xerr+ yval yerr- yerr+ -2.261000e+00 1.170000e-01 1.170000e-01 1.428628e-02 2.233514e-03 2.233514e-03 -3.897000e+00 1.430000e-01 1.430000e-01 1.460563e-02 1.847901e-03 1.847901e-03 -5.398000e+00 1.730000e-01 1.730000e-01 1.477874e-02 1.831193e-03 1.831193e-03 -6.724000e+00 2.070000e-01 2.070000e-01 1.556897e-02 1.915336e-03 1.915336e-03 -8.447000e+00 2.540000e-01 2.540000e-01 1.579411e-02 1.934233e-03 1.934233e-03 -1.079300e+01 3.240000e-01 3.240000e-01 1.673327e-02 1.990169e-03 1.990169e-03 -1.345700e+01 4.040000e-01 4.040000e-01 1.700366e-02 2.109334e-03 2.109334e-03 -1.651300e+01 4.950000e-01 4.950000e-01 1.663973e-02 2.031010e-03 2.031010e-03 -2.129400e+01 6.390000e-01 6.390000e-01 1.691547e-02 2.237534e-03 2.237534e-03 -END YODA_SCATTER2D_V2 diff --git a/analyses/pluginALICE/ALICE_HM_FLOW.cc b/analyses/pluginALICE/ALICE_HM_FLOW.cc deleted file mode 100644 --- a/analyses/pluginALICE/ALICE_HM_FLOW.cc +++ /dev/null @@ -1,271 +0,0 @@ -// -*- C++ -*- -#include "Rivet/Analysis.hh" -#include "Rivet/Projections/ChargedFinalState.hh" -#include "Rivet/Tools/Correlators.hh" -#include "Rivet/Tools/AliceCommon.hh" -#include "Rivet/Projections/AliceCommon.hh" -#include "YODA/Scatter2D.h" - -namespace Rivet { - - - /// @brief Add a short analysis description here - class ALICE_HM_FLOW : public CumulantAnalysis { - public: - - /// Constructor - ALICE_HM_FLOW() : CumulantAnalysis("ALICE_HM_FLOW"){}; - - - /// @name Analysis methods - //@{ - - /// Book histograms and initialise projections before the run - void init() { - // Initialise and register projections - // Declare the trigger projection. - declare(ALICE::V0AndTrigger(),"V0-AND"); - declare(ALICE::V0MMultiplicity(),"V0M"); - - // The full central charged final state. - const ChargedFinalState& cfs = ChargedFinalState(Cuts::abseta < 0.8 && - Cuts::pT > 0.2*GeV && Cuts::pT < 5.0*GeV); - declare(cfs, "CFS"); - - // The positive eta side used for rapidity gap. - const ChargedFinalState& cfsp = ChargedFinalState(Cuts::eta > 0.5 && - Cuts::eta < 0.8 && Cuts::pT > 0.2*GeV && Cuts::pT < 3.0*GeV); - declare(cfsp, "CFSP"); - // ..negative ditto. - const ChargedFinalState& cfsn = ChargedFinalState(Cuts::eta < -0.5 && - Cuts::eta > -0.8 && Cuts::pT > 0.2*GeV && Cuts::pT < 3.0*GeV); - declare(cfsn, "CFSN"); - - // We need a larger gap for v2. - const ChargedFinalState& cfspl = ChargedFinalState(Cuts::eta > 0.7 && - Cuts::eta < 0.8 && Cuts::pT > 0.2*GeV && Cuts::pT < 3.0*GeV); - declare(cfsp, "CFSPL"); - // ..negative ditto. - const ChargedFinalState& cfsnl = ChargedFinalState(Cuts::eta < -0.7 && - Cuts::eta > -0.8 && Cuts::pT > 0.2*GeV && Cuts::pT < 3.0*GeV); - declare(cfsn, "CFSNL"); - - h_v22 = bookScatter2D("v22",120,0,120); - h_v32 = bookScatter2D("v32",120,0,120); - h_v42 = bookScatter2D("v42",120,0,120); - - h_v22gap = bookScatter2D("v22gap",120,0,120); - h_v32gap = bookScatter2D("v32gap",120,0,120); - h_v42gap = bookScatter2D("v42gap",120,0,120); - - vector pTbins = {0.2,0.5,1.0,2.0,3.0,4.0,5.0,7.0,9.0,11.0,15.0}; - - h_v22gappT = bookScatter2D("v22gappT",pTbins,"v22gappT","pT","v2"); - h_v32gappT = bookScatter2D("v32gappT",pTbins,"v32gappT","pT","v3"); - h_v42gappT = bookScatter2D("v42gappT",pTbins,"v42gappT","pT","v4"); - - h_v24 = bookScatter2D("v24",120,0,120); - - h_c22 = bookScatter2D("c22",120,0,120); - h_c32 = bookScatter2D("c32",120,0,120); - h_c42 = bookScatter2D("c42",120,0,120); - - h_c22gap = bookScatter2D("c22gap",120,0,120); - h_c32gap = bookScatter2D("c32gap",120,0,120); - h_c42gap = bookScatter2D("c42gap",120,0,120); - - h_c24 = bookScatter2D("c24",120,0,120); - - h_ec22 = bookScatter2D("ec22",120,0,120); - h_ec24 = bookScatter2D("ec24",120,0,120); - - h_ec22fm = bookScatter2D("ec22fm",120,0,120); - h_ec24fm = bookScatter2D("ec24fm",120,0,120); - - h_c22fm = bookScatter2D("c22fm",120,0,120); - h_c24fm = bookScatter2D("c24fm",120,0,120); - - // Corresponding event averaged correlators. - // Integrated, no gap. - ec22 = bookECorrelator<2,2>("ec22",h_v22); - ec32 = bookECorrelator<3,2>("ec32",h_v32); - ec42 = bookECorrelator<4,2>("ec42",h_v42); - ec24 = bookECorrelator<2,4>("ec24",h_v24); - - ec22fm = bookECorrelator<2,2>("ec22fm",h_c22fm); - ec24fm = bookECorrelator<2,2>("ec24fm",h_c24fm); - - // ... with gap. - ec22gap = bookECorrelatorGap<2,2>("ec22gap",h_v22gap); - ec32gap = bookECorrelatorGap<3,2>("ec32gap",h_v32gap); - ec42gap = bookECorrelatorGap<4,2>("ec42gap",h_v42gap); - // ... pT binned - ec22gappT = bookECorrelatorGap<2,2>("ec22gappT",h_v22gappT); - ec32gappT = bookECorrelatorGap<3,2>("ec32gappT",h_v32gappT); - ec42gappT = bookECorrelatorGap<4,2>("ec42gappT",h_v42gappT); - - - // Maximal N and P. - pair max = getMaxValues(); - // Declare correlator projections. - declare(Correlators(cfs, max.first, max.second),"Correlators"); - declare(Correlators(cfsp, max.first, max.second, pTbins),"CorrelatorsPos"); - declare(Correlators(cfsn, max.first, max.second, pTbins),"CorrelatorsNeg"); - declare(Correlators(cfspl, max.first, max.second),"CorrelatorsPosL"); - declare(Correlators(cfsnl, max.first, max.second),"CorrelatorsNegL"); - } - - - /// Perform the per-event analysis - void analyze(const Event& event) { - // Event trigger. - if (!apply(event, "V0-AND")() ) vetoEvent; - - // High multiplicity trigger. - if (apply(event, "V0M")() < 160.) vetoEvent; - - - const double w = event.weight(); - - // Projections. - const ChargedFinalState& cfs = applyProjection(event,"CFS"); - - // The cumulants projections. - const Correlators& c = applyProjection(event,"Correlators"); - const Correlators& cn = applyProjection(event,"CorrelatorsPos"); - const Correlators& cp = applyProjection(event,"CorrelatorsNeg"); - const Correlators& cnl = applyProjection(event,"CorrelatorsPosL"); - const Correlators& cpl = applyProjection(event,"CorrelatorsNegL"); - - // Number of charged particles. - const double nch = cfs.particles().size(); - - const double nchFwd = apply(event,"V0M")(); - - ec22->fill(nch, c, w); - ec32->fill(nch, c, w); - ec42->fill(nch, c, w); - ec24->fill(nch, c, w); - - ec22fm->fill(nchFwd, c, w); - ec24fm->fill(nchFwd, c, w); - - ec22gap->fill(nch, cpl, cnl, w); - ec32gap->fill(nch, cp, cn, w); - ec42gap->fill(nch, cp, cn, w); - - ec22gappT->fill(cp, cn, w); - ec32gappT->fill(cp, cn, w); - ec42gappT->fill(cp, cn, w); - } - - - - /// Normalise histograms etc., after the run - void finalize() { - - cnTwoInt(h_c22,ec22); - cnTwoInt(h_c32,ec32); - cnTwoInt(h_c42,ec42); - - cnTwoInt(h_c22gap,ec22gap); - cnTwoInt(h_c32gap,ec32gap); - cnTwoInt(h_c42gap,ec42gap); - - cnFourInt(h_c24, ec22, ec24); - - vnTwoInt(h_v22,ec22); - vnTwoInt(h_v32,ec32); - vnTwoInt(h_v42,ec42); - - vnTwoInt(h_v22gap,ec22gap); - vnTwoInt(h_v32gap,ec32gap); - vnTwoInt(h_v42gap,ec42gap); - - vnFourInt(h_v24, ec22, ec24); - - vnTwoDiff(h_v22gappT, ec22gappT); - vnTwoDiff(h_v32gappT, ec32gappT); - vnTwoDiff(h_v42gappT, ec42gappT); - - corrPlot(h_ec22, ec22); - corrPlot(h_ec24, ec24); - - cnTwoInt(h_c22fm, ec22fm); - cnFourInt(h_c24fm, ec22fm, ec24fm); - corrPlot(h_ec22fm, ec22fm); - corrPlot(h_ec24fm, ec24fm); - - - - } - //@} - - - /// @name Histograms - //@{ - Scatter2DPtr h_v22; - Scatter2DPtr h_v32; - Scatter2DPtr h_v42; - - Scatter2DPtr h_v22gap; - Scatter2DPtr h_v32gap; - Scatter2DPtr h_v42gap; - - Scatter2DPtr h_v22gappT; - Scatter2DPtr h_v32gappT; - Scatter2DPtr h_v42gappT; - - Scatter2DPtr h_v24; - - Scatter2DPtr h_c22; - Scatter2DPtr h_c32; - Scatter2DPtr h_c42; - - Scatter2DPtr h_c22fm; - Scatter2DPtr h_c24fm; - - Scatter2DPtr h_c22gap; - Scatter2DPtr h_c32gap; - Scatter2DPtr h_c42gap; - - Scatter2DPtr h_c24; - - // Test histos - Scatter2DPtr h_ec22; - Scatter2DPtr h_ec24; - Scatter2DPtr h_ec22fm; - Scatter2DPtr h_ec24fm; - - // Event correlators. - ECorrPtr ec22; - ECorrPtr ec32; - ECorrPtr ec42; - // ... with gap. - ECorrPtr ec22gap; - ECorrPtr ec32gap; - ECorrPtr ec42gap; - - ECorrPtr ec22gappT; - ECorrPtr ec32gappT; - ECorrPtr ec42gappT; - // Four particle. - ECorrPtr ec24; - - ECorrPtr ec22fm; - ECorrPtr ec24fm; - - - - - //@} - - - }; - - - // The hook for the plugin system - DECLARE_RIVET_PLUGIN(ALICE_HM_FLOW); - - -} diff --git a/analyses/pluginALICE/ALICE_PP_FLOW.cc b/analyses/pluginALICE/ALICE_PP_FLOW.cc deleted file mode 100644 --- a/analyses/pluginALICE/ALICE_PP_FLOW.cc +++ /dev/null @@ -1,267 +0,0 @@ -// -*- C++ -*- -#include "Rivet/Analysis.hh" -#include "Rivet/Projections/ChargedFinalState.hh" -#include "Rivet/Tools/Correlators.hh" -#include "Rivet/Tools/AliceCommon.hh" -#include "Rivet/Projections/AliceCommon.hh" -#include "YODA/Scatter2D.h" - -namespace Rivet { - - - /// @brief Add a short analysis description here - class ALICE_PP_FLOW : public CumulantAnalysis { - public: - - /// Constructor - ALICE_PP_FLOW() : CumulantAnalysis("ALICE_PP_FLOW"){}; - - - /// @name Analysis methods - //@{ - - /// Book histograms and initialise projections before the run - void init() { - // Initialise and register projections - // Declare the trigger projection. - declare(ALICE::V0AndTrigger(),"V0-AND"); - declare(ALICE::V0MMultiplicity(),"V0M"); - - // The full central charged final state. - const ChargedFinalState& cfs = ChargedFinalState(Cuts::abseta < 0.8 && - Cuts::pT > 0.2*GeV && Cuts::pT < 5.0*GeV); - declare(cfs, "CFS"); - - // The positive eta side used for rapidity gap. - const ChargedFinalState& cfsp = ChargedFinalState(Cuts::eta > 0.5 && - Cuts::eta < 0.8 && Cuts::pT > 0.2*GeV && Cuts::pT < 3.0*GeV); - declare(cfsp, "CFSP"); - // ..negative ditto. - const ChargedFinalState& cfsn = ChargedFinalState(Cuts::eta < -0.5 && - Cuts::eta > -0.8 && Cuts::pT > 0.2*GeV && Cuts::pT < 3.0*GeV); - declare(cfsn, "CFSN"); - - // We need a larger gap for v2. - const ChargedFinalState& cfspl = ChargedFinalState(Cuts::eta > 0.7 && - Cuts::eta < 0.8 && Cuts::pT > 0.2*GeV && Cuts::pT < 3.0*GeV); - declare(cfsp, "CFSPL"); - // ..negative ditto. - const ChargedFinalState& cfsnl = ChargedFinalState(Cuts::eta < -0.7 && - Cuts::eta > -0.8 && Cuts::pT > 0.2*GeV && Cuts::pT < 3.0*GeV); - declare(cfsn, "CFSNL"); - - h_v22 = bookScatter2D("v22",120,0,120); - h_v32 = bookScatter2D("v32",120,0,120); - h_v42 = bookScatter2D("v42",120,0,120); - - h_v22gap = bookScatter2D("v22gap",120,0,120); - h_v32gap = bookScatter2D("v32gap",120,0,120); - h_v42gap = bookScatter2D("v42gap",120,0,120); - - vector pTbins = {0.2,0.5,1.0,2.0,3.0,4.0,5.0,7.0,9.0,11.0,15.0}; - - h_v22gappT = bookScatter2D("v22gappT",pTbins,"v22gappT","pT","v2"); - h_v32gappT = bookScatter2D("v32gappT",pTbins,"v32gappT","pT","v3"); - h_v42gappT = bookScatter2D("v42gappT",pTbins,"v42gappT","pT","v4"); - - h_v24 = bookScatter2D("v24",120,0,120); - - h_c22 = bookScatter2D("c22",120,0,120); - h_c32 = bookScatter2D("c32",120,0,120); - h_c42 = bookScatter2D("c42",120,0,120); - - h_c22gap = bookScatter2D("c22gap",120,0,120); - h_c32gap = bookScatter2D("c32gap",120,0,120); - h_c42gap = bookScatter2D("c42gap",120,0,120); - - h_c24 = bookScatter2D("c24",120,0,120); - - h_ec22 = bookScatter2D("ec22",120,0,120); - h_ec24 = bookScatter2D("ec24",120,0,120); - - h_ec22fm = bookScatter2D("ec22fm",120,0,120); - h_ec24fm = bookScatter2D("ec24fm",120,0,120); - - h_c22fm = bookScatter2D("c22fm",120,0,120); - h_c24fm = bookScatter2D("c24fm",120,0,120); - - // Corresponding event averaged correlators. - // Integrated, no gap. - ec22 = bookECorrelator<2,2>("ec22",h_v22); - ec32 = bookECorrelator<3,2>("ec32",h_v32); - ec42 = bookECorrelator<4,2>("ec42",h_v42); - ec24 = bookECorrelator<2,4>("ec24",h_v24); - - ec22fm = bookECorrelator<2,2>("ec22fm",h_c22fm); - ec24fm = bookECorrelator<2,2>("ec24fm",h_c24fm); - - // ... with gap. - ec22gap = bookECorrelatorGap<2,2>("ec22gap",h_v22gap); - ec32gap = bookECorrelatorGap<3,2>("ec32gap",h_v32gap); - ec42gap = bookECorrelatorGap<4,2>("ec42gap",h_v42gap); - // ... pT binned - ec22gappT = bookECorrelatorGap<2,2>("ec22gappT",h_v22gappT); - ec32gappT = bookECorrelatorGap<3,2>("ec32gappT",h_v32gappT); - ec42gappT = bookECorrelatorGap<4,2>("ec42gappT",h_v42gappT); - - - // Maximal N and P. - pair max = getMaxValues(); - // Declare correlator projections. - declare(Correlators(cfs, max.first, max.second),"Correlators"); - declare(Correlators(cfsp, max.first, max.second, pTbins),"CorrelatorsPos"); - declare(Correlators(cfsn, max.first, max.second, pTbins),"CorrelatorsNeg"); - declare(Correlators(cfspl, max.first, max.second),"CorrelatorsPosL"); - declare(Correlators(cfsnl, max.first, max.second),"CorrelatorsNegL"); - } - - - /// Perform the per-event analysis - void analyze(const Event& event) { - // Event trigger. - if (!apply(event, "V0-AND")() ) vetoEvent; - - const double w = event.weight(); - - // Projections. - const ChargedFinalState& cfs = applyProjection(event,"CFS"); - - // The cumulants projections. - const Correlators& c = applyProjection(event,"Correlators"); - const Correlators& cn = applyProjection(event,"CorrelatorsPos"); - const Correlators& cp = applyProjection(event,"CorrelatorsNeg"); - const Correlators& cnl = applyProjection(event,"CorrelatorsPosL"); - const Correlators& cpl = applyProjection(event,"CorrelatorsNegL"); - - // Number of charged particles. - const double nch = cfs.particles().size(); - - const double nchFwd = apply(event,"V0M")(); - - ec22->fill(nch, c, w); - ec32->fill(nch, c, w); - ec42->fill(nch, c, w); - ec24->fill(nch, c, w); - - ec22fm->fill(nchFwd, c, w); - ec24fm->fill(nchFwd, c, w); - - ec22gap->fill(nch, cpl, cnl, w); - ec32gap->fill(nch, cp, cn, w); - ec42gap->fill(nch, cp, cn, w); - - ec22gappT->fill(cp, cn, w); - ec32gappT->fill(cp, cn, w); - ec42gappT->fill(cp, cn, w); - } - - - - /// Normalise histograms etc., after the run - void finalize() { - - cnTwoInt(h_c22,ec22); - cnTwoInt(h_c32,ec32); - cnTwoInt(h_c42,ec42); - - cnTwoInt(h_c22gap,ec22gap); - cnTwoInt(h_c32gap,ec32gap); - cnTwoInt(h_c42gap,ec42gap); - - cnFourInt(h_c24, ec22, ec24); - - vnTwoInt(h_v22,ec22); - vnTwoInt(h_v32,ec32); - vnTwoInt(h_v42,ec42); - - vnTwoInt(h_v22gap,ec22gap); - vnTwoInt(h_v32gap,ec32gap); - vnTwoInt(h_v42gap,ec42gap); - - vnFourInt(h_v24, ec22, ec24); - - vnTwoDiff(h_v22gappT, ec22gappT); - vnTwoDiff(h_v32gappT, ec32gappT); - vnTwoDiff(h_v42gappT, ec42gappT); - - corrPlot(h_ec22, ec22); - corrPlot(h_ec24, ec24); - - cnTwoInt(h_c22fm, ec22fm); - cnFourInt(h_c24fm, ec22fm, ec24fm); - corrPlot(h_ec22fm, ec22fm); - corrPlot(h_ec24fm, ec24fm); - - - - } - //@} - - - /// @name Histograms - //@{ - Scatter2DPtr h_v22; - Scatter2DPtr h_v32; - Scatter2DPtr h_v42; - - Scatter2DPtr h_v22gap; - Scatter2DPtr h_v32gap; - Scatter2DPtr h_v42gap; - - Scatter2DPtr h_v22gappT; - Scatter2DPtr h_v32gappT; - Scatter2DPtr h_v42gappT; - - Scatter2DPtr h_v24; - - Scatter2DPtr h_c22; - Scatter2DPtr h_c32; - Scatter2DPtr h_c42; - - Scatter2DPtr h_c22fm; - Scatter2DPtr h_c24fm; - - Scatter2DPtr h_c22gap; - Scatter2DPtr h_c32gap; - Scatter2DPtr h_c42gap; - - Scatter2DPtr h_c24; - - // Test histos - Scatter2DPtr h_ec22; - Scatter2DPtr h_ec24; - Scatter2DPtr h_ec22fm; - Scatter2DPtr h_ec24fm; - - // Event correlators. - ECorrPtr ec22; - ECorrPtr ec32; - ECorrPtr ec42; - // ... with gap. - ECorrPtr ec22gap; - ECorrPtr ec32gap; - ECorrPtr ec42gap; - - ECorrPtr ec22gappT; - ECorrPtr ec32gappT; - ECorrPtr ec42gappT; - // Four particle. - ECorrPtr ec24; - - ECorrPtr ec22fm; - ECorrPtr ec24fm; - - - - - //@} - - - }; - - - // The hook for the plugin system - DECLARE_RIVET_PLUGIN(ALICE_PP_FLOW); - - -}