From 6b079439910f7b5ad315a6963f2995648a9cbdb4 Mon Sep 17 00:00:00 2001 From: cholm Date: Thu, 19 Sep 2013 12:53:22 +0000 Subject: [PATCH] Small fixes --- .../FORWARD/analysis2/scripts/OtherPNchData.C | 529 +++++++++--------- 1 file changed, 271 insertions(+), 258 deletions(-) diff --git a/PWGLF/FORWARD/analysis2/scripts/OtherPNchData.C b/PWGLF/FORWARD/analysis2/scripts/OtherPNchData.C index 51a1f3d8f04..69c184cb7dc 100644 --- a/PWGLF/FORWARD/analysis2/scripts/OtherPNchData.C +++ b/PWGLF/FORWARD/analysis2/scripts/OtherPNchData.C @@ -14,6 +14,7 @@ struct OtherPNch */ TGraphAsymmErrors* CMSsqrts900eta05() { + //printf("CMS sqrt(s)=900GeV |eta|<05"); double x[] = { 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, @@ -106,6 +107,7 @@ struct OtherPNch */ TGraphAsymmErrors* CMSsqrts900eta10() { + //printf("CMS sqrt(s)=900GeV |eta|<10"); double x[] = { 0.0, 1.0, 2.0, 3.0, @@ -247,58 +249,59 @@ struct OtherPNch */ TGraphAsymmErrors* CMSsqrts900eta15() { - double x[] = { 0.0, 1.0, - 2.0, 3.0, - 4.0, 5.0, - 6.0, 7.0, - 8.0, 9.0, - 10.0, 11.0, - 12.0, 13.0, - 14.0, 15.0, - 16.0, 17.0, - 18.0, 19.0, - 20.0, 21.0, - 22.0, 23.0, - 24.0, 25.0, - 26.0, 27.0, - 28.0, 29.0, - 30.0, 31.0, - 32.0, 33.0, - 34.0, 35.0, - 36.0, 37.0, - 38.0, 39.0, - 40.0, 41.0, - 42.0, 43.0, - 44.0, 45.0, - 46.5, 49.0, - 53.0, 58.0, - 63.0, 68.0 }; - double xem[] = { 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 1.0, 1.5, - 2.5, 2.5, - 2.5, 2.5 }; + //printf("CMS sqrt(s)=900GeV |eta|<15"); + double x[] = { 0.0, 1.0, // 2 + 2.0, 3.0, // 4 + 4.0, 5.0, // 6 + 6.0, 7.0, // 8 + 8.0, 9.0, // 10 + 10.0, 11.0, // 12 + 12.0, 13.0, // 14 + 14.0, 15.0, // 16 + 16.0, 17.0, // 18 + 18.0, 19.0, // 20 + 20.0, 21.0, // 22 + 22.0, 23.0, // 24 + 24.0, 25.0, // 26 + 26.0, 27.0, // 28 + 28.0, 29.0, // 30 + 30.0, 31.0, // 32 + 32.0, 33.0, // 34 + 34.0, 35.0, // 36 + 36.0, 37.0, // 38 + 38.0, 39.0, // 40 + 40.0, 41.0, // 42 + 42.0, 43.0, // 44 + 44.0, 45.0, // 46 + 46.5, 49.0, // 48 + 53.0, 58.0, // 50 + 63.0, 68.0 }; // 52 + double xem[] = { 0.5, 0.5, // 2 + 0.5, 0.5, // 4 + 0.5, 0.5, // 6 + 0.5, 0.5, // 8 + 0.5, 0.5, // 10 + 0.5, 0.5, // 12 + 0.5, 0.5, // 14 + 0.5, 0.5, // 16 + 0.5, 0.5, // 18 + 0.5, 0.5, // 20 + 0.5, 0.5, // 22 + 0.5, 0.5, // 24 + 0.5, 0.5, // 26 + 0.5, 0.5, // 28 + 0.5, 0.5, // 30 + 0.5, 0.5, // 32 + 0.5, 0.5, // 34 + 0.5, 0.5, // 36 + 0.5, 0.5, // 38 + 0.5, 0.5, // 40 + 0.5, 0.5, // 42 + 0.5, 0.5, // 44 + 0.5, 0.5, // 46 + 1.0, 1.5, // 48 + 2.5, 2.5, // 50 + 2.5, 2.5 }; // 52 double xep[] = { 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, @@ -415,212 +418,6 @@ struct OtherPNch g->SetTitle("/HepData/8068/d4x1y1"); return g; } - /** - * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\lt2.4@f$ - * http://hepdata.cedar.ac.uk/view/p8068/d5 - * - * @return p8068_d5x1y1 - */ - TGraphAsymmErrors* CMSsqrts900eta24() - { - double x[] = { 0.0, 1.0, - 2.0, 3.0, - 4.0, 5.0, - 6.0, 7.0, - 8.0, 9.0, - 10.0, 11.0, - 12.0, 13.0, - 14.0, 15.0, - 16.0, 17.0, - 18.0, 19.0, - 20.0, 21.0, - 22.0, 23.0, - 24.0, 25.0, - 26.0, 27.0, - 28.0, 29.0, - 30.0, 31.0, - 32.0, 33.0, - 34.0, 35.0, - 36.0, 37.0, - 38.0, 39.0, - 40.0, 41.0, - 42.0, 43.0, - 44.0, 45.0, - 46.0, 47.0, - 48.0, 49.0, - 50.0, 51.0, - 52.0, 53.0, - 54.0, 55.0, - 56.0, 57.5, - 60.0, 66.5, - 76.5, 86.5 }; - double xem[] = { 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 1.0, - 1.5, 5.0, - 5.0, 5.0 }; - double xep[] = { 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 0.5, - 0.5, 1.0, - 1.5, 5.0, - 5.0, 5.0 }; - double y[] = { 0.0561, 0.02321, - 0.02645, 0.031, - 0.0362, 0.04081, - 0.04402, 0.04548, - 0.04592, 0.04515, - 0.04342, 0.04101, - 0.03836, 0.03573, - 0.03316, 0.03079, - 0.02868, 0.02679, - 0.02506, 0.02343, - 0.02187, 0.02039, - 0.01901, 0.0178, - 0.01668, 0.01556, - 0.01445, 0.01334, - 0.01229, 0.01132, - 0.01045, 0.00967, - 0.00896, 0.00829, - 0.00763, 0.00699, - 0.00637, 0.00578, - 0.00527, 0.00482, - 0.00441, 0.00412, - 0.00363, 0.00328, - 0.00292, 0.00261, - 0.00233, 0.00209, - 0.00185, 0.00164, - 0.00146, 0.00131, - 0.00117, 0.00106, - 9.7E-4, 8.9E-4, - 8.2E-4, 6.63E-4, - 4.53E-4, 2.12E-4, - 7.7E-5, 1.12E-5 }; - double yem[] = { 0.01767766952966369, 0.003660737630587584, - 0.005079094407470687, 0.006069810540700591, - 0.006147235150862541, 0.005458580401532985, - 0.004473712105176193, 0.0035970126494078384, - 0.003456848275524976, 0.003228761372415125, - 0.0029797483115189447, 0.002739361239413305, - 0.0025401181074902798, 0.002313114783144148, - 0.0020573040611440983, 0.001820109886792553, - 0.0016324827717314506, 0.001484318025222358, - 0.0013727709204379296, 0.001287866452703851, - 0.001217538500417954, 0.0011279184367674819, - 0.0010412012293500234, 9.741663102366044E-4, - 9.360555539069249E-4, 9.139474820797965E-4, - 9.013878188659973E-4, 8.697700845625814E-4, - 8.190848551890091E-4, 7.433034373659253E-4, - 6.741661516273269E-4, 6.239390995922599E-4, - 5.961543424315552E-4, 5.8309518948453E-4, - 5.738466694161429E-4, 5.515432893255072E-4, - 5.292447448959697E-4, 4.976946855251721E-4, - 4.609772228646444E-4, 4.295346318982906E-4, - 4.0249223594996216E-4, 4.386342439892262E-4, - 3.712142238654117E-4, 3.6674241641784496E-4, - 3.6249137920783716E-4, 3.5341194094144583E-4, - 3.3105890714493693E-4, 3.130495168499705E-4, - 2.8178005607210744E-4, 2.5553864678361276E-4, - 2.3323807579381201E-4, 2.1633307652783938E-4, - 1.8601075237738274E-4, 1.6401219466856724E-4, - 1.562049935181331E-4, 1.562049935181331E-4, - 1.5480633061990717E-4, 1.395886814895821E-4, - 1.0756393447619885E-4, 6.676076692189808E-5, - 2.9966648127543395E-5, 6.551335741663679E-6 }; - double yep[] = { 0.09323352401362935, 0.004710169848317574, - 0.005736898116578331, 0.006618164398078972, - 0.006715772479767314, 0.006096630216767292, - 0.00519042387479096, 0.004362247127341595, - 0.004172217635742412, 0.003913693907295255, - 0.003644283194264683, 0.0033638668225719043, - 0.0031145144083789367, 0.0028765604460883488, - 0.002629087294100369, 0.0023903974564912843, - 0.002181604913819182, 0.0020026232796010335, - 0.0018514858897652987, 0.0017280335644888382, - 0.00160822883943797, 0.0014787156589419076, - 0.0013612494260788505, 0.0012539936203984452, - 0.0011763077828527702, 0.0011157060544785082, - 0.001084158659975559, 0.001033247308247159, - 9.630160954002794E-4, 9.244457799135652E-4, - 8.35224520712844E-4, 6.989277502002622E-4, - 6.519202405202648E-4, 6.203224967708329E-4, - 6.10982814815605E-4, 5.88727441181401E-4, - 5.664803615307418E-4, 5.348831648126533E-4, - 4.976946855251721E-4, 4.477722635447622E-4, - 4.204759208325728E-4, 3.9357337308308855E-4, - 3.981205847478877E-4, 3.8483762809787716E-4, - 3.807886552931954E-4, 3.716180835212409E-4, - 3.584689665786984E-4, 3.4928498393145964E-4, - 3.1780497164141406E-4, 2.9068883707497264E-4, - 2.505992817228334E-4, 2.2472205054244234E-4, - 2.0248456731316588E-4, 1.8027756377319944E-4, - 1.6401219466856724E-4, 1.6401219466856724E-4, - 1.5640012787718557E-4, 1.4223923509355636E-4, - 1.0288342918079667E-4, 6.86804193347711E-5, - 3.361547262794322E-5, 7.071067811865475E-6 }; - int np = 62; - TGraphAsymmErrors* g = new TGraphAsymmErrors(np, - x, - y, - xem, - xep, - yem, - yep); - g->SetName("/HepData/8068/d5x1y1"); - g->SetTitle("/HepData/8068/d5x1y1"); - return g; - } /** * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\lt2.0@f$ * http://hepdata.cedar.ac.uk/view/p8068/d6 @@ -629,6 +426,7 @@ struct OtherPNch */ TGraphAsymmErrors* CMSsqrts900eta20() { + //printf("CMS sqrt(s)=900GeV |eta|<20"); double x[] = { 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, @@ -845,6 +643,213 @@ struct OtherPNch g->SetTitle("/HepData/8068/d6x1y1"); return g; } + /** + * CMS @f$ P(N_{ch})@f$ at @f$\sqrt{s}=900GeV@f$ for @f$|\eta|\lt2.4@f$ + * http://hepdata.cedar.ac.uk/view/p8068/d5 + * + * @return p8068_d5x1y1 + */ + TGraphAsymmErrors* CMSsqrts900eta24() + { + //printf("CMS sqrt(s)=900GeV |eta|<24"); + double x[] = { 0.0, 1.0, + 2.0, 3.0, + 4.0, 5.0, + 6.0, 7.0, + 8.0, 9.0, + 10.0, 11.0, + 12.0, 13.0, + 14.0, 15.0, + 16.0, 17.0, + 18.0, 19.0, + 20.0, 21.0, + 22.0, 23.0, + 24.0, 25.0, + 26.0, 27.0, + 28.0, 29.0, + 30.0, 31.0, + 32.0, 33.0, + 34.0, 35.0, + 36.0, 37.0, + 38.0, 39.0, + 40.0, 41.0, + 42.0, 43.0, + 44.0, 45.0, + 46.0, 47.0, + 48.0, 49.0, + 50.0, 51.0, + 52.0, 53.0, + 54.0, 55.0, + 56.0, 57.5, + 60.0, 66.5, + 76.5, 86.5 }; + double xem[] = { 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 1.0, + 1.5, 5.0, + 5.0, 5.0 }; + double xep[] = { 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 0.5, + 0.5, 1.0, + 1.5, 5.0, + 5.0, 5.0 }; + double y[] = { 0.0561, 0.02321, + 0.02645, 0.031, + 0.0362, 0.04081, + 0.04402, 0.04548, + 0.04592, 0.04515, + 0.04342, 0.04101, + 0.03836, 0.03573, + 0.03316, 0.03079, + 0.02868, 0.02679, + 0.02506, 0.02343, + 0.02187, 0.02039, + 0.01901, 0.0178, + 0.01668, 0.01556, + 0.01445, 0.01334, + 0.01229, 0.01132, + 0.01045, 0.00967, + 0.00896, 0.00829, + 0.00763, 0.00699, + 0.00637, 0.00578, + 0.00527, 0.00482, + 0.00441, 0.00412, + 0.00363, 0.00328, + 0.00292, 0.00261, + 0.00233, 0.00209, + 0.00185, 0.00164, + 0.00146, 0.00131, + 0.00117, 0.00106, + 9.7E-4, 8.9E-4, + 8.2E-4, 6.63E-4, + 4.53E-4, 2.12E-4, + 7.7E-5, 1.12E-5 }; + double yem[] = { 0.01767766952966369, 0.003660737630587584, + 0.005079094407470687, 0.006069810540700591, + 0.006147235150862541, 0.005458580401532985, + 0.004473712105176193, 0.0035970126494078384, + 0.003456848275524976, 0.003228761372415125, + 0.0029797483115189447, 0.002739361239413305, + 0.0025401181074902798, 0.002313114783144148, + 0.0020573040611440983, 0.001820109886792553, + 0.0016324827717314506, 0.001484318025222358, + 0.0013727709204379296, 0.001287866452703851, + 0.001217538500417954, 0.0011279184367674819, + 0.0010412012293500234, 9.741663102366044E-4, + 9.360555539069249E-4, 9.139474820797965E-4, + 9.013878188659973E-4, 8.697700845625814E-4, + 8.190848551890091E-4, 7.433034373659253E-4, + 6.741661516273269E-4, 6.239390995922599E-4, + 5.961543424315552E-4, 5.8309518948453E-4, + 5.738466694161429E-4, 5.515432893255072E-4, + 5.292447448959697E-4, 4.976946855251721E-4, + 4.609772228646444E-4, 4.295346318982906E-4, + 4.0249223594996216E-4, 4.386342439892262E-4, + 3.712142238654117E-4, 3.6674241641784496E-4, + 3.6249137920783716E-4, 3.5341194094144583E-4, + 3.3105890714493693E-4, 3.130495168499705E-4, + 2.8178005607210744E-4, 2.5553864678361276E-4, + 2.3323807579381201E-4, 2.1633307652783938E-4, + 1.8601075237738274E-4, 1.6401219466856724E-4, + 1.562049935181331E-4, 1.562049935181331E-4, + 1.5480633061990717E-4, 1.395886814895821E-4, + 1.0756393447619885E-4, 6.676076692189808E-5, + 2.9966648127543395E-5, 6.551335741663679E-6 }; + double yep[] = { 0.09323352401362935, 0.004710169848317574, + 0.005736898116578331, 0.006618164398078972, + 0.006715772479767314, 0.006096630216767292, + 0.00519042387479096, 0.004362247127341595, + 0.004172217635742412, 0.003913693907295255, + 0.003644283194264683, 0.0033638668225719043, + 0.0031145144083789367, 0.0028765604460883488, + 0.002629087294100369, 0.0023903974564912843, + 0.002181604913819182, 0.0020026232796010335, + 0.0018514858897652987, 0.0017280335644888382, + 0.00160822883943797, 0.0014787156589419076, + 0.0013612494260788505, 0.0012539936203984452, + 0.0011763077828527702, 0.0011157060544785082, + 0.001084158659975559, 0.001033247308247159, + 9.630160954002794E-4, 9.244457799135652E-4, + 8.35224520712844E-4, 6.989277502002622E-4, + 6.519202405202648E-4, 6.203224967708329E-4, + 6.10982814815605E-4, 5.88727441181401E-4, + 5.664803615307418E-4, 5.348831648126533E-4, + 4.976946855251721E-4, 4.477722635447622E-4, + 4.204759208325728E-4, 3.9357337308308855E-4, + 3.981205847478877E-4, 3.8483762809787716E-4, + 3.807886552931954E-4, 3.716180835212409E-4, + 3.584689665786984E-4, 3.4928498393145964E-4, + 3.1780497164141406E-4, 2.9068883707497264E-4, + 2.505992817228334E-4, 2.2472205054244234E-4, + 2.0248456731316588E-4, 1.8027756377319944E-4, + 1.6401219466856724E-4, 1.6401219466856724E-4, + 1.5640012787718557E-4, 1.4223923509355636E-4, + 1.0288342918079667E-4, 6.86804193347711E-5, + 3.361547262794322E-5, 7.071067811865475E-6 }; + int np = 62; + TGraphAsymmErrors* g = new TGraphAsymmErrors(np, + x, + y, + xem, + xep, + yem, + yep); + g->SetName("/HepData/8068/d5x1y1"); + g->SetTitle("/HepData/8068/d5x1y1"); + return g; + } /** * Get CMS data * @@ -857,6 +862,7 @@ struct OtherPNch { Double_t aEta = TMath::Abs(eta); TGraphAsymmErrors* g = 0; + //printf("CMS %dGeV |eta|<%f ->", sNN, aEta); switch (sNN) { case 900: if (aEta <= 0.51) g = CMSsqrts900eta05(); @@ -883,6 +889,7 @@ struct OtherPNch */ TGraphAsymmErrors* ALICEsqrts900eta05() { + //printf("ALICE sqrt(s)=900GeV |eta|<05"); double x[] = { 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, @@ -975,6 +982,7 @@ struct OtherPNch */ TGraphAsymmErrors* ALICEsqrts900eta10() { + //printf("ALICE sqrt(s)=900GeV |eta|<10"); double x[] = { 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, @@ -1122,6 +1130,7 @@ struct OtherPNch */ TGraphAsymmErrors* ALICEsqrts900eta13() { + //printf("ALICE sqrt(s)=900GeV |eta|<13"); double x[] = { 0.0, 1.0, 2.0, 3.0, 4.0, 5.0, @@ -1290,6 +1299,7 @@ struct OtherPNch { Double_t aEta = TMath::Abs(eta); TGraphAsymmErrors* g = 0; + //printf("CMS %dGeV |eta|<%f ->", sNN, aEta); switch (sNN) { case 900: if (aEta <= 0.51) g = ALICEsqrts900eta05(); @@ -1308,6 +1318,8 @@ struct OtherPNch { OtherPNch d; TGraphAsymmErrors* g = 0; + //printf("type=%d (%s), eta=%f, sqrt(s)=%dGeV ->", + // type, (type==0 ? "CMS" : "ALICE"), eta, sNN); switch (type) { case 0: g = d.GetCMS(eta, sNN); break; case 1: g = d.GetALICE(eta, sNN); break; @@ -1317,6 +1329,7 @@ struct OtherPNch // type, eta, sNN); return g; } + //printf("%p (%s)\n", g, (g ? g->GetName() : "")); return g; } }; -- 2.39.3