PWGCF/FLOW/papers/PRL107_032301/fig3abc.C

   1 // For the lists of output files see the bottom of this macro.
   2
   3 static  int      myDarkRed     = TColor::GetColor(128,0,0);
   4 static  int      myDarkGreen   = TColor::GetColor(0,128,0);
   5 static  int      myDarkBlue    = TColor::GetColor(0,0,128);
   6
   7 void fig3abc(Int_t rWrite = 0, Int_t rPerformance = 0)
   8 {
   9  myOptions();
  10  gROOT->ForceStyle();
  11  TDatime now;
  12  int iDate = now.GetDate();
  13  int iYear=iDate/10000;
  14  int iMonth=(iDate%10000)/100;
  15  int iDay=iDate%100;
  16  char* cMonth[12]={"Jan","Feb","Mar","Apr","May","Jun","Jul","Aug","Sep","Oct","Nov","Dec"};
  17  char cStamp1[25],cStamp2[25];
  18  sprintf(cStamp1,"%i %s %i",iDay,cMonth[iMonth-1],iYear);
  19  sprintf(cStamp2,"%i/%.2d/%i",iDay,iMonth,iYear);
  20
  21     //===================================================================================================================
  22     //  v2{SP}(pt) for 30-40%, rapidity gap = 1.0:
  23     const Int_t nPointsSP_3040ALICE_v2_etaGap10 = 22;
  24     Double_t xSP_3040ALICE_v2_etaGap10[nPointsSP_3040ALICE_v2_etaGap10] = {0.250000,0.350000,0.450000,0.550000,0.650000,
  25         0.750000,0.850000,0.950000,1.100000,1.300000,1.500000,1.700000,1.900000,2.100000,2.300000,2.500000,2.800000,3.200000,
  26         3.600000,4.000000,4.400000,4.800000};
  27     Double_t ySP_3040ALICE_v2_etaGap10[nPointsSP_3040ALICE_v2_etaGap10] = {0.039107,0.053955,0.068788,0.083037,0.095749,
  28         0.107080,0.119547,0.127823,0.143002,0.160203,0.180128,0.191935,0.201014,0.210955,0.220726,0.227580,0.244640,0.230650,
  29         0.234024,0.220566,0.217360,0.208241};
  30     Double_t xErrSP_3040ALICE_v2_etaGap10[nPointsSP_3040ALICE_v2_etaGap10] = {0.};
  31     Double_t yErrSP_3040ALICE_v2_etaGap10[nPointsSP_3040ALICE_v2_etaGap10] = {0.000746,0.000744,0.000790,0.000867,0.000960,
  32         0.001069,0.001187,0.001321,0.001091,0.001339,0.001642,0.002016,0.002483,0.003050,0.003736,0.004558,0.004267,0.006108,
  33         0.008517,0.011491,0.017248,0.017496};
  34     TGraphErrors *GrSP_3040ALICE_v2_etaGap10 = new TGraphErrors(nPointsSP_3040ALICE_v2_etaGap10,xSP_3040ALICE_v2_etaGap10,
  35                                                               ySP_3040ALICE_v2_etaGap10,xErrSP_3040ALICE_v2_etaGap10,
  36                                                               yErrSP_3040ALICE_v2_etaGap10);
  37     myGraphSetUp(GrSP_3040ALICE_v2_etaGap10,1.0,kOpenCircle,kBlue,1,kBlack);
  38
  39     //  v3{SP}(pt) for 30-40%, rapidity gap = 1.0:
  40     const Int_t nPointsSP_3040ALICE_v3_etaGap10 = 18;
  41     Double_t xSP_3040ALICE_v3_etaGap10[nPointsSP_3040ALICE_v3_etaGap10] = {0.250000,0.350000,0.450000,0.550000,0.650000,
  42         0.750000,0.850000,0.950000,1.100000,1.300000,1.500000,1.700000,1.900000,2.100000,2.400000,2.800000,3.500000,4.500000};
  43     Double_t ySP_3040ALICE_v3_etaGap10[nPointsSP_3040ALICE_v3_etaGap10] = {0.009546,0.013907,0.013673,0.021384,0.029748,
  44         0.035769,0.037629,0.044245,0.053656,0.061235,0.063584,0.080152,0.089683,0.089203,0.094393,0.115172,0.107330,0.062049};
  45     Double_t xErrSP_3040ALICE_v3_etaGap10[nPointsSP_3040ALICE_v3_etaGap10] = {0.};
  46     Double_t yErrSP_3040ALICE_v3_etaGap10[nPointsSP_3040ALICE_v3_etaGap10] = {0.001753,0.001759,0.001898,0.002071,0.002291,
  47         0.002545,0.002838,0.003136,0.002585,0.003170,0.003875,0.004736,0.005828,0.007141,0.006755,0.009900,0.010855,0.023521};
  48     TGraphErrors *GrSP_3040ALICE_v3_etaGap10 = new TGraphErrors(nPointsSP_3040ALICE_v3_etaGap10,xSP_3040ALICE_v3_etaGap10,ySP_3040ALICE_v3_etaGap10,
  49                                                               xErrSP_3040ALICE_v3_etaGap10,yErrSP_3040ALICE_v3_etaGap10);
  50     myGraphSetUp(GrSP_3040ALICE_v3_etaGap10,1.2,kOpenTriangleUp,kGreen+2,1,kBlack);
  51
  52     //  v4{SP}(pt) for 30-40%, rapidity gap = 1.0:
  53     const Int_t nPointsSP_3040ALICE_v4_etaGap10 = 12;
  54     Double_t xSP_3040ALICE_v4_etaGap10[nPointsSP_3040ALICE_v4_etaGap10] = {0.300000,0.500000,0.700000,0.900000,1.150000,1.450000,1.750000,2.050000,2.400000,2.800000,3.500000,4.500000};
  55     Double_t ySP_3040ALICE_v4_etaGap10[nPointsSP_3040ALICE_v4_etaGap10] = {0.002125,0.005223,0.017491,0.020758,0.026932,0.043126,0.055596,0.055244,0.081049,0.097779,0.074055,0.158711};
  56     Double_t xErrSP_3040ALICE_v4_etaGap10[nPointsSP_3040ALICE_v4_etaGap10] = {0.};
  57     Double_t yErrSP_3040ALICE_v4_etaGap10[nPointsSP_3040ALICE_v4_etaGap10] = {0.002391,0.002677,0.003227,0.003998,0.004216,0.005692,0.007729,0.010472,0.012836,0.018727,0.020749,0.044065};
  58     TGraphErrors *GrSP_3040ALICE_v4_etaGap10 = new TGraphErrors(nPointsSP_3040ALICE_v4_etaGap10,xSP_3040ALICE_v4_etaGap10,ySP_3040ALICE_v4_etaGap10,
  59                                                               xErrSP_3040ALICE_v4_etaGap10,yErrSP_3040ALICE_v4_etaGap10);
  60     myGraphSetUp(GrSP_3040ALICE_v4_etaGap10,1.,kOpenSquare,kRed,1,kBlack);
  61
  62     //  v5{SP}(pt) for 30-40%, rapidity gap = 1.0:
  63     const Int_t nPointsSP_3040ALICE_v5_etaGap10 = 12;
  64     Double_t xSP_3040ALICE_v5_etaGap10[nPointsSP_3040ALICE_v5_etaGap10] = {0.300000,0.500000,0.700000,0.900000,1.150000,1.450000,1.750000,2.050000,2.400000,2.800000,3.500000,4.500000};
  65     Double_t ySP_3040ALICE_v5_etaGap10[nPointsSP_3040ALICE_v5_etaGap10] = {-0.001660,-0.000953,0.009396,0.020289,-0.003699,0.016010,0.045845,0.044840,0.055518,0.125920,0.014951,0.017336};
  66     Double_t xErrSP_3040ALICE_v5_etaGap10[nPointsSP_3040ALICE_v5_etaGap10] = {0.};
  67     Double_t yErrSP_3040ALICE_v5_etaGap10[nPointsSP_3040ALICE_v5_etaGap10] = {0.006239,0.007340,0.008126,0.010441,0.010926,0.014386,0.020467,0.027239,0.033738,0.051262,0.053235,0.113824};
  68     TGraphErrors *GrSP_3040ALICE_v5_etaGap10 = new TGraphErrors(10,xSP_3040ALICE_v5_etaGap10,ySP_3040ALICE_v5_etaGap10,
  69                                                               xErrSP_3040ALICE_v5_etaGap10,yErrSP_3040ALICE_v5_etaGap10);
  70     myGraphSetUp(GrSP_3040ALICE_v5_etaGap10,1.4,kOpenDiamond,kBlack,1,kBlack);
  71     //===================================================================================================================
  72
  73     //===================================================================================================================
  74     //  v2{SP}(pt) for 30-40%, rapidity gap = 0.2:
  75     const Int_t nPointsSP_3040ALICE_v2_etaGap02 = 23;
  76     Double_t xSP_3040ALICE_v2_etaGap02[nPointsSP_3040ALICE_v2_etaGap02] = {0.250000,0.350000,0.450000,0.550000,0.650000,
  77         0.750000,0.850000,0.950000,1.100000,1.300000,1.500000,1.700000,1.900000,2.100000,2.300000,2.500000,2.700000,2.900000,
  78         3.200000,3.600000,4.000000,4.400000,4.800000};
  79     Double_t ySP_3040ALICE_v2_etaGap02[nPointsSP_3040ALICE_v2_etaGap02] = {0.039506,0.054606,0.068971,0.084048,0.096809,
  80         0.109131,0.119861,0.130735,0.145594,0.161632,0.180863,0.194553,0.208946,0.218347,0.223919,0.230006,0.237168,0.243976,
  81         0.239478,0.237796,0.229257,0.211401,0.199559};
  82     Double_t xErrSP_3040ALICE_v2_etaGap02[nPointsSP_3040ALICE_v2_etaGap02] = {0.};
  83     Double_t yErrSP_3040ALICE_v2_etaGap02[nPointsSP_3040ALICE_v2_etaGap02] = {0.000453,0.000449,0.000471,0.000511,0.000564,
  84         0.000626,0.000697,0.000777,0.000643,0.000799,0.000986,0.001220,0.001512,0.001867,0.002282,0.002804,0.003396,0.004077,
  85         0.003726,0.005269,0.007165,0.010680,0.010646};
  86     TGraphErrors *GrSP_3040ALICE_v2_etaGap02 = new TGraphErrors(nPointsSP_3040ALICE_v2_etaGap02,xSP_3040ALICE_v2_etaGap02,
  87                                                               ySP_3040ALICE_v2_etaGap02,xErrSP_3040ALICE_v2_etaGap02,yErrSP_3040ALICE_v2_etaGap02);
  88     myGraphSetUp(GrSP_3040ALICE_v2_etaGap02,1.,kFullCircle,kBlue,1,kBlack);
  89
  90     //  v3{SP}(pt) for 30-40%, rapidity gap = 0.2:
  91     const Int_t nPointsSP_3040ALICE_v3_etaGap02 = 18;
  92     Double_t xSP_3040ALICE_v3_etaGap02[nPointsSP_3040ALICE_v3_etaGap02] = {0.250000,0.350000,0.450000,0.550000,0.650000,
  93         0.750000,0.850000,0.950000,1.100000,1.300000,1.500000,1.700000,1.900000,2.100000,2.400000,2.800000,3.500000,4.500000};
  94     Double_t ySP_3040ALICE_v3_etaGap02[nPointsSP_3040ALICE_v3_etaGap02] = {0.008459,0.013237,0.017909,0.024466,0.029466,
  95         0.036073,0.040105,0.046120,0.051994,0.064426,0.070241,0.081317,0.092806,0.091925,0.108180,0.114513,0.116033,0.113098};
  96     Double_t xErrSP_3040ALICE_v3_etaGap02[nPointsSP_3040ALICE_v3_etaGap02] = {0.};
  97     Double_t yErrSP_3040ALICE_v3_etaGap02[nPointsSP_3040ALICE_v3_etaGap02] = {0.000868,0.000874,0.000933,0.001020,0.001132,
  98         0.001255,0.001395,0.001548,0.001273,0.001560,0.001911,0.002346,0.002890,0.003561,0.003350,0.004920,0.005420,0.011451};
  99     TGraphErrors *GrSP_3040ALICE_v3_etaGap02 = new TGraphErrors(nPointsSP_3040ALICE_v3_etaGap02,xSP_3040ALICE_v3_etaGap02,ySP_3040ALICE_v3_etaGap02,
 100                                                               xErrSP_3040ALICE_v3_etaGap02,yErrSP_3040ALICE_v3_etaGap02);
 101     myGraphSetUp(GrSP_3040ALICE_v3_etaGap02,1.2,kFullTriangleUp,kGreen+2,1,kBlack);
 102
 103     //  v4{SP}(pt) for 30-40%, rapidity gap = 0.2:
 104     const Int_t nPointsSP_3040ALICE_v4_etaGap02 = 12;
 105     Double_t xSP_3040ALICE_v4_etaGap02[nPointsSP_3040ALICE_v4_etaGap02] = {0.300000,0.500000,0.700000,0.900000,1.150000,1.450000,1.750000,2.050000,2.400000,2.800000,3.500000,4.500000};
 106     Double_t ySP_3040ALICE_v4_etaGap02[nPointsSP_3040ALICE_v4_etaGap02] = {0.002607,0.007077,0.014585,0.020707,0.032150,0.044317,0.052834,0.062102,0.085004,0.087401,0.082642,0.094469};
 107     Double_t xErrSP_3040ALICE_v4_etaGap02[nPointsSP_3040ALICE_v4_etaGap02] = {0.};
 108     Double_t yErrSP_3040ALICE_v4_etaGap02[nPointsSP_3040ALICE_v4_etaGap02] = {0.001134,0.001267,0.001541,0.001899,0.001991,0.002707,0.003668,0.005000,0.006097,0.008925,0.009853,0.020906};
 109     TGraphErrors *GrSP_3040ALICE_v4_etaGap02 = new TGraphErrors(nPointsSP_3040ALICE_v4_etaGap02,xSP_3040ALICE_v4_etaGap02,ySP_3040ALICE_v4_etaGap02,
 110                                                               xErrSP_3040ALICE_v4_etaGap02,yErrSP_3040ALICE_v4_etaGap02);
 111     myGraphSetUp(GrSP_3040ALICE_v4_etaGap02,1.,kFullSquare,kRed,1,kBlack);
 112
 113     //  v5{SP}(pt) for 30-40%, rapidity gap = 0.2:
 114     const Int_t nPointsSP_3040ALICE_v5_etaGap02 = 12;
 115     Double_t xSP_3040ALICE_v5_etaGap02[nPointsSP_3040ALICE_v5_etaGap02] = {0.300000,0.500000,0.700000,0.900000,1.150000,1.450000,1.750000,2.050000,2.400000,2.800000,3.500000,4.500000};
 116     Double_t ySP_3040ALICE_v5_etaGap02[nPointsSP_3040ALICE_v5_etaGap02] = {0.005428,0.002807,0.003034,0.008406,0.011691,0.016794,0.033101,0.025018,0.030943,0.045353,0.071972,0.064883};
 117     Double_t xErrSP_3040ALICE_v5_etaGap02[nPointsSP_3040ALICE_v5_etaGap02] = {0.};
 118     Double_t yErrSP_3040ALICE_v5_etaGap02[nPointsSP_3040ALICE_v5_etaGap02] = {0.002227,0.002588,0.003172,0.003872,0.004061,0.005520,0.007585,0.010182,0.012424,0.018294,0.020202,0.042575};
 119     TGraphErrors *GrSP_3040ALICE_v5_etaGap02 = new TGraphErrors(nPointsSP_3040ALICE_v5_etaGap02,xSP_3040ALICE_v5_etaGap02,ySP_3040ALICE_v5_etaGap02,
 120                                                               xErrSP_3040ALICE_v5_etaGap02,yErrSP_3040ALICE_v5_etaGap02);
 121     myGraphSetUp(GrSP_3040ALICE_v5_etaGap02,1.4,33,kBlack,1,kBlack);
 122     //===================================================================================================================
 123
 124     //===================================================================================================================
 125     // Schenke:
 126     // ideal v3(pt) for 30-40%:
 127     Double_t xIdealSchenke3040_v3[] = {0.007776,0.041058,0.101300,0.189163,0.305635,0.452111,0.630493,0.843352,
 128         1.094190,1.387867,1.731373,2.135325,2.617293,3.210890,4.002174};
 129     Double_t yIdealSchenke3040_v3[] = {1.479520e-04,1.124939e-04,6.584102e-04,3.220716e-03,8.219176e-03,1.547538e-02,2.635454e-02,4.111273e-02,
 130         6.115409e-02,8.484866e-02,1.116099e-01,1.404768e-01,1.692029e-01,1.963983e-01,2.213624e-01};
 131     Double_t xErrIdealSchenke3040_v3[15] = {0.};
 132     Double_t yErrIdealSchenke3040_v3[] = {4.047338e-05,1.796386e-05,5.504791e-05,1.993759e-04,4.913364e-04,8.945152e-04,1.473105e-03,2.234133e-03,
 133         3.251702e-03,4.456732e-03,5.843234e-03,7.403563e-03,9.146037e-03,1.122705e-02,1.408298e-02};
 134     Int_t nPointsIdealSchenke3040_v3 = sizeof(xIdealSchenke3040_v3)/sizeof(Double_t);
 135     TGraphErrors *GrIdealSchenke3040_v3 = new TGraphErrors(nPointsIdealSchenke3040_v3,xIdealSchenke3040_v3,yIdealSchenke3040_v3,
 136                                                          xErrIdealSchenke3040_v3,yErrIdealSchenke3040_v3);
 137     myGraphSetUp(GrIdealSchenke3040_v3,1.,33,kBlack,kDashed,kGreen+2,kGreen-10,1001);
 138
 139     // viscous v3(pt) for 30-40%:
 140     Double_t xViscousSchenke3040_v3[] = {0.007776,0.041058,0.101300,0.189163,0.305635,0.452111,0.630493,0.843352,
 141         1.094190,1.387867,1.731373,2.135325,2.617293,3.210890,4.002174};
 142     Double_t yViscousSchenke3040_v3[] = {1.810421e-04,1.088683e-04,4.874871e-04,2.624849e-03,6.650429e-03,1.214362e-02,2.019854e-02,3.120146e-02,
 143         4.615661e-02,6.383687e-02,8.328145e-02,1.029784e-01,1.200315e-01,1.310700e-01,1.294474e-01};
 144     Double_t xErrViscousSchenke3040_v3[15] = {0.};
 145     Double_t yErrViscousSchenke3040_v3[] = {2.895564e-05,9.863071e-06,4.089573e-05,1.529722e-04,3.555346e-04,6.109573e-04,9.704488e-04,1.454140e-03,
 146         2.115684e-03,2.906832e-03,3.816101e-03,4.841144e-03,6.011145e-03,7.466822e-03,9.715151e-03};
 147     Int_t nPointsViscousSchenke3040_v3 = sizeof(xViscousSchenke3040_v3)/sizeof(Double_t);
 148     TGraphErrors *GrViscousSchenke3040_v3 = new TGraphErrors(nPointsViscousSchenke3040_v3,xViscousSchenke3040_v3,yViscousSchenke3040_v3,
 149                                                            xErrViscousSchenke3040_v3,yErrViscousSchenke3040_v3);
 150     myGraphSetUp(GrViscousSchenke3040_v3,1.,33,kBlack,1,kGreen+2,kGreen-10,1001);
 151     //===================================================================================================================
 152
 153     //===================================================================================================================
 154     // Schenke:
 155     // ideal v2(pt) for 30-40%:
 156     Double_t xIdealSchenke3040_v2[] = {0.007776,0.041058,0.101300,0.189163,0.305635,0.452111,0.630493,0.843352,
 157         1.094190,1.387867,1.731373,2.135325,2.617293,3.210890,4.002174};
 158     Double_t yIdealSchenke3040_v2[] = {2.558414e-05,9.867626e-04,5.661819e-03,1.649046e-02,3.377584e-02,5.560800e-02,8.121079e-02,1.086796e-01,
 159         1.375457e-01,1.667749e-01,1.964467e-01,2.260053e-01,2.560743e-01,2.863471e-01,3.177713e-01};
 160     Double_t xErrIdealSchenke3040_v2[15] = {0.};
 161     Double_t yErrIdealSchenke3040_v2[] = {2.497460e-05,4.656790e-05,2.650216e-04,7.813507e-04,1.647879e-03,2.800632e-03,4.085668e-03,5.417398e-03,
 162         6.708719e-03,7.947069e-03,9.168744e-03,1.043873e-02,1.204135e-02,1.416794e-02,1.715206e-02};
 163     Int_t nPointsIdealSchenke3040_v2 = sizeof(xIdealSchenke3040_v2)/sizeof(Double_t);
 164     TGraphErrors *GrIdealSchenke3040_v2 = new TGraphErrors(nPointsIdealSchenke3040_v2,xIdealSchenke3040_v2,yIdealSchenke3040_v2,
 165                                                          xErrIdealSchenke3040_v2,yErrIdealSchenke3040_v2);
 166     myGraphSetUp(GrIdealSchenke3040_v2,1.,33,kBlack,kDashed,kBlue,kBlue-10,1001);
 167
 168     // viscous v2(pt) for 30-40%:
 169     Double_t xViscousSchenke3040_v2[] = {0.007776,0.041058,0.101300,0.189163,0.305635,0.452111,0.630493,0.843352,
 170         1.094190,1.387867,1.731373,2.135325,2.617293,3.210890,4.002174};
 171     Double_t yViscousSchenke3040_v2[] = {-2.625537e-05,9.441659e-04,5.499670e-03,1.580113e-02,3.170383e-02,5.125177e-02,7.392717e-02,9.807519e-02,
 172         1.229014e-01,1.475803e-01,1.719447e-01,1.952542e-01,2.174093e-01,2.368600e-01,2.511733e-01};
 173     Double_t xErrViscousSchenke3040_v2[15] = {0.};
 174     Double_t yErrViscousSchenke3040_v2[] = {2.033230e-05,3.826003e-05,2.169308e-04,6.176627e-04,1.242795e-03,2.027983e-03,2.914887e-03,3.833102e-03,
 175         4.703648e-03,5.516849e-03,6.311937e-03,7.179270e-03,8.381866e-03,1.020568e-02,1.325437e-02};
 176     Int_t nPointsViscousSchenke3040_v2 = sizeof(xViscousSchenke3040_v2)/sizeof(Double_t);
 177     TGraphErrors *GrViscousSchenke3040_v2 = new TGraphErrors(nPointsViscousSchenke3040_v2,xViscousSchenke3040_v2,yViscousSchenke3040_v2,
 178                                                            xErrViscousSchenke3040_v2,yErrViscousSchenke3040_v2);
 179
 180     myGraphSetUp(GrViscousSchenke3040_v2,1.,33,kBlack,1,kBlue,kBlue-10,1001);
 181     //===================================================================================================================
 182
 183     //===================================================================================================================
 184     // Figure 3b:
 185     //===================================================================================================================
 186     //  v2{SP}(pt) for 00-05%:
 187     const Int_t nPointsSP_0005ALICE_v2_etaGap10 = 17;
 188     Double_t xSP_0005ALICE_v2_etaGap10[nPointsSP_0005ALICE_v2_etaGap10] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.900000,1.100000,
 189         1.300000,1.500000,1.700000,1.900000,2.200000,2.600000,3.000000,3.650000,4.550000};
 190     Double_t ySP_0005ALICE_v2_etaGap10[nPointsSP_0005ALICE_v2_etaGap10] = {0.013672,0.015745,0.019944,0.024169,0.026921,0.030296,0.034916,0.038472,
 191         0.043103,0.046626,0.052386,0.057132,0.060070,0.068419,0.061459,0.056816,0.050311};
 192     Double_t xErrSP_0005ALICE_v2_etaGap10[nPointsSP_0005ALICE_v2_etaGap10] = {0.};
 193     Double_t yErrSP_0005ALICE_v2_etaGap10[nPointsSP_0005ALICE_v2_etaGap10] = {0.000813,0.000804,0.000852,0.000930,0.001029,0.001140,0.000939,0.001155,
 194         0.001410,0.001736,0.002131,0.002620,0.002502,0.003759,0.005615,0.006617,0.014242};
 195     TGraphErrors *GrSP_0005ALICE_v2_etaGap10 = new TGraphErrors(nPointsSP_0005ALICE_v2_etaGap10,xSP_0005ALICE_v2_etaGap10,ySP_0005ALICE_v2_etaGap10,
 196                                                               xErrSP_0005ALICE_v2_etaGap10,yErrSP_0005ALICE_v2_etaGap10);
 197     myGraphSetUp(GrSP_0005ALICE_v2_etaGap10,1.,kOpenCircle,kBlue,1,kBlack,1.,kBlack);
 198
 199     //  v3{SP}(pt) for 00-05%:
 200     const Int_t nPointsSP_0005ALICE_v3_etaGap10 = 16;
 201     Double_t xSP_0005ALICE_v3_etaGap10[nPointsSP_0005ALICE_v3_etaGap10] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.900000,1.100000,
 202         1.300000,1.500000,1.700000,1.900000,2.300000,2.900000,3.650000,4.550000};
 203     Double_t ySP_0005ALICE_v3_etaGap10[nPointsSP_0005ALICE_v3_etaGap10] = {0.006303,0.009800,0.011143,0.014246,0.017628,0.019437,0.028412,0.030580,
 204         0.038730,0.045653,0.052469,0.062303,0.071522,0.082223,0.083373,0.076951};
 205     Double_t xErrSP_0005ALICE_v3_etaGap10[nPointsSP_0005ALICE_v3_etaGap10] = {0.};
 206     Double_t yErrSP_0005ALICE_v3_etaGap10[nPointsSP_0005ALICE_v3_etaGap10] = {0.001012,0.000996,0.001062,0.001158,0.001277,0.001415,0.001158,0.001422,
 207         0.001734,0.002124,0.002610,0.003206,0.002724,0.004977,0.008123,0.017180};
 208     TGraphErrors *GrSP_0005ALICE_v3_etaGap10 = new TGraphErrors(nPointsSP_0005ALICE_v3_etaGap10,xSP_0005ALICE_v3_etaGap10,ySP_0005ALICE_v3_etaGap10,
 209                                                               xErrSP_0005ALICE_v3_etaGap10,yErrSP_0005ALICE_v3_etaGap10);
 210     //ShiftAlongXaxis(SP_0005ALICE_v3_etaGap10,-0.025);
 211     myGraphSetUp(GrSP_0005ALICE_v3_etaGap10,1.2,kOpenTriangleUp,kGreen+2,1.,kBlack);
 212
 213     //  v4{SP}(pt) for 00-05%:
 214     const Int_t nPointsSP_0005ALICE_v4_etaGap10 = 11;
 215     Double_t xSP_0005ALICE_v4_etaGap10[nPointsSP_0005ALICE_v4_etaGap10] = {0.300000,0.500000,0.700000,0.950000,1.250000,1.550000,1.850000,2.300000,2.900000,3.650000,4.550000};
 216     Double_t ySP_0005ALICE_v4_etaGap10[nPointsSP_0005ALICE_v4_etaGap10] = {0.002042,0.002556,0.009693,0.013286,0.016780,0.027865,0.031797,0.051101,0.060164,
 217         0.095985,0.094607};
 218     Double_t xErrSP_0005ALICE_v4_etaGap10[nPointsSP_0005ALICE_v4_etaGap10] = {0.};
 219     Double_t yErrSP_0005ALICE_v4_etaGap10[nPointsSP_0005ALICE_v4_etaGap10] = {0.001460,0.001624,0.001930,0.002021,0.002737,0.003717,0.005042,0.005564,0.010160,
 220         0.016472,0.035083};
 221     TGraphErrors *GrSP_0005ALICE_v4_etaGap10 = new TGraphErrors(nPointsSP_0005ALICE_v4_etaGap10,xSP_0005ALICE_v4_etaGap10,ySP_0005ALICE_v4_etaGap10,
 222                                                               xErrSP_0005ALICE_v4_etaGap10,yErrSP_0005ALICE_v4_etaGap10);
 223     //ShiftAlongXaxis(SP_0005ALICE_v4_etaGap10,0.025);
 224     myGraphSetUp(GrSP_0005ALICE_v4_etaGap10,1.,kOpenSquare,kRed,1.,kBlack);
 225
 226     //  v5{SP}(pt) for 00-05%:
 227     const Int_t nPointsSP_0005ALICE_v5_etaGap10 = 12;
 228     Double_t xSP_0005ALICE_v5_etaGap10[nPointsSP_0005ALICE_v5_etaGap10] = {0.300000,0.500000,0.700000,0.900000,1.100000,1.300000,1.600000,2.000000,2.400000,
 229         2.800000,3.500000,4.500000};
 230     Double_t ySP_0005ALICE_v5_etaGap10[nPointsSP_0005ALICE_v5_etaGap10] = {0.002016,0.003409,0.004029,0.002665,0.002765,0.003042,0.013241,0.015430,0.031845,
 231         0.031373,0.068504,0.017964};
 232     Double_t xErrSP_0005ALICE_v5_etaGap10[nPointsSP_0005ALICE_v5_etaGap10] = {0.};
 233     Double_t yErrSP_0005ALICE_v5_etaGap10[nPointsSP_0005ALICE_v5_etaGap10] = {0.001260,0.001386,0.001696,0.002101,0.002560,0.003119,0.002970,0.004472,0.006802,
 234         0.010073,0.011899,0.027756};
 235     TGraphErrors *GrSP_0005ALICE_v5_etaGap10 = new TGraphErrors(nPointsSP_0005ALICE_v5_etaGap10,xSP_0005ALICE_v5_etaGap10,ySP_0005ALICE_v5_etaGap10,
 236                                                               xErrSP_0005ALICE_v5_etaGap10,yErrSP_0005ALICE_v5_etaGap10);
 237     myGraphSetUp(GrSP_0005ALICE_v5_etaGap10,1.4,kOpenDiamond,kBlack,1.,kBlack);
 238     //===================================================================================================================
 239
 240     //===================================================================================================================
 241     //  v2{SP}(pt) for 00-05%:
 242     const Int_t nPointsSP_0005ALICE_v2_etaGap02 = 26;
 243     Double_t xSP_0005ALICE_v2_etaGap02[nPointsSP_0005ALICE_v2_etaGap02] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.850000,0.950000,
 244         1.050000,1.150000,1.250000,1.350000,1.450000,1.550000,1.650000,1.750000,1.850000,1.950000,
 245         2.100000,2.300000,2.500000,2.700000,2.900000,3.250000,3.750000,4.500000};
 246     Double_t ySP_0005ALICE_v2_etaGap02[nPointsSP_0005ALICE_v2_etaGap02] = {0.012030,0.016538,0.020312,0.024888,0.028227,0.032597,0.035418,0.039010,
 247         0.040737,0.044450,0.045448,0.048067,0.050659,0.055360,0.056007,0.057680,0.059254,0.062111,
 248         0.062703,0.064158,0.066930,0.074260,0.073086,0.073927,0.080837,0.064785};
 249     Double_t xErrSP_0005ALICE_v2_etaGap02[nPointsSP_0005ALICE_v2_etaGap02] = {0.};
 250     Double_t yErrSP_0005ALICE_v2_etaGap02[nPointsSP_0005ALICE_v2_etaGap02] = {0.000436,0.000430,0.000453,0.000492,0.000543,0.000596,0.000668,0.000737,
 251         0.000823,0.000912,0.001015,0.001121,0.001251,0.001386,0.001545,0.001714,0.001907,0.002121,
 252         0.001744,0.002149,0.002632,0.003226,0.003930,0.003436,0.005465,0.006891};
 253     TGraphErrors *GrSP_0005ALICE_v2_etaGap02 = new TGraphErrors(nPointsSP_0005ALICE_v2_etaGap02,xSP_0005ALICE_v2_etaGap02,ySP_0005ALICE_v2_etaGap02,
 254                                                               xErrSP_0005ALICE_v2_etaGap02,yErrSP_0005ALICE_v2_etaGap02);
 255     myGraphSetUp(GrSP_0005ALICE_v2_etaGap02,1.,kFullCircle,kBlue,1.,kBlack);
 256
 257     //  v3{SP}(pt) for 00-05%:
 258     const Int_t nPointsSP_0005ALICE_v3_etaGap02 = 20;
 259     Double_t xSP_0005ALICE_v3_etaGap02[nPointsSP_0005ALICE_v3_etaGap02] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.850000,0.950000,
 260         1.100000,1.300000,1.500000,1.700000,1.900000,2.100000,2.300000,2.600000,3.000000,3.500000,
 261         4.100000,4.700000};
 262     Double_t ySP_0005ALICE_v3_etaGap02[nPointsSP_0005ALICE_v3_etaGap02] = {0.005264,0.008365,0.011069,0.014565,0.015858,0.021911,0.026581,0.029673,
 263         0.034343,0.039329,0.044124,0.054249,0.059637,0.066966,0.076106,0.080674,0.091843,0.089706,
 264         0.092376,0.054630};
 265     Double_t xErrSP_0005ALICE_v3_etaGap02[nPointsSP_0005ALICE_v3_etaGap02] = {0.};
 266     Double_t yErrSP_0005ALICE_v3_etaGap02[nPointsSP_0005ALICE_v3_etaGap02] = {0.000596,0.000590,0.000626,0.000678,0.000751,0.000823,0.000909,0.001005,
 267         0.000831,0.001014,0.001251,0.001539,0.001895,0.002331,0.002882,0.002720,0.004047,0.005268,
 268         0.008874,0.014420};
 269     TGraphErrors *GrSP_0005ALICE_v3_etaGap02 = new TGraphErrors(nPointsSP_0005ALICE_v3_etaGap02,xSP_0005ALICE_v3_etaGap02,ySP_0005ALICE_v3_etaGap02,
 270                                                               xErrSP_0005ALICE_v3_etaGap02,yErrSP_0005ALICE_v3_etaGap02);
 271     myGraphSetUp(GrSP_0005ALICE_v3_etaGap02,1.2,kFullTriangleUp,kGreen+2,1.,kBlack);
 272
 273     //  v4{SP}(pt) for 00-05%:
 274     const Int_t nPointsSP_0005ALICE_v4_etaGap02 = 20;
 275     Double_t xSP_0005ALICE_v4_etaGap02[nPointsSP_0005ALICE_v4_etaGap02] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.850000,0.950000,
 276         1.100000,1.300000,1.500000,1.700000,1.900000,2.100000,2.300000,2.600000,3.000000,3.500000,
 277         4.100000,4.700000};
 278     Double_t ySP_0005ALICE_v4_etaGap02[nPointsSP_0005ALICE_v4_etaGap02] = {0.001858,0.002950,0.003658,0.005777,0.008362,0.012289,0.012218,0.012171,
 279         0.017361,0.020719,0.024185,0.030187,0.037104,0.042695,0.050235,0.056076,0.075308,0.080679,
 280         0.090205,0.077760};
 281     Double_t xErrSP_0005ALICE_v4_etaGap02[nPointsSP_0005ALICE_v4_etaGap02] = {0.};
 282     Double_t yErrSP_0005ALICE_v4_etaGap02[nPointsSP_0005ALICE_v4_etaGap02] = {0.000955,0.000951,0.001013,0.001098,0.001204,0.001326,0.001472,0.001634,
 283         0.001341,0.001641,0.002016,0.002476,0.003059,0.003750,0.004616,0.004365,0.006467,0.008390,
 284         0.014154,0.022657};
 285     TGraphErrors *GrSP_0005ALICE_v4_etaGap02 = new TGraphErrors(nPointsSP_0005ALICE_v4_etaGap02,xSP_0005ALICE_v4_etaGap02,ySP_0005ALICE_v4_etaGap02,
 286                                                               xErrSP_0005ALICE_v4_etaGap02,yErrSP_0005ALICE_v4_etaGap02);
 287     myGraphSetUp(GrSP_0005ALICE_v4_etaGap02,1.,kFullSquare,kRed,1.,kBlack);
 288
 289     //  v5{SP}(pt) for 00-05%:
 290     const Int_t nPointsSP_0005ALICE_v5_etaGap02 = 13;
 291     Double_t xSP_0005ALICE_v5_etaGap02[nPointsSP_0005ALICE_v5_etaGap02] = {0.300000,0.500000,0.700000,0.900000,1.200000,1.600000,2.000000,2.400000,2.800000,
 292         3.300000,3.900000,4.500000,4.900000};
 293     Double_t ySP_0005ALICE_v5_etaGap02[nPointsSP_0005ALICE_v5_etaGap02] = {0.000227,-0.000991,0.006407,0.010032,0.002788,0.008782,0.015454,0.012575,0.044690,
 294         0.061348,0.089657,-0.023651,-0.184366};
 295     Double_t xErrSP_0005ALICE_v5_etaGap02[nPointsSP_0005ALICE_v5_etaGap02] = {0.};
 296     Double_t yErrSP_0005ALICE_v5_etaGap02[nPointsSP_0005ALICE_v5_etaGap02] = {0.001875,0.002131,0.002423,0.003018,0.002881,0.004339,0.006567,0.009819,0.014850,
 297         0.019363,0.034072,0.053514,0.131967};
 298     TGraphErrors *GrSP_0005ALICE_v5_etaGap02 = new TGraphErrors(nPointsSP_0005ALICE_v5_etaGap02,xSP_0005ALICE_v5_etaGap02,ySP_0005ALICE_v5_etaGap02,
 299                                                               xErrSP_0005ALICE_v5_etaGap02,yErrSP_0005ALICE_v5_etaGap02);
 300     myGraphSetUp(GrSP_0005ALICE_v5_etaGap02,1.4,33,kBlack,1.,kBlack);
 301     //===================================================================================================================
 302
 303
 304     //===================================================================================================================
 305     // Figure 3c:
 306  //===================================================================================================================
 307  //  v2{SP}(pt) for 00-02%, eta gap = 0.2:
 308  const Int_t nPointsSP_0002_v2_etaGap02 = 27;
 309  Double_t xSP_0002_v2_etaGap02[nPointsSP_0002_v2_etaGap02] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,
 310 0.850000,0.950000,1.050000,1.150000,1.250000,1.350000,1.450000,1.550000,1.650000,1.750000,1.900000,2.100000,2.300000,
 311 2.500000,2.700000,2.900000,3.200000,3.600000,4.000000,4.400000,4.800000};
 312  Double_t ySP_0002_v2_etaGap02[nPointsSP_0002_v2_etaGap02] = {0.010458,0.013621,0.017307,0.020609,0.024900,0.027688,
 313 0.029414,0.031976,0.034382,0.035864,0.038812,0.040049,0.043550,0.044058,0.045320,0.046360,0.047605,0.053355,0.056907,
 314 0.057136,0.060536,0.058016,0.060921,0.054174,0.042666,0.046099,0.038495};
 315  Double_t xErrSP_0002_v2_etaGap02[nPointsSP_0002_v2_etaGap02] = {0.};
 316  Double_t yErrSP_0002_v2_etaGap02[nPointsSP_0002_v2_etaGap02] = {0.000310,0.000306,0.000323,0.000352,0.000386,0.000428,
 317 0.000476,0.000529,0.000588,0.000653,0.000722,0.000802,0.000892,0.000991,0.001103,0.001224,0.001009,0.001243,0.001524,
 318 0.001876,0.002292,0.002802,0.002624,0.003827,0.005442,0.008369,0.008633};
 319  TGraphErrors *GrSP_0002_v2_etaGap02 = new TGraphErrors(nPointsSP_0002_v2_etaGap02,xSP_0002_v2_etaGap02,ySP_0002_v2_etaGap02,
 320                                                       xErrSP_0002_v2_etaGap02,yErrSP_0002_v2_etaGap02);
 321  myGraphSetUp(GrSP_0002_v2_etaGap02,1.,kFullCircle,kBlue,1.,kBlack);
 322
 323  //  v2{SP}(pt) for 00-02%, eta gap = 1.0:
 324  const Int_t nPointsSP_0002_v2_etaGap10 = 15;
 325  Double_t xSP_0002_v2_etaGap10[nPointsSP_0002_v2_etaGap10] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,
 326 0.900000,1.100000,1.350000,1.650000,1.950000,2.250000,2.700000,3.500000,4.500000};
 327  Double_t ySP_0002_v2_etaGap10[nPointsSP_0002_v2_etaGap10] = {0.010171,0.013190,0.017342,0.020629,0.022617,0.026549,
 328 0.027423,0.032261,0.037467,0.041001,0.045763,0.049327,0.049688,0.051480,0.038527};
 329  Double_t xErrSP_0002_v2_etaGap10[nPointsSP_0002_v2_etaGap10] = {0.};
 330  Double_t yErrSP_0002_v2_etaGap10[nPointsSP_0002_v2_etaGap10] = {0.000600,0.000590,0.000625,0.000683,0.000757,0.000839,
 331 0.000692,0.000848,0.000888,0.001209,0.001653,0.002252,0.002465,0.003968,0.009391};
 332  TGraphErrors *GrSP_0002_v2_etaGap10 = new TGraphErrors(nPointsSP_0002_v2_etaGap10,xSP_0002_v2_etaGap10,ySP_0002_v2_etaGap10,
 333                                                   xErrSP_0002_v2_etaGap10,yErrSP_0002_v2_etaGap10);
 334  myGraphSetUp(GrSP_0002_v2_etaGap10,1.,kOpenCircle,kBlue,1.,kBlack);
 335
 336  //  v3{SP}(pt) for 00-02%, eta gap = 0.2:
 337  const Int_t nPointsSP_0002_v3_etaGap02 = 25;
 338  Double_t xSP_0002_v3_etaGap02[nPointsSP_0002_v3_etaGap02] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,
 339 0.850000,0.950000,1.050000,1.150000,1.250000,1.350000,1.450000,1.550000,1.650000,1.750000,1.900000,2.100000,2.300000,
 340 2.500000,2.800000,3.200000,3.600000,4.000000,4.600000};
 341  Double_t ySP_0002_v3_etaGap02[nPointsSP_0002_v3_etaGap02] = {0.005351,0.007675,0.010356,0.013931,0.017761,0.019780,
 342 0.023915,0.026834,0.030480,0.032369,0.037061,0.037480,0.041826,0.044985,0.048539,0.053263,0.058862,0.063724,0.070478,
 343 0.079161,0.080104,0.088543,0.091927,0.081413,0.075981};
 344  Double_t xErrSP_0002_v3_etaGap02[nPointsSP_0002_v3_etaGap02] = {0.};
 345  Double_t yErrSP_0002_v3_etaGap02[nPointsSP_0002_v3_etaGap02] = {0.000353,0.000350,0.000370,0.000402,0.000441,0.000487,
 346 0.000540,0.000597,0.000663,0.000734,0.000810,0.000899,0.001002,0.001112,0.001234,0.001368,0.001130,0.001391,0.001705,
 347 0.002089,0.001980,0.002936,0.004254,0.006087,0.006718};
 348  TGraphErrors *GrSP_0002_v3_etaGap02 = new TGraphErrors(nPointsSP_0002_v3_etaGap02,xSP_0002_v3_etaGap02,ySP_0002_v3_etaGap02,
 349                                                       xErrSP_0002_v3_etaGap02,yErrSP_0002_v3_etaGap02);
 350  myGraphSetUp(GrSP_0002_v3_etaGap02,1.2,kFullTriangleUp,kGreen+2,1.,kBlack);
 351
 352  //  v3{SP}(pt) for 00-02%, eta gap = 1.0:
 353  const Int_t nPointsSP_0002_v3_etaGap10 = 15;
 354  Double_t xSP_0002_v3_etaGap10[nPointsSP_0002_v3_etaGap10] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,
 355 0.900000,1.100000,1.350000,1.650000,1.950000,2.250000,2.700000,3.500000,4.500000};
 356  Double_t ySP_0002_v3_etaGap10[nPointsSP_0002_v3_etaGap10] = {0.006592,0.007286,0.012180,0.012242,0.017416,0.018393,
 357 0.024716,0.030980,0.037703,0.046558,0.051285,0.064613,0.074831,0.077093,0.082442};
 358  Double_t xErrSP_0002_v3_etaGap10[nPointsSP_0002_v3_etaGap10] = {0.};
 359  Double_t yErrSP_0002_v3_etaGap10[nPointsSP_0002_v3_etaGap10] = {0.000682,0.000676,0.000713,0.000782,0.000860,0.000953,
 360 0.000782,0.000957,0.001002,0.001361,0.001862,0.002541,0.002767,0.004466,0.010586};
 361  TGraphErrors *GrSP_0002_v3_etaGap10 = new TGraphErrors(nPointsSP_0002_v3_etaGap10,xSP_0002_v3_etaGap10,ySP_0002_v3_etaGap10,
 362                                                       xErrSP_0002_v3_etaGap10,yErrSP_0002_v3_etaGap10);
 363  myGraphSetUp(GrSP_0002_v3_etaGap10,1.,kOpenTriangleUp,kGreen+2,1.,kBlack);
 364
 365  //  v4{SP}(pt) for 00-02%, eta gap = 0.2:
 366  const Int_t nPointsSP_0002_v4_etaGap02 = 20;
 367  Double_t xSP_0002_v4_etaGap02[nPointsSP_0002_v4_etaGap02] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,
 368 0.850000,0.950000,1.100000,1.300000,1.500000,1.700000,1.900000,2.200000,2.600000,3.000000,3.400000,3.800000,4.250000,4.750000};
 369  Double_t ySP_0002_v4_etaGap02[nPointsSP_0002_v4_etaGap02] = {0.001803,0.003590,0.003214,0.005683,0.007166,0.007113,
 370 0.010206,0.012110,0.016560,0.019941,0.026784,0.034169,0.036075,0.046846,0.058077,0.060063,0.072744,0.076006,
 371 0.067621,0.044664};
 372  Double_t xErrSP_0002_v4_etaGap02[nPointsSP_0002_v4_etaGap02] = {0.};
 373  Double_t yErrSP_0002_v4_etaGap02[nPointsSP_0002_v4_etaGap02] = {0.000630,0.000623,0.000666,0.000721,0.000794,0.000881,
 374 0.000973,0.001076,0.000882,0.001080,0.001331,0.001640,0.002016,0.001922,0.002879,0.004289,0.006267,0.009064,
 375 0.011705,0.017076};
 376  TGraphErrors *GrSP_0002_v4_etaGap02 = new TGraphErrors(nPointsSP_0002_v4_etaGap02,xSP_0002_v4_etaGap02,ySP_0002_v4_etaGap02,
 377                                                       xErrSP_0002_v4_etaGap02,yErrSP_0002_v4_etaGap02);
 378  myGraphSetUp(GrSP_0002_v4_etaGap02,1.,kFullSquare,kRed,1.,kBlack);
 379
 380  //  v4{SP}(pt) for 00-02%, eta gap = 1.0:
 381  const Int_t nPointsSP_0002_v4_etaGap10 = 15;
 382  Double_t xSP_0002_v4_etaGap10[nPointsSP_0002_v4_etaGap10] = {0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,
 383 0.900000,1.100000,1.350000,1.650000,1.950000,2.250000,2.700000,3.500000,4.500000};
 384  Double_t ySP_0002_v4_etaGap10[nPointsSP_0002_v4_etaGap10] = {-0.000533,0.001167,0.002081,0.005218,0.006826,0.008440,
 385 0.013009,0.014812,0.017125,0.030106,0.038279,0.050488,0.067640,0.071637,0.084239};
 386  Double_t xErrSP_0002_v4_etaGap10[nPointsSP_0002_v4_etaGap10] = {0.};
 387  Double_t yErrSP_0002_v4_etaGap10[nPointsSP_0002_v4_etaGap10] = {0.001427,0.001398,0.001482,0.001594,0.001758,0.001945,
 388 0.001593,0.001951,0.002046,0.002787,0.003802,0.005182,0.005663,0.009064,0.021449};
 389  TGraphErrors *GrSP_0002_v4_etaGap10 = new TGraphErrors(nPointsSP_0002_v4_etaGap10,xSP_0002_v4_etaGap10,ySP_0002_v4_etaGap10,
 390                                                       xErrSP_0002_v4_etaGap10,yErrSP_0002_v4_etaGap10);
 391  myGraphSetUp(GrSP_0002_v4_etaGap10,1.,kOpenSquare,kRed,1.,kBlack);
 392
 393  //  v5{SP}(pt) for 00-02%, eta gap = 0.2:
 394  const Int_t nPointsSP_0002_v5_etaGap02 = 13;
 395  Double_t xSP_0002_v5_etaGap02[nPointsSP_0002_v5_etaGap02] = {0.300000,0.500000,0.700000,0.900000,1.100000,1.300000,1.500000,
 396 1.700000,2.000000,2.550000,3.250000,3.950000,4.650000};
 397  Double_t ySP_0002_v5_etaGap02[nPointsSP_0002_v5_etaGap02] = {0.000570,0.002922,0.002151,0.005256,0.006287,0.005849,0.009399,
 398 0.011420,0.012455,0.032134,0.057009,0.020607,0.013551};
 399  Double_t xErrSP_0002_v5_etaGap02[nPointsSP_0002_v5_etaGap02] = {0.};
 400  Double_t yErrSP_0002_v5_etaGap02[nPointsSP_0002_v5_etaGap02] = {0.001074,0.001155,0.001433,0.001725,0.002123,0.002608,0.003196,
 401 0.003930,0.003755,0.004869,0.009719,0.018353,0.031814};
 402  TGraphErrors *GrSP_0002_v5_etaGap02 = new TGraphErrors(11,xSP_0002_v5_etaGap02,ySP_0002_v5_etaGap02,
 403                                                       xErrSP_0002_v5_etaGap02,yErrSP_0002_v5_etaGap02);
 404  myGraphSetUp(GrSP_0002_v5_etaGap02,1.4,33,kBlack,1.,kBlack);
 405
 406  //  v5{SP}(pt) for 00-02%, eta gap = 1.0:
 407  const Int_t nPointsSP_0002_v5_etaGap10 = 13;
 408  Double_t xSP_0002_v5_etaGap10[nPointsSP_0002_v5_etaGap10] = {0.300000,0.500000,0.700000,0.900000,1.100000,1.300000,1.500000,
 409 1.700000,2.000000,2.550000,3.250000,3.950000,4.650000};
 410  Double_t ySP_0002_v5_etaGap10[nPointsSP_0002_v5_etaGap10] = {0.004709,0.002816,0.004686,0.003099,-0.000160,0.006887,0.007711,
 411 0.014473,0.009113,0.033936,0.056363,-0.013242,0.023228};
 412  Double_t xErrSP_0002_v5_etaGap10[nPointsSP_0002_v5_etaGap10] = {0.};
 413  Double_t yErrSP_0002_v5_etaGap10[nPointsSP_0002_v5_etaGap10] = {0.001725,0.001974,0.002383,0.002966,0.003646,0.004383,0.005409,
 414 0.006631,0.006321,0.008218,0.016372,0.030968,0.053822};
 415  TGraphErrors *GrSP_0002_v5_etaGap10 = new TGraphErrors(11,xSP_0002_v5_etaGap10,ySP_0002_v5_etaGap10,
 416                                                       xErrSP_0002_v5_etaGap10,yErrSP_0002_v5_etaGap10);
 417  myGraphSetUp(GrSP_0002_v5_etaGap10,1.4,kOpenDiamond,kBlack,1.,kBlack);
 418  //===================================================================================================================
 419
 420
 421  // Style histogramm:
 422  TH1F *myBlankHisto = new TH1F("myBlankHisto","Blank Histogram",100,0.,10.);
 423  //myBlankHisto->SetXTitle("p_{t} (GeV/#font[72]{c})");
 424  myBlankHisto->SetYTitle("v_{n}");
 425  myBlankHisto->SetNdivisions(505,"y");
 426  myBlankHisto->GetYaxis()->SetTitleOffset(0.6);
 427  //myBlankHisto->GetYaxis()->SetLabelSize(0.05);
 428  //myBlankHisto->GetYaxis()->SetTitleSize(0.07);
 429  //myBlankHisto->GetYaxis()->SetRangeUser(-0.02101,0.399);
 430  myBlankHisto->GetXaxis()->SetRangeUser(0.,5.);
 431
 432  // Main canvas:
 433  TCanvas *myCan = new TCanvas("myCan",cStamp1,500,1000);
 434  myCan->Divide(1,3,0.,0.);
 435  myCan->SetTopMargin(0.0002);
 436  myCan->Draw();
 437
 438  //=================================================================================================
 439
 440  // Figure 3a:
 441  myCan->cd(1);
 442  gPad->SetTopMargin(0.001);
 443  gPad->SetBottomMargin(0.);
 444  gPad->SetRightMargin(0.001);
 445  TH1D* styleHistPad1 = myBlankHisto->Clone("styleHistPad1");
 446  styleHistPad1->GetYaxis()->SetRangeUser(-0.02101,0.3799);
 447  styleHistPad1->Draw();
 448  //TLatex *system = new TLatex(0.2,0.93,"ALICE Preliminary, Pb-Pb events at  #sqrt{s_{NN}} = 2.76 TeV");
 449  //system->SetNDC();
 450  //  system->SetTextFont(42);
 451  //system->SetTextSize(0.0475);
 452  //system->SetTextColor(myDarkRed);
 453  //system->Draw();
 454  TLatex *a = new TLatex(0.95,0.93,"(a)");
 455  a->SetNDC();
 456  a->SetTextSize(0.05);
 457  a->Draw();
 458  //GrIdealSchenke3040_v3->Draw("same3");
 459  //GrIdealSchenke3040_v3->Draw("lsameX");
 460  GrViscousSchenke3040_v3->Draw("same3");
 461  GrViscousSchenke3040_v3->Draw("lsameX");
 462  GrIdealSchenke3040_v3->Draw("same3");
 463  GrIdealSchenke3040_v3->Draw("lsameX");
 464  //GrIdealSchenke3040_v2->Draw("same3");
 465  //GrIdealSchenke3040_v2->Draw("lsameX");
 466  GrViscousSchenke3040_v2->Draw("same3");
 467  GrViscousSchenke3040_v2->Draw("lsameX");
 468  GrIdealSchenke3040_v2->Draw("same3");
 469  GrIdealSchenke3040_v2->Draw("lsameX");
 470  GrSP_3040ALICE_v2_etaGap02->Draw("psame");
 471  GrSP_3040ALICE_v2_etaGap10->Draw("psame");
 472  GrSP_3040ALICE_v3_etaGap02->Draw("psame");
 473  GrSP_3040ALICE_v3_etaGap10->Draw("psame");
 474  GrSP_3040ALICE_v4_etaGap02->Draw("psame");
 475  GrSP_3040ALICE_v4_etaGap10->Draw("psame");
 476  GrSP_3040ALICE_v5_etaGap02->Draw("psame");
 477  GrSP_3040ALICE_v5_etaGap10->Draw("psame");
 478  // Common legend:
 479  TLegend *myLegendPad1 = new TLegend(0.12,0.50,0.48,0.98,"Centrality 30-40%");
 480  myLegendSetUp(myLegendPad1,0.04);
 481  //myLegendPad1->AddEntry(GrSP_3040ALICE_v2_etaGap10,"v_{2}{SP} (#Delta#eta=1.0)","p");
 482  //myLegendPad1->AddEntry(GrSP_3040ALICE_v2_etaGap02,"v_{2}{SP} (#Delta#eta=0.2)","p");
 483  myLegendPad1->AddEntry(GrSP_3040ALICE_v2_etaGap02,"v_{2}{2}","p");
 484  //myLegendPad1->AddEntry(GrSP_3040ALICE_v3_etaGap10,"v_{3}{SP} (#Delta#eta=1.0)","p");
 485  //myLegendPad1->AddEntry(GrSP_3040ALICE_v3_etaGap02,"v_{3}{SP} (#Delta#eta=0.2)","p");
 486  myLegendPad1->AddEntry(GrSP_3040ALICE_v3_etaGap02,"v_{3}{2}","p");
 487  //myLegendPad1->AddEntry(GrSP_3040ALICE_v4_etaGap10,"v_{4}{SP} (#Delta#eta=1.0)","p");
 488  //myLegendPad1->AddEntry(GrSP_3040ALICE_v4_etaGap02,"v_{4}{SP} (#Delta#eta=0.2)","p");
 489  myLegendPad1->AddEntry(GrSP_3040ALICE_v4_etaGap02,"v_{4}{2}","p");
 490  //myLegendPad1->AddEntry(GrSP_3040ALICE_v5_etaGap10,"v_{5}{SP} (#Delta#eta=1.0)","p");
 491  //myLegendPad1->AddEntry(GrSP_3040ALICE_v5_etaGap02,"v_{5}{SP} (#Delta#eta=0.2)","p");
 492  myLegendPad1->AddEntry(GrSP_3040ALICE_v5_etaGap02,"v_{5}{2}","p");
 493  myLegendPad1->AddEntry(GrIdealSchenke3040_v2,"v_{2} (#eta/s = 0.0)","fl");
 494  myLegendPad1->AddEntry(GrViscousSchenke3040_v2,"v_{2} (#eta/s = 0.08)","fl");
 495  myLegendPad1->AddEntry(GrIdealSchenke3040_v3,"v_{3} (#eta/s = 0.0)","fl");
 496  myLegendPad1->AddEntry(GrViscousSchenke3040_v3,"v_{3} (#eta/s = 0.08)","fl");
 497  //myLegendPad1->SetTextSize(0.024);
 498  myLegendPad1->SetMargin(0.2);
 499  myLegendPad1->SetHeader("Centrality 30-40%");
 500  //TLegendEntry *header = (TLegendEntry*)myLegendPad1->GetListOfPrimitives()->First();
 501  //header->SetTextAlign(22); // center header
 502  //header->SetTextSize(0.027);
 503  //myLegendPad1->SetTextFont(22); // 22 = Times New Roman (bold)
 504  myLegendPad1->Draw("same");
 505
 506  //=================================================================================================
 507
 508  // Figure 3b:
 509  myCan->cd(2);
 510  gPad->SetTopMargin(0.);
 511  gPad->SetBottomMargin(0.);
 512  gPad->SetRightMargin(0.001);
 513  TH1D* styleHistPad2 = myBlankHisto->Clone("styleHistPad2");
 514  styleHistPad2->GetYaxis()->SetRangeUser(-0.01,0.1299);
 515  styleHistPad2->Draw();
 516  TLatex *b = new TLatex(0.95,0.93,"(b)");
 517  b->SetNDC();
 518  b->SetTextSize(0.05);
 519  b->Draw();
 520  TLegend *myLegendPad2 = new TLegend(0.12,0.6,0.45,0.98,"Centrality 0-5%");
 521  myLegendSetUp(myLegendPad2,0.04);
 522  //myLegendPad2->AddEntry(SP_0005ALICE_v2_etaGap10,"v_{2}{SP} (#Delta#eta=1.0)","p");
 523  //myLegendPad2->AddEntry(SP_0005ALICE_v2_etaGap02,"v_{2}{SP} (#Delta#eta=0.2)","p");
 524  //myLegendPad2->AddEntry(SP_0005ALICE_v3_etaGap10,"v_{3}{SP} (#Delta#eta=1.0)","p");
 525  //myLegendPad2->AddEntry(SP_0005ALICE_v3_etaGap02,"v_{3}{SP} (#Delta#eta=0.2)","p");
 526  //myLegendPad2->AddEntry(SP_0005ALICE_v4_etaGap10,"v_{4}{SP} (#Delta#eta=1.0)","p");
 527  //myLegendPad2->AddEntry(SP_0005ALICE_v4_etaGap02,"v_{4}{SP} (#Delta#eta=0.2)","p");
 528  //myLegendPad2->AddEntry(SP_0005ALICE_v5_etaGap10,"v_{5}{SP} (#Delta#eta=1.0)","p");
 529  //myLegendPad2->AddEntry(SP_0005ALICE_v5_etaGap02,"v_{5}{SP} (#Delta#eta=0.2)","p");
 530  myLegendPad2->AddEntry(GrSP_0005ALICE_v2_etaGap02,"v_{2}{2}","p");
 531  myLegendPad2->AddEntry(GrSP_0005ALICE_v3_etaGap02,"v_{3}{2}","p");
 532  myLegendPad2->AddEntry(GrSP_0005ALICE_v4_etaGap02,"v_{4}{2}","p");
 533  myLegendPad2->AddEntry(GrSP_0005ALICE_v5_etaGap02,"v_{5}{2}","p");
 534  myLegendPad2->Draw("same");
 535  GrSP_0005ALICE_v2_etaGap02->Draw("psame");
 536  GrSP_0005ALICE_v3_etaGap02->Draw("psame");
 537  GrSP_0005ALICE_v4_etaGap02->Draw("psame");
 538  GrSP_0005ALICE_v5_etaGap02->Draw("psame");
 539  GrSP_0005ALICE_v2_etaGap10->Draw("psame");
 540  GrSP_0005ALICE_v3_etaGap10->Draw("psame");
 541  GrSP_0005ALICE_v4_etaGap10->Draw("psame");
 542  GrSP_0005ALICE_v5_etaGap10->Draw("psame");
 543
 544  //=================================================================================================
 545
 546  // Figure 3c:
 547  myCan->cd(3);
 548  gPad->SetTopMargin(0.);
 549  gPad->SetBottomMargin(0.15);
 550  gPad->SetRightMargin(0.001);
 551  TH1D* styleHistPad3 = myBlankHisto->Clone("styleHistPad3");
 552  styleHistPad3->SetXTitle("p_{t} (GeV/#font[72]{c})");
 553  styleHistPad3->GetYaxis()->SetRangeUser(-0.01,0.1099);
 554  styleHistPad3->Draw();
 555  TLatex *c = new TLatex(0.95,0.93,"(c)");
 556  c->SetNDC();
 557  c->SetTextSize(0.05);
 558  c->Draw();
 559  TLegend *myLegendPad3 = new TLegend(0.12,0.65,0.45,0.98,"Centrality 0-2%");
 560  myLegendSetUp(myLegendPad3,0.04);
 561  // myLegendPad3->AddEntry(SP_0002_v2_etaGap02,"v_{2}{SP} (#Delta#eta = 0.2)","p");
 562  myLegendPad3->AddEntry(GrSP_0002_v2_etaGap02,"v_{2}{2}","p");
 563  // myLegendPad3->AddEntry(SP_0002_v2_etaGap10,"v_{2}{SP} (#Delta#eta = 1.0)","p");
 564  // myLegendPad3->AddEntry(SP_0002_v3_etaGap02,"v_{3}{SP} (#Delta#eta = 0.2)","p");
 565  myLegendPad3->AddEntry(GrSP_0002_v3_etaGap02,"v_{3}{2}","p");
 566  // myLegendPad3->AddEntry(SP_0002_v3_etaGap10,"v_{3}{SP} (#Delta#eta = 1.0)","p");
 567  myLegendPad3->AddEntry(GrSP_0002_v4_etaGap02,"v_{4}{2}","p");
 568  // myLegendPad3->AddEntry(SP_0002_v4_etaGap02,"v_{4}{SP} (#Delta#eta = 0.2)","p");
 569  // myLegendPad3->AddEntry(SP_0002_v4_etaGap10,"v_{4}{SP} (#Delta#eta = 1.0)","p");
 570  // myLegendPad3->AddEntry(SP_0002_v5_etaGap02,"v_{5}{SP} (#Delta#eta = 0.2)","p");
 571  myLegendPad3->AddEntry(GrSP_0002_v5_etaGap02,"v_{5}{2}","p");
 572  myLegendPad3->Draw("same");
 573  // myLegendPad3->AddEntry(SP_0002_v5_etaGap10,"v_{5}{SP} (#Delta#eta = 1.0)","p");
 574  GrSP_0002_v2_etaGap02->Draw("psame");
 575  GrSP_0002_v2_etaGap10->Draw("psame");
 576  GrSP_0002_v3_etaGap02->Draw("psame");
 577  GrSP_0002_v3_etaGap10->Draw("psame");
 578  GrSP_0002_v4_etaGap02->Draw("psame");
 579  GrSP_0002_v4_etaGap10->Draw("psame");
 580  GrSP_0002_v5_etaGap02->Draw("psame");
 581  GrSP_0002_v5_etaGap10->Draw("psame");
 582
 583
 584  return;
 585
 586
 587     /*
 588     TPad *myPadLogo = new TPad("myPadLogo", "Pad for ALICE Logo",0.47,0.37,0.62,0.60);
 589     //myPadLogo->SetFillColor(2); // color to first figure out where is the pad then comment !
 590     myPadSetUp(myPadLogo,0,0,0,0);
 591     myPadLogo->Draw();
 592     myPadLogo->cd();
 593     TASImage *myAliceLogo = new TASImage("alice_logo_transparent.png");
 594     //.TASImage *myAliceLogo = new TASImage("alice_logo.pdf");
 595     myAliceLogo->Draw();
 596     */
 597
 598     if (rPerformance){
 599         TLatex *alice = new TLatex(0.75,0.34,"Performance");
 600         alice->SetNDC();
 601         alice->SetTextColor(myDarkRed);
 602         //    alice->SetTextFont(42);
 603         alice->SetTextSize(0.05);
 604         alice->SetLineWidth(2);
 605         alice->Draw();
 606
 607         TText *date = new TText(0.78,0.28,cStamp2);
 608         date->SetNDC();
 609         date->SetTextFont(42);
 610         date->SetTextSize(0.04);
 611         date->Draw();
 612
 613         TPad *myPadLogo = new TPad("myPadLogo", "Pad for ALICE Logo",0.77,0.37,0.92,0.60);
 614         myPadLogo->SetFillColor(2); // color to first figure out where is the pad then comment !
 615         myPadSetUp(myPadLogo,0,0,0,0);
 616         myPadLogo->Draw();
 617         myPadLogo->cd();
 618         TASImage *myAliceLogo = new TASImage("alice_logo_transparent.png");
 619         //.TASImage *myAliceLogo = new TASImage("alice_logo.pdf");
 620
 621         myAliceLogo->Draw();
 622     }
 623     if (rWrite == 1)  myCan->SaveAs("fig_template.pdf");
 624     if (rWrite == 2)  myCan->SaveAs("fig_template.png");
 625 }
 626
 627 void myLegendSetUp(TLegend *currentLegend=0,float currentTextSize=0.07){
 628     currentLegend->SetTextFont(42);
 629     currentLegend->SetBorderSize(0);
 630     currentLegend->SetFillStyle(0);
 631     currentLegend->SetFillColor(0);
 632     currentLegend->SetMargin(0.25);
 633     currentLegend->SetTextSize(currentTextSize);
 634     currentLegend->SetEntrySeparation(0.5);
 635     return;
 636 }
 637
 638 void myPadSetUp(TPad *currentPad, float currentLeft=0.11, float currentTop=0.04, float currentRight=0.04, float currentBottom=0.15){
 639     currentPad->SetLeftMargin(currentLeft);
 640     currentPad->SetTopMargin(currentTop);
 641     currentPad->SetRightMargin(currentRight);
 642     currentPad->SetBottomMargin(currentBottom);
 643     return;
 644 }
 645
 646 void myGraphSetUp(TGraphErrors *currentGraph=0, Float_t currentMarkerSize = 1.0,
 647                   int currentMarkerStyle=20, int currentMarkerColor=0,
 648                   int currentLineStyle=1, int currentLineColor=0, int currentFillColor=0, int currentFillStyle=0)
 649 {
 650     currentGraph->SetMarkerSize(currentMarkerSize);
 651     currentGraph->SetMarkerStyle(currentMarkerStyle);
 652     currentGraph->SetMarkerColor(currentMarkerColor);
 653     currentGraph->SetLineStyle(currentLineStyle);
 654     currentGraph->SetLineColor(currentLineColor);
 655     currentGraph->SetLineWidth(2.);
 656     currentGraph->SetFillColor(currentFillColor);
 657     currentGraph->SetFillStyle(currentFillStyle);
 658     return;
 659 }
 660
 661 void ShiftAlongXaxis(TGraphErrors *ge, Double_t shift)
 662 {
 663  // Shift original TGraphErrors along x-axis by Double_t shift.
 664
 665  if(!ge){cout<<"!ge"<<endl;exit(0);}
 666
 667  Int_t nPoints = ge->GetN();
 668  Double_t x = 0.;
 669  Double_t y = 0.;
 670  for(Int_t p=0;p<nPoints;p++)
 671  {
 672   ge->GetPoint(p,x,y);
 673   x+=shift;
 674   ge->SetPoint(p,x,y);
 675  } // end of for(Int_t p=0;p<nPoints;p++)
 676
 677 } // end of void ShiftAlongXaxis(TGraphErrors *ge, Double_t shift)
 678
 679 void myOptions(Int_t lStat=0){
 680     // Set gStyle
 681     int font = 42;
 682     // From plain
 683     gStyle->SetFrameBorderMode(0);
 684     gStyle->SetFrameFillColor(0);
 685     gStyle->SetCanvasBorderMode(0);
 686     gStyle->SetPadBorderMode(0);
 687     gStyle->SetPadColor(10);
 688     gStyle->SetCanvasColor(10);
 689     gStyle->SetTitleFillColor(10);
 690     gStyle->SetTitleBorderSize(1);
 691     gStyle->SetStatColor(10);
 692     gStyle->SetStatBorderSize(1);
 693     gStyle->SetLegendBorderSize(1);
 694     //
 695     gStyle->SetDrawBorder(0);
 696     gStyle->SetTextFont(font);
 697     gStyle->SetStatFont(font);
 698     gStyle->SetStatFontSize(0.05);
 699     gStyle->SetStatX(0.97);
 700     gStyle->SetStatY(0.98);
 701     gStyle->SetStatH(0.03);
 702     gStyle->SetStatW(0.3);
 703     gStyle->SetTickLength(0.02,"y");
 704     gStyle->SetEndErrorSize(3);
 705     gStyle->SetLabelSize(0.05,"xyz");
 706     gStyle->SetLabelFont(font,"xyz");
 707     gStyle->SetLabelOffset(0.01,"xyz");
 708     gStyle->SetTitleFont(font,"xyz");
 709     gStyle->SetTitleOffset(1.0,"xyz");
 710     gStyle->SetTitleSize(0.06,"xyz");
 711     gStyle->SetMarkerSize(1);
 712     gStyle->SetPalette(1,0);
 713     if (lStat){
 714         gStyle->SetOptTitle(1);
 715         gStyle->SetOptStat(1111);
 716         gStyle->SetOptFit(1111);
 717     }
 718     else {
 719         gStyle->SetOptTitle(0);
 720         gStyle->SetOptStat(0);
 721         gStyle->SetOptFit(0);
 722     }
 723 }
 724
 725 TH1D* StyleHistogram()
 726 {
 727     // Style histogram for top panel:
 728
 729     Int_t nPtBins = 100;
 730     Double_t ptMin = 0.;
 731     Double_t ptMax = 10.0;
 732     TString xTitle  = "p_{t} (GeV/#font[72]{c})";
 733     TString yTitle  = "v_{n}";
 734
 735     TH1D* hist = new TH1D("","",nPtBins,ptMin,ptMax);
 736     hist->SetXTitle(xTitle.Data());
 737     hist->SetYTitle(yTitle.Data());
 738     hist->GetYaxis()->SetLabelSize(0.05);
 739     hist->GetYaxis()->SetTitleSize(0.07);
 740     hist->GetYaxis()->SetTitleOffset(1.10);
 741     hist->SetMinimum(-0.02101); // minimum y-value
 742     hist->SetMaximum(0.3799); // maximum y-value
 743     hist->GetXaxis()->SetRangeUser(0.,plotUpTo);
 744
 745     return hist;
 746
 747 } // end of TH1D* StyleHistogram()
 748
 749 //============================================================================================================
 750
 751 // OUTPUT FILES:
 752
 753 // Figure 3a in "Other harmonics" paper:
 754 // (pt dependence of harmonics v2, v3, v4 and v5 in centrality 30-40%).
 755 // Output files (all harmonics):
 756 //  dEta = 1.0: /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/QisQb/etaGap10/merged/outputCentrality.root
 757 //  dEta = 0.2: /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/QisQb/etaGap02/merged/outputCentrality.root
 758 // Objects holding measurements are following TGraphErrors:
 759 //  1.) SP_3040ALICE_v2_etaGap10: v2{SP}(pt) in centrality 30-40%, rapidity gap = 1.0;
 760 //  2.) SP_3040ALICE_v3_etaGap10: v3{SP}(pt) in centrality 30-40%, rapidity gap = 1.0;
 761 //  3.) SP_3040ALICE_v4_etaGap10: v4{SP}(pt) in centrality 30-40%, rapidity gap = 1.0;
 762 //  4.) SP_3040ALICE_v5_etaGap10: v5{SP}(pt) in centrality 30-40%, rapidity gap = 1.0;
 763 //  5.) SP_3040ALICE_v2_etaGap02: v2{SP}(pt) in centrality 30-40%, rapidity gap = 0.2;
 764 //  6.) SP_3040ALICE_v3_etaGap02: v3{SP}(pt) in centrality 30-40%, rapidity gap = 0.2;
 765 //  7.) SP_3040ALICE_v4_etaGap02: v4{SP}(pt) in centrality 30-40%, rapidity gap = 0.2;
 766 //  8.) SP_3040ALICE_v5_etaGap02: v5{SP}(pt) in centrality 30-40%, rapidity gap = 0.2.
 767
 768 // Figure 3b in "Other harmonics" paper:
 769 // (pt dependence of harmonics v2, v3, v4 and v5 in centrality 00-05%).
 770 // Output files (all harmonics):
 771 //  dEta = 1.0: /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/QisQb/etaGap10/merged/outputCentrality.root
 772 //  dEta = 0.2: /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/QisQb/etaGap02/merged/outputCentrality.root
 773 // Objects holding measurements are following TGraphErrors:
 774 //  1.) SP_0005ALICE_v2_etaGap10: v2{SP}(pt) in centrality 00-05%, rapidity gap = 1.0;
 775 //  2.) SP_0005ALICE_v3_etaGap10: v3{SP}(pt) in centrality 00-05%, rapidity gap = 1.0;
 776 //  3.) SP_0005ALICE_v4_etaGap10: v4{SP}(pt) in centrality 00-05%, rapidity gap = 1.0;
 777 //  4.) SP_0005ALICE_v5_etaGap10: v5{SP}(pt) in centrality 00-05%, rapidity gap = 1.0;
 778 //  5.) SP_0005ALICE_v2_etaGap02: v2{SP}(pt) in centrality 00-05%, rapidity gap = 0.2;
 779 //  6.) SP_0005ALICE_v3_etaGap02: v3{SP}(pt) in centrality 00-05%, rapidity gap = 0.2;
 780 //  7.) SP_0005ALICE_v4_etaGap02: v4{SP}(pt) in centrality 00-05%, rapidity gap = 0.2;
 781 //  8.) SP_0005ALICE_v5_etaGap02: v5{SP}(pt) in centrality 00-05%, rapidity gap = 0.2.
 782
 783 // Figure 3c in "Other harmonics" paper:
 784 // (pt dependence of harmonics v2, v3, v4, v5 in centrality 00-02%).
 785 // (Output files - see at the end of the macro).
 786 // Objects holding measurements are following TGraphErrors:
 787 //  1.) SP_0002_v2_etaGap02: v2{SP}(pt) in 0-2%, eta gap = 0.2;
 788 //  2.) SP_0002_v3_etaGap02: v3{SP}(pt) in 0-2%, eta gap = 0.2;
 789 //  3.) SP_0002_v4_etaGap02: v4{SP}(pt) in 0-2%, eta gap = 0.2;
 790 //  4.) SP_0002_v5_etaGap02: v5{SP}(pt) in 0-2%, eta gap = 0.2;
 791 //  5.) SP_0002_v2_etaGap10: v2{SP}(pt) in 0-2%, eta gap = 1.0;
 792 //  6.) SP_0002_v3_etaGap10: v3{SP}(pt) in 0-2%, eta gap = 1.0;
 793 //  7.) SP_0002_v4_etaGap10: v4{SP}(pt) in 0-2%, eta gap = 1.0;
 794 //  8.) SP_0002_v5_etaGap10: v5{SP}(pt) in 0-2%, eta gap = 1.0.
 795 //  eta gap = 0.2, i.e. Qa = [-0.8,-0.1], Qb = [0.1,0.8], POIs in [-0.8,-0.1], Q = Qb
 796 //                  OR  Qa = [-0.8,-0.1], Qb = [0.1,0.8], POIs in [0.1,0.8], Q = Qa
 797 //            /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/cent02/QisQb/etaGap02/merged/outputCentrality.root
 798 //            /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/cent02/QisQa/etaGap02/merged/outputCentrality.root
 799 //            /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/cent02/mergedQisQa_QisQb/etaGap02/outputCentrality.root
 800 //  eta gap = 1.0, i.e. Qa = [-0.8,-0.5], Qb = [0.5,0.8], POIs in [-0.8,-0.5], Q = Qb
 801 //                  OR  Qa = [-0.8,-0.5], Qb = [0.5,0.8], POIs in [0.5,0.8], Q = Qa
 802 //            /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/cent02/QisQb/etaGap10/merged/outputCentrality.root
 803 //            /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/cent02/QisQa/etaGap10/merged/outputCentrality.root
 804 //            /data/alice2/ab/grid/21-19/trainV0/TPConly/all/sp/cent02/mergedQisQa_QisQb/etaGap10/outputCentrality.root