|
[Rivet-svn] r2164 - in trunk: . data/anainfo data/plotinfo data/refdata include/Rivet src/Analysesblackhole at projects.hepforge.org blackhole at projects.hepforge.orgThu Dec 10 13:37:27 GMT 2009
Author: hoeth Date: Thu Dec 10 13:37:27 2009 New Revision: 2164 Log: STAR_2006_S6860818: add <pT> vs mass plot, validation Modified: trunk/ChangeLog trunk/data/anainfo/STAR_2006_S6860818.info trunk/data/plotinfo/STAR_2006_S6860818.plot trunk/data/refdata/STAR_2006_S6860818.aida trunk/data/refdata/STAR_2006_S6860818.dat trunk/include/Rivet/ParticleName.hh trunk/src/Analyses/STAR_2006_S6860818.cc Modified: trunk/ChangeLog ============================================================================== --- trunk/ChangeLog Thu Dec 10 13:32:00 2009 (r2163) +++ trunk/ChangeLog Thu Dec 10 13:37:27 2009 (r2164) @@ -2,6 +2,7 @@ * propagate SPECIAL and HISTOGRAM sections from .plot files through compare-histos + * STAR_2006_S6860818: <pT> vs particle mass, validate analysis 2009-12-04 Andy Buckley <andy at insectnation.org> Modified: trunk/data/anainfo/STAR_2006_S6860818.info ============================================================================== --- trunk/data/anainfo/STAR_2006_S6860818.info Thu Dec 10 13:32:00 2009 (r2163) +++ trunk/data/anainfo/STAR_2006_S6860818.info Thu Dec 10 13:37:27 2009 (r2164) @@ -4,7 +4,7 @@ Experiment: STAR Collider: RHIC pp 200 GeV SpiresID: 6860818 -Status: UNVALIDATED, under DEVELOPMENT, NOT WORKING +Status: VALIDATED Authors: - Hendrik Hoeth <hendrik.hoeth at cern.ch> References: @@ -17,4 +17,5 @@ Description: pT distributions of identified strange particles in pp collisions at $\sqrt{s} = 200$ GeV, measured by the STAR experiment at RHIC in - non-single-diffractive minbias events. + non-single-diffractive minbias events. WARNING The <pT> vs. particle + mass plot is not validated yet and might be wrong. Modified: trunk/data/plotinfo/STAR_2006_S6860818.plot ============================================================================== --- trunk/data/plotinfo/STAR_2006_S6860818.plot Thu Dec 10 13:32:00 2009 (r2163) +++ trunk/data/plotinfo/STAR_2006_S6860818.plot Thu Dec 10 13:37:27 2009 (r2164) @@ -66,6 +66,7 @@ # BEGIN PLOT /STAR_2006_S6860818/d02-x02-y01 Title=Ratio of $\bar \Lambda / \Lambda$ as function of $p_\perp$ XLabel=$p_\perp$ [GeV] +YLabel=$\bar \Lambda / \Lambda$ LogY=0 YMin=0 YMax=2 @@ -74,7 +75,38 @@ # BEGIN PLOT /STAR_2006_S6860818/d02-x03-y01 Title=Ratio of $\Xi^+ / \Xi^-$ as function of $p_\perp$ XLabel=$p_\perp$ [GeV] +YLabel=$\Xi^+ / \Xi^-$ LogY=0 YMin=0 YMax=2 # END PLOT + +# BEGIN PLOT /STAR_2006_S6860818/d03-x01-y01 +Title=Mean $p_\perp$ vs particle mass +XLabel=mass [GeV] +YLabel=$\langle p_\perp \rangle$ [GeV] +LogY=0 +XMin=0 +XMax=1.8 +# END PLOT +# BEGIN HISTOGRAM /STAR_2006_S6860818/d03-x01-y01 +ErrorBars=1 +PolyMarker=o +# END HISTOGRAM +# BEGIN SPECIAL /STAR_2006_S6860818/d03-x01-y01 +\rput[B]\physicscoor(0.145, 0.14){\small $\pi^-$} +\rput[B]\physicscoor(0.480, 0.14){\small $K^-$} +\rput[B]\physicscoor(0.510, 0.25){\small $K^0$} +\rput[B]\physicscoor(0.775, 0.14){\small $\rho^0$} +\rput[B]\physicscoor(0.896, 0.14){\small $K^{*0}$} +\rput[B]\physicscoor(0.938, 0.25){\small $\bar p$} +\rput[B]\physicscoor(1.019, 0.14){\small $\phi$} +\rput[B]\physicscoor(1.105, 0.14){\small $\Lambda$} +\rput[B]\physicscoor(1.125, 0.25){\small $\bar \Lambda$} +\rput[B]\physicscoor(1.312, 0.14){\small $\Xi^-$} +\rput[B]\physicscoor(1.332, 0.25){\small $\bar \Xi^+$} +\rput[B]\physicscoor(1.384, 0.14){\small $\Sigma$} +\rput[B]\physicscoor(1.520, 0.14){\small $\Lambda$} +\rput[B]\physicscoor(1.672, 0.14){\small $(\Omega + \bar\Omega)$} +# END SPECIAL + Modified: trunk/data/refdata/STAR_2006_S6860818.aida ============================================================================== --- trunk/data/refdata/STAR_2006_S6860818.aida Thu Dec 10 13:32:00 2009 (r2163) +++ trunk/data/refdata/STAR_2006_S6860818.aida Thu Dec 10 13:37:27 2009 (r2164) @@ -666,4 +666,68 @@ <measurement errorPlus="4.497560e-01" value="1.006650e+00" errorMinus="4.497560e-01"/> </dataPoint> </dataPointSet> + <dataPointSet name="d03-x01-y01" dimension="2" + path="/REF/STAR_2006_S6860818" title="Mean $p_\perp$ vs particle mass"> + <annotation> + <item key="Title" value="Mean $p_\perp$ vs particle mass" sticky="true"/> + <item key="AidaPath" value="/REF/STAR_2006_S6860818/d03-x01-y01" sticky="true"/> + <item key="FullPath" value="/STAR_2006_S6860818.aida/REF/STAR_2006_S6860818/d03-x01-y01" sticky="true"/> + </annotation> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="1.396000e-01" errorMinus="1.000000e-02"/> + <measurement errorPlus="1.800000e-02" value="3.480000e-01" errorMinus="1.800000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="4.856000e-01" errorMinus="1.000000e-02"/> + <measurement errorPlus="7.300000e-02" value="6.050000e-01" errorMinus="7.300000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="5.056000e-01" errorMinus="1.000000e-02"/> + <measurement errorPlus="3.000000e-02" value="6.000000e-01" errorMinus="3.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="7.755000e-01" errorMinus="1.000000e-02"/> + <measurement errorPlus="6.000000e-02" value="6.400000e-01" errorMinus="6.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="8.960000e-01" errorMinus="1.000000e-02"/> + <measurement errorPlus="1.400000e-01" value="8.100000e-01" errorMinus="1.400000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="9.383000e-01" errorMinus="1.000000e-02"/> + <measurement errorPlus="4.100000e-02" value="6.830000e-01" errorMinus="4.100000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.160000e-02" value="1.020600e+00" errorMinus="1.160000e-02"/> + <measurement errorPlus="5.000000e-02" value="8.200000e-01" errorMinus="5.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="1.105000e+00" errorMinus="1.000000e-02"/> + <measurement errorPlus="4.000000e-02" value="7.800000e-01" errorMinus="4.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="1.125000e+00" errorMinus="1.000000e-02"/> + <measurement errorPlus="4.000000e-02" value="7.600000e-01" errorMinus="4.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="1.312000e+00" errorMinus="1.000000e-02"/> + <measurement errorPlus="1.300000e-01" value="9.240000e-01" errorMinus="1.300000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="1.332000e+00" errorMinus="1.000000e-02"/> + <measurement errorPlus="1.300000e-01" value="8.810000e-01" errorMinus="1.300000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="1.384000e+00" errorMinus="1.000000e-02"/> + <measurement errorPlus="8.000000e-02" value="1.020000e+00" errorMinus="8.000000e-02"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="1.520000e+00" errorMinus="1.000000e-02"/> + <measurement errorPlus="1.100000e-01" value="1.080000e+00" errorMinus="1.100000e-01"/> + </dataPoint> + <dataPoint> + <measurement errorPlus="1.000000e-02" value="1.672000e+00" errorMinus="1.000000e-02"/> + <measurement errorPlus="3.000000e-01" value="1.080000e+00" errorMinus="3.000000e-01"/> + </dataPoint> + </dataPointSet> </aida> Modified: trunk/data/refdata/STAR_2006_S6860818.dat ============================================================================== --- trunk/data/refdata/STAR_2006_S6860818.dat Thu Dec 10 13:32:00 2009 (r2163) +++ trunk/data/refdata/STAR_2006_S6860818.dat Thu Dec 10 13:37:27 2009 (r2164) @@ -219,44 +219,49 @@ # END HISTOGRAM -# -------------------------------------------------------------------------- -# -# Figure 11: Mean-Pt vs particle mass -# -# +--------------------------------------------------+ -# | 200 GeV p+p data (red) | -# |--------------------------------------------------| -# |particle name |mean-pT (GeV/c)|error (stat only) | -# |---------------+---------------+------------------| -# |pi- |0.348000 |0.018000 | -# |---------------+---------------+------------------| -# |K- |0.605000 |0.073000 | -# |---------------+---------------+------------------| -# |K0 |0.600000 |0.030000 | -# |---------------+---------------+------------------| -# |rho^{0} |0.640000 |0.060000 | -# |---------------+---------------+------------------| -# |K*{0} |0.810000 |0.140000 | -# |---------------+---------------+------------------| -# |p-bar |0.683000 |0.041000 | -# |---------------+---------------+------------------| -# |phi |0.820000 |0.050000 | -# |---------------+---------------+------------------| -# |lambda |0.780000 |0.040000 | -# |---------------+---------------+------------------| -# |lambda-bar |0.760000 |0.040000 | -# |---------------+---------------+------------------| -# |Xi{-} |0.924000 |0.130000 | -# |---------------+---------------+------------------| -# |Xi-bar{+} |0.881000 |0.130000 | -# |---------------+---------------+------------------| -# |sigma(1385) |1.020000 |0.080000 | -# |---------------+---------------+------------------| -# |lambda(1520) |1.080000 |0.110000 | -# |---------------+---------------+------------------| -# |omega+Omega-bar|1.080000 |0.300000 | -# |--------------------------------------------------| -# ------------------------------------------------------------------------ +## Figure 11: Mean-Pt vs particle mass +## +## +--------------------------------------------------+ +## | 200 GeV p+p data (red) | +## |--------------------------------------------------+----------------------------------+ +## |particle name |mean-pT (GeV/c)|error (stat only) | mass for filling the histograms | +## |---------------+---------------+------------------+----------------------------------+ +## |pi- |0.348000 |0.018000 | 0.1396 | +## |K- |0.605000 |0.073000 | 0.4856 | +## |K0 |0.600000 |0.030000 | 0.5056 | +## |rho^{0} |0.640000 |0.060000 | 0.7755 | +## |K*{0} |0.810000 |0.140000 | 0.8960 | +## |p-bar |0.683000 |0.041000 | 0.9383 | +## |phi |0.820000 |0.050000 | 1.0190 | +## |lambda |0.780000 |0.040000 | 1.1050 | +## |lambda-bar |0.760000 |0.040000 | 1.1250 | +## |Xi{-} |0.924000 |0.130000 | 1.3120 | +## |Xi-bar{+} |0.881000 |0.130000 | 1.3320 | +## |sigma(1385) |1.020000 |0.080000 | 1.3840 | +## |lambda(1520) |1.080000 |0.110000 | 1.5200 | +## |omega+Omega-bar|1.080000 |0.300000 | 1.6720 | +## |-------------------------------------------------------------------------------------| + +# BEGIN HISTOGRAM mean_pt_vs_mass +AidaPath=/REF/STAR_2006_S6860818/d03-x01-y01 +Title=Mean $p_\perp$ vs particle mass +0.1296 0.1496 0.348000 0.018000 +0.4756 0.4956 0.605000 0.073000 +0.4956 0.5156 0.600000 0.030000 +0.7655 0.7855 0.640000 0.060000 +0.8860 0.9060 0.810000 0.140000 +0.9283 0.9483 0.683000 0.041000 +1.0090 1.0322 0.820000 0.050000 +1.0950 1.1150 0.780000 0.040000 +1.1150 1.1350 0.760000 0.040000 +1.3020 1.3220 0.924000 0.130000 +1.3220 1.3420 0.881000 0.130000 +1.3740 1.3940 1.020000 0.080000 +1.5100 1.5300 1.080000 0.110000 +1.6620 1.6820 1.080000 0.300000 +# END HISTOGRAM + + # # Figure 20: Particle yield ratios vs THERMUS model calculations # Modified: trunk/include/Rivet/ParticleName.hh ============================================================================== --- trunk/include/Rivet/ParticleName.hh Thu Dec 10 13:32:00 2009 (r2163) +++ trunk/include/Rivet/ParticleName.hh Thu Dec 10 13:37:27 2009 (r2164) @@ -25,9 +25,7 @@ PIPLUS = 211, PIMINUS = -211, K0L = 130, - K0LBAR = -130, K0S = 310, - K0SBAR = -310, KPLUS = 321, KMINUS = -321, LAMBDA = 3122, Modified: trunk/src/Analyses/STAR_2006_S6860818.cc ============================================================================== --- trunk/src/Analyses/STAR_2006_S6860818.cc Thu Dec 10 13:32:00 2009 (r2163) +++ trunk/src/Analyses/STAR_2006_S6860818.cc Thu Dec 10 13:37:27 2009 (r2164) @@ -33,7 +33,7 @@ addProjection(bbc1, "BBC1"); addProjection(bbc2, "BBC2"); - UnstableFinalState ufs(-2.5, 2.5, 0.2*GeV); + UnstableFinalState ufs(-2.5, 2.5, 0.0*GeV); addProjection(ufs, "UFS"); _h_pT_k0s = bookHistogram1D(1, 1, 1); @@ -45,6 +45,7 @@ _h_pT_xiplus = bookHistogram1D(1, 7, 1); //_h_pT_omega = bookHistogram1D(1, 8, 1); _h_antibaryon_baryon_ratio = bookDataPointSet(2, 1, 1); + _h_pT_vs_mass = bookProfile1D(3, 1, 1); } @@ -65,7 +66,11 @@ const PdgId pid = p.pdgId(); const double pT = p.momentum().pT() / GeV; switch (abs(pid)) { + case PIPLUS: + if (pid < 0) _h_pT_vs_mass->fill(0.1396, pT, weight); + break; case PROTON: + if (pid < 0) _h_pT_vs_mass->fill(0.9383, pT, weight); if (pT > 0.4) { pid > 0 ? _nBaryon[0]++ : _nAntiBaryon[0]++; pid > 0 ? _nWeightedBaryon[0]+=weight : _nWeightedAntiBaryon[0]+=weight; @@ -75,13 +80,34 @@ if (pT > 0.2) { _h_pT_k0s->fill(pT, weight/pT); } + _h_pT_vs_mass->fill(0.5056, pT, weight); + break; + case K0L: + _h_pT_vs_mass->fill(0.5056, pT, weight); + break; + case 113: // rho0(770) + _h_pT_vs_mass->fill(0.7755, pT, weight); + break; + case 313: // K0*(892) + _h_pT_vs_mass->fill(0.8960, pT, weight); + break; + case 333: // phi(1020) + _h_pT_vs_mass->fill(1.0190, pT, weight); + break; + case 3214: // Sigma(1385) + _h_pT_vs_mass->fill(1.3840, pT, weight); + break; + case 3124: // Lambda(1520) + _h_pT_vs_mass->fill(1.5200, pT, weight); break; case KPLUS: + if (pid < 0) _h_pT_vs_mass->fill(0.4856, pT, weight); if (pT > 0.2) { pid > 0 ? _h_pT_kplus->fill(pT, weight/pT) : _h_pT_kminus->fill(pT, weight/pT); } break; case LAMBDA: + pid > 0 ? _h_pT_vs_mass->fill(1.1050, pT, weight) : _h_pT_vs_mass->fill(1.1250, pT, weight); if (pT > 0.3) { pid > 0 ? _h_pT_lambda->fill(pT, weight/pT) : _h_pT_lambdabar->fill(pT, weight/pT); pid > 0 ? _nBaryon[1]++ : _nAntiBaryon[1]++; @@ -89,6 +115,7 @@ } break; case XIMINUS: + pid > 0 ? _h_pT_vs_mass->fill(1.3120, pT, weight) : _h_pT_vs_mass->fill(1.3320, pT, weight); if (pT > 0.5) { pid > 0 ? _h_pT_ximinus->fill(pT, weight/pT) : _h_pT_xiplus->fill(pT, weight/pT); pid > 0 ? _nBaryon[2]++ : _nAntiBaryon[2]++; @@ -96,6 +123,7 @@ } break; case OMEGAMINUS: + _h_pT_vs_mass->fill(1.6720, pT, weight); if (pT > 0.5) { //_h_pT_omega->fill(pT, weight/pT); pid > 0 ? _nBaryon[3]++ : _nAntiBaryon[3]++; @@ -167,8 +195,8 @@ AIDA::IHistogram1D * _h_pT_ximinus; AIDA::IHistogram1D * _h_pT_xiplus; //AIDA::IHistogram1D * _h_pT_omega; - AIDA::IDataPointSet* _h_antibaryon_baryon_ratio; + AIDA::IProfile1D* _h_pT_vs_mass; };
More information about the Rivet-svn mailing list |