Index: trunk/code/CMS_HIN_12_004.cc =================================================================== --- trunk/code/CMS_HIN_12_004.cc (revision 461) +++ trunk/code/CMS_HIN_12_004.cc (revision 462) @@ -1,109 +1,109 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Tools/Logging.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/ChargedFinalState.hh" #include "Rivet/Projections/FastJets.hh" /// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder... namespace Rivet { class CMS_HIN_12_004 : public Analysis { public: /// @name Constructors etc. //@{ /// Constructor CMS_HIN_12_004() : Analysis("CMS_HIN_12_004") { - setNeedsCrossSection(true); + //setNeedsCrossSection(true); } //@} public: /// @name Analysis methods //@{ /// Book histograms and initialise projections before the run void init() { _centedges += 0.0, 0.05, 0.1, 0.3, 0.5, 0.7, 0.9; _ncentbins = _centedges.size()-1; FinalState fs(-5.0, 5.0, 0.*GeV); addProjection(fs, "FS"); FastJets fj(fs, FastJets::ANTIKT, 0.3); fj.useInvisibles(); addProjection(fj, "jets03"); for (size_t i = 0; i < _ncentbins; ++i) { _sumwtcentbins[i]=0.; stringstream ss; ss << "centbin_" << i; const string pname = ss.str(); _histos[i] = bookHisto1D(pname,refData(1, 1, i+1)); _ratios[i] = bookHisto1D(1, 1, i+1); } } /// Perform the per-event analysis void analyze(const Event& event) { const double weight = event.weight(); const double cent = (event.genEvent()->heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.); Cut cuts = Cuts::abseta < 2.0 && Cuts::pT > 100.*GeV && Cuts::pT < 300.*GeV; const Jets jets = applyProjection(event, "jets03").jetsByPt(cuts); for (size_t i = 0; i < _ncentbins; ++i) { if (cent < 0. || (cent >= _centedges[i] && cent < _centedges[i+1])) { _sumwtcentbins[i]+=weight; // std::cout<fill(jet.momentum().pT(),weight); _ratios[i]->fill(jet.momentum().pT(),weight); } } } } /// Normalise histograms etc., after the run void finalize() { for (size_t i = 0; i < _ncentbins; ++i) { scale(_histos[i],crossSection()/(_sumwtcentbins[i]>0.?_sumwtcentbins[i]:1.)); scale(_ratios[i],crossSection()/(_sumwtcentbins[i]>0.?_sumwtcentbins[i]:1.)); } } //@} private: vector _centedges; double _sumwtcentbins[6]; size_t _ncentbins; private: Histo1DPtr _histos[6]; Histo1DPtr _ratios[6]; }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(CMS_HIN_12_004); } Index: trunk/code/ALICE_GIRTH_4MOMSUB.plot =================================================================== --- trunk/code/ALICE_GIRTH_4MOMSUB.plot (revision 461) +++ trunk/code/ALICE_GIRTH_4MOMSUB.plot (revision 462) @@ -1,7 +1,10 @@ # BEGIN PLOT /ALICE_GIRTH_4MOMSUB/d01-x01-y01 LogY=0 +Title=Pb+Pb $(0-10\%)$ (2.76 TeV) +XLabel=$g$ +YLabel=$(1/N_{\text{jets}}) \ \mathrm{d}\,N/\mathrm{d}\,g$ # END PLOT # BEGIN PLOT /ALICE_GIRTH_4MOMSUB/pp_g_* LogY=0 # END PLOT Index: trunk/code/ATLAS_2013_I1228693.plot =================================================================== --- trunk/code/ATLAS_2013_I1228693.plot (revision 461) +++ trunk/code/ATLAS_2013_I1228693.plot (revision 462) @@ -1,8 +1,84 @@ # BEGIN PLOT /ATLAS_2013_I1228693/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=$|y| < 0.3$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d02-x01-y01 +Title=$0.3 < |y| < 0.8$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d03-x01-y01 +Title=$0.8 < |y| < 1.2$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d04-x01-y01 +Title=$1.2 < |y| < 2.1$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d05-x01-y01 +Title=$2.1 < |y| < 2.8$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d06-x01-y01 +Title=$2.8 < |y| < 3.6$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d07-x01-y01 +Title=$3.6 < |y| < 4.4$, $R=0.4$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d08-x01-y01 +Title=$|y| < 0.3$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d09-x01-y01 +Title=$0.3 < |y| < 0.8$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d10-x01-y01 +Title=$0.8 < |y| < 1.2$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d11-x01-y01 +Title=$1.2 < |y| < 2.1$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d12-x01-y01 +Title=$2.1 < |y| < 2.8$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d13-x01-y01 +Title=$2.8 < |y| < 3.6$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1228693/d14-x01-y01 +Title=$3.6 < |y| < 4.4$, $R=0.6$, $\sqrt(s)=2.76\,\mathrm{TeV}$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$\mathrm{d}^2 \sigma/\mathrm{d} p_\perp\, \mathrm{d}y\ [\mathrm{nb/GeV}]$ # END PLOT -# ... add more histograms as you need them ... Index: trunk/code/hadronspectra.cc =================================================================== --- trunk/code/hadronspectra.cc (revision 461) +++ trunk/code/hadronspectra.cc (revision 462) @@ -1,324 +1,354 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Tools/Logging.hh" #include "Rivet/Projections/ChargedFinalState.hh" #include "Rivet/Projections/IdentifiedFinalState.hh" namespace Rivet { /// @brief hadron spectra in pp class HADRONSPECTRA : public Analysis { public: /// Constructor HADRONSPECTRA() : Analysis("HADRONSPECTRA"), _sumWeightSelected(0.0) { setNeedsCrossSection(true); } /// Book projections and histograms void init() { IdentifiedFinalState pi0fs(-10., 10., 0.*GeV); pi0fs.acceptId(111); addProjection(pi0fs, "PI0FS"); IdentifiedFinalState pifs(-10., 10., 0.*GeV); pifs.acceptId(211); addProjection(pifs, "PIFS"); IdentifiedFinalState kfs(-10., 10., 0.*GeV); kfs.acceptId(130); kfs.acceptId(310); kfs.acceptId(311); kfs.acceptId(321); addProjection(kfs, "KFS"); IdentifiedFinalState pfs(-10., 10., 0.*GeV); pfs.acceptId(2212); addProjection(pfs, "PFS"); ChargedFinalState cfs(-10., 10., 0.*GeV); addProjection(cfs, "CFS"); FinalState fs(-10., 10., 0.*GeV); addProjection(fs, "FS"); vector ptbins; ptbins += 12.0, 14.4,19.2, 24.0, 28.8, 35.2, 41.6, 48.0, 60.8,73.6,86.4,103.6,120.8,140,165,250,400; _h_pT_pi0_035 = bookHisto1D("pi0ptspec", 50, 0.0, 200.0, "pi0 $p_\\perp$ spectrum $|\\eta|<0.35$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_pT_ch_08 = bookHisto1D("chptspec08", 50, 0.0, 200.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<0.8$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_pT_ch_1 = bookHisto1D("chptspec1", 50, 0.0, 200.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<1.$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_pT_ch_2 = bookHisto1D("chptspec2", 50, 0.0, 200.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<2.$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_pT_ch_3 = bookHisto1D("chptspec3", 50, 0.0, 200.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<3.$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_pT_ch_10 = bookHisto1D("chptspec10", 50, 0.0, 200.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<10.$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_E_ch_1 = bookHisto1D("chespec1", 50, 0.0, 1000.0, "charged hadron energy spectrum $|\\eta|<1.$", "$E\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}E$"); _h_E_ch_2 = bookHisto1D("chespec2", 50, 0.0, 1000.0, "charged hadron energy spectrum $|\\eta|<2.$", "$E\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}E$"); _h_E_ch_3 = bookHisto1D("chespec3", 50, 0.0, 1000.0, "charged hadron energy spectrum $|\\eta|<3.$", "$E\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}E$"); _h_E_ch_10 = bookHisto1D("chespec10", 50, 0.0, 1000.0, "charged hadron energy spectrum $|\\eta|<10.$", "$E\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}E$"); _h_pT_ch_1h = bookHisto1D("chptspec1high", 40, 0.0, 1000.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<1.$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_pT_ch_1h2 = bookHisto1D("chptspec1high2", 20, 112.2, 1000.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<1.$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_pT_ch_2h = bookHisto1D("chptspec2high", 40, 0.0, 1000.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<2.$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_pT_ch_3h = bookHisto1D("chptspec3high", 40, 0.0, 1000.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<3.$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_pT_ch_10h = bookHisto1D("chptspec10high", 40, 0.0, 1000.0, "charged hadron $p_\\perp$ spectrum $|\\eta|<10.$", "$p_\\perp\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); _h_E_ch_1h = bookHisto1D("chespec1high", 100, 0.0, 2000.0, "charged hadron energy spectrum $|\\eta|<1.$", "$E\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}E$"); _h_E_ch_2h = bookHisto1D("chespec2high", 100, 0.0, 2000.0, "charged hadron energy spectrum $|\\eta|<2.$", "$E\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}E$"); _h_E_ch_3h = bookHisto1D("chespec3high", 100, 0.0, 2000.0, "charged hadron energy spectrum $|\\eta|<3.$", "$E\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}E$"); _h_E_ch_10h = bookHisto1D("chespec10high", 100, 0.0, 2000.0, "charged hadron energy spectrum $|\\eta|<10.$", "$E\\ [\\mathrm{GeV}]$", "$\\mathrm{d}\\sigma/\\mathrm{d}E$"); + _h_abseta_pt1 = bookHisto1D("abseta_pt1", 30, 0.0, 3.0, + "hadrons $\\p_\\perp > 1\\ \\mathrm{GeV}$", + "$|\\eta|$", + "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$"); + + _h_abseta_pt10 = bookHisto1D("abseta_pt10", 30, 0.0, 3.0, + "hadrons $\\p_\\perp > 10\\ \\mathrm{GeV}$", + "$|\\eta|$", + "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$"); + + _h_abseta_pt50 = bookHisto1D("abseta_pt50", 30, 0.0, 3.0, + "hadrons $\\p_\\perp > 50\\ \\mathrm{GeV}$", + "$|\\eta|$", + "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$"); + + _h_abseta_pt100 = bookHisto1D("abseta_pt100", 30, 0.0, 3.0, + "hadrons $\\p_\\perp > 100\\ \\mathrm{GeV}$", + "$|\\eta|$", + "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$"); + _h_chpT1 = bookHisto1D("chpT1", ptbins); _h_chpT10 = bookHisto1D("chpT10", ptbins); _h_allpT1 = bookHisto1D("allpT1", ptbins); _h_allpT10 = bookHisto1D("allpT10", ptbins); _h_pi0pT1 = bookHisto1D("pi0pT1", ptbins); _h_pi0pT10 = bookHisto1D("pi0pT10", ptbins); _h_pipT1 = bookHisto1D("pipT1", ptbins); _h_pipT10 = bookHisto1D("pipT10", ptbins); _h_KpT1 = bookHisto1D("KpT1", ptbins); _h_KpT10 = bookHisto1D("KpT10", ptbins); _h_ppT1 = bookHisto1D("ppT1", ptbins); _h_ppT10 = bookHisto1D("ppT10", ptbins); } /// Do the analysis void analyze(const Event& event) { const double weight = event.weight(); const IdentifiedFinalState& pi0fs = applyProjection(event, "PI0FS"); foreach (const Particle& p, pi0fs.particles()) { const double pT = p.momentum().pT() / GeV; if (fabs(p.momentum().eta()) < 0.35) _h_pT_pi0_035->fill(pT, weight); if (fabs(p.momentum().eta()) < 1.) _h_pi0pT1->fill(pT, weight); _h_pi0pT10->fill(pT, weight); } const IdentifiedFinalState& pifs = applyProjection(event, "PIFS"); foreach (const Particle& p, pifs.particles()) { const double pT = p.momentum().pT() / GeV; if (fabs(p.momentum().eta()) < 1.) _h_pipT1->fill(pT, weight); _h_pipT10->fill(pT, weight); } const IdentifiedFinalState& kfs = applyProjection(event, "KFS"); foreach (const Particle& p, kfs.particles()) { const double pT = p.momentum().pT() / GeV; if (fabs(p.momentum().eta()) < 1.) _h_KpT1->fill(pT, weight); _h_KpT10->fill(pT, weight); } const IdentifiedFinalState& pfs = applyProjection(event, "PFS"); foreach (const Particle& p, pfs.particles()) { const double pT = p.momentum().pT() / GeV; if (fabs(p.momentum().eta()) < 1.) _h_ppT1->fill(pT, weight); _h_ppT10->fill(pT, weight); } const ChargedFinalState& cfs = applyProjection(event, "CFS"); foreach (const Particle& p, cfs.particles()) { if (fabs(p.momentum().eta()) < 0.8) { const double pT = p.momentum().pT() / GeV; _h_pT_ch_08->fill(pT, weight); } if (fabs(p.momentum().eta()) < 1.) { const double pT = p.momentum().pT() / GeV; const double E = p.momentum().E() / GeV; _h_pT_ch_1->fill(pT, weight); _h_E_ch_1->fill(E, weight); _h_pT_ch_1h->fill(pT, weight); _h_pT_ch_1h2->fill(pT, weight/pT); _h_E_ch_1h->fill(E, weight); _h_chpT1->fill(pT,weight); } if (fabs(p.momentum().eta()) < 2.) { const double pT = p.momentum().pT() / GeV; const double E = p.momentum().E() / GeV; _h_pT_ch_2->fill(pT, weight); _h_E_ch_2->fill(E, weight); _h_pT_ch_2h->fill(pT, weight); _h_E_ch_2h->fill(E, weight); } if (fabs(p.momentum().eta()) < 3.) { const double pT = p.momentum().pT() / GeV; const double E = p.momentum().E() / GeV; _h_pT_ch_3->fill(pT, weight); _h_E_ch_3->fill(E, weight); _h_pT_ch_3h->fill(pT, weight); _h_E_ch_3h->fill(E, weight); } if (fabs(p.momentum().eta()) < 10.) { const double pT = p.momentum().pT() / GeV; const double E = p.momentum().E() / GeV; _h_pT_ch_10->fill(pT, weight); _h_E_ch_10->fill(E, weight); _h_pT_ch_10h->fill(pT, weight); _h_E_ch_10h->fill(E, weight); _h_chpT10->fill(pT,weight); } } const FinalState& fs = applyProjection(event, "FS"); foreach (const Particle& p, fs.particles()) { const double pT = p.momentum().pT() / GeV; - if (fabs(p.momentum().eta()) < 1.) _h_allpT1->fill(pT,weight); - if (fabs(p.momentum().eta()) < 10.) _h_allpT10->fill(pT,weight); + const double eta = p.momentum().abseta(); + if (eta < 1.) _h_allpT1->fill(pT,weight); + if (eta < 10.) _h_allpT10->fill(pT,weight); + if (pT > 1.) _h_abseta_pt1->fill(eta,weight); + if (pT > 10.) _h_abseta_pt10->fill(eta,weight); + if (pT > 50.) _h_abseta_pt50->fill(eta,weight); + if (pT > 100.) _h_abseta_pt100->fill(eta,weight); } _sumWeightSelected += event.weight(); } /// Finalize void finalize() { scale(_h_pT_pi0_035,crossSection()/(_sumWeightSelected)); scale(_h_pT_ch_08,crossSection()/(_sumWeightSelected)); scale(_h_pT_ch_1,crossSection()/(_sumWeightSelected)); scale(_h_pT_ch_2,crossSection()/(_sumWeightSelected)); scale(_h_pT_ch_3,crossSection()/(_sumWeightSelected)); scale(_h_pT_ch_10,crossSection()/(_sumWeightSelected)); scale(_h_E_ch_1,crossSection()/(_sumWeightSelected)); scale(_h_E_ch_2,crossSection()/(_sumWeightSelected)); scale(_h_E_ch_3,crossSection()/(_sumWeightSelected)); scale(_h_E_ch_10,crossSection()/(_sumWeightSelected)); scale(_h_pT_ch_1h,crossSection()/(_sumWeightSelected)); scale(_h_pT_ch_1h2,1./(2.*M_PI*_sumWeightSelected)); scale(_h_pT_ch_2h,crossSection()/(_sumWeightSelected)); scale(_h_pT_ch_3h,crossSection()/(_sumWeightSelected)); scale(_h_pT_ch_10h,crossSection()/(_sumWeightSelected)); scale(_h_E_ch_1h,crossSection()/(_sumWeightSelected)); scale(_h_E_ch_2h,crossSection()/(_sumWeightSelected)); scale(_h_E_ch_3h,crossSection()/(_sumWeightSelected)); scale(_h_E_ch_10h,crossSection()/(_sumWeightSelected)); scale(_h_chpT1,crossSection()/(_sumWeightSelected)); scale(_h_chpT10,crossSection()/(_sumWeightSelected)); scale(_h_allpT1,crossSection()/(_sumWeightSelected)); scale(_h_allpT10,crossSection()/(_sumWeightSelected)); scale(_h_pi0pT1,crossSection()/(_sumWeightSelected)); scale(_h_pi0pT10,crossSection()/(_sumWeightSelected)); scale(_h_pipT1,crossSection()/(_sumWeightSelected)); scale(_h_pipT10,crossSection()/(_sumWeightSelected)); scale(_h_KpT1,crossSection()/(_sumWeightSelected)); scale(_h_KpT10,crossSection()/(_sumWeightSelected)); scale(_h_ppT1,crossSection()/(_sumWeightSelected)); scale(_h_ppT10,crossSection()/(_sumWeightSelected)); + scale(_h_abseta_pt1,crossSection()/(_sumWeightSelected)); + scale(_h_abseta_pt10,crossSection()/(_sumWeightSelected)); + scale(_h_abseta_pt50,crossSection()/(_sumWeightSelected)); + scale(_h_abseta_pt100,crossSection()/(_sumWeightSelected)); getLog() << Log::DEBUG << "sumOfWeights() = " << sumOfWeights() << std::endl; getLog() << Log::DEBUG << "_sumWeightSelected = " << _sumWeightSelected << std::endl; } private: double _sumWeightSelected; Histo1DPtr _h_pT_pi0_035; Histo1DPtr _h_pT_ch_08; Histo1DPtr _h_pT_ch_1; Histo1DPtr _h_pT_ch_2; Histo1DPtr _h_pT_ch_3; Histo1DPtr _h_pT_ch_10; Histo1DPtr _h_E_ch_1; Histo1DPtr _h_E_ch_2; Histo1DPtr _h_E_ch_3; Histo1DPtr _h_E_ch_10; Histo1DPtr _h_pT_ch_1h; Histo1DPtr _h_pT_ch_1h2; Histo1DPtr _h_pT_ch_2h; Histo1DPtr _h_pT_ch_3h; Histo1DPtr _h_pT_ch_10h; Histo1DPtr _h_E_ch_1h; Histo1DPtr _h_E_ch_2h; Histo1DPtr _h_E_ch_3h; Histo1DPtr _h_E_ch_10h; Histo1DPtr _h_chpT1; Histo1DPtr _h_chpT10; Histo1DPtr _h_allpT1; Histo1DPtr _h_allpT10; Histo1DPtr _h_pi0pT1; Histo1DPtr _h_pi0pT10; Histo1DPtr _h_pipT1; Histo1DPtr _h_pipT10; Histo1DPtr _h_KpT1; Histo1DPtr _h_KpT10; Histo1DPtr _h_ppT1; Histo1DPtr _h_ppT10; + Histo1DPtr _h_abseta_pt1, _h_abseta_pt10, _h_abseta_pt50, _h_abseta_pt100; }; // This global object acts as a hook for the plugin system AnalysisBuilder plugin_HADRONSPECTRA; } Index: trunk/code/CMS_2012_I1090064.plot =================================================================== --- trunk/code/CMS_2012_I1090064.plot (revision 461) +++ trunk/code/CMS_2012_I1090064.plot (revision 462) @@ -1,155 +1,215 @@ -# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y0* -#NormalizeToIntegral=1 -Scale=0.10472 -# END PLOT - -# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y0* -#NormalizeToIntegral=1 -Scale=0.06 -# END PLOT - -# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y0* -#NormalizeToIntegral=1 -Scale=0.06 +# BEGIN PLOT /CMS_2012_I1090064/* +FullRange=1 # END PLOT -# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y0* -#NormalizeToIntegral=1 -Scale=0.1 +# BEGIN PLOT /CMS_2012_I1090064/d01-x01-y0* +Scale=0.10472 # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d01-x01-y01 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y01 +Title=PbPb 0-20, $120\ \mathrm{GeV} < p_{\perp,1} < 150\ \mathrm{GeV}$ +XLabel=$\Delta \phi$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d01-x01-y02 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y02 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y02 +Title=PbPb 0-20, $150\ \mathrm{GeV} < p_{\perp,1} < 180\ \mathrm{GeV}$ +XLabel=$\Delta \phi$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d01-x01-y03 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y03 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y03 +Title=PbPb 0-20, $180\ \mathrm{GeV} < p_{\perp,1} < 220\ \mathrm{GeV}$ +XLabel=$\Delta \phi$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d01-x01-y04 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y04 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y04 +Title=PbPb 0-20, $220\ \mathrm{GeV} < p_{\perp,1} < 260\ \mathrm{GeV}$ +XLabel=$\Delta \phi$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d01-x01-y05 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y05 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y05 -# END PLOT +#Title=PbPb 0-20, $260\ \mathrm{GeV} < p_{\perp,1} < 300\ \mathrm{GeV}$ +XLabel=$\Delta \phi$ +YLabel=event fraction + END PLOT # BEGIN PLOT /CMS_2012_I1090064/d01-x01-y06 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d01-x01-y06 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d01-x01-y06 +Title=PbPb 0-20, $300\ \mathrm{GeV} < p_{\perp,1} < 500\ \mathrm{GeV}$ +XLabel=$\Delta \phi$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d02-x01-y01 LogY=0 XMax=220. YMin=0.7 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d02-x01-y01 /home/jewel/trunk/data/jewel-2.0.0_vac_fixed_5-350.aida/CMS_2012_I1090064/d02-x01-y01 +Title=di-jet fraction in pp +XLabel=$p_{\perp,1}$ [GeV] +YLabel=$N_\mathrm{dijet}/N_\mathrm{leading jet}$ # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d02-x01-y02 LogY=0 XMax=370. YMin=0.7 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d02-x01-y02 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d02-x01-y02 +Title=di-jet fraction in PbPb 0-20 +XLabel=$p_{\perp,1}$ [GeV] +YLabel=$N_\mathrm{dijet}/N_\mathrm{leading jet}$ # END PLOT -# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y01 +# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y0* +Scale=0.06 LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y01 +# END PLOT + +# BEGIN PLOT /CMS_2012_I1090064/d03-x01-y01 +Title=di-jet asymmetry in PbPb 0-20, $120\ \mathrm{GeV} < p_{\perp,1} < 150\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d03-x01-y02 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y02 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y02 +Title=di-jet asymmetry in PbPb 0-20, $150\ \mathrm{GeV} < p_{\perp,1} < 180\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d03-x01-y03 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y03 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y03 +Title=di-jet asymmetry in PbPb 0-20, $180\ \mathrm{GeV} < p_{\perp,1} < 220\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d03-x01-y04 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y04 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y04 +Title=di-jet asymmetry in PbPb 0-20, $220\ \mathrm{GeV} < p_{\perp,1} < 260\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d03-x01-y05 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y05 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y05 +Title=di-jet asymmetry in PbPb 0-20, $260\ \mathrm{GeV} < p_{\perp,1} < 300\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d03-x01-y06 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d03-x01-y06 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d03-x01-y06 +Title=di-jet asymmetry in PbPb 0-20, $300\ \mathrm{GeV} < p_{\perp,1} < 500\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d04-x01-y0* +Scale=0.06 LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d04-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-10_wo_1000-1199.aida/CMS_2012_I1090064/d04-x01-y01 # END PLOT -# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y01 +# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y01 +Title=di-jet asymmetry in PbPb 0-10, $p_{\perp,1} > 120\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction +# END PLOT + +# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y02 +Title=di-jet asymmetry in PbPb 10-20, $p_{\perp,1} > 120\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction +# END PLOT + +# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y03 +Title=di-jet asymmetry in PbPb 20-30, $p_{\perp,1} > 120\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction +# END PLOT + +# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y04 +Title=di-jet asymmetry in PbPb 30-50, $p_{\perp,1} > 120\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction +# END PLOT + +# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y05 +Title=di-jet asymmetry in PbPb 50-70, $p_{\perp,1} > 120\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction +# END PLOT + +# BEGIN PLOT /CMS_2012_I1090064/d04-x01-y06 +Title=di-jet asymmetry in PbPb 70-100, $p_{\perp,1} > 120\ \mathrm{GeV}$ +XLabel=$A_J$ +YLabel=event fraction +# END PLOT + +# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y0* +Scale=0.1 LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y01 +# END PLOT + +# BEGIN PLOT /CMS_2012_I1090064/d05-x01-y01 +Title=di-jet $p_\perp$ ratio in PbPb 0-20, $120\ \mathrm{GeV} < p_{\perp,1} < 150\ \mathrm{GeV}$ +XLabel=$p_{\perp,2}/p_{\perp,1}$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d05-x01-y02 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y02 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y02 +Title=di-jet $p_\perp$ ratio in PbPb 0-20, $150\ \mathrm{GeV} < p_{\perp,1} < 180\ \mathrm{GeV}$ +XLabel=$p_{\perp,2}/p_{\perp,1}$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d05-x01-y03 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y03 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y03 +Title=di-jet $p_\perp$ ratio in PbPb 0-20, $180\ \mathrm{GeV} < p_{\perp,1} < 220\ \mathrm{GeV}$ +XLabel=$p_{\perp,2}/p_{\perp,1}$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d05-x01-y04 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y04 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y04 +Title=di-jet $p_\perp$ ratio in PbPb 0-20, $220\ \mathrm{GeV} < p_{\perp,1} < 260\ \mathrm{GeV}$ +XLabel=$p_{\perp,2}/p_{\perp,1}$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d05-x01-y05 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y05 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y05 +Title=di-jet $p_\perp$ ratio in PbPb 0-20, $260\ \mathrm{GeV} < p_{\perp,1} < 300\ \mathrm{GeV}$ +XLabel=$p_{\perp,2}/p_{\perp,1}$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d05-x01-y06 -LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d05-x01-y06 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-20_wo_1000-1199.aida/CMS_2012_I1090064/d05-x01-y06 +Title=di-jet $p_\perp$ ratio in PbPb 0-20, $300\ \mathrm{GeV} < p_{\perp,1} < 500\ \mathrm{GeV}$ +XLabel=$p_{\perp,2}/p_{\perp,1}$ +YLabel=event fraction # END PLOT # BEGIN PLOT /CMS_2012_I1090064/d06-x01-y0* YMin=0.5 XMax=300. LogY=0 -FullRange=1 -#DrawOnly=/home/jewel/trunk/CMS_2012_I1090064.aida/REF/CMS_2012_I1090064/d06-x01-y01 /home/jewel/trunk/data/jewel-2.0.0_vac_fixed_5-350.aida/CMS_2012_I1090064/d06-x01-y01 # END PLOT +# BEGIN PLOT /CMS_2012_I1090064/d06-x01-y01 +Title=mean di-jet $p_\perp$ ratio in pp +XLabel=$p_{\perp,1}$ [GeV] +YLabel=$\langle p_{\perp,2}/p_{\perp,1}\rangle$ +# END PLOT +# BEGIN PLOT /CMS_2012_I1090064/d06-x01-y02 +Title=mean di-jet $p_\perp$ ratio in PbPb 0-20 +XLabel=$p_{\perp,1}$ [GeV] +YLabel=$\langle p_{\perp,2}/p_{\perp,1}\rangle$ +# END PLOT + +# BEGIN PLOT /CMS_2012_I1090064/d06-x01-y03 +Title=mean di-jet $p_\perp$ ratio in PbPb 20-50 +XLabel=$p_{\perp,1}$ [GeV] +YLabel=$\langle p_{\perp,2}/p_{\perp,1}\rangle$ +# END PLOT +# BEGIN PLOT /CMS_2012_I1090064/d06-x01-y04 +Title=mean di-jet $p_\perp$ ratio in PbPb 50-100 +XLabel=$p_{\perp,1}$ [GeV] +YLabel=$\langle p_{\perp,2}/p_{\perp,1}\rangle$ +# END PLOT Index: trunk/code/ALICE_2015_I1343112.cc =================================================================== --- trunk/code/ALICE_2015_I1343112.cc (revision 461) +++ trunk/code/ALICE_2015_I1343112.cc (revision 462) @@ -1,181 +1,181 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Tools/Logging.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/FastJets.hh" #include "Rivet/Tools/ParticleIdUtils.hh" #include #include "Rivet/Tools/Logging.hh" /// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder... namespace Rivet { class ALICE_2015_I1343112 : public Analysis { public: /// @name Constructors etc. //@{ /// Constructor ALICE_2015_I1343112() : Analysis("ALICE_2015_I1343112"), _hadronetamax(1.), _jetradius(0.2), _jetetamax(0.5), _jetptmin(20.), _jetptmax(120.) { } //@} public: /// @name Analysis methods //@{ /// Book histograms and initialise projections before the run void init() { _centedges += 0.0, 0.1, 0.3; _biases += 0., 3., 5., 7., 10.; _ncentbins = 2; _nbiasbins = 5; FinalState fs(-_hadronetamax, _hadronetamax, 0.*GeV); addProjection(fs, "FS"); FastJets fj(fs, FastJets::ANTIKT, _jetradius); fj.useInvisibles(); addProjection(fj, "Jets"); _sumwtcentbins[0] = 0.; _sumwtcentbins[1] = 0.; _sumwtcentbins[2] = 0.; _h_jetpt[0] = bookHisto1D(1, 1, 1); _h_jetpt[1] = bookHisto1D(1, 1, 2); _h_jetpt[2] = bookHisto1D(1, 1, 3); _h_jetpt_pp_incl = bookHisto1D("jetpt_pp_incl",refData(2, 1, 1)); _h_jetpt_pp_bias = bookHisto1D("jetpt_pp_bias",refData(2, 1, 1)); _s_jetpt_pp_ratio = bookScatter2D(2, 1, 1); _h_jetpt_AA_bias[0] = bookHisto1D("jetpt_AA_bias_0",refData(3, 1, 1)); _h_jetpt_AA_bias[1] = bookHisto1D("jetpt_AA_bias_3",refData(3, 1, 2)); _h_jetpt_AA_bias[2] = bookHisto1D("jetpt_AA_bias_5",refData(3, 1, 1)); _h_jetpt_AA_bias[3] = bookHisto1D("jetpt_AA_bias_7",refData(3, 1, 3)); _h_jetpt_AA_bias[4] = bookHisto1D("jetpt_AA_bias_10",refData(3, 1, 4)); _s_jetpt_AA_ratios[0] = bookScatter2D(3, 1, 1); _s_jetpt_AA_ratios[1] = bookScatter2D(3, 1, 2); _s_jetpt_AA_ratios[2] = bookScatter2D(3, 1, 3); _s_jetpt_AA_ratios[3] = bookScatter2D(3, 1, 4); _h_raa[0] = bookHisto1D(4, 1, 1); _h_raa[1] = bookHisto1D(4, 1, 2); } /// Perform the per-event analysis void analyze(const Event& event) { const double weight = event.weight(); const double cent = (event.genEvent()->heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.); Cut cuts = Cuts::abseta < _jetetamax && Cuts::pT > _jetptmin*GeV && Cuts::pT < _jetptmax*GeV; const Jets jets = applyProjection(event, "Jets").jetsByPt(cuts); if (cent < 0.) _sumwtcentbins[0] += weight; if (cent < 0. || (cent >= _centedges[0] && cent < _centedges[1])) _sumwtcentbins[1] += weight; if (cent < 0. || (cent >= _centedges[1] && cent < _centedges[2])) _sumwtcentbins[2] += weight; for (size_t k = 0; k < _nbiasbins; ++k) { foreach (Jet jet, jets) { // track bias foreach (Particle part, jet.particles()) { if (PID::charge(part.pdgId()) != 0. && part.momentum().pT() > _biases[k]*GeV) { double jetpt(jet.momentum().pT()); if (cent < 0.) { if (k == 0) _h_jetpt_pp_incl->fill(jetpt,weight); if (k == 2) { _h_jetpt[0]->fill(jetpt,weight); _h_jetpt_pp_bias->fill(jetpt,weight); _h_raa[0]->fill(jetpt,weight); _h_raa[1]->fill(jetpt,weight); } } if (cent >= _centedges[0] && cent < _centedges[1]) { _h_jetpt_AA_bias[k]->fill(jetpt,weight); if (k == 2) { _h_jetpt[1]->fill(jetpt,weight); _h_raa[0]->fill(jetpt,weight); } } if (cent >= _centedges[0] && cent < _centedges[2]) { if (k == 2) { _h_jetpt[2]->fill(jetpt,weight); _h_raa[1]->fill(jetpt,weight); } } break; } } } } } /// Normalise histograms etc., after the run void finalize() { scale(_h_jetpt[0],1./(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.)); scale(_h_jetpt[1],1./(_sumwtcentbins[1]>0.?_sumwtcentbins[1]:1.)); scale(_h_jetpt[2],1./(_sumwtcentbins[2]>0.?_sumwtcentbins[2]:1.)); scale(_h_raa[0],crossSection()/(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.)); scale(_h_raa[1],crossSection()/(_sumwtcentbins[1]>0.?_sumwtcentbins[1]:1.)); scale(_h_jetpt_pp_incl,crossSection()/(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.)); scale(_h_jetpt_pp_bias,crossSection()/(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.)); for (size_t k = 0; k < _nbiasbins; ++k) { scale(_h_jetpt_AA_bias[k],crossSection()/(_sumwtcentbins[0]>0.?_sumwtcentbins[0]:1.)); } - divide(_h_jetpt_pp_bias, _h_jetpt_pp_bias, _s_jetpt_pp_ratio); + divide(_h_jetpt_pp_bias, _h_jetpt_pp_incl, _s_jetpt_pp_ratio); divide(_h_jetpt_AA_bias[0], _h_jetpt_AA_bias[2], _s_jetpt_AA_ratios[0]); divide(_h_jetpt_AA_bias[1], _h_jetpt_AA_bias[2], _s_jetpt_AA_ratios[1]); divide(_h_jetpt_AA_bias[3], _h_jetpt_AA_bias[2], _s_jetpt_AA_ratios[2]); divide(_h_jetpt_AA_bias[4], _h_jetpt_AA_bias[2], _s_jetpt_AA_ratios[3]); getLog() << Log::DEBUG << "sumOfWeights() = " << sumOfWeights() << std::endl; } //@} private: double _hadronetamax, _jetradius, _jetetamax, _jetptmin, _jetptmax; double _sumwtcentbins[3]; vector _centedges, _biases; size_t _ncentbins, _nbiasbins; private: /// @name Histograms //@{ Histo1DPtr _h_jetpt[3]; Histo1DPtr _h_jetpt_pp_incl; Histo1DPtr _h_jetpt_pp_bias; Histo1DPtr _h_jetpt_AA_bias[5]; Histo1DPtr _h_raa[2]; Scatter2DPtr _s_jetpt_pp_ratio; Scatter2DPtr _s_jetpt_AA_ratios[4]; //@} }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(ALICE_2015_I1343112); } Index: trunk/code/ALICE_2014_I1263194.cc =================================================================== --- trunk/code/ALICE_2014_I1263194.cc (revision 461) +++ trunk/code/ALICE_2014_I1263194.cc (revision 462) @@ -1,156 +1,192 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Tools/Logging.hh" #include "Rivet/Projections/ChargedFinalState.hh" #include "Rivet/Projections/FastJets.hh" /// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder... namespace Rivet { class ALICE_2014_I1263194 : public Analysis { public: /// @name Constructors etc. //@{ /// Constructor ALICE_2014_I1263194() : Analysis("ALICE_2014_I1263194") { } //@} public: /// @name Analysis methods //@{ /// Book histograms and initialise projections before the run void init() { _centedges += 0.0, 0.1, 0.3, 0.5, 0.8; _jetparams += 0.2, 0.3; _biases += 0., 5., 10.; _ncentbins = 4; _njetbins = 2; _nbiasbins = 3; ChargedFinalState cfs(-0.9, 0.9, 0.15*GeV); addProjection(cfs, "CFS"); for (size_t j = 0; j < _njetbins; ++j) { stringstream ss; ss << "Jets" << j; const string pname = ss.str(); _names[j] = pname; FastJets fj(cfs, FastJets::ANTIKT, _jetparams[j]); addProjection(fj, pname); } for (size_t i = 0; i < _ncentbins; ++i) { _sumwtcentbins[i] = 0.; if (i < _ncentbins-1) _rsumwtcentbins[i] = 0.; for (size_t j = 0; j < _njetbins; ++j) { for (size_t k = 0; k < _nbiasbins; ++k) { _h_jetspec[i][j][k] = bookHisto1D(2+4*k+i, 1, j+1); if (i < _ncentbins-1) { _h_rcp[i][j][k] = bookHisto1D(18+3*k+i, 1, j+1); } } } - } + } + + for (size_t i = 0; i < 2; i++) { + for (size_t j = 0; j < _njetbins; ++j) { + for (size_t k = 0; k < _nbiasbins-1; ++k) { + _s_rbias[i][j][k] = bookScatter2D(14+2*k+i, 1, j+1); + stringstream ss1; ss1 << "num_"<heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.); for (size_t i = 0; i < _ncentbins; ++i) { if (cent < 0. || (cent >= _centedges[i] && cent < _centedges[i+1])) { _sumwtcentbins[i] += weight; if (i < _ncentbins-1) _rsumwtcentbins[i] += weight; else { for (size_t l = 0; l < _ncentbins-1; ++l) { _rsumwtcentbins[l] += weight; } } } } for (size_t j = 0; j < _njetbins; ++j) { Cut cuts = Cuts::abseta < 0.5 && Cuts::pT > 10*GeV && Cuts::pT < 120.*GeV; const Jets jets = applyProjection(event, _names[j]).jetsByPt(cuts); for (size_t i = 0; i < _ncentbins; ++i) { if (cent < 0. || (cent >= _centedges[i] && cent < _centedges[i+1])) { // track bias for (size_t k = 0; k < _nbiasbins; ++k) { foreach (Jet jet, jets) { vector< Particle > & particles(jet.particles()); foreach (Particle part, particles) { if (part.momentum().pT() > _biases[k]*GeV) { _h_jetspec[i][j][k]->fill(jet.momentum().pT(),weight); if (i < _ncentbins-1) _h_rcp[i][j][k]->fill(jet.momentum().pT(),weight); else { for (size_t l = 0; l < _ncentbins-1; ++l) { _h_rcp[l][j][k]->fill(jet.momentum().pT(),weight); } - } + } + if (i == 0 && k == 0) _h_rbiasnum[0][j][0]->fill(jet.momentum().pT(),weight); + if (i == _ncentbins-1 && k == 0) _h_rbiasnum[1][j][0]->fill(jet.momentum().pT(),weight); + if (i == 0 && k == 1) { + _h_rbiasden[0][j][0]->fill(jet.momentum().pT(),weight); + _h_rbiasden[0][j][1]->fill(jet.momentum().pT(),weight); + } + if (i == _ncentbins-1 && k == 0) { + _h_rbiasden[1][j][0]->fill(jet.momentum().pT(),weight); + _h_rbiasden[1][j][1]->fill(jet.momentum().pT(),weight); + } + if (i == 0 && k == 2) _h_rbiasnum[0][j][1]->fill(jet.momentum().pT(),weight); + if (i == _ncentbins-1 && k == 2) _h_rbiasnum[1][j][1]->fill(jet.momentum().pT(),weight); break; } } } } } } } } /// Normalise histograms etc., after the run void finalize() { for (size_t i = 0; i < _ncentbins; ++i) { for (size_t j = 0; j < _njetbins; ++j) { for (size_t k = 0; k < _nbiasbins; ++k) { scale(_h_jetspec[i][j][k],1.e6/(_sumwtcentbins[i]>0.?_sumwtcentbins[i]:1.)); if (i < _ncentbins-1) scale(_h_rcp[i][j][k],crossSection()/(_rsumwtcentbins[i]>0.?_rsumwtcentbins[i]:1.)); } } - } - } + } + + for (size_t i = 0; i < 2; ++i) { + for (size_t j = 0; j < _njetbins; ++j) { + for (size_t k = 0; k < _nbiasbins-1; ++k) { + divide(_h_rbiasnum[i][j][k], _h_rbiasden[i][j][k], _s_rbias[i][j][k]); + } + } + } + + } //@} private: vector _jetparams, _centedges, _biases; double _sumwtcentbins[4],_rsumwtcentbins[3]; size_t _ncentbins, _njetbins, _nbiasbins; string _names[2]; private: /// @name Histograms //@{ Histo1DPtr _h_jetspec[4][2][3]; - Histo1DPtr _h_rcp[3][2][3]; + Histo1DPtr _h_rcp[3][2][3]; + Scatter2DPtr _s_rbias[2][2][2]; + Histo1DPtr _h_rbiasnum[2][2][2]; + Histo1DPtr _h_rbiasden[2][2][2]; //@} }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(ALICE_2014_I1263194); } Index: trunk/code/ALICE_2012_I1127497.plot =================================================================== --- trunk/code/ALICE_2012_I1127497.plot (revision 461) +++ trunk/code/ALICE_2012_I1127497.plot (revision 462) @@ -1,20 +1,126 @@ +# BEGIN PLOT /ALICE_2012_I1127497/d01-x01-y01 +Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 00-05 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d02-x01-y01 +Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 05-10 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d03-x01-y01 +Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 10-20 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d04-x01-y01 +Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 20-30 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d05-x01-y01 +Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 30-40 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d06-x01-y01 +Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 40-50 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d07-x01-y01 +Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 50-60 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d08-x01-y01 +Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 60-70 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d09-x01-y01 +Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 70-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ +# END PLOT + # BEGIN PLOT /ALICE_2012_I1127497/d10-x01-y01 Title=charged hadrons ($|\eta|<0.8$) in PbPb at $\sqrt{s}=2.76\ \mathrm{TeV}$, 00-10 XLabel=$p_\perp\ [\mathrm{GeV}]$ YLabel=$E \mathrm{d}^3N/\mathrm{d}^3p\ [\mathrm{GeV}^{-2}]$ LogY=1 FullRange=1 NormalizeToIntegral=1 -#DrawOnly=/home/jewel/trunk/ALICE_2012_I1127497.aida/REF/ALICE_2012_I1127497/d10-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-10_wo_1000-1099.aida/ALICE_2012_I1127497/d10-x01-y01 /home/jewel/trunk/data/jewel-2.0.0_vac_fixed2_5-350.aida/ALICE_2012_I1127497/d10-x01-y01 +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d16-x01-y01 +Title=charged hadrons ($|\eta|<1$), 00-05/pp +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{AA}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d17-x01-y01 +Title=charged hadrons ($|\eta|<1$), 05-10/pp +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{AA}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d18-x01-y01 +Title=charged hadrons ($|\eta|<1$), 10-20/pp +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{AA}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d19-x01-y01 +Title=charged hadrons ($|\eta|<1$), 20-30/pp +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{AA}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d20-x01-y01 +Title=charged hadrons ($|\eta|<1$), 30-40/pp +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{AA}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d21-x01-y01 +Title=charged hadrons ($|\eta|<1$), 40-50/pp +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{AA}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d22-x01-y01 +Title=charged hadrons ($|\eta|<1$), 50-60/pp +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{AA}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d23-x01-y01 +Title=charged hadrons ($|\eta|<1$), 60-70/pp +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{AA}$ +# END PLOT + +# BEGIN PLOT /ALICE_2012_I1127497/d24-x01-y01 +Title=charged hadrons ($|\eta|<1$), 70-80/pp +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{AA}$ # END PLOT # BEGIN PLOT /ALICE_2012_I1127497/d25-x01-y01 Title=charged hadrons ($|\eta|<1$), 00-10/pp XLabel=$p_\perp\ [\mathrm{GeV}]$ YLabel=$R_{AA}$ LogY=0 FullRange=1 YMax=1. -#DrawOnly=/home/jewel/trunk/ALICE_2012_I1127497.aida/REF/ALICE_2012_I1127497/d25-x01-y01 /home/jewel/trunk/eventfiles/jewel-2.0.0_360MeV-boost_md09_0-10_yres_wo_1000-1099.aida/ALICE_2012_I1127497/d25-x01-y01 # END PLOT Index: trunk/code/ATLAS_2012_I1126965.plot =================================================================== --- trunk/code/ATLAS_2012_I1126965.plot (revision 461) +++ trunk/code/ATLAS_2012_I1126965.plot (revision 462) @@ -1,118 +1,154 @@ -# BEGIN PLOT /ATLAS_2012_I1126965/d0* +# BEGIN PLOT /ATLAS_2012_I1126965/* LogY=0 FullRange=1 YMax=1. # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y01 Title=R=0.2, 00-10/60-80 XLabel=$p_\perp\ [\mathrm{GeV}]$ YLabel=$R_{CP}$ -LogY=0 -FullRange=1 -YMax=1. # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y02 Title=R=0.3, 00-10/60-80 XLabel=$p_\perp\ [\mathrm{GeV}]$ YLabel=$R_{CP}$ -LogY=0 -FullRange=1 -YMax=1. # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y03 Title=R=0.4, 00-10/60-80 XLabel=$p_\perp\ [\mathrm{GeV}]$ YLabel=$R_{CP}$ -LogY=0 -FullRange=1 -YMax=1. # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y04 Title=R=0.5, 00-10/60-80 XLabel=$p_\perp\ [\mathrm{GeV}]$ YLabel=$R_{CP}$ -LogY=0 -FullRange=1 -YMax=1. # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/jetpt0* #DrawOnly= # END PLOT -# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y02 -#DrawOnly= -# END PLOT - -# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y03 -#DrawOnly= -# END PLOT - -# BEGIN PLOT /ATLAS_2012_I1126965/d02-x01-y04 -#DrawOnly= +# BEGIN PLOT /ATLAS_2012_I1126965/d03-x01-y01 +Title=R=0.2, 10-20/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d03-x01-y02 -#DrawOnly= +Title=R=0.3, 10-20/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d03-x01-y03 -#DrawOnly= +Title=R=0.4, 10-20/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d03-x01-y04 -#DrawOnly= +Title=R=0.5, 10-20/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ +# END PLOT + +# BEGIN PLOT /ATLAS_2012_I1126965/d04-x01-y01 +Title=R=0.2, 20-30/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d04-x01-y02 -#DrawOnly= +Title=R=0.3, 20-30/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d04-x01-y03 -#DrawOnly= +Title=R=0.4, 20-30/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d04-x01-y04 -#DrawOnly= +Title=R=0.5, 20-30/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ +# END PLOT + +# BEGIN PLOT /ATLAS_2012_I1126965/d05-x01-y01 +Title=R=0.2, 30-40/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d05-x01-y02 -#DrawOnly= +Title=R=0.3, 30-40/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d05-x01-y03 -#DrawOnly= +Title=R=0.4, 30-40/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d05-x01-y04 -#DrawOnly= +Title=R=0.5, 30-40/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ +# END PLOT + +# BEGIN PLOT /ATLAS_2012_I1126965/d06-x01-y01 +Title=R=0.2, 40-50/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d06-x01-y02 -#DrawOnly= +Title=R=0.3, 40-50/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d06-x01-y03 -#DrawOnly= +Title=R=0.4, 40-50/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d06-x01-y04 -#DrawOnly= +Title=R=0.5, 40-50/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ +# END PLOT + +# BEGIN PLOT /ATLAS_2012_I1126965/d07-x01-y01 +Title=R=0.2, 50-60/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d07-x01-y02 -#DrawOnly= +Title=R=0.3, 50-60/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d07-x01-y03 -#DrawOnly= +Title=R=0.4, 50-60/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT # BEGIN PLOT /ATLAS_2012_I1126965/d07-x01-y04 -#DrawOnly= +Title=R=0.5, 50-60/60-80 +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_{CP}$ # END PLOT Index: trunk/code/CMS_2013_I1256590.yoda =================================================================== --- trunk/code/CMS_2013_I1256590.yoda (revision 461) +++ trunk/code/CMS_2013_I1256590.yoda (revision 462) @@ -1,206 +1,206 @@ -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y03 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y03 Path=/REF/CMS_2013_I1256590/d01-x01-y03 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-02 2.500000e-02 2.500000e-02 9.996849e-01 1.534981e-02 1.496040e-02 7.500000e-02 2.500000e-02 2.500000e-02 9.914790e-01 4.180134e-02 4.529186e-02 1.250000e-01 2.500000e-02 2.500000e-02 9.037203e-01 4.181550e-02 3.772316e-02 1.750000e-01 2.500000e-02 2.500000e-02 1.016727e+00 5.312961e-02 5.671926e-02 2.250000e-01 2.500000e-02 2.500000e-02 1.273673e+00 1.439611e-01 1.400741e-01 2.750000e-01 2.500000e-02 2.500000e-02 1.382885e+00 1.060964e-01 1.097498e-01 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y02 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y02 Path=/REF/CMS_2013_I1256590/d01-x01-y02 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-02 2.500000e-02 2.500000e-02 1.266784e+01 0.000000e+00 0.000000e+00 7.500000e-02 2.500000e-02 2.500000e-02 4.565196e+00 0.000000e+00 0.000000e+00 1.250000e-01 2.500000e-02 2.500000e-02 1.470847e+00 0.000000e+00 0.000000e+00 1.750000e-01 2.500000e-02 2.500000e-02 6.714706e-01 0.000000e+00 0.000000e+00 2.250000e-01 2.500000e-02 2.500000e-02 3.694601e-01 0.000000e+00 0.000000e+00 2.750000e-01 2.500000e-02 2.500000e-02 2.007718e-01 0.000000e+00 0.000000e+00 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y01 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d01-x01-y01 Path=/REF/CMS_2013_I1256590/d01-x01-y01 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-02 2.500000e-02 2.500000e-02 1.266784e+01 0.000000e+00 0.000000e+00 7.500000e-02 2.500000e-02 2.500000e-02 4.565196e+00 0.000000e+00 0.000000e+00 1.250000e-01 2.500000e-02 2.500000e-02 1.298714e+00 0.000000e+00 0.000000e+00 1.750000e-01 2.500000e-02 2.500000e-02 6.883952e-01 0.000000e+00 0.000000e+00 2.250000e-01 2.500000e-02 2.500000e-02 4.738880e-01 0.000000e+00 0.000000e+00 2.750000e-01 2.500000e-02 2.500000e-02 2.774874e-01 0.000000e+00 0.000000e+00 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y01 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y01 Path=/REF/CMS_2013_I1256590/d03-x01-y01 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 2.500000e-01 2.500000e-01 4.589493e-02 2.640026e-02 2.607364e-02 7.500000e-01 2.500000e-01 2.500000e-01 3.302992e-01 1.372295e-01 1.308755e-01 1.250000e+00 2.500000e-01 2.500000e-01 7.553553e-01 2.064520e-01 2.113501e-01 1.750000e+00 2.500000e-01 2.500000e-01 1.219367e+00 1.948854e-01 1.904400e-01 2.250000e+00 2.500000e-01 2.500000e-01 1.560549e+00 1.710532e-01 1.668595e-01 2.750000e+00 2.500000e-01 2.500000e-01 1.804274e+00 2.437255e-01 2.817904e-01 3.250000e+00 2.500000e-01 2.500000e-01 2.086064e+00 4.080724e-01 3.258002e-01 3.750000e+00 2.500000e-01 2.500000e-01 2.447121e+00 4.787017e-01 4.655116e-01 4.250000e+00 2.500000e-01 2.500000e-01 2.556013e+00 3.771229e-01 3.566198e-01 4.750000e+00 2.500000e-01 2.500000e-01 1.273626e+00 4.731263e-01 4.790329e-01 5.250000e+00 2.500000e-01 2.500000e-01 1.213483e-01 5.143852e-02 1.293548e-01 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y03 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y03 Path=/REF/CMS_2013_I1256590/d03-x01-y03 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 2.500000e-01 2.500000e-01 9.511278e-01 2.506266e-01 2.631579e-01 7.500000e-01 2.500000e-01 2.500000e-01 9.949875e-01 1.817043e-01 1.879699e-01 1.250000e+00 2.500000e-01 2.500000e-01 8.508772e-01 1.190476e-01 1.190476e-01 1.750000e+00 2.500000e-01 2.500000e-01 8.132832e-01 1.002506e-01 9.398496e-02 2.250000e+00 2.500000e-01 2.500000e-01 7.380952e-01 8.145363e-02 8.145363e-02 2.750000e+00 2.500000e-01 2.500000e-01 7.506266e-01 8.145363e-02 8.145363e-02 3.250000e+00 2.500000e-01 2.500000e-01 8.508772e-01 8.771930e-02 8.771930e-02 3.750000e+00 2.500000e-01 2.500000e-01 1.095238e+00 1.190476e-01 1.127820e-01 4.250000e+00 2.500000e-01 2.500000e-01 1.508772e+00 1.879699e-01 1.879699e-01 4.750000e+00 2.500000e-01 2.500000e-01 1.765664e+00 3.320802e-01 3.320802e-01 5.250000e+00 2.500000e-01 2.500000e-01 2.223058e+00 6.892231e-01 6.892231e-01 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y02 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d03-x01-y02 Path=/REF/CMS_2013_I1256590/d03-x01-y02 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-01 2.500000e-01 2.500000e-01 4.793716e-02 0.000000e+00 0.000000e+00 7.500000e-01 2.500000e-01 2.500000e-01 3.302992e-01 0.000000e+00 0.000000e+00 1.250000e+00 2.500000e-01 2.500000e-01 8.990452e-01 0.000000e+00 0.000000e+00 1.750000e+00 2.500000e-01 2.500000e-01 1.494066e+00 0.000000e+00 0.000000e+00 2.250000e+00 2.500000e-01 2.500000e-01 2.147498e+00 0.000000e+00 0.000000e+00 2.750000e+00 2.500000e-01 2.500000e-01 2.447121e+00 0.000000e+00 0.000000e+00 3.250000e+00 2.500000e-01 2.500000e-01 2.447121e+00 0.000000e+00 0.000000e+00 3.750000e+00 2.500000e-01 2.500000e-01 2.243057e+00 0.000000e+00 0.000000e+00 4.250000e+00 2.500000e-01 2.500000e-01 1.702521e+00 0.000000e+00 0.000000e+00 4.750000e+00 2.500000e-01 2.500000e-01 7.231755e-01 0.000000e+00 0.000000e+00 5.250000e+00 2.500000e-01 2.500000e-01 5.462544e-02 0.000000e+00 0.000000e+00 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y02 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y02 Path=/REF/CMS_2013_I1256590/d04-x01-y02 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 1.097304e+00 9.730405e-02 9.730405e-02 9.843526e-01 0.000000e+00 0.000000e+00 1.310850e+00 1.162410e-01 1.162410e-01 1.032045e+00 0.000000e+00 0.000000e+00 1.565954e+00 1.388630e-01 1.388630e-01 1.032045e+00 0.000000e+00 0.000000e+00 1.870703e+00 1.658865e-01 1.658865e-01 9.689500e-01 0.000000e+00 0.000000e+00 2.234760e+00 1.981700e-01 1.981700e-01 9.537884e-01 0.000000e+00 0.000000e+00 2.669665e+00 2.367350e-01 2.367350e-01 8.407319e-01 0.000000e+00 0.000000e+00 3.189208e+00 2.828065e-01 2.828065e-01 7.648242e-01 0.000000e+00 0.000000e+00 3.809856e+00 3.378430e-01 3.378430e-01 6.132976e-01 0.000000e+00 0.000000e+00 4.551290e+00 4.035905e-01 4.035905e-01 5.238152e-01 0.000000e+00 0.000000e+00 5.437013e+00 4.821335e-01 4.821335e-01 4.473886e-01 0.000000e+00 0.000000e+00 6.495107e+00 5.759600e-01 5.759600e-01 3.644549e-01 0.000000e+00 0.000000e+00 7.759116e+00 6.880480e-01 6.880480e-01 2.922493e-01 0.000000e+00 0.000000e+00 9.269112e+00 8.219480e-01 8.219480e-01 2.418587e-01 0.000000e+00 0.000000e+00 1.107298e+01 9.819200e-01 9.819200e-01 1.764312e-01 0.000000e+00 0.000000e+00 1.322790e+01 1.173000e+00 1.173000e+00 1.370839e-01 0.000000e+00 0.000000e+00 1.580216e+01 1.401255e+00 1.401255e+00 1.015896e-01 0.000000e+00 0.000000e+00 1.887738e+01 1.673970e+00 1.673970e+00 6.957701e-02 0.000000e+00 0.000000e+00 2.255110e+01 1.999745e+00 1.999745e+00 4.917913e-02 0.000000e+00 0.000000e+00 2.693975e+01 2.388910e+00 2.388910e+00 3.368195e-02 0.000000e+00 0.000000e+00 3.218248e+01 2.853820e+00 2.853820e+00 2.306820e-02 0.000000e+00 0.000000e+00 3.844549e+01 3.409190e+00 3.409190e+00 1.287031e-02 0.000000e+00 0.000000e+00 4.592735e+01 4.072655e+00 4.072655e+00 6.532336e-03 0.000000e+00 0.000000e+00 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y03 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y03 Path=/REF/CMS_2013_I1256590/d04-x01-y03 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 1.097304e+00 9.730405e-02 9.730405e-02 1.050239e+00 5.598086e-01 5.645933e-01 1.310850e+00 1.162410e-01 1.162410e-01 6.339713e-01 4.354067e-01 4.354067e-01 1.565954e+00 1.388630e-01 1.388630e-01 6.818182e-01 3.636364e-01 3.684211e-01 1.870703e+00 1.658865e-01 1.658865e-01 4.043062e-01 2.918660e-01 2.918660e-01 2.234760e+00 1.981700e-01 1.981700e-01 2.224880e-01 2.344498e-01 2.392344e-01 2.669665e+00 2.367350e-01 2.367350e-01 5.502392e-02 1.818182e-01 1.866029e-01 3.189208e+00 2.828065e-01 2.828065e-01 -5.502392e-02 1.387560e-01 1.435407e-01 3.809856e+00 3.378430e-01 3.378430e-01 -6.937799e-02 1.052632e-01 1.004785e-01 4.551290e+00 4.035905e-01 4.035905e-01 -9.808612e-02 7.177033e-02 7.177033e-02 5.437013e+00 4.821335e-01 4.821335e-01 -1.028708e-01 5.263158e-02 5.263158e-02 6.495107e+00 5.759600e-01 5.759600e-01 -8.851675e-02 3.827751e-02 3.349282e-02 7.759116e+00 6.880480e-01 6.880480e-01 -8.373206e-02 2.870813e-02 2.392344e-02 9.269112e+00 8.219480e-01 8.219480e-01 -6.937799e-02 1.435407e-02 2.392344e-02 1.107298e+01 9.819200e-01 9.819200e-01 -4.066986e-02 1.913876e-02 1.435407e-02 1.322790e+01 1.173000e+00 1.173000e+00 -3.588517e-02 1.435407e-02 1.435407e-02 1.580216e+01 1.401255e+00 1.401255e+00 -2.631579e-02 9.569378e-03 1.435407e-02 1.887738e+01 1.673970e+00 1.673970e+00 -7.177033e-03 1.435407e-02 4.784689e-03 2.255110e+01 1.999745e+00 1.999745e+00 -7.177033e-03 9.569378e-03 4.784689e-03 2.693975e+01 2.388910e+00 2.388910e+00 -2.392344e-03 0.000000e+00 0.000000e+00 3.218248e+01 2.853820e+00 2.853820e+00 2.392344e-03 0.000000e+00 0.000000e+00 3.844549e+01 3.409190e+00 3.409190e+00 2.392344e-03 0.000000e+00 0.000000e+00 4.592735e+01 4.072655e+00 4.072655e+00 -2.392344e-03 0.000000e+00 0.000000e+00 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y01 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d04-x01-y01 Path=/REF/CMS_2013_I1256590/d04-x01-y01 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 1.097304e+00 9.730405e-02 9.730405e-02 2.033383e+00 4.782024e-01 4.627072e-01 1.310850e+00 1.162410e-01 1.162410e-01 1.656449e+00 3.489591e-01 3.451167e-01 1.565954e+00 1.388630e-01 1.388630e-01 1.709530e+00 2.947629e-01 3.238530e-01 1.870703e+00 1.658865e-01 1.658865e-01 1.370839e+00 2.541162e-01 2.341778e-01 2.234760e+00 1.981700e-01 1.981700e-01 1.189440e+00 2.204898e-01 1.813987e-01 2.669665e+00 2.367350e-01 2.367350e-01 8.954777e-01 1.774117e-01 1.696391e-01 3.189208e+00 2.828065e-01 2.828065e-01 7.068302e-01 1.400369e-01 1.339018e-01 3.809856e+00 3.378430e-01 3.378430e-01 5.491943e-01 1.088062e-01 1.040393e-01 4.551290e+00 4.035905e-01 4.035905e-01 4.267141e-01 8.454043e-02 8.890473e-02 5.437013e+00 4.821335e-01 4.821335e-01 3.421737e-01 6.779131e-02 5.845289e-02 6.495107e+00 5.759600e-01 5.759600e-01 2.700890e-01 4.656965e-02 4.613880e-02 7.759116e+00 6.880480e-01 6.880480e-01 2.065706e-01 3.561760e-02 3.528807e-02 9.269112e+00 8.219480e-01 8.219480e-01 1.709530e-01 2.947629e-02 2.607165e-02 1.107298e+01 9.819200e-01 9.819200e-01 1.328274e-01 2.115517e-02 2.269066e-02 1.322790e+01 1.173000e+00 1.173000e+00 1.015896e-01 1.883195e-02 1.924512e-02 1.580216e+01 1.401255e+00 1.401255e+00 7.648242e-02 1.611232e-02 1.593491e-02 1.887738e+01 1.673970e+00 1.673970e+00 5.667932e-02 1.332960e-02 1.073732e-02 2.255110e+01 1.999745e+00 1.999745e+00 4.006265e-02 1.083773e-02 1.149923e-02 2.693975e+01 2.388910e+00 2.388910e+00 2.831749e-02 8.615028e-03 9.295896e-03 3.218248e+01 2.853820e+00 2.853820e+00 1.909070e-02 6.620015e-03 6.266981e-03 3.844549e+01 3.409190e+00 3.409190e+00 1.170828e-02 4.750580e-03 4.856209e-03 4.592735e+01 4.072655e+00 4.072655e+00 6.636175e-03 2.933691e-03 3.207351e-03 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d02-x01-y01 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d02-x01-y01 Path=/REF/CMS_2013_I1256590/d02-x01-y01 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-02 2.500000e-02 2.500000e-02 6.375449e-01 0.000000e+00 0.000000e+00 7.500000e-02 2.500000e-02 2.500000e-02 8.643001e-01 0.000000e+00 0.000000e+00 1.250000e-01 2.500000e-02 2.500000e-02 9.278651e-01 0.000000e+00 0.000000e+00 1.750000e-01 2.500000e-02 2.500000e-02 9.612714e-01 0.000000e+00 0.000000e+00 2.250000e-01 2.500000e-02 2.500000e-02 9.858084e-01 0.000000e+00 0.000000e+00 2.750000e-01 2.500000e-02 2.500000e-02 9.979320e-01 0.000000e+00 0.000000e+00 -# END YODA_SCATTER2D +END YODA_SCATTER2D -# BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d02-x01-y02 +BEGIN YODA_SCATTER2D /REF/CMS_2013_I1256590/d02-x01-y02 Path=/REF/CMS_2013_I1256590/d02-x01-y02 Title= Type=Scatter2D # xval xerr- xerr+ yval yerr- yerr+ 2.500000e-02 2.500000e-02 2.500000e-02 6.357711e-01 0.000000e+00 0.000000e+00 7.500000e-02 2.500000e-02 2.500000e-02 8.678478e-01 0.000000e+00 0.000000e+00 1.250000e-01 2.500000e-02 2.500000e-02 9.366820e-01 0.000000e+00 0.000000e+00 1.750000e-01 2.500000e-02 2.500000e-02 9.700882e-01 0.000000e+00 0.000000e+00 2.250000e-01 2.500000e-02 2.500000e-02 9.892897e-01 0.000000e+00 0.000000e+00 2.750000e-01 2.500000e-02 2.500000e-02 9.978656e-01 0.000000e+00 0.000000e+00 -# END YODA_SCATTER2D +END YODA_SCATTER2D Index: trunk/code/ALICE_2015_I1343112.plot =================================================================== --- trunk/code/ALICE_2015_I1343112.plot (revision 461) +++ trunk/code/ALICE_2015_I1343112.plot (revision 462) @@ -1,8 +1,67 @@ # BEGIN PLOT /ALICE_2015_I1343112/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=p+p, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$1/N_\text{coll}\, 1/N_\text{evt} \mathrm{d}^2N/\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2015_I1343112/d01-x01-y02 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$1/N_\text{coll}\, 1/N_\text{evt} \mathrm{d}^2N/\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2015_I1343112/d01-x01-y03 +Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$1/N_\text{coll}\, 1/N_\text{evt} \mathrm{d}^2N/\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2015_I1343112/d02-x01-y01 +Title=p+p, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$N_\text{jets}^\text{5 GeV bias}/N_\text{jets}^\text{0 GeV bias}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2015_I1343112/d03-x01-y01 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$N_\text{jets}^\text{0 GeV bias}/N_\text{jets}^\text{5 GeV bias}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2015_I1343112/d03-x01-y02 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 3\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$N_\text{jets}^\text{3 GeV bias}/N_\text{jets}^\text{5 GeV bias}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2015_I1343112/d03-x01-y03 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 7\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$N_\text{jets}^\text{7 GeV bias}/N_\text{jets}^\text{5 GeV bias}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2015_I1343112/d03-x01-y04 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$N_\text{jets}^\text{10 GeV bias}/N_\text{jets}^\text{5 GeV bias}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2015_I1343112/d04-x01-y01 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{AA}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2015_I1343112/d04-x01-y02 +Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{AA}$ +LogY=0 # END PLOT -# ... add more histograms as you need them ... Index: trunk/code/jetspectrum.cc =================================================================== --- trunk/code/jetspectrum.cc (revision 461) +++ trunk/code/jetspectrum.cc (revision 462) @@ -1,148 +1,197 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Projections/FastJets.hh" #include "Rivet/Projections/ChargedFinalState.hh" #include "Rivet/Tools/Logging.hh" #include "Rivet/Tools/ParticleIdUtils.hh" #include namespace Rivet { /// @brief jet pt spectrum class JETSPECTRUM : public Analysis { public: /// Constructor JETSPECTRUM() : Analysis("JETSPECTRUM"), _sumWeightSelected1(0.0),_sumWeightSelected2(0.0) { setNeedsCrossSection(true); } /// Book projections and histograms void init() { ChargedFinalState cfs(-0.9, 0.9, 0.150*GeV); + FinalState fs(-3., 3., 0.); addProjection(cfs, "CFS"); addProjection(FastJets(cfs, FastJets::ANTIKT, 0.2), "Jets1"); addProjection(FastJets(cfs, FastJets::ANTIKT, 0.3), "Jets2"); + addProjection(FastJets(fs, FastJets::ANTIKT, 0.4), "Jets3"); pt_02 = bookHisto1D("jetspec1", 24, 30.0, 150.0, "jet spectrum R=0.2", "$p_\\perp\\ [GeV]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); pt_03 = bookHisto1D("jetspec2", 24, 30.0, 150.0, "jet spectrum R=0.3", "$p_\\perp\\ [GeV]$", "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); n_02 = bookHisto1D("njet1", 20, 0.0, 20.0, "number of jet2 R=0.2", "$N_\\mathrm{jet}$", "$\\mathrm{d}\\sigma/\\mathrm{d}N_\\mathrm{jet}$"); n_03 = bookHisto1D("njet2", 20, 0.0, 20.0, "number of jet2 R=0.3", "$N_\\mathrm{jet}$", "$\\mathrm{d}\\sigma/\\mathrm{d}N_\\mathrm{jet}$"); + _h_jetpt_05 = bookHisto1D("jetpt_eta05",45 , 50.0, 500.0, + "jet spectrum R=0.4 |\\eta|<0.5", + "$p_\\perp\\ [GeV]$", + "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); + _h_jetpt_1 = bookHisto1D("jetpt_eta1", 45, 50.0, 500.0, + "jet spectrum R=0.4 |\\eta|<1", + "$p_\\perp\\ [GeV]$", + "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); + _h_jetpt_2 = bookHisto1D("jetpt_eta2", 45, 50.0, 500.0, + "jet spectrum R=0.4 |\\eta|<0.5", + "$p_\\perp\\ [GeV]$", + "$\\mathrm{d}\\sigma/\\mathrm{d}p_\\perp$"); + _h_jeteta_50 = bookHisto1D("jeteta_pt50",30 , 0.0, 3.0, + "jet eta R=0.4 $p_\\perp > 50$ GeV", + "$|\\eta|$", + "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$"); + _h_jeteta_100 = bookHisto1D("jeteta_pt100",30 , 0.0, 3.0, + "jet eta R=0.4 $p_\\perp > 100$ GeV", + "$|\\eta|$", + "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$"); + _h_jeteta_250 = bookHisto1D("jeteta_pt250",30 , 0.0, 3.0, + "jet eta R=0.4 $p_\\perp > 250$ GeV", + "$|\\eta|$", + "$\\mathrm{d}\\sigma/\\mathrm{d}|\\eta|$"); } /// Do the analysis void analyze(const Event& event) { const double weight = event.weight(); // const double sigma = 15.; // std::cout<(event, "Jets1"); /* std::list jetlist1; foreach (Jet jet , fastjets1.jets()){ // std::cout<<"Et before: "<30.) jets1.push_back(jet); } const FastJets& fastjets2 = applyProjection(event, "Jets2"); /* std::list jetlist2; foreach (Jet jet , fastjets2.jets()){ // std::cout<<"Et before: "<30.) jets2.push_back(jet); } n_02->fill(jets1.size(),weight); if (jets1.size() > 0) { _sumWeightSelected1 += event.weight(); foreach (Jet jet, jets1){ pt_02->fill(jet.momentum().pT(),weight); } } n_03->fill(jets2.size(),weight); if (jets2.size() > 0) { _sumWeightSelected2 += event.weight(); foreach (Jet jet, jets2){ pt_03->fill(jet.momentum().pT(),weight); } } + + const FastJets& fastjets3 = applyProjection(event, "Jets3"); + const Jets alljets3 = fastjets3.jetsByPt(); + foreach (Jet jet, alljets3) { + double pt = jet.momentum().pT(); + double eta = jet.momentum().abseta(); + if (eta < 0.5) _h_jetpt_05->fill(pt, weight); + if (eta < 1.0) _h_jetpt_1->fill(pt, weight); + if (eta < 2.0) _h_jetpt_2->fill(pt, weight); + if (pt > 50.) _h_jeteta_50->fill(eta, weight); + if (pt > 100.) _h_jeteta_100->fill(eta, weight); + if (pt > 250.) _h_jeteta_250->fill(eta, weight); + + + } } /// Finalize void finalize() { /* scale(pt_02,1./_sumWeightSelected1); scale(pt_03,1./_sumWeightSelected2);*/ scale(pt_02,crossSection()/sumOfWeights()); scale(pt_03,crossSection()/sumOfWeights()); scale(n_02,1./sumOfWeights()); scale(n_03,1./sumOfWeights()); + scale(_h_jetpt_05,crossSection()/sumOfWeights()); + scale(_h_jetpt_1,crossSection()/sumOfWeights()); + scale(_h_jetpt_2,crossSection()/sumOfWeights()); + scale(_h_jeteta_50,crossSection()/sumOfWeights()); + scale(_h_jeteta_100,crossSection()/sumOfWeights()); + scale(_h_jeteta_250,crossSection()/sumOfWeights()); getLog() << Log::DEBUG << "sumOfWeights() = " << sumOfWeights() << std::endl; } private: double _sumWeightSelected1; double _sumWeightSelected2; Histo1DPtr pt_02; Histo1DPtr pt_03; Histo1DPtr n_02; Histo1DPtr n_03; + Histo1DPtr _h_jetpt_05, _h_jetpt_1, _h_jetpt_2; + Histo1DPtr _h_jeteta_50, _h_jeteta_100, _h_jeteta_250; }; // This global object acts as a hook for the plugin system AnalysisBuilder plugin_JETSPECTRUM; } Index: trunk/code/ALICE_2017_I1512107_4MOMSUB.plot =================================================================== --- trunk/code/ALICE_2017_I1512107_4MOMSUB.plot (revision 461) +++ trunk/code/ALICE_2017_I1512107_4MOMSUB.plot (revision 462) @@ -1,19 +1,28 @@ # BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/chfrac_rat LogY=0 # END PLOT # BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/jetmass LogY=0 # END PLOT # BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/d01-x01-y01 +XLabel=$M_\text{ch jet}\ [\text{GeV}^2]$ +YLabel=$1/N_\text{jets}\ \mathrm{d}N/\mathrm{d}M_\text{ch jet}\ [\mathrm{GeV}^{-2}]$ +Title=$60\,\mathrm{GeV} < p_{\perp,\,\text{ch jet}} < 80\,\mathrm{GeV}$ - 0-10 Pb+Pb LogY=0 # END PLOT # BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/d02-x01-y01 +XLabel=$M_\text{ch jet}\ [\text{GeV}^2]$ +YLabel=$1/N_\text{jets}\ \mathrm{d}N/\mathrm{d}M_\text{ch jet}\ [\mathrm{GeV}^{-2}]$ +Title=$80\,\mathrm{GeV} < p_{\perp,\,\text{ch jet}} < 100\,\mathrm{GeV}$ - 0-10 Pb+Pb LogY=0 # END PLOT # BEGIN PLOT /ALICE_2017_I1512107_4MOMSUB/d03-x01-y01 +XLabel=$M_\text{ch jet}\ [\text{GeV}^2]$ +YLabel=$1/N_\text{jets}\ \mathrm{d}N/\mathrm{d}M_\text{ch jet}\ [\mathrm{GeV}^{-2}]$ +Title=$100\,\mathrm{GeV} < p_{\perp,\,\text{ch jet}} < 120\,\mathrm{GeV}$ - 0-10 Pb+Pb LogY=0 # END PLOT Index: trunk/code/ATLAS_CONF_2014_028.plot =================================================================== --- trunk/code/ATLAS_CONF_2014_028.plot (revision 461) +++ trunk/code/ATLAS_CONF_2014_028.plot (revision 462) @@ -1,8 +1,546 @@ # BEGIN PLOT /ATLAS_CONF_2014_028/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=$R=0.2$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ # END PLOT -# ... add more histograms as you need them ... +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x01-y02 +Title=$R=0.2$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x01-y03 +Title=$R=0.2$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x01-y04 +Title=$R=0.2$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x02-y01 +Title=$R=0.2$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x02-y02 +Title=$R=0.2$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x02-y03 +Title=$R=0.2$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x02-y04 +Title=$R=0.2$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x03-y01 +Title=$R=0.2$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x03-y02 +Title=$R=0.2$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x03-y03 +Title=$R=0.2$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d01-x03-y04 +Title=$R=0.2$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x01-y01 +Title=$R=0.3$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x01-y02 +Title=$R=0.3$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x01-y03 +Title=$R=0.3$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x01-y04 +Title=$R=0.3$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x02-y01 +Title=$R=0.3$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x02-y02 +Title=$R=0.3$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x02-y03 +Title=$R=0.3$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x02-y04 +Title=$R=0.3$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x03-y01 +Title=$R=0.3$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x03-y02 +Title=$R=0.3$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x03-y03 +Title=$R=0.3$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d02-x03-y04 +Title=$R=0.3$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x01-y01 +Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x01-y02 +Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x01-y03 +Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x01-y04 +Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x02-y01 +Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x02-y02 +Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x02-y03 +Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x02-y04 +Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x03-y01 +Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x03-y02 +Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x03-y03 +Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d03-x03-y04 +Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$R_{\Delta R}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x01-y01 +Title=$R=0.2$, $E_\perp^\text{test} > 80\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x01-y02 +Title=$R=0.2$, $E_\perp^\text{test} > 80\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x01-y03 +Title=$R=0.2$, $E_\perp^\text{test} > 80\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x01-y04 +Title=$R=0.2$, $E_\perp^\text{test} > 80\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x02-y01 +Title=$R=0.2$, $E_\perp^\text{test} > 90\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x02-y02 +Title=$R=0.2$, $E_\perp^\text{test} > 90\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x02-y03 +Title=$R=0.2$, $E_\perp^\text{test} > 90\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x02-y04 +Title=$R=0.2$, $E_\perp^\text{test} > 90\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x03-y01 +Title=$R=0.2$, $E_\perp^\text{test} > 110\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x03-y02 +Title=$R=0.2$, $E_\perp^\text{test} > 110\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x03-y03 +Title=$R=0.2$, $E_\perp^\text{test} > 110\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d11-x03-y04 +Title=$R=0.2$, $E_\perp^\text{test} > 110\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x01-y01 +Title=$R=0.3$, $E_\perp^\text{test} > 80\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x01-y02 +Title=$R=0.3$, $E_\perp^\text{test} > 80\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x01-y03 +Title=$R=0.3$, $E_\perp^\text{test} > 80\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x01-y04 +Title=$R=0.3$, $E_\perp^\text{test} > 80\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x02-y01 +Title=$R=0.3$, $E_\perp^\text{test} > 90\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x02-y02 +Title=$R=0.3$, $E_\perp^\text{test} > 90\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x02-y03 +Title=$R=0.3$, $E_\perp^\text{test} > 90\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x02-y04 +Title=$R=0.3$, $E_\perp^\text{test} > 90\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x03-y01 +Title=$R=0.3$, $E_\perp^\text{test} > 110\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x03-y02 +Title=$R=0.3$, $E_\perp^\text{test} > 110\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x03-y03 +Title=$R=0.3$, $E_\perp^\text{test} > 110\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d12-x03-y04 +Title=$R=0.3$, $E_\perp^\text{test} > 110\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x01-y01 +Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x01-y02 +Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x01-y03 +Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x01-y04 +Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x02-y01 +Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x02-y02 +Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x02-y03 +Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x02-y04 +Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x03-y01 +Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 0-10 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x03-y02 +Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 10-20 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x03-y03 +Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 20-40 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d13-x03-y04 +Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\text{d}\,R_{\Delta R}/\text{d}\,E_\perp^\text{nbr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x01-y01 +Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 0-10 / 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x01-y02 +Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 10-20 / 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x01-y03 +Title=$R=0.4$, $E_\perp^\text{nbr} > 30\,\text{GeV}$, 20-40 / 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x02-y01 +Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 0-10 / 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x02-y02 +Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 10-20 / 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x02-y03 +Title=$R=0.4$, $E_\perp^\text{nbr} > 45\,\text{GeV}$, 20-40 / 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x03-y01 +Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 0-10 / 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x03-y02 +Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 10-20 / 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d21-x03-y03 +Title=$R=0.4$, $E_\perp^\text{nbr} > 60\,\text{GeV}$, 20-40 / 40-80 +XLabel=$E_\perp^\text{test}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x01-y01 +Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 0-10 / 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x01-y02 +Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 10-20 / 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x01-y03 +Title=$R=0.4$, $E_\perp^\text{test} > 80\,\text{GeV}$, 20-40 / 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x02-y01 +Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 0-10 / 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x02-y02 +Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 10-20 / 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x02-y03 +Title=$R=0.4$, $E_\perp^\text{test} > 90\,\text{GeV}$, 20-40 / 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x03-y01 +Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 0-10 / 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x03-y02 +Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 10-20 / 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +# BEGIN PLOT /ATLAS_CONF_2014_028/d22-x03-y03 +Title=$R=0.4$, $E_\perp^\text{test} > 110\,\text{GeV}$, 20-40 / 40-80 +XLabel=$E_\perp^\text{nbr}\ [\text{GeV}]$ +YLabel=$\rho_{R_{\Delta R}}$ +# END PLOT + +## BEGIN PLOT /ATLAS_CONF_2014_028/d9* +#Title=dummy histogram +## END PLOT + + + Index: trunk/code/ALICE_2014_I1263194.plot =================================================================== --- trunk/code/ALICE_2014_I1263194.plot (revision 461) +++ trunk/code/ALICE_2014_I1263194.plot (revision 462) @@ -1,6 +1,327 @@ -## BEGIN PLOT /ALICE_2014_I1263194/* -#RatioPlotYMin=0 -#RatioPlotYMin=1 -#RatioPlotYMax=10 -## END PLOT +# BEGIN PLOT /ALICE_2014_I1263194/d02-x01-y01 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d02-x01-y02 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d03-x01-y01 +Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d03-x01-y02 +Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d04-x01-y01 +Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d04-x01-y02 +Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d05-x01-y01 +Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d05-x01-y02 +Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d06-x01-y01 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d06-x01-y02 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d07-x01-y01 +Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d07-x01-y02 +Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + + +# BEGIN PLOT /ALICE_2014_I1263194/d08-x01-y01 +Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d08-x01-y02 +Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d09-x01-y01 +Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d09-x01-y02 +Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d10-x01-y01 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d10-x01-y02 +Title=Pb+Pb 00-10, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d11-x01-y01 +Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d11-x01-y02 +Title=Pb+Pb 10-30, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d12-x01-y01 +Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d12-x01-y02 +Title=Pb+Pb 30-50, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d13-x01-y01 +Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d13-x01-y02 +Title=Pb+Pb 50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$10^6\,1/N_\text{coll}\,1/N_\text{evt}\,\mathrm{d}^2 N /\mathrm{d}p_\perp\,\mathrm{d}\eta\ [\mathrm{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d14-x01-y01 +Title=Pb+Pb 0-10, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=N( PT(leading track) > 0.15 GeV ) / N( PT(leading track) > 5 GeV ) +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d14-x01-y02 +Title=Pb+Pb 0-10, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=N( PT(leading track) > 0.15 GeV ) / N( PT(leading track) > 5 GeV ) +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d15-x01-y01 +Title=Pb+Pb 50-80, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=N( PT(leading track) > 0.15 GeV ) / N( PT(leading track) > 5 GeV ) +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d15-x01-y02 +Title=Pb+Pb 50-80, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=N( PT(leading track) > 0.15 GeV ) / N( PT(leading track) > 5 GeV ) +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d16-x01-y01 +Title=Pb+Pb 0-10, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=N( PT(leading track) > 10 GeV ) / N( PT(leading track) > 5 GeV ) +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d16-x01-y02 +Title=Pb+Pb 0-10, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=N( PT(leading track) > 10 GeV ) / N( PT(leading track) > 5 GeV ) +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d17-x01-y01 +Title=Pb+Pb 50-80, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=N( PT(leading track) > 10 GeV ) / N( PT(leading track) > 5 GeV ) +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d17-x01-y02 +Title=Pb+Pb 50-80, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=N( PT(leading track) > 10 GeV ) / N( PT(leading track) > 5 GeV ) +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d18-x01-y01 +Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d18-x01-y02 +Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d19-x01-y01 +Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d19-x01-y02 +Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d20-x01-y01 +Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d20-x01-y02 +Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 0.15\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d21-x01-y01 +Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d21-x01-y02 +Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d22-x01-y01 +Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d22-x01-y02 +Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d23-x01-y01 +Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d23-x01-y02 +Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 5\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d24-x01-y01 +Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d24-x01-y02 +Title=Pb+Pb 00-10/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d25-x01-y01 +Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d25-x01-y02 +Title=Pb+Pb 10-30/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d26-x01-y01 +Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.2$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT + +# BEGIN PLOT /ALICE_2014_I1263194/d26-x01-y02 +Title=Pb+Pb 30-50/50-80, $|\eta| < 0.5$, $p_\perp^\text{lead,ch} > 10\,\mathrm{GeV}$, $R=0.3$ +XLabel=$p_\perp\ [\mathrm{GeV}]$ +YLabel=$R_\text{cp}$ +LogY=0 +# END PLOT Index: trunk/code/ALICE_2017_I1512107.cc =================================================================== --- trunk/code/ALICE_2017_I1512107.cc (revision 461) +++ trunk/code/ALICE_2017_I1512107.cc (revision 462) @@ -1,515 +1,520 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/ChargedFinalState.hh" #include "Rivet/Projections/FastJets.hh" #include "fastjet/ClusterSequence.hh" #include "fastjet/ClusterSequence.hh" +#ifndef MODE +#warning MODE not defined, using default +#define MODE 0 +#endif + namespace Rivet { using namespace fastjet; #if MODE==0 class ALICE_2017_I1512107_4MOMSUB : public Analysis { #elif MODE==1 class ALICE_2017_I1512107_GRIDSUB1 : public Analysis { #elif MODE==2 class ALICE_2017_I1512107_GRIDSUB2 : public Analysis { #elif MODE==3 class ALICE_2017_I1512107_CHARGED : public Analysis { #endif public: /// Constructor #if MODE==0 ALICE_2017_I1512107_4MOMSUB() : Analysis("ALICE_2017_I1512107_4MOMSUB") #elif MODE==1 ALICE_2017_I1512107_GRIDSUB1() : Analysis("ALICE_2017_I1512107_GRIDSUB1") #elif MODE==2 ALICE_2017_I1512107_GRIDSUB2() : Analysis("ALICE_2017_I1512107_GRIDSUB2") #elif MODE==3 ALICE_2017_I1512107_CHARGED() : Analysis("ALICE_2017_I1512107_CHARGED") #endif { _jetR = 0.4; _ptmin = 50.; _etamax = 0.5; _etaloose = 0.9; _nphibins = 120; _netabins = 160; //_nphibins = 60; _netabins = 80; _ptscalefac = 0.75; _massscalefac = 0.67; } /// @name Analysis methods //@{ /// Book histograms and initialise projections before the run void init() { _ptedges += 60.0, 80.0, 100.0, 120.0; //FinalState fs(-_etaloose, _etaloose, 0.*GeV); FinalState fs(-3., 3., 0.*GeV); addProjection(fs,"FS"); ChargedFinalState cfs(-_etaloose, _etaloose, 0.*GeV); addProjection(cfs,"CFS"); FastJets tfj(cfs, FastJets::ANTIKT, _jetR); addProjection(tfj, "TrackJets"); FastJets ffj(fs, FastJets::ANTIKT, _jetR); addProjection(ffj, "FullJets"); //_h_scalefac = bookHisto1D("scalefac", 50, 0., 1.); _h_jetmass[0] = bookHisto1D(1, 1, 1); _h_jetmass[1] = bookHisto1D(2, 1, 1); _h_jetmass[2] = bookHisto1D(3, 1, 1); _h_jetpt = bookHisto1D("jetpT", 40, 100., 400.); _h_jetpt_resc = bookHisto1D("jetpT_resc", 40, 100., 400.); _h_JetMass = bookHisto1D("JetMass", 50, 0., 100.); _h_JetpT = bookHisto1D("JetpT", 30, 100, 400); } PseudoJets extractThermalMomenta(const Event & event) { PseudoJets thermom; foreach (const HepMC::GenParticle* p, particles(event.genEvent())) { FourMomentum mom(p->momentum()); if (p->status() == 3) { PseudoJet pjet(mom.px(),mom.py(),mom.pz(),mom.E()); thermom.push_back(pjet); } } return thermom; } /// 4-momentum subtraction FourMomentum SubtractJetMom(Jet jet, PseudoJets * thermom, double fac = 1.) { if (thermom->empty()) return jet.momentum(); FourMomentum sub(0.,0.,0.,0.); foreach (Particle part, jet.constituents()) { if (part.E() < 1e-5 && part.E() > 1e-7) { FourMomentum test(part.momentum()); test *= 10000.; //cout<<"Found dummy :"< > fillGrid(const Event & event, std::set > & gridevent) { // initialise grid vector > grid; grid.resize(_nphibins); for (size_t i = 0; i < _nphibins; ++i) { grid[i].resize(_netabins); for (size_t k = 0; k < _nphibins; ++k) { grid[i][k] = FourMomentum(0., 0., 0., 0.); } } // fill grid foreach (const HepMC::GenParticle* p, particles(event.genEvent())) { FourMomentum mom(p->momentum()); if (fabs(mom.eta()) < 4. && mom.E() > 1e-5) { size_t phibin, etabin; phibin = int(mom.phi(ZERO_2PI)*_nphibins/(2.*M_PI)); etabin = int((mom.eta() + 4.)*_netabins/8.); if (phibin < 0 || phibin > _nphibins-1) {std::cout<<"Error: "< "< _netabins-1) {std::cout<<"Error: "< "<status() == 3) { grid[phibin][etabin] -= mom; } else { grid[phibin][etabin] += mom; }*/ double cellphi((phibin+0.5)*2*M_PI/_nphibins), celleta((etabin+0.5)*8./_netabins-4.); double celltheta(2.*atan(exp(-celleta))); if (p->status() == 3) { double rho(grid[phibin][etabin].E() - mom.E()); grid[phibin][etabin].setE(rho); grid[phibin][etabin].setPx(fabs(rho)*sin(celltheta)*cos(cellphi)); grid[phibin][etabin].setPy(fabs(rho)*sin(celltheta)*sin(cellphi)); grid[phibin][etabin].setPz(fabs(rho)*cos(celltheta)); } else { double rho(grid[phibin][etabin].E() + mom.E()); grid[phibin][etabin].setE(rho); grid[phibin][etabin].setPx(rho*sin(celltheta)*cos(cellphi)); grid[phibin][etabin].setPy(rho*sin(celltheta)*sin(cellphi)); grid[phibin][etabin].setPz(rho*cos(celltheta)); } gridevent.insert(std::pair(phibin,etabin)); } } return grid; } /// determine fraction of pt carried by charged particles in jet double getScaleFac(const Jet jet, double ptcut, const FinalState * fs = NULL) { Particles parts; if (fs) parts = fs->particles(); else parts = jet.constituents(); double scalefac; double ntracks(0.), nparts(0.); FourMomentum chmom(0.,0.,0.,0.), allmom(0.,0.,0.,0.); foreach (Particle part, parts) { if (deltaR(part,jet) < _jetR && part.momentum().E() > 1e-5) { nparts += 1.; allmom += part.momentum(); if (PID::charge(part.pdgId()) != 0 && part.pT() > ptcut) { ntracks += 1.; chmom += part.momentum(); } } } scalefac = (nparts>0.?ntracks/nparts:1); //scalefac = (allmom.pT()>0.?chmom.pT()/allmom.pT():1.); return scalefac; } /// initialise and fill grid for grid subtration version 1 void fillJetGrid(Jet jet, PseudoJets * thermom, vector > & grid, vector > & thgrid, std::set > & gridevent) { // fill grid foreach (Particle part, jet.constituents()) { FourMomentum mom(part.momentum()); size_t phibin, etabin; if (mom.E() < 1e-5 && mom.E() > 1e-7) { /*FourMomentum test(mom); test *= 10000.; //cout<<"Found dummy :"< _nphibins-1) {std::cout<<"Error: "< "< _netabins-1) {std::cout<<"Error: "< "<(phibin,etabin)); } else cout<<"Error: did not find scattering centre matching dummy.\n";*/ } else { phibin = int(mom.phi(ZERO_2PI)*_nphibins/(2.*M_PI)); etabin = int((mom.eta() + 4.)*_netabins/8.); if (phibin < 0 || phibin > _nphibins-1) {std::cout<<"Error: "< "< _netabins-1) {std::cout<<"Error: "< "<(phibin,etabin)); } } foreach (PseudoJet pjet, *thermom) { FourMomentum mom = Get4Mom(pjet); double dR(deltaR(mom, jet)); if (dR < _jetR) { size_t phibin, etabin; phibin = int(mom.phi(ZERO_2PI)*_nphibins/(2.*M_PI)); etabin = int((mom.eta() + 4.)*_netabins/8.); if (phibin < 0 || phibin > _nphibins-1) {std::cout<<"Error: "< "< _netabins-1) {std::cout<<"Error: "< "<(phibin,etabin)); } } return; } /// 4-momentum adapter FourMomentum Get4Mom(PseudoJet pjet){ FourMomentum mom(pjet.E(), pjet.px(), pjet.py(), pjet.pz()); return mom; } /// Perform the per-event analysis void analyze(const Event& event) { const double weight = event.weight(); const double cent = (event.genEvent()->heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.); if (cent > 0.1) vetoEvent; #if MODE==0 // cout<<"Karo: event number "<event_number()< _ptmin*GeV; Cut cuts = Cuts::abseta < 2.0 && Cuts::pT > _ptmin*GeV; const Jets jets = applyProjection(event, "FullJets").jetsByPt(cuts); PseudoJets thermom = extractThermalMomenta(event); foreach (Jet jet, jets) { //double scalefac; //scalefac = getScaleFac(jet, 0.); //_h_scalefac->fill(scalefac, weight); FourMomentum jetsub = SubtractJetMom(jet, &thermom, 1.); // if (jet.momentum().pT() > 100.) { // cout<<"Karo: "< "<0.?sqrt(mass2):-1.); //mass *= scalefac; mass *= _massscalefac; //double jetpt(scalefac*jetsub.pT()); double jetpt(jetsub.pT()); _h_jetpt->fill(jetpt,weight); jetpt *= _ptscalefac; _h_jetpt_resc->fill(jetpt,weight); for (size_t i=0; i<3; ++i) { if (jetpt > _ptedges[i] && jetpt <= _ptedges[i+1]) _h_jetmass[i]->fill(mass,weight); } } if (jetsub.pT() > 100.*GeV) { // cout<<"Karo: "<0.?sqrt(mass2):-1.); _h_JetMass->fill(mass,weight); _h_JetpT->fill(jetsub.pT(),weight); } } #elif MODE==1 //Cut cuts = Cuts::abseta < _etaloose && Cuts::pT > _ptmin*GeV; Cut cuts = Cuts::abseta < 2.0 && Cuts::pT > _ptmin*GeV; const Jets jets = applyProjection(event, "FullJets").jetsByPt(cuts); PseudoJets thermom = extractThermalMomenta(event); // initialise grid vector > grid, thgrid; grid.resize(_nphibins); thgrid.resize(_nphibins); for (size_t i = 0; i < _nphibins; ++i) { grid[i].resize(_netabins); thgrid[i].resize(_netabins); for (size_t k = 0; k < _nphibins; ++k) { grid[i][k] = FourMomentum(0., 0., 0., 0.); thgrid[i][k] = FourMomentum(0., 0., 0., 0.); } } std::set > gridevent; foreach (Jet jet, jets) { fillJetGrid(jet, &thermom, grid, thgrid, gridevent); } double negenergy(0.); PseudoJets pevent; for (set >::iterator it=gridevent.begin(); it!=gridevent.end(); ++it) { size_t i,k; i = it->first; k = it->second; FourMomentum diff(grid[i][k] - thgrid[i][k]); if (grid[i][k].pT() > thgrid[i][k].pT() && diff.E() > 0.) { PseudoJet pjet(diff.px(),diff.py(),diff.pz(),diff.E()); pevent.push_back(pjet); } else negenergy += diff.E(); } JetDefinition jet_def(antikt_algorithm, _jetR); ClusterSequence cs(pevent, jet_def); PseudoJets pjets = sorted_by_pt(cs.inclusive_jets(_ptmin)); //cout<<"size comparison: "<0.?sqrt(mass2):-1.); mass *= _massscalefac; double jetpt(jetsub.pT()); _h_jetpt->fill(jetpt,weight); jetpt *= _ptscalefac; _h_jetpt_resc->fill(jetpt,weight); for (size_t i=0; i<3; ++i) { if (jetpt > _ptedges[i] && jetpt <= _ptedges[i+1]){ _h_jetmass[i]->fill(mass,weight); } } } if (jetsub.pT() > 100.*GeV) { double mass2(jetsub.mass2()); double mass(mass2>0.?sqrt(mass2):-1.); _h_JetMass->fill(mass,weight); _h_JetpT->fill(jetsub.pT(),weight); } } #elif MODE==2 const FinalState parts = applyProjection(event, "FS"); PseudoJets pevent, pjets; std::set > gridevent; vector > grid = fillGrid(event, gridevent); for (set >::iterator it=gridevent.begin(); it!=gridevent.end(); ++it) { size_t i,k; i = it->first; k = it->second; if (grid[i][k].E() > 0. ) { PseudoJet part(grid[i][k].px(),grid[i][k].py(),grid[i][k].pz(),grid[i][k].E()); pevent.push_back(part); } } JetDefinition jet_def(antikt_algorithm, _jetR); ClusterSequence cs(pevent, jet_def); pjets = sorted_by_pt(cs.inclusive_jets(_ptmin)); foreach (PseudoJet pjet, pjets) { FourMomentum jetsub = Get4Mom(pjet); if (abs(jetsub.eta()) < _etamax) { //double scalefac; //scalefac = getScaleFac(Jet(pjet),0.,&parts); //_h_scalefac->fill(scalefac,weight); double mass2(jetsub.mass2()); double mass(mass2>0.?sqrt(mass2):-1.); //mass *= scalefac; mass *= _massscalefac; //double jetpt(scalefac*jetsub.pT()); double jetpt(jetsub.pT()); _h_jetpt->fill(jetpt,weight); jetpt *= _ptscalefac; _h_jetpt_resc->fill(jetpt,weight); for (size_t i=0; i<3; ++i) { if (jetpt > _ptedges[i] && jetpt <= _ptedges[i+1]){ _h_jetmass[i]->fill(mass,weight); } } } if (jetsub.pT() > 100.*GeV) { double mass2(jetsub.mass2()); double mass(mass2>0.?sqrt(mass2):-1.); _h_JetMass->fill(mass,weight); _h_JetpT->fill(jetsub.pT(),weight); } } #elif MODE==3 const FinalState parts = applyProjection(event, "FS"); Cut cuts = Cuts::abseta < _etaloose && Cuts::pT > _ptmin*GeV; const Jets jets = applyProjection(event, "TrackJets").jetsByPt(cuts); foreach (Jet jet, jets) { //double scalefac; //scalefac = getScaleFac(jet, 0., &parts); //_h_scalefac->fill(scalefac, weight); PseudoJets thermom = extractThermalMomenta(event); //FourMomentum jetsub = SubtractJetMom(jet, &thermom, scalefac); FourMomentum jetsub = SubtractJetMom(jet, &thermom, _ptscalefac); if (abs(jetsub.eta()) < _etamax) { double mass2(jetsub.mass2()); double mass(mass2>0.?sqrt(mass2):-1.); double jetpt(jetsub.pT()); _h_jetpt_resc->fill(jetpt,weight); for (size_t i=0; i<3; ++i) { if (jetpt > _ptedges[i] && jetpt <= _ptedges[i+1]) _h_jetmass[i]->fill(mass,weight); } } } #endif } /// Normalise histograms etc., after the run void finalize() { //normalize(_h_scalefac); for (size_t i=0; i<3; ++i) { normalize(_h_jetmass[i]); } scale(_h_jetpt, crossSection()/sumOfWeights()); scale(_h_jetpt_resc, crossSection()/sumOfWeights()); normalize(_h_JetMass); scale(_h_JetpT, crossSection()/picobarn/sumOfWeights()); } private: double _jetR, _ptmin, _etamax, _etaloose; double _massscalefac, _ptscalefac; size_t _netabins, _nphibins; vector _ptedges; //Histo1DPtr _h_scalefac; Histo1DPtr _h_jetmass[3]; Histo1DPtr _h_jetpt, _h_jetpt_resc; Histo1DPtr _h_JetMass, _h_JetpT; }; // The hook for the plugin system #if MODE==0 DECLARE_RIVET_PLUGIN(ALICE_2017_I1512107_4MOMSUB); #elif MODE==1 DECLARE_RIVET_PLUGIN(ALICE_2017_I1512107_GRIDSUB1); #elif MODE==2 DECLARE_RIVET_PLUGIN(ALICE_2017_I1512107_GRIDSUB2); #elif MODE==3 DECLARE_RIVET_PLUGIN(ALICE_2017_I1512107_CHARGED); #endif } Index: trunk/code/ATLAS_2013_I1240088.plot =================================================================== --- trunk/code/ATLAS_2013_I1240088.plot (revision 461) +++ trunk/code/ATLAS_2013_I1240088.plot (revision 462) @@ -1,13 +1,295 @@ # BEGIN PLOT /ATLAS_2013_I1240088/* -#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] +YMin=0.5 +YMax=0.8 NormalizeToIntegral=1 -LogY=0 # END PLOT -# BEGIN PLOT /ATLAS_2013_I1240088/* -YMin=0.5 -YMax=0.8 +# BEGIN PLOT /ATLAS_2013_I1240088/d01-x01-y01 +Title=$5-10$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d01-x02-y01 +Title=$5-10$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d01-x01-y02 +Title=$5-10$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d01-x02-y02 +Title=$5-10$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d01-x01-y03 +Title=$5-10$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d01-x02-y03 +Title=$5-10$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d01-x01-y04 +Title=$5-10$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d01-x02-y04 +Title=$5-10$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d02-x01-y01 +Title=$10-20$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d02-x02-y01 +Title=$10-20$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d02-x01-y02 +Title=$10-20$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d02-x02-y02 +Title=$10-20$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d02-x01-y03 +Title=$10-20$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d02-x02-y03 +Title=$10-20$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d02-x01-y04 +Title=$10-20$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d02-x02-y04 +Title=$10-20$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d03-x01-y01 +Title=$20-30$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +PolyMarker=* +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d03-x02-y01 +Title=$20-30$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d03-x01-y02 +Title=$20-30$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d03-x02-y02 +Title=$20-30$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d03-x01-y03 +Title=$20-30$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d03-x02-y03 +Title=$20-30$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d03-x01-y04 +Title=$20-30$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d03-x02-y04 +Title=$20-30$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d04-x01-y01 +Title=$30-40$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d04-x02-y01 +Title=$30-40$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d04-x01-y02 +Title=$30-40$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d04-x02-y02 +Title=$30-40$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d04-x01-y03 +Title=$30-40$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d04-x02-y03 +Title=$30-40$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d04-x01-y04 +Title=$30-40$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d04-x02-y04 +Title=$30-40$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d05-x01-y01 +Title=$40-50$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d05-x02-y01 +Title=$40-50$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d05-x01-y02 +Title=$40-50$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d05-x02-y02 +Title=$40-50$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d05-x01-y03 +Title=$40-50$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d05-x02-y03 +Title=$40-50$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d05-x01-y04 +Title=$40-50$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d05-x02-y04 +Title=$40-50$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d06-x01-y01 +Title=$50-60$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d06-x02-y01 +Title=$50-60$, $45 < p_\perp < 60\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d06-x01-y02 +Title=$50-60$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d06-x02-y02 +Title=$50-60$, $60 < p_\perp < 80\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d06-x01-y03 +Title=$50-60$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d06-x02-y03 +Title=$50-60$, $80 < p_\perp < 110\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d06-x01-y04 +Title=$50-60$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{meas}\ [\text{GeV}^{-1}]$ +# END PLOT + +# BEGIN PLOT /ATLAS_2013_I1240088/d06-x02-y04 +Title=$50-60$, $110 < p_\perp < 160\,\text{GeV}$ +XLabel=$\Delta \phi$ +YLabel=$1/N_\text{jet}\ \left.\frac{\text{d}^2 N_\text{jet}}{\text{d}p_\perp\,\text{d}\Delta\phi}\right|_\text{corr}\ [\text{GeV}^{-1}]$ # END PLOT Index: trunk/code/ATLAS_2013_I1240088.cc =================================================================== --- trunk/code/ATLAS_2013_I1240088.cc (revision 461) +++ trunk/code/ATLAS_2013_I1240088.cc (revision 462) @@ -1,124 +1,125 @@ // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Tools/Logging.hh" #include "Rivet/Projections/FinalState.hh" #include "Rivet/Projections/FastJets.hh" #include /// @todo Include more projections as required, e.g. ChargedFinalState, FastJets, ZFinder... namespace Rivet { class ATLAS_2013_I1240088 : public Analysis { public: /// @name Constructors etc. //@{ /// Constructor ATLAS_2013_I1240088() : Analysis("ATLAS_2013_I1240088") { - normalise = false; } //@} public: /// @name Analysis methods //@{ /// Book histograms and initialise projections before the run void init() { FinalState fs(-4.0, 4.0, 0.*GeV); addProjection(fs, "FS"); FastJets fj02(fs, FastJets::ANTIKT, 0.2); fj02.useInvisibles(); addProjection(fj02, "Jets02"); - _centedges += 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6; + //_centedges += 0.05, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6; + _centedges += 0.05, 0.1; _ptedges += 45.0, 60.0, 80.0, 110.0, 160.0; _ncentbins = _centedges.size()-1; _nptbins = _ptedges.size()-1; for (size_t i = 0; i < _ncentbins; ++i) { for (size_t j = 0; j < _nptbins; ++j) { + normalise[i][j] = false; _histos[i][j] = bookHisto1D(i+1, 2, j+1); } } } /// Perform the per-event analysis void analyze(const Event& event) { const double weight = event.weight(); const double cent = (event.genEvent()->heavy_ion()?event.genEvent()->heavy_ion()->impact_parameter():-1.); Cut cuts = Cuts::abseta < 2.1 && Cuts::pT > _ptedges.front()*GeV && Cuts::pT < _ptedges.back()*GeV; Jets alljets = applyProjection(event, "Jets02").jetsByPt(cuts); for (size_t i = 0; i < _ncentbins; ++i) { if (cent < 0. || (cent >= _centedges[i] && cent < _centedges[i+1])) { for (size_t j = 0; j < _nptbins; ++j) { foreach (Jet jet, alljets) { double pt(jet.momentum().pT()); if (pt > _ptedges[j]*GeV && pt < _ptedges[j+1]*GeV) { - normalise = true; + normalise[i][j] = true; double phi(jet.phi()); if (phi > M_PI/2. && phi < M_PI) phi = M_PI - phi; if (phi > M_PI && phi < 3.*M_PI/2.) phi = phi - M_PI; if (phi > 3.*M_PI/2.) phi = 2.*M_PI - phi; _histos[i][j]->fill(phi,weight); } } } } } } /// Normalise histograms etc., after the run void finalize() { - if (!normalise) return; + //if (!normalise) return; for (size_t i = 0; i < _ncentbins; ++i) { for (size_t j = 0; j < _nptbins; ++j) { - normalize(_histos[i][j]); + if (normalise[i][j]) normalize(_histos[i][j]); } } } //@} private: // Data members like post-cuts event weight counters go here vector _ptedges; vector _centedges; size_t _ncentbins, _nptbins; - bool normalise; + bool normalise[1][4]; private: /// @name Histograms //@{ - Histo1DPtr _histos[6][4]; + Histo1DPtr _histos[1][4]; //@} }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(ATLAS_2013_I1240088); }