PWG0/multVScentPbPb/DrawMult.C

   1 #if !defined(__CINT__) || defined(__MAKECINT__)
   2 #include "TMath.h"
   3 #include "TH1F.h"
   4 #include "TStyle.h"
   5 #include "TCanvas.h"
   6 #include "TLegend.h"
   7 #include "TGraphErrors.h"
   8 #endif
   9
  10 void    GetRatE(double x,double xe, double y,double ye, double &rat, double &rate);
  11 TGraphErrors* CreateGraph(const char* nm, int np, const double *x,const double *xe, const double* y,const double *ye,Bool_t rat=kFALSE,int col=kBlack,int mrk=20,double sz=1.2,int lst=1);
  12
  13 double Npart_Geom[11] =  {383.851,331.921,262.756,188.639,130.721 ,86.5384,54.0822,31.1483,16.0198,7.75875,3.8324};
  14 double Npart_rms_Geom[11] = {16.2332,17.7241,26.9529,22.2198,19.1009,15.9985,13.4516,10.7482,7.93099,4.90297,2.52505};
  15 //
  16 double Npart_Geom_4b[4]     = {396.624,371.028,344.477,319.018};
  17 double Npart_rms_Geom_4b[4] = {8.84252,11.016,12.1732,12.4974};
  18 //
  19
  20 //double v0SelNov21[11] = 0,79,247,577,1185,2155,3565,5527,8203,12167,15073,19889
  21 double Npart_V0CSelNov21[11] ={380.692,328.323,259.018,185.05,128.121,84.9204,52.6411,29.79,15.1951,6.73934,2.95315};
  22 double Npart2_V0CSelNov21[11]={379.112,323.171,253.417,183.054,129.116,86.9088,54.2474,30.2884,15.1553,6.7094,2.95315};
  23 double Npart_rms_V0CSelNov21[11] = { 18.4857,19.5242,27.0422,21.1484,16.5675,12.8486,9.64889,6.81264,4.5665,2.72584,1.2575};
  24
  25 //double spdSelNov21[11] = 0,21,69,165,345,637,1069,1683,2511,3721,4585,6213}
  26 double Npart_spdSelNov21[11] ={378.888,324.787,254.892,183.689,128.448,85.0047,52.3777,28.8656,14.3408,6.70305,3.16337};
  27 double Npart2_spdSelNov21[11]={379.973,328.414,259.244,185.094,128.309,84.9333,52.7138,29.9245,15.3195,6.98915,3.16337};
  28 double Npart_rms_spdSelNov21[11] = { 21.8302, 25.7908,30.7147, 24.5131,19.7898,15.9204,12.1924,8.81696,6.18204,3.78075,1.76665 };
  29
  30 //double spdSel[11] = {0,29,85,191,385,687,1133,1743,2567,3765,4611,6213};
  31 double Npart_spdSel[11] ={379.393,328.651,259.304,185.004,128.265,84.8281,52.7315,29.0842,15.0519,6.876,3.14467};
  32 double Npart2_spdSel[11]={378.367,322.458,248.503,176.048,121.769,80.3259,49.7426,28.1977,14.5535,6.88646,3.19064};
  33 double Npart_rms_spdSel[11] = { 21.8302, 25.7908,30.7147, 24.5131,19.7898,15.9204,12.1924,8.81696,6.18204,3.78075,1.76665 };
  34 //
  35 //double v0Sel[11]  = {0,107,297,659,1305,2301,3747,5715,8361,12307,15153,19889 };
  36 double Npart_V0CSel[11]  ={381.040,328.781,259.314,185.197,128.160,84.8531,52.6495,29.7634,15.0721,6.67109,2.93958};
  37 double Npart2_V0CSel[11] ={375.596,313.697,240.935,170.22,118.96,79.1409,49.3726,28.2845,14.6765,6.62758,2.91953};
  38 double Npart_rms_V0CSel[11] = {18.4857,19.5242,27.0422,21.1484,16.5675,12.8486,9.64889,6.81264,4.5665,2.72584,1.2575};
  39 //
  40 //double v0rSel[11]  = {0,32.5,89.5,198.5,396.5,700.5,1140.5,1744.5,2562.5,3767.5,4647.5,6367.5};
  41 double Npart_V0CRSel[11]  ={379.114,329.655,260.318,185.996,130.454,87.0497,54.7531,31.4739,15.9883,7.28226,3.23659};
  42 double Npart2_V0CRSel[11] ={376.795,319.268,244.685,172.492,120.471,80.0706,49.9794,28.5653,14.7274,7.13181,3.27423};
  43 double Npart_rms_V0CRSel[11] = {18.4857,19.5242,27.0422,21.1484,16.5675,12.8486,9.64889,6.81264,4.5665,2.72584,1.2575};
  44 //
  45 //double tpcSel[11]  = {0,13,37,85,171,307,507,783,1157,1685,2055,2787}; // these are  tracks
  46 double Npart_TPCSel[11]  ={378.667,329.234,259.17,185.129,128.221,84.7075,52.6338,29.7764,15.1032,6.98597,3.22604}; // these are  tracks
  47 double Npart2_TPCSel[11] ={378.513,325.472,251.739,178.55,123.373,81.187,50.2066,28.4005,14.466,6.92418,3.30427};
  48 double Npart_rms_TPCSel[11] = {18.4857,19.5242,27.0422,21.1484,16.5675,12.8486,9.64889,6.81264,4.5665,2.72584,1.2575}; // wrong
  49 //
  50
  51 //double spdSelNov21[11] = 0,21,69,165,345,637,1069,1683,2511,3721,4585,6213}
  52 double Npart_spdSelNov21_4b[4] ={396.624,371.02,344.477,319.018};
  53 double Npart_rms_spdSelNov21_4b[11] = {8.84252,11.016,12.1732,12.4974};
  54 //
  55
  56 //----------------------------------------------------------------------------------
  57 double mltV0Ch11Nov27[] = {1590.09,1272.03,951.40,640.01,419.10,256.36,146.50,75.01,34.15,13.10,3.37};
  58 double mltV0Ch11eNov27[] = {2.29,1.98,1.17,0.94,0.75,0.56,0.41,0.28,0.19,0.11,0.05};
  59
  60 // new baseline with ZDC TDC cleanup
  61 double mltV0Ch11Nov29[]  = {1600.98,1293.63,966.21,649.30,426.09,260.65,149.44,76.23,34.75,13.28,3.38};
  62 double mltV0Ch11eNov29[] = {2.34,1.96,1.17,0.93,0.74,0.55,0.41,0.28,0.19,0.11,0.05};
  63
  64 // new baseline with ZDC TDC cleanup
  65 double mltV0Ch11Nov29_4b[4]  = {1693.65,1512.59,1365.12,1222.52};
  66 double mltV0Ch11Nov29_4be[4] = {3.69,3.04,2.91,2.66};
  67
  68
  69 //ZDC vs V0C, LHC10h_000137161_ZDCV0_24Nov10RS
  70 // result from setting ZDC&V0 to 80% in the same V0 as V01D
  71 ////double mltZDCV0Ch8Nov24RS[11]={1599.53,1278.17,941.59,625.66,402.07,240.92,130.04,56.54,29.69,15.04,8.12};
  72 ////double mltZDCV0Ch8eNov24RS[11]={2.31,1.89,1.13,0.89,0.69,0.51,0.35,0.22,0.16,0.12,0.11};
  73 // result from setting ZDC&V0 to have the same excess at 0 as V01D
  74 double mltZDCV0Ch8Nov24RS[11]={1606.24,1294.40,964.34,651.15,427.81,264.23,149.39,72.52,34.75,18.79,8.92};
  75 double mltZDCV0Ch8eNov24RS[11]={2.34,1.93,1.16,0.92,0.74,0.55,0.40,0.27,0.18,0.14,0.12};
  76
  77 //----------------------------------------------------------------------------------
  78 //V0C, LHC10h8_000137161_V0C
  79 double mltV0Ch8Nov21[11]={1606.14,1293.78,963.88,650.18,427.15,263.20,151.27,77.95,36.07,13.72,4.19};
  80 double mltV0Ch8eNov21[11]={2.34,1.93,1.16,0.92,0.74,0.55,0.40,0.28,0.19,0.11,0.06};
  81
  82
  83 //-----------------------------------
  84 //SPD, LHC10h8_000137161_SPD
  85 double mltSPDh8Nov21[11]={1608.81,1293.03,962.20,647.64,424.17,259.92,149.15,75.83,34.56,12.68,3.40};
  86 double mltSPDh8eNov21[11]={2.31,1.96,1.15,0.91,0.72,0.54,0.40,0.28,0.19,0.11,0.06};
  87
  88
  89 //SPD, LHC10h6_000137161_SPD (off-vertex primaries removed)
  90 double mltSPDh6[11]={1603.99,1298.08,972.26,661.06,441.62,276.34,163.10,86.02,40.75,16.28,4.13};
  91 double mltSPDh6e[11]={2.54,2.11,1.26,1.01,0.80,0.61,0.46,0.33,0.23,0.14,0.06};
  92
  93 //-----------------------------------
  94 //SPD, LHC10h8_000137161_SPD
  95 double mltSPDh8[11]={1613.30,1304.45,977.46,665.32,442.60,277.46,163.26,86.15,40.97,16.36,4.17};
  96 double mltSPDh8e[11]={2.35,1.99,1.18,0.94,0.76,0.57,0.43,0.31,0.21,0.13,0.06};
  97
  98 //-----------------------
  99 //V0C, LHC10h8_000137161_V0C
 100 double mltV0Ch8[11]={1610.21,1304.62,977.99,667.33,444.72,279.27,164.76,87.79,42.34,17.53,4.89};
 101 double mltV0Ch8e[11] = {2.38,1.97,1.18,0.95,0.76,0.57,0.43,0.31,0.21,0.13,0.06};
 102
 103 //-----------------------
 104 //V0CR, LHC10h8_000137161_V0CR
 105 double mltV0CRh8[11]={1611.13,1304.76,994.16,678.31,449.72,282.53,166.06,87.97,42.40,17.49,4.99};
 106 double mltV0CRh8e[11]={2.38,1.98,1.19,0.97,0.77,0.58,0.43,0.31,0.21,0.13,0.06};
 107
 108 //-----------------------
 109 //TPC, LHC10h8_000137161_TPC
 110 double mltTPCh8[11]={1616.65,1304.32,957.72,656.67,435.59,273.39,160.58,86.09,40.96,16.58,4.16};
 111 double mltTPCh8e[11]={2.44,2.01,1.15,0.92,0.74,0.56,0.42,0.31,0.21,0.13,0.06};
 112
 113 //--------------------------
 114 //137366, SPD, LHC10h8_000137366_SPD
 115 double mlt366SPDh8[11]={1574.05,1278.00,955.98,651.66,432.23,272.50,158.64,84.24,39.71,16.16,4.01};
 116 double mlt366SPDh8e[11]={1.96,1.69,0.98,0.81,0.62,0.48,0.36,0.25,0.17,0.10,0.05};
 117
 118 //--------------------------
 119 //SPD, LHC10h2_000137161_SPD  DPMJET
 120 double  mltSPDdpmjh2[11]={1624.52,1304.84,976.87,664.33,443.03,277.41,163.68,85.95,40.76,16.58,4.25};
 121 double  mltSPDdpmjh2e[11]={2.12,2.46,1.53,1.29,1.06,0.81,0.62,0.43,0.29,0.19,0.09};
 122
 123 //-------------------------------------------------------------------------
 124 // top 4x2.5% bins in SPD
 125 double mltSPDh8_4b[4] = {1702.89,1516.20,1362.55,1221.55}; // w.means: 1587.4, 1259.0
 126 double mltSPDh8_4be[4] = {3.55,3.06,2.91,2.65};
 127 //double mltSPDh8_4b[4] = {1681.25,1480.93,1322.75,1188.11}; // w.means: 1587.4, 1259.0
 128 //double mltSPDh8_4be[4] = {3.51,3.01,2.83,2.59};
 129
 130 //-------------------------------------------------------------------------
 131 // top 4x2.5% bins in V0
 132 double mltV0Ch8_4b[4]  = {1692.60,1514.45,1363.53,1221.21};
 133 double mltV0Ch8_4be[4] = {3.69,3.02,2.85,2.61};
 134 //double mltSPDh8_4b[4] = {1681.25,1480.93,1322.75,1188.11}; // w.means: 1587.4, 1259.0
 135 //double mltSPDh8_4be[4] = {3.51,3.01,2.83,2.59};
 136
 137
 138
 139 TGraphErrors *MV0Ch11Nov27A=0,*RV0Ch11Nov27A=0;
 140
 141 TGraphErrors *MV0Ch11Nov29=0,*RV0Ch11Nov29=0;
 142
 143 TGraphErrors *MzdcV0CNov24=0,*RzdcV0CNov24=0;
 144 //
 145 TGraphErrors *Mspdh8Nov21=0,*Mspdh8aNov21=0,*Rspdh8Nov21=0,*Rspdh8aNov21=0;
 146 TGraphErrors *MV0Ch8Nov21=0,*MV0Ch8aNov21=0,*RV0Ch8Nov21=0,*RV0Ch8aNov21=0;
 147
 148 TGraphErrors *Mspdh6=0,*Mspdh6a=0,*Rspdh6=0,*Rspdh6a=0;
 149 TGraphErrors *Mspdh8=0,*Mspdh8a=0,*Rspdh8=0,*Rspdh8a=0;
 150 TGraphErrors *MV0Ch8=0,*MV0Ch8a=0,*RV0Ch8=0,*RV0Ch8a=0;
 151 TGraphErrors *MV0CRh8=0,*MV0CRh8a=0,*RV0CRh8=0,*RV0CRh8a=0;
 152 TGraphErrors *Mtpch8=0,*Mtpch8a=0,*Rtpch8=0,*Rtpch8a=0;
 153 TGraphErrors *Mspd366h8=0,*Mspd366h8a=0,*Rspd366h8=0,*Rspd366h8a=0;
 154 TGraphErrors *MspdDPMJh2=0,*MspdDPMJh2a=0,*RspdDPMJh2=0,*RspdDPMJh2a=0;
 155
 156 TGraphErrors *Mspd4b=0,*Rspd4b=0;
 157 TGraphErrors *MV0C4b=0,*RV0C4b=0;
 158
 159 TGraphErrors *MV0C29Nov4b=0,*RV0C29Nov4b=0;
 160
 161 TCanvas *cnvMult=0;
 162 TCanvas *cnvMultR=0;
 163 TH1F *frMlt,*frMltR=0;
 164
 165 void DrawMult()
 166 {
 167   gStyle->SetOptStat(0);
 168
 169   /*
 170   Mspd4b       =  CreateGraph("Mspd4b" ,4,Npart_spdSelNov21_4b,Npart_rms_spdSelNov21_4b, mltSPDh8_4b,mltSPDh8_4be ,kFALSE, kGreen, 20,1.2,1);
 171   Rspd4b       =  CreateGraph("Rspd4b" ,4,Npart_spdSelNov21_4b,Npart_rms_spdSelNov21_4b, mltSPDh8_4b,mltSPDh8_4be ,kTRUE,  kGreen, 20,1.2,1);
 172   //
 173   Mspdh8Nov21  = CreateGraph("Mspdh8Nov21" ,11,Npart_spdSelNov21 ,Npart_rms_spdSelNov21,mltSPDh8Nov21,mltSPDh8eNov21, kFALSE, kRed, 20,1.2,1);
 174   Mspdh8aNov21 = CreateGraph("Mspdh8aNov21",11,Npart2_spdSelNov21,Npart_rms_spdSelNov21,mltSPDh8Nov21,mltSPDh8eNov21, kFALSE, kRed, 24,1.42,2);
 175   Rspdh8Nov21  = CreateGraph("Rspdh8Nov21" ,11,Npart_spdSelNov21 ,Npart_rms_spdSelNov21,mltSPDh8Nov21,mltSPDh8eNov21, kTRUE,  kRed, 20,1.2,1);
 176   Rspdh8aNov21 = CreateGraph("Rspdh8aNov21",11,Npart2_spdSelNov21,Npart_rms_spdSelNov21,mltSPDh8Nov21,mltSPDh8eNov21, kTRUE,  kRed, 24,1.42,2);
 177   //
 178   MV0Ch8Nov21  = CreateGraph("MV0Ch8Nov21" ,11,Npart_V0CSelNov21 ,Npart_rms_V0CSelNov21,mltV0Ch8Nov21,mltV0Ch8eNov21, kFALSE, kBlue, 20,1.2, 1);
 179   MV0Ch8aNov21 = CreateGraph("MV0Ch8aNov21",11,Npart2_V0CSelNov21,Npart_rms_V0CSelNov21,mltV0Ch8Nov21,mltV0Ch8eNov21, kFALSE, kBlue, 24,1.42, 2);
 180   RV0Ch8Nov21  = CreateGraph("RV0Ch8Nov21" ,11,Npart_V0CSelNov21 ,Npart_rms_V0CSelNov21,mltV0Ch8Nov21,mltV0Ch8eNov21, kTRUE,  kBlue, 20,1.2, 1);
 181   RV0Ch8aNov21 = CreateGraph("RV0Ch8aNov21",11,Npart2_V0CSelNov21,Npart_rms_V0CSelNov21,mltV0Ch8Nov21,mltV0Ch8eNov21, kTRUE,  kBlue, 24,1.42, 2);
 182   //
 183   Mtpch8  = CreateGraph("Mtpch8" ,11,Npart_TPCSel ,Npart_rms_TPCSel,mltTPCh8,mltTPCh8e, kFALSE, kBlack, 23, 1.2, 1);
 184   Mtpch8a = CreateGraph("Mtpch8a",11,Npart2_TPCSel,Npart_rms_TPCSel,mltTPCh8,mltTPCh8e, kFALSE, kBlack, 32, 1.2, 2);
 185   Rtpch8  = CreateGraph("Rtpch8" ,11,Npart_TPCSel ,Npart_rms_TPCSel,mltTPCh8,mltTPCh8e, kTRUE,  kBlack, 23, 1.2, 1);
 186   Rtpch8a = CreateGraph("Rtpch8a",11,Npart2_TPCSel,Npart_rms_TPCSel,mltTPCh8,mltTPCh8e, kTRUE,  kBlack, 32, 1.2, 2);
 187   //  //
 188   */
 189   Mspd4b       =  CreateGraph("Mspd4b" ,4,Npart_Geom_4b,Npart_rms_Geom_4b, mltSPDh8_4b,mltSPDh8_4be ,kFALSE, kGreen, 20,1.2,1);
 190   Rspd4b       =  CreateGraph("Rspd4b" ,4,Npart_Geom_4b,Npart_rms_Geom_4b, mltSPDh8_4b,mltSPDh8_4be ,kTRUE,  kGreen, 20,1.2,1);
 191   //
 192   MV0C4b       =  CreateGraph("MV0C4b" ,4,Npart_Geom_4b,Npart_rms_Geom_4b, mltV0Ch8_4b,mltV0Ch8_4be ,kFALSE, kMagenta, 21,1.2,1);
 193   RV0C4b       =  CreateGraph("RV0C4b" ,4,Npart_Geom_4b,Npart_rms_Geom_4b, mltV0Ch8_4b,mltV0Ch8_4be ,kTRUE,  kMagenta, 21,1.2,1);
 194   //
 195   MV0C29Nov4b  =  CreateGraph("MV0C29Nov4b" ,4,Npart_Geom_4b,Npart_rms_Geom_4b, mltV0Ch11Nov29_4b,mltV0Ch11Nov29_4be ,kFALSE, kGreen+1, 26,1.2,1);
 196   RV0C29Nov4b  =  CreateGraph("RV0C29Nov4b" ,4,Npart_Geom_4b,Npart_rms_Geom_4b, mltV0Ch11Nov29_4b,mltV0Ch11Nov29_4be ,kTRUE,  kGreen+1, 26,1.2,1);
 197   //
 198   MzdcV0CNov24 = CreateGraph("MzdcV0CNov24" ,11,Npart_Geom ,Npart_rms_Geom,mltZDCV0Ch8Nov24RS,mltZDCV0Ch8eNov24RS, kFALSE, kBlack, 25,1.4,1);
 199   RzdcV0CNov24 = CreateGraph("RzdcV0CNov24" ,11,Npart_Geom ,Npart_rms_Geom,mltZDCV0Ch8Nov24RS,mltZDCV0Ch8eNov24RS, kTRUE,  kBlack, 25,1.4,1);
 200   //
 201   Mspdh8Nov21  = CreateGraph("Mspdh8Nov21" ,11,Npart_Geom ,Npart_rms_Geom,mltSPDh8Nov21,mltSPDh8eNov21, kFALSE, kRed, 20,1.4,1);
 202   Rspdh8Nov21  = CreateGraph("Rspdh8Nov21" ,11,Npart_Geom ,Npart_rms_Geom,mltSPDh8Nov21,mltSPDh8eNov21, kTRUE,  kRed, 20,1.4,1);
 203   //
 204   MV0Ch8Nov21  = CreateGraph("MV0Ch8Nov21" ,11,Npart_Geom ,Npart_rms_Geom,mltV0Ch8Nov21,mltV0Ch8eNov21, kFALSE, kBlue, 20,1.2, 1);
 205   RV0Ch8Nov21  = CreateGraph("RV0Ch8Nov21" ,11,Npart_Geom ,Npart_rms_Geom,mltV0Ch8Nov21,mltV0Ch8eNov21, kTRUE,  kBlue, 20,1.2, 1);
 206   //
 207   //
 208   MV0Ch11Nov27A  = CreateGraph("MV0Ch8Nov27A" ,11,Npart_Geom ,Npart_rms_Geom,mltV0Ch11Nov27,mltV0Ch11eNov27, kFALSE, kCyan, 20,1.2, 1);
 209   RV0Ch11Nov27A  = CreateGraph("RV0Ch8Nov27A" ,11,Npart_Geom ,Npart_rms_Geom,mltV0Ch11Nov27,mltV0Ch11eNov27, kTRUE,  kCyan, 20,1.2, 1);
 210   //
 211   MV0Ch11Nov29  = CreateGraph("MV0Ch8Nov29" ,11,Npart_Geom ,Npart_rms_Geom,mltV0Ch11Nov29,mltV0Ch11eNov29, kFALSE, kGreen+1, 22,1., 1);
 212   RV0Ch11Nov29  = CreateGraph("RV0Ch8Nov29" ,11,Npart_Geom ,Npart_rms_Geom,mltV0Ch11Nov29,mltV0Ch11eNov29, kTRUE,  kGreen+1, 22,1., 1);
 213   //
 214   Mtpch8  = CreateGraph("Mtpch8" ,11,Npart_Geom ,Npart_rms_Geom,mltTPCh8,mltTPCh8e, kFALSE, kBlack, 23, 1.2, 1);
 215   Rtpch8  = CreateGraph("Rtpch8" ,11,Npart_Geom ,Npart_rms_Geom,mltTPCh8,mltTPCh8e, kTRUE,  kBlack, 23, 1.2, 1);
 216   //
 217   Mspdh6  = CreateGraph("Mspdh6" ,11,Npart_spdSel ,Npart_rms_spdSel,mltSPDh6,mltSPDh6e, kFALSE, kRed, 29, 1.0,1);
 218   Mspdh6a = CreateGraph("Mspdh6a",11,Npart2_spdSel,Npart_rms_spdSel,mltSPDh6,mltSPDh6e, kFALSE, kRed, 30, 1.0,2);
 219   Rspdh6  = CreateGraph("Rspdh6" ,11,Npart_spdSel ,Npart_rms_spdSel,mltSPDh6,mltSPDh6e, kTRUE,  kRed, 29, 1.0,1);
 220   Rspdh6a = CreateGraph("Rspdh6a",11,Npart2_spdSel,Npart_rms_spdSel,mltSPDh6,mltSPDh6e, kTRUE,  kRed, 30, 1.0,2);
 221   //
 222   Mspdh8  = CreateGraph("Mspdh8" ,11,Npart_spdSel ,Npart_rms_spdSel,mltSPDh8,mltSPDh8e, kFALSE, kMagenta, 20,1.02,1);
 223   Mspdh8a = CreateGraph("Mspdh8a",11,Npart2_spdSel,Npart_rms_spdSel,mltSPDh8,mltSPDh8e, kFALSE, kMagenta, 24,1.02,2);
 224   Rspdh8  = CreateGraph("Rspdh8" ,11,Npart_spdSel ,Npart_rms_spdSel,mltSPDh8,mltSPDh8e, kTRUE,  kMagenta, 20,1.02,1);
 225   Rspdh8a = CreateGraph("Rspdh8a",11,Npart2_spdSel,Npart_rms_spdSel,mltSPDh8,mltSPDh8e, kTRUE,  kMagenta, 24,1.02,2);
 226   //
 227   MV0Ch8  = CreateGraph("MV0Ch8" ,11,Npart_V0CSel ,Npart_rms_V0CSel,mltV0Ch8,mltV0Ch8e, kFALSE, kCyan, 21,1.02, 1);
 228   MV0Ch8a = CreateGraph("MV0Ch8a",11,Npart2_V0CSel,Npart_rms_V0CSel,mltV0Ch8,mltV0Ch8e, kFALSE, kCyan, 25,1.02, 2);
 229   RV0Ch8  = CreateGraph("RV0Ch8" ,11,Npart_V0CSel ,Npart_rms_V0CSel,mltV0Ch8,mltV0Ch8e, kTRUE,  kCyan, 21,1.02, 1);
 230   RV0Ch8a = CreateGraph("RV0Ch8a",11,Npart2_V0CSel,Npart_rms_V0CSel,mltV0Ch8,mltV0Ch8e, kTRUE,  kCyan, 25,1.02, 2);
 231   //
 232   MV0CRh8  = CreateGraph("MV0CRh8" ,11,Npart_V0CRSel ,Npart_rms_V0CRSel,mltV0CRh8,mltV0CRh8e, kFALSE, kGreen, 22,1.2, 1);
 233   MV0CRh8a = CreateGraph("MV0CRh8a",11,Npart2_V0CRSel,Npart_rms_V0CRSel,mltV0CRh8,mltV0CRh8e, kFALSE, kGreen, 26,1.2, 2);
 234   RV0CRh8  = CreateGraph("RV0CRh8" ,11,Npart_V0CRSel ,Npart_rms_V0CRSel,mltV0CRh8,mltV0CRh8e, kTRUE,  kGreen, 22,1.2, 1);
 235   RV0CRh8a = CreateGraph("RV0CRh8a",11,Npart2_V0CRSel,Npart_rms_V0CRSel,mltV0CRh8,mltV0CRh8e, kTRUE,  kGreen, 26,1.2, 2);
 236   //
 237   Mspd366h8  = CreateGraph("Mspd366h8" ,11,Npart_spdSel ,Npart_rms_spdSel,mlt366SPDh8,mlt366SPDh8e, kFALSE, kRed, 27, 1.0, 1);
 238   Mspd366h8a = CreateGraph("Mspd366h8a",11,Npart2_spdSel,Npart_rms_spdSel,mlt366SPDh8,mlt366SPDh8e, kFALSE, kRed, 33, 1.0, 2);
 239   Rspd366h8  = CreateGraph("Rspd366h8" ,11,Npart_spdSel ,Npart_rms_spdSel,mlt366SPDh8,mlt366SPDh8e, kTRUE,  kRed, 27, 1.0, 1);
 240   Rspd366h8a = CreateGraph("Rspd366h8a",11,Npart2_spdSel,Npart_rms_spdSel,mlt366SPDh8,mlt366SPDh8e, kTRUE,  kRed, 33, 1.0, 2);
 241   //
 242   MspdDPMJh2  = CreateGraph("MspdDPMJh2" ,11,Npart_spdSel ,Npart_rms_spdSel,mltSPDdpmjh2,mltSPDdpmjh2e, kFALSE, kRed, 28, 1.0, 1);
 243   MspdDPMJh2a = CreateGraph("MspdDPMJh2a",11,Npart2_spdSel,Npart_rms_spdSel,mltSPDdpmjh2,mltSPDdpmjh2e, kFALSE, kRed, 34, 1.0, 2);
 244   RspdDPMJh2  = CreateGraph("RspdDPMJh2" ,11,Npart_spdSel ,Npart_rms_spdSel,mltSPDdpmjh2,mltSPDdpmjh2e, kTRUE,  kRed, 28, 1.0, 1);
 245   RspdDPMJh2a = CreateGraph("RspdDPMJh2a",11,Npart2_spdSel,Npart_rms_spdSel,mltSPDdpmjh2,mltSPDdpmjh2e, kTRUE,  kRed, 34, 1.0, 2);
 246   //
 247   TLegend *lg = new TLegend(0.4,0.2,0.65,0.5);
 248   lg->SetFillColor(kWhite);
 249   lg->SetHeader("N_{part} Glauber");
 250   //  lg->AddEntry(Rspdh8,"SPD corr","pl");
 251   //  lg->AddEntry(RV0Ch8,"V0 corr","pl");
 252   //  lg->AddEntry(RV0CRh8,"V0 corr resc","pl");
 253   //  lg->AddEntry(Rtpch8,"TPC","pl");
 254   //
 255   //  lg->AddEntry(Rspdh6,"SPD corr, old hijing","pl");
 256   //  lg->AddEntry(Rspd366h8,"SPD corr, 137366","pl");
 257   //  lg->AddEntry(RspdDPMJh2,"SPD corr, dpmjet","pl");
 258   //
 259   //  lg->AddEntry(Rspdh8Nov21,"SPD corr, 21/11/10","pl");
 260   lg->AddEntry(RV0Ch8Nov21,"V0  corr, 21/11/10","pl");
 261   //  lg->AddEntry(RzdcV0CNov24,"ZDC & V0 corr","pl");
 262   //  lg->AddEntry(Rspd4b,"SPD corr, 4#times2.5%","pl");
 263   lg->AddEntry(RV0C4b,"V0 corr, 4#times2.5%","pl");
 264   //  lg->AddEntry(RV0Ch11Nov27A,"V0 corr, 5+ 27/11/10","pl");
 265   lg->AddEntry(RV0Ch11Nov29,"V0 corr, 5+ 29/11/10","pl");
 266   lg->AddEntry(RV0C29Nov4b,"V0 corr, 5+ 29/11/10, 4#times2.5%","pl");
 267   //
 268   TLegend *lga = new TLegend(0.65,0.2,0.8,0.5);
 269   lga->SetFillColor(kWhite);
 270   lga->SetHeader("N_{part} data");
 271   lga->AddEntry(Rspdh8a," ","pl");
 272   lga->AddEntry(RV0Ch8a," ","pl");
 273   //  lga->AddEntry(RV0CRh8a, " ","pl");
 274   //  lga->AddEntry(Rtpch8a," ","pl");
 275   //
 276   //  lga->AddEntry(Rspdh6a," ","pl");
 277   //  lga->AddEntry(Rspd366h8a," ","pl");
 278   //  lga->AddEntry(RspdDPMJh2a," ","pl");
 279   //
 280   lga->AddEntry(Rspdh8aNov21,"SPD corr, 21/11/10","pl");
 281   lga->AddEntry(RV0Ch8aNov21,"V0  corr, 21/11/10","pl");
 282   //
 283
 284
 285   frMlt = new TH1F("dndeta","dN/d#eta",86,1,430);
 286   frMlt->SetMinimum(1e-3);
 287   frMlt->SetMaximum(2000);
 288   frMlt->GetXaxis()->SetTitle("N_{part}");
 289   frMlt->GetYaxis()->SetTitle("dN/d#eta_{|#eta|<0.5}");
 290
 291   frMltR = new TH1F("dndetaperpart","dN/d#eta / part.pair",86,1,430);
 292   frMltR->SetMinimum(1e-3);
 293   frMltR->SetMaximum(10);
 294   frMltR->GetXaxis()->SetTitle("N_{part}");
 295   frMltR->GetYaxis()->SetTitle("dN/d#eta_{|#eta|<0.5} / part.pair");
 296
 297   cnvMult  = new TCanvas("cnvMult","",10,10,1000,700);
 298   cnvMultR = new TCanvas("cnvMultR","",20,20,1000,700);
 299   //
 300   cnvMult->cd();
 301   gPad->Clear();
 302   gPad->SetLeftMargin(0.15);
 303   gPad->SetBottomMargin(0.15);
 304   frMlt->Draw();
 305   frMlt->GetYaxis()->SetTitleOffset(1.6);
 306   frMlt->GetYaxis()->SetLabelSize(0.03);
 307   frMlt->GetXaxis()->SetTitleOffset(1.6);
 308   frMlt->GetXaxis()->SetLabelSize(0.03);
 309   gPad->SetGrid(1,1);
 310   //
 311   //  Mspdh6->Draw("p");
 312   //  Mspdh6a->Draw("p");
 313   //
 314   //  Mspd366h8->Draw("p");
 315   //  Mspd366h8a->Draw("p");
 316   //
 317   //  MspdDPMJh2->Draw("p");
 318   //  MspdDPMJh2a->Draw("p");
 319   //
 320   //  Mspdh8->Draw("p");
 321   //  Mspdh8a->Draw("p");
 322   //
 323   //  Mspd4b->Draw("p");
 324   MV0C4b->Draw("p");
 325   MV0C29Nov4b->Draw("p");
 326   //
 327   //  MV0Ch8->Draw("p");
 328   //  MV0Ch8a->Draw("p");
 329   //
 330   //  MV0CRh8->Draw("p");
 331   //  MV0CRh8a->Draw("p");
 332   //
 333   //  Mtpch8->Draw("p");
 334   //  Mtpch8a->Draw("p");
 335   //
 336   //
 337   //  Mspdh8Nov21->Draw("p");
 338   //  Mspdh8aNov21->Draw("p");
 339   //
 340   MV0Ch8Nov21->Draw("p");
 341   //  MV0Ch8aNov21->Draw("p");
 342   //
 343   //  MzdcV0CNov24->Draw("p");
 344   //
 345
 346   //  MV0Ch11Nov27A->Draw("p");
 347   MV0Ch11Nov29->Draw("p");
 348
 349   lg->Draw();
 350   //  lga->Draw();
 351   //----------------------------------------
 352   //
 353   cnvMultR->cd();
 354   gPad->Clear();
 355   gPad->SetLeftMargin(0.15);
 356   gPad->SetBottomMargin(0.15);
 357   frMltR->Draw();
 358   frMltR->GetYaxis()->SetTitleOffset(1.6);
 359   frMltR->GetYaxis()->SetLabelSize(0.03);
 360   frMltR->GetXaxis()->SetTitleOffset(1.6);
 361   frMltR->GetXaxis()->SetLabelSize(0.03);
 362   gPad->SetGrid(1,1);
 363   //
 364   //  /*
 365   //  Rspdh6->Draw("p");
 366   //  Rspdh6a->Draw("p");
 367   //
 368   //  Rspd366h8->Draw("p");
 369   //  Rspd366h8a->Draw("p");
 370   //
 371   //  RspdDPMJh2->Draw("p");
 372   //  RspdDPMJh2a->Draw("p");
 373   //
 374   //  */
 375   //  Rspdh8->Draw("p");
 376   //  Rspdh8a->Draw("p");
 377   //
 378   //  Rspd4b->Draw("p");
 379   RV0C4b->Draw("p");
 380   //
 381   //  RV0Ch8->Draw("p");
 382   //  RV0Ch8a->Draw("p");
 383   //
 384   //  RV0CRh8->Draw("p");
 385   //  RV0CRh8a->Draw("p");
 386   //
 387   //  Rtpch8->Draw("p");
 388   //  Rtpch8a->Draw("p");
 389   //
 390   //  Rspdh8Nov21->Draw("p");
 391   //  Rspdh8aNov21->Draw("p");
 392   //
 393   RV0Ch8Nov21->Draw("p");
 394   //  RV0Ch8aNov21->Draw("p");
 395   //  RzdcV0CNov24->Draw("p");
 396
 397   //  RV0Ch11Nov27A->Draw("p");
 398   RV0Ch11Nov29->Draw("p");
 399   RV0C29Nov4b->Draw("p");
 400   //
 401   lg->Draw();
 402   //  lga->Draw();
 403   //
 404 }
 405
 406
 407
 408 TGraphErrors* CreateGraph(const char* nm, int np, const double *x,const double *xe, const double* y,const double *ye,Bool_t rat,int col,int mrk,double sz,int lst)
 409 {
 410   TGraphErrors* gr = new TGraphErrors(np);
 411   gr->SetName(nm);
 412   gr->SetTitle(nm);
 413   gr->SetMarkerStyle(mrk);
 414   gr->SetMarkerSize(sz);
 415   gr->SetMarkerColor(col);
 416   gr->SetLineColor(col);
 417   gr->SetLineStyle(lst);
 418   int ip=0;
 419   for (int i=0;i<np;i++) {
 420     double px  = x[i];
 421     double pxe = xe[i];
 422     double py = y[i];
 423     double pye = ye[i];
 424     if (rat) {
 425       double rt,rte;
 426       GetRatE(py,pye, px/2, 0, rt,rte);
 427       py = rt;
 428       pye = rte;
 429     }
 430     gr->SetPoint(ip, px,py);
 431     gr->SetPointError(ip, pxe,pye);
 432     ip++;
 433   }
 434   //
 435   return gr;
 436 }
 437
 438
 439
 440
 441 //____________________________________________________________________________
 442 void GetRatE(double x,double xe, double y,double ye, double &rat, double &rate)
 443 {
 444   rat = 0; rate = 0;
 445   if (TMath::Abs(y)<1e-16 || TMath::Abs(x)<1e-16) return;
 446   rat = x/y;
 447   rate = rat*TMath::Sqrt( xe*xe/(x*x) + ye*ye/(y*y));
 448   //  printf("RAT %e %e / %e %e -> %e %e\n",x,xe,y,ye,rat,rate);
 449 }