| fee95d0f |
1 | // Merged original figures 2a and 2b into one figure. |
| 2 | |
| 3 | void figure2() |
| 4 | { |
| 5 | // Set style: |
| 6 | SetFlowStyle(); |
| 7 | |
| 8 | |
| 9 | Int_t kDarkGreen = TColor::GetColor(34,139,34); |
| 10 | Int_t kDarkCyan = TColor::GetColor(0,180,180); |
| 11 | |
| 12 | |
| 13 | |
| 14 | const int nbins = 17; |
| 15 | |
| 16 | //-----70-80% --------------------------------------------------------------------------- |
| 17 | |
| 18 | double pt_7080[nbins] = { |
| 19 | 0.29644,0.475364,0.700708,0.895582,1.07612, |
| 20 | 1.30166,1.4961,1.67731,1.90243,2.18302,2.56019,2.96531, |
| 21 | 3.39006,3.93562,5.10687,7.32286, 10.2651 |
| 22 | }; |
| 23 | double errPtLow_7080[nbins] = { |
| 24 | 0.0964398,0.075364,0.100708,0.0955816, |
| 25 | 0.076121,0.101662,0.0960967,0.0773092,0.102431, |
| 26 | 0.183019,0.160193,0.165308,0.19006,0.335623,0.606867,0.822858,0.765051 |
| 27 | }; |
| 28 | double errPtHigh_7080[nbins] = { |
| 29 | 0.10356,0.124636,0.0992923,0.104418, |
| 30 | 0.123879,0.0983377,0.103903,0.122691,0.097569,0.216981, |
| 31 | 0.239806,0.234692,0.209939,0.564377,1.39313,2.17714,1.73495 |
| 32 | }; |
| 33 | double v2_7080[nbins] ={ |
| 34 | 6.18034,7.78103,9.58698,16.708,16.1315,14.3096, |
| 35 | 11.8754,18.648,5.1453,6.53926,-6.88641,-13.9267,-52.5485, |
| 36 | -191.067,-415.917,291.021,45.6615 |
| 37 | }; |
| 38 | double errV2_7080[nbins] = { |
| 39 | 0.566557,0.582667,0.830618,1.38276,1.49568, |
| 40 | 2.06431,3.40244,3.51732,4.83874,5.55657,8.97981, |
| 41 | 13.9386,22.9194,31.4821,76.8651,-9999,-9999 |
| 42 | }; |
| 43 | |
| 44 | //-----60-70% ---------------------------------------------------------------------------------- |
| 45 | |
| 46 | double pt_6070[nbins] = { |
| 47 | 0.296676,0.47597,0.701133,0.895742,1.07645, |
| 48 | 1.30194,1.49616,1.67738,1.90256,2.18364,2.56013,2.96414, |
| 49 | 3.38945,3.93604,5.09071,7.3686,10.683 |
| 50 | }; |
| 51 | double errPtLow_6070[nbins] = { |
| 52 | 0.0966761,0.0759702,0.101133,0.0957418, |
| 53 | 0.0764486,0.101935,0.0961564,0.0773785,0.102563,0.183642, |
| 54 | 0.160132,0.164139,0.189453,0.33604,0.59071,0.868602,1.18298 |
| 55 | }; |
| 56 | double errPtHigh_6070[nbins] = { |
| 57 | 0.103324,0.12403,0.0988673,0.104258, |
| 58 | 0.123551,0.0980649,0.103844,0.122621,0.0974373, |
| 59 | 0.216358,0.239868,0.235861,0.210547,0.563961,1.40929,2.1314,1.31702 |
| 60 | }; |
| 61 | double v2_6070[nbins] = { |
| 62 | 3.49974,5.88458,8.12386,10.0714,12.5145,11.6577, |
| 63 | 11.5439,13.5126,11.4789,10.7059,11.6547,9.49895, |
| 64 | 1.37392,5.47184,-1.3349,-45.6053,156.746 |
| 65 | }; |
| 66 | double errV2_6070[nbins] = { |
| 67 | 0.177159,0.179905,0.250419,0.406175,0.431074, |
| 68 | 0.595281,0.940248,0.987307,1.32458,1.50529,2.40035, |
| 69 | 3.48782,6.0306,7.22813,13.5038,56.712,-9999 |
| 70 | }; |
| 71 | |
| 72 | |
| 73 | //-----50-60% ----------------------------------------------------------------------------------- |
| 74 | |
| 75 | double pt_5060[nbins] = { |
| 76 | 0.296932,0.47652,0.701555,0.895881,1.07678,1.30221, |
| 77 | 1.49621,1.67759,1.90265,2.18367,2.55994,2.96436, |
| 78 | 3.38993,3.93277,5.08075,7.40325,10.5265 |
| 79 | }; |
| 80 | double errPtLow_5060[nbins] = { |
| 81 | 0.0969318,0.0765199,0.101555,0.0958806, |
| 82 | 0.0767766,0.102209,0.0962076,0.0775871,0.102653,0.183672, |
| 83 | 0.159938,0.164359,0.189935,0.332773,0.580747,0.903252,1.02646 |
| 84 | }; |
| 85 | double errPtHigh_5060[nbins] = { |
| 86 | 0.103068,0.12348,0.0984446,0.104119, |
| 87 | 0.123223,0.097791,0.103792,0.122413,0.0973468,0.216328, |
| 88 | 0.240062,0.235642,0.210065,0.567227,1.41925,2.09675,1.47354 |
| 89 | }; |
| 90 | double v2_5060[nbins] = { |
| 91 | 3.44798,5.90898,8.4642,10.5149,11.9192,13.3181, |
| 92 | 14.0252,14.3745,15.0161,14.8689,13.5232,13.4977,11.6071, |
| 93 | 9.47076,7.71149,17.5352,-5.54787 |
| 94 | }; |
| 95 | double errV2_5060[nbins] = { |
| 96 | 0.0606734,0.0606551,0.0829413,0.132536, |
| 97 | 0.139612,0.190289,0.299361,0.317186,0.414861,0.475373, |
| 98 | 0.769568,1.05829,1.91803,2.09282,3.91973,13.8477,40.6673 |
| 99 | }; |
| 100 | |
| 101 | //-----40-50% ---------------------------------------------------------------------------------- |
| 102 | |
| 103 | double pt_4050[nbins] = { |
| 104 | 0.297194,0.47704,0.701961,0.896023,1.07714,1.30245, |
| 105 | 1.4963,1.6777,1.90274,2.18372,2.55965,2.96361,3.38975, |
| 106 | 3.93175,5.08329,7.35969,10.4992 |
| 107 | }; |
| 108 | double errPtLow_4050[nbins] = { |
| 109 | 0.0971937,0.0770399,0.101961,0.0960231,0.0771427, |
| 110 | 0.102452,0.0963024,0.077695,0.10274,0.183724, |
| 111 | 0.159645,0.16361,0.189751,0.331754,0.583295,0.859692,0.999224 |
| 112 | }; |
| 113 | double errPtHigh_4050[nbins] = { |
| 114 | 0.102806,0.12296,0.0980395,0.103977,0.122857,0.0975482, |
| 115 | 0.103698,0.122305,0.0972599,0.216276,0.240355, |
| 116 | 0.23639,0.210248,0.568246,1.41671,2.14031,1.50078 |
| 117 | }; |
| 118 | double v2_4050[nbins] = { |
| 119 | 3.64534,5.97555,8.71339,10.8948,12.3118,14.0208,15.1841, |
| 120 | 16.2186,17.0478,17.4933,17.837,17.5645,16.0996,14.2722, |
| 121 | 11.2284,14.8254,39.8013 |
| 122 | }; |
| 123 | double errV2_4050[nbins] = { |
| 124 | 0.0280054,0.0275449,0.0370108,0.0584979,0.0610596, |
| 125 | 0.0821168,0.129148,0.133787,0.179269,0.201004,0.319255, |
| 126 | 0.465269,0.805499,0.917085,1.69463,5.79271,-9999 |
| 127 | }; |
| 128 | |
| 129 | //-----30-40% ----------------------------------------------------------------------------------- |
| 130 | |
| 131 | double pt_3040[nbins] = { |
| 132 | 0.297459,0.477509,0.702306,0.896152,1.07742,1.30267, |
| 133 | 1.49633,1.67775,1.90278,2.18368,2.55957,2.96322,3.38921, |
| 134 | 3.92869,5.07593,7.41037,10.408 |
| 135 | }; |
| 136 | double errPtLow_3040[nbins] = { |
| 137 | 0.0974586,0.0775094,0.102306,0.0961525, |
| 138 | 0.0774162,0.102673,0.096333,0.077749,0.10278,0.183675, |
| 139 | 0.159572,0.163219,0.189209,0.32869,0.575929,0.910369,0.907955 |
| 140 | }; |
| 141 | double errPtHigh_3040[nbins] = { |
| 142 | 0.102541,0.122491,0.0976941,0.103848, |
| 143 | 0.122584,0.0973271,0.103667,0.122251,0.0972203,0.216325, |
| 144 | 0.240428,0.236781,0.210791,0.57131,1.42407,2.08963,1.59205 |
| 145 | }; |
| 146 | double v2_3040[nbins] = { |
| 147 | 3.32144,5.60995,8.20984,10.1864,11.7526,13.5713, |
| 148 | 14.8011,15.8389,16.7527,17.5895,18.1343,18.2392,18.4816,16.8442, |
| 149 | 14.3587,7.07376,-26.4391 |
| 150 | }; |
| 151 | double errV2_3040[nbins] = { |
| 152 | 0.0171702,0.0166421,0.0220506,0.034594,0.0358166, |
| 153 | 0.0478107,0.0746667,0.0773614,0.103065,0.115688, |
| 154 | 0.184685,0.270488,0.46637,0.53962,1.01023,3.51725,12.6064 |
| 155 | }; |
| 156 | |
| 157 | |
| 158 | //-----20-30% ----------------------------------------------------------------------------------- |
| 159 | |
| 160 | double pt_2030[nbins] = { |
| 161 | 0.297722,0.477941,0.702611,0.896268,1.07768,1.30284, |
| 162 | 1.49638,1.67786,1.9028,2.18362,2.55924,2.96217,3.38881, |
| 163 | 3.92699,5.07304,7.42509,10.5134 |
| 164 | }; |
| 165 | double errPtLow_2030[nbins] = { |
| 166 | 0.0977221,0.0779407,0.102611,0.0962682, |
| 167 | 0.0776817,0.102839,0.0963758,0.0778584,0.102797, |
| 168 | 0.183618,0.159242,0.162174,0.188808,0.326986,0.573036,0.925089,1.01339 |
| 169 | }; |
| 170 | double errPtHigh_2030[nbins] = { |
| 171 | 0.102278,0.122059,0.0973894,0.103732, |
| 172 | 0.122318,0.0971605,0.103624,0.122142,0.0972035,0.216382, |
| 173 | 0.240758,0.237826,0.211192,0.573015,1.42696,2.07491,1.48661 |
| 174 | }; |
| 175 | double v2_2030[nbins] = { |
| 176 | 2.95888,4.90692,7.19851,8.89996,10.2714,11.926,13.1472, |
| 177 | 14.1925,15.2307,16.2584,16.8146,16.351,16.1069,14.8237,13.7406, |
| 178 | 6.26249,-2.67304 |
| 179 | }; |
| 180 | double errV2_2030[nbins] = { |
| 181 | 0.0134584,0.0128658,0.016865,0.0262622, |
| 182 | 0.0270303,0.03587,0.0557706,0.0576958,0.0765681,0.0859164, |
| 183 | 0.13768,0.203798,0.357643,0.419151,0.826191,2.70724,8.88004 |
| 184 | }; |
| 185 | |
| 186 | //-----10-20% ----------------------------------------------------------------------------------- |
| 187 | |
| 188 | double pt_1020[nbins] = { |
| 189 | 0.298001,0.478337,0.702859,0.896349,1.07787,1.30299, |
| 190 | 1.49642,1.6779,1.90278,2.18349,2.55915,2.96176,3.38798, |
| 191 | 3.92357,5.07047,7.40607,10.4205 |
| 192 | }; |
| 193 | double errPtLow_1020[nbins] = { |
| 194 | 0.0980014,0.0783374,0.102859,0.0963492, |
| 195 | 0.0778747,0.102987,0.0964151,0.0779025,0.102782,0.183491, |
| 196 | 0.159151,0.161757,0.187976,0.323571,0.57047,0.906065,0.920542 |
| 197 | }; |
| 198 | double errPtHigh_1020[nbins] = { |
| 199 | 0.101999,0.121663,0.0971407,0.103651, |
| 200 | 0.122125,0.097013,0.103585,0.122097,0.0972182,0.216509, |
| 201 | 0.240849,0.238243,0.212024,0.576429,1.42953,2.09393,1.57946 |
| 202 | }; |
| 203 | double v2_1020[nbins] = { |
| 204 | 2.17225,3.66955,5.32934,6.65091,7.67923,8.95035,9.92419, |
| 205 | 10.6285,11.4726,12.4663,13.2892,13.1702,13.383,12.3723,8.93792, |
| 206 | 12.5296,15.0085 |
| 207 | }; |
| 208 | double errV2_1020[nbins] = { |
| 209 | 0.0162515,0.0153472,0.0199456,0.0308611, |
| 210 | 0.031683,0.0419599,0.0651578,0.0672946,0.0898034,0.100566, |
| 211 | 0.16173,0.243491,0.434311,0.50903,1.02401,3.4382,9.82681 |
| 212 | }; |
| 213 | |
| 214 | |
| 215 | //-----5-10% ----------------------------------------------------------------------------------- |
| 216 | |
| 217 | double pt_510[nbins] = { |
| 218 | 0.298249,0.478648,0.703017,0.896405,1.078,1.30306,1.49641, |
| 219 | 1.6779,1.90274,2.18332,2.55882,2.96126,3.38753,3.9217,5.065,7.43055,10.4663 |
| 220 | }; |
| 221 | double errPtLow_510[nbins] = { |
| 222 | 0.0982494,0.0786479,0.103017,0.0964049, |
| 223 | 0.0779983,0.10306,0.0964067,0.0778999,0.102742,0.183315, |
| 224 | 0.158818,0.161262,0.187531,0.321702,0.565,0.930555,0.966266 |
| 225 | }; |
| 226 | double errPtHigh_510[nbins] = { |
| 227 | 0.101751,0.121352,0.0969831,0.103595, |
| 228 | 0.122002,0.0969397,0.103593,0.1221,0.0972583,0.216685, |
| 229 | 0.241182,0.238738,0.212469,0.578298,1.435,2.06945,1.53373 |
| 230 | }; |
| 231 | double v2_510[nbins] = { |
| 232 | 1.336,2.31457,3.34265,4.2049,4.76444,5.67546,5.92831, |
| 233 | 6.72256,6.95237,7.25684,8.04808,8.10241,5.13828,9.00966, |
| 234 | 1.0595,-7.48113,34.1494 |
| 235 | }; |
| 236 | double errV2_510[nbins] = { |
| 237 | 0.0490133,0.0459249,0.0593962,0.0917575, |
| 238 | 0.0941489,0.124114,0.193714,0.198987,0.266057,0.298391, |
| 239 | 0.482814,0.729011,1.32001,1.56362,3.12691,11.0502,56.8878 |
| 240 | }; |
| 241 | |
| 242 | |
| 243 | |
| 244 | double errPt[nbins] = {0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.}; |
| 245 | |
| 246 | double v2_7080_ecc[nbins], errV2_7080_ecc[nbins]; |
| 247 | double v2_6070_ecc[nbins], errV2_6070_ecc[nbins]; |
| 248 | double v2_5060_ecc[nbins], errV2_5060_ecc[nbins]; |
| 249 | double v2_4050_ecc[nbins], errV2_4050_ecc[nbins]; |
| 250 | double v2_3040_ecc[nbins], errV2_3040_ecc[nbins]; |
| 251 | double v2_2030_ecc[nbins], errV2_2030_ecc[nbins]; |
| 252 | double v2_1020_ecc[nbins], errV2_1020_ecc[nbins]; |
| 253 | double v2_510_ecc[nbins], errV2_510_ecc[nbins]; |
| 254 | |
| 255 | |
| 256 | for (int i= 0; i< nbins; i++) { |
| 257 | v2_7080_ecc[i] = 0.01*v2_7080[i]; |
| 258 | errV2_7080_ecc[i] = 0.01*errV2_7080[i]; |
| 259 | |
| 260 | v2_6070_ecc[i] = 0.01*v2_6070[i]; |
| 261 | errV2_6070_ecc[i] = 0.01*errV2_6070[i]; |
| 262 | |
| 263 | v2_5060_ecc[i] = 0.01*v2_5060[i]; |
| 264 | errV2_5060_ecc[i] = 0.01*errV2_5060[i]; |
| 265 | |
| 266 | v2_4050_ecc[i] = 0.01*v2_4050[i]; |
| 267 | errV2_4050_ecc[i] = 0.01*errV2_4050[i]; |
| 268 | |
| 269 | v2_3040_ecc[i] = 0.01*v2_3040[i]; |
| 270 | errV2_3040_ecc[i] = 0.01*errV2_3040[i]; |
| 271 | |
| 272 | v2_2030_ecc[i] = 0.01*v2_2030[i]; |
| 273 | errV2_2030_ecc[i] = 0.01*errV2_2030[i]; |
| 274 | |
| 275 | v2_1020_ecc[i] = 0.01*v2_1020[i]; |
| 276 | errV2_1020_ecc[i] = 0.01*errV2_1020[i]; |
| 277 | |
| 278 | v2_510_ecc[i] = 0.01*v2_510[i]; |
| 279 | errV2_510_ecc[i] = 0.01*errV2_510[i]; |
| 280 | } |
| 281 | |
| 282 | //=================================================================================================================== |
| 283 | // centrality 10-20% |
| 284 | Double_t xCumulant4th1020ALICE[] = {0.050000,0.150000,0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.850000,0.950000, |
| 285 | 1.100000,1.300000,1.500000,1.700000,1.900000,2.250000,2.750000,3.250000,3.750000,4.500000, |
| 286 | 5.500000,7.000000,9.000000}; |
| 287 | Double_t yCumulant4th1020ALICE[] = {0.000000,0.000000,0.024075,0.031505,0.040413,0.044981,0.055358,0.060563,0.063378,0.070030, |
| 288 | 0.082692,0.091611,0.099641,0.107223,0.122376,0.131240,0.137425,0.146050,0.131365,0.124708, |
| 289 | 0.000000,0.000000,0.000000}; |
| 290 | Double_t xErrCumulant4th1020ALICE[23] = {0.}; |
| 291 | Double_t yErrCumulant4th1020ALICE[] = {0.000000,0.000000,0.002413,0.002931,0.003444,0.003950,0.004338,0.004835,0.005059,0.005586, |
| 292 | 0.004521,0.005278,0.005999,0.007072,0.008260,0.007279,0.011897,0.017409,0.023995,0.025701, |
| 293 | 0.000000,0.000000,0.000000}; |
| 294 | Int_t nPointsCumulant4th1020ALICE = sizeof(xCumulant4th1020ALICE)/sizeof(Double_t); |
| 295 | TGraphErrors *Cumulant4th1020ALICE = new TGraphErrors(nPointsCumulant4th1020ALICE,xCumulant4th1020ALICE,yCumulant4th1020ALICE, |
| 296 | xErrCumulant4th1020ALICE,yErrCumulant4th1020ALICE); |
| 297 | Cumulant4th1020ALICE->SetMarkerStyle(kFullCircle); |
| 298 | Cumulant4th1020ALICE->SetMarkerColor(kBlue); |
| 299 | Cumulant4th1020ALICE->SetMarkerSize(1.2); |
| 300 | Cumulant4th1020ALICE->SetFillStyle(1001); |
| 301 | Cumulant4th1020ALICE->SetFillColor(kBlue-10); |
| 302 | //=================================================================================================================== |
| 303 | |
| 304 | //=================================================================================================================== |
| 305 | // centrality 20-30% |
| 306 | Double_t xCumulant4th2030ALICE[] = {0.050000,0.150000,0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.850000,0.950000, |
| 307 | 1.100000,1.300000,1.500000,1.700000,1.900000,2.250000,2.750000,3.250000,3.750000,4.500000, |
| 308 | 5.500000,7.000000,9.000000}; |
| 309 | Double_t yCumulant4th2030ALICE[] = {0.000000,0.000000,0.030926,0.041076,0.052063,0.059429,0.070500,0.084461,0.086745,0.099254, |
| 310 | 0.109691,0.116398,0.130831,0.141959,0.158932,0.169680,0.171387,0.178858,0.171475,0.140358, |
| 311 | 0.000000,0.000000,0.000000}; |
| 312 | Double_t xErrCumulant4th2030ALICE[23] = {0.}; |
| 313 | Double_t yErrCumulant4th2030ALICE[] = {0.000000,0.000000,0.002857,0.003451,0.003567,0.003859,0.004609,0.004976,0.005412,0.006277, |
| 314 | 0.004748,0.005808,0.006896,0.007987,0.008683,0.008080,0.013278,0.018413,0.024873,0.026057, |
| 315 | 0.000000,0.000000,0.000000}; |
| 316 | Int_t nPointsCumulant4th2030ALICE = sizeof(xCumulant4th2030ALICE)/sizeof(Double_t); |
| 317 | TGraphErrors *Cumulant4th2030ALICE = new TGraphErrors(nPointsCumulant4th2030ALICE,xCumulant4th2030ALICE,yCumulant4th2030ALICE, |
| 318 | xErrCumulant4th2030ALICE,yErrCumulant4th2030ALICE); |
| 319 | Cumulant4th2030ALICE->SetMarkerStyle(kFullSquare); |
| 320 | Cumulant4th2030ALICE->SetMarkerColor(kRed); |
| 321 | Cumulant4th2030ALICE->SetMarkerSize(1.2); |
| 322 | Cumulant4th2030ALICE->SetFillStyle(1001); |
| 323 | Cumulant4th2030ALICE->SetFillColor(kRed-10); |
| 324 | //ShiftAlongXaxis(Cumulant4th2030ALICE,0.03); |
| 325 | //=================================================================================================================== |
| 326 | |
| 327 | //=================================================================================================================== |
| 328 | // centrality 30-40% |
| 329 | Double_t xCumulant4th3040ALICE[] = {0.050000,0.150000,0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.850000,0.950000, |
| 330 | 1.100000,1.300000,1.500000,1.700000,1.900000,2.250000,2.750000,3.250000,3.750000,4.500000, |
| 331 | 5.500000,7.000000,9.000000}; |
| 332 | Double_t yCumulant4th3040ALICE[] = {0.000000,0.000000,0.037071,0.048566,0.061083,0.070910,0.078831,0.091396,0.102026,0.109691, |
| 333 | 0.124449,0.139819,0.155561,0.165701,0.173678,0.191149,0.202015,0.204540,0.212560,0.195885, |
| 334 | 0.000000,0.000000,0.000000}; |
| 335 | Double_t xErrCumulant4th3040ALICE[23] = {0.}; |
| 336 | Double_t yErrCumulant4th3040ALICE[] = {0.000000,0.000000,0.002992,0.003364,0.003669,0.003931,0.004698,0.005261,0.005446,0.006151, |
| 337 | 0.004980,0.005741,0.007198,0.008576,0.010868,0.009926,0.015269,0.020691,0.027601,0.031834, |
| 338 | 0.000000,0.000000,0.000000}; |
| 339 | Int_t nPointsCumulant4th3040ALICE = sizeof(xCumulant4th3040ALICE)/sizeof(Double_t); |
| 340 | TGraphErrors *Cumulant4th3040ALICE = new TGraphErrors(nPointsCumulant4th3040ALICE,xCumulant4th3040ALICE,yCumulant4th3040ALICE, |
| 341 | xErrCumulant4th3040ALICE,yErrCumulant4th3040ALICE); |
| 342 | Cumulant4th3040ALICE->SetMarkerStyle(kFullTriangleUp); |
| 343 | Cumulant4th3040ALICE->SetMarkerColor(kGreen+2); |
| 344 | Cumulant4th3040ALICE->SetMarkerSize(1.2); |
| 345 | Cumulant4th3040ALICE->SetFillStyle(1001); |
| 346 | Cumulant4th3040ALICE->SetFillColor(kGreen+2); |
| 347 | //ShiftAlongXaxis(Cumulant4th3040ALICE,0.0); |
| 348 | //=================================================================================================================== |
| 349 | |
| 350 | //=================================================================================================================== |
| 351 | // centrality 40-50% |
| 352 | // v2{2} |
| 353 | Double_t xCumulant2nd4050ALICE[] = {0.050000,0.150000,0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.850000,0.950000, |
| 354 | 1.100000,1.300000,1.500000,1.700000,1.900000,2.250000,2.750000,3.250000,3.750000,4.500000, |
| 355 | 5.500000,7.000000,9.000000}; |
| 356 | Double_t yCumulant2nd4050ALICE[] = {0.000000,0.000000,0.043400,0.059911,0.073516,0.089756,0.105486,0.117391,0.128199,0.138013, |
| 357 | 0.158271,0.177726,0.196383,0.208277,0.216648,0.242954,0.249961,0.240131,0.269006,0.207796, |
| 358 | 0.000000,0.000000,0.000000}; |
| 359 | Double_t xErrCumulant2nd4050ALICE[23] = {0.}; |
| 360 | Double_t yErrCumulant2nd4050ALICE[] = {0.000000,0.000000,0.001496,0.001611,0.001700,0.001895,0.002081,0.002335,0.002598,0.002839, |
| 361 | 0.002385,0.002908,0.003639,0.004457,0.005427,0.004813,0.007465,0.011109,0.016086,0.018067, |
| 362 | 0.000000,0.000000,0.000000}; |
| 363 | Int_t nPointsCumulant2nd4050ALICE = sizeof(xCumulant2nd4050ALICE)/sizeof(Double_t); |
| 364 | TGraphErrors *Cumulant2nd4050ALICE = new TGraphErrors(nPointsCumulant2nd4050ALICE,xCumulant2nd4050ALICE,yCumulant2nd4050ALICE, |
| 365 | xErrCumulant2nd4050ALICE,yErrCumulant2nd4050ALICE); |
| 366 | Cumulant2nd4050ALICE->SetMarkerStyle(kStar); |
| 367 | Cumulant2nd4050ALICE->SetMarkerColor(kBlue); |
| 368 | Cumulant2nd4050ALICE->SetMarkerSize(1.2); |
| 369 | Cumulant2nd4050ALICE->SetFillStyle(1001); |
| 370 | Cumulant2nd4050ALICE->SetFillColor(kBlue-10); |
| 371 | //ShiftAlongXaxis(Cumulant2nd4050ALICE,0.044); |
| 372 | |
| 373 | // v2{4} |
| 374 | Double_t xCumulant4th4050ALICE[] = {0.050000,0.150000,0.250000,0.350000,0.450000,0.550000,0.650000,0.750000,0.850000,0.950000, |
| 375 | 1.100000,1.300000,1.500000,1.700000,1.900000,2.250000,2.750000,3.250000,3.750000,4.500000, |
| 376 | 5.500000,7.000000,9.000000}; |
| 377 | Double_t yCumulant4th4050ALICE[] = {0.000000,0.000000,0.038646,0.049824,0.066662,0.075856,0.081583,0.099778,0.104674,0.118545, |
| 378 | 0.131874,0.152959,0.155348,0.169751,0.179052,0.178532,0.198851,0.185737,0.239901,0.186098, |
| 379 | 0.000000,0.000000,0.000000}; |
| 380 | Double_t xErrCumulant4th4050ALICE[23] = {0.}; |
| 381 | Double_t yErrCumulant4th4050ALICE[] = {0.000000,0.000000,0.004853,0.006008,0.006304,0.006887,0.007797,0.008258,0.009678,0.010500, |
| 382 | 0.008706,0.010322,0.012939,0.015934,0.019074,0.016231,0.023799,0.033236,0.046821,0.044138, |
| 383 | 0.000000,0.000000,0.000000}; |
| 384 | Int_t nPointsCumulant4th4050ALICE = sizeof(xCumulant4th4050ALICE)/sizeof(Double_t); |
| 385 | TGraphErrors *Cumulant4th4050ALICE = new TGraphErrors(nPointsCumulant4th4050ALICE,xCumulant4th4050ALICE,yCumulant4th4050ALICE, |
| 386 | xErrCumulant4th4050ALICE,yErrCumulant4th4050ALICE); |
| 387 | Cumulant4th4050ALICE->SetMarkerStyle(kFullTriangleDown); |
| 388 | Cumulant4th4050ALICE->SetMarkerColor(kRed); |
| 389 | Cumulant4th4050ALICE->SetMarkerSize(1.2); |
| 390 | Cumulant4th4050ALICE->SetFillStyle(1001); |
| 391 | Cumulant4th4050ALICE->SetFillColor(kRed-10); |
| 392 | //ShiftAlongXaxis(Cumulant4th4050ALICE,0.044); |
| 393 | //=================================================================================================================== |
| 394 | |
| 395 | v2_1020_gr = new TGraphErrors(nbins,pt_1020, v2_1020_ecc, errPt, errV2_1020_ecc); |
| 396 | v2_1020_gr->SetMarkerStyle(22); |
| 397 | v2_1020_gr->SetMarkerSize(1.2); |
| 398 | v2_1020_gr->SetLineColor(kBlack); |
| 399 | v2_1020_gr->SetMarkerColor(4); |
| 400 | v2_1020_gr->SetLineWidth(2); |
| 401 | v2_1020_gr->SetFillStyle(3244); |
| 402 | v2_1020_gr->SetFillColor(kBlack); |
| 403 | |
| 404 | v2_2030_gr = new TGraphErrors(nbins,pt_2030, v2_2030_ecc, errPt, errV2_2030_ecc); |
| 405 | v2_2030_gr->SetMarkerStyle(22); |
| 406 | v2_2030_gr->SetMarkerSize(1.2); |
| 407 | v2_2030_gr->SetLineColor(kBlack); |
| 408 | v2_2030_gr->SetMarkerColor(4); |
| 409 | v2_2030_gr->SetLineWidth(2); |
| 410 | v2_2030_gr->SetFillStyle(3144); |
| 411 | v2_2030_gr->SetFillColor(kBlack); |
| 412 | |
| 413 | v2_3040_gr = new TGraphErrors(nbins,pt_3040, v2_3040_ecc, errPt, errV2_3040_ecc); |
| 414 | v2_3040_gr->SetMarkerStyle(22); |
| 415 | v2_3040_gr->SetMarkerSize(1.2); |
| 416 | v2_3040_gr->SetLineColor(kBlack); |
| 417 | v2_3040_gr->SetMarkerColor(4); |
| 418 | v2_3040_gr->SetLineWidth(2); |
| 419 | v2_3040_gr->SetFillStyle(1001); |
| 420 | v2_3040_gr->SetFillColor(kGray+1); |
| 421 | |
| 422 | v2_4050_gr = new TGraphErrors(nbins,pt_4050, v2_4050_ecc, errPt, errV2_4050_ecc); |
| 423 | v2_4050_gr->SetMarkerStyle(22); |
| 424 | v2_4050_gr->SetMarkerSize(1.2); |
| 425 | v2_4050_gr->SetLineColor(kBlack); |
| 426 | v2_4050_gr->SetMarkerColor(4); |
| 427 | v2_4050_gr->SetLineWidth(2); |
| 428 | v2_4050_gr->SetFillStyle(1001); |
| 429 | v2_4050_gr->SetFillColor(kGray+1); |
| 430 | |
| 431 | // Legend in uppper panel: |
| 432 | TLegend *legendUp = new TLegend(0.20,0.57,0.37,0.94); |
| 433 | legendUp->SetFillStyle(0); // white legendUp background |
| 434 | legendUp->AddEntry(Cumulant2nd4050ALICE,"v_{2}{2}","p"); |
| 435 | legendUp->AddEntry(Cumulant4th4050ALICE,"v_{2}{4}","p"); |
| 436 | legendUp->AddEntry(v2_4050_gr,"v_{2}{4} (STAR)","f"); |
| 437 | |
| 438 | legendUp->SetTextSize(0.050); |
| 439 | legendUp->SetTextFont(22); // 22 = Times New Roman (bold) |
| 440 | |
| 441 | // Text in uppper panel: |
| 442 | TPaveText *paveTextUp = new TPaveText(2.49,0.04,2.74,0.07,"br"); |
| 443 | paveTextUp->AddText("centrality 40-50%"); |
| 444 | paveTextUp->SetTextFont(22); // 22 = Times New Roman (bold) |
| 445 | paveTextUp->SetTextSize(0.10); |
| 446 | paveTextUp->SetFillColor(0); // white pave text background |
| 447 | |
| 448 | // Final drawing u=in upper panel (order is important!): |
| 449 | TCanvas *v2pt = new TCanvas("v2pt","v2pt",10,10,750,700); |
| 450 | v2pt->Divide(1,2,0,0); |
| 451 | v2pt->cd(1); |
| 452 | StyleHistogramVsPtUp()->Draw(); |
| 453 | v2_4050_gr->Draw("same3"); // mesh |
| 454 | v2_4050_gr->Draw("lsameX"); // line |
| 455 | //gPad->SetTopMargin(0); |
| 456 | gPad->SetBottomMargin(0.0); |
| 457 | Cumulant2nd4050ALICE->Draw("psame"); |
| 458 | Cumulant4th4050ALICE->Draw("psame"); |
| 459 | legendUp->Draw("same"); |
| 460 | paveTextUp->Draw("same"); |
| 461 | v2pt->cd(2); |
| 462 | StyleHistogramVsPtUp()->Draw(); |
| 463 | v2_4050_gr->Draw("same3"); // mesh |
| 464 | v2_4050_gr->Draw("lsameX"); // line |
| 465 | gPad->SetTopMargin(0); |
| 466 | Cumulant2nd4050ALICE->Draw("psame"); |
| 467 | Cumulant4th4050ALICE->Draw("psame"); |
| 468 | |
| 469 | // Legend in lower panel: |
| 470 | v2pt->cd(2); |
| 471 | StyleHistogramVsPtDown()->Draw(); |
| 472 | TLegend *legendDown = new TLegend(0.20,0.59,0.39,0.96); |
| 473 | legendDown->SetFillStyle(0); // white legend background |
| 474 | legendDown->AddEntry(Cumulant4th1020ALICE,"10-20%","p"); |
| 475 | legendDown->AddEntry(Cumulant4th2030ALICE,"20-30%","p"); |
| 476 | legendDown->AddEntry(Cumulant4th3040ALICE,"30-40%","p"); |
| 477 | legendDown->AddEntry(v2_1020_gr,"10-20% (STAR)","f"); |
| 478 | legendDown->AddEntry(v2_2030_gr,"20-30% (STAR)","f"); |
| 479 | legendDown->AddEntry(v2_3040_gr,"30-40% (STAR)","f"); |
| 480 | v2_3040_gr->Draw("same3"); // mash |
| 481 | v2_3040_gr->Draw("lsameX"); // line |
| 482 | v2_2030_gr->Draw("same3"); // mash |
| 483 | v2_2030_gr->Draw("lsameX"); // line |
| 484 | v2_1020_gr->Draw("same3"); // mash |
| 485 | v2_1020_gr->Draw("lsameX"); // line |
| 486 | Cumulant4th1020ALICE->Draw("psame"); |
| 487 | Cumulant4th2030ALICE->Draw("psame"); |
| 488 | Cumulant4th3040ALICE->Draw("psame"); |
| 489 | legendDown->Draw("same"); |
| 490 | |
| 491 | legendDown->SetTextSize(0.038); |
| 492 | legendDown->SetTextFont(22); // 22 = Times New Roman (bold) |
| 493 | |
| 494 | //v2pt->Print("fig2.pdf"); |
| 495 | |
| 496 | } // end of void figure2b() |
| 497 | |
| 498 | //================================================================================== |
| 499 | |
| 500 | TH1D* StyleHistogramVsPtUp() |
| 501 | { |
| 502 | // Style histogram: |
| 503 | |
| 504 | Int_t nPtBins = 100; // to be improved - hardwired 100 |
| 505 | Double_t ptMin = 0.; // in [GeV/c] // to be improved - hardwired 0 |
| 506 | Double_t ptMax = 5.044; // in [GeV/c] // to be improved - hardwired 5.75 |
| 507 | TString xTitle = "p_{t} (GeV/#font[72]{c})"; |
| 508 | TString yTitle = "v_{2}"; |
| 509 | |
| 510 | TH1D* hist = new TH1D("","",nPtBins,ptMin,ptMax); |
| 511 | hist->SetXTitle(xTitle.Data()); |
| 512 | hist->SetYTitle(yTitle.Data()); |
| 513 | hist->GetYaxis()->SetLabelSize(0.06); |
| 514 | hist->GetYaxis()->SetTitleSize(0.07); |
| 515 | hist->GetYaxis()->SetTitleOffset(1.10); |
| 516 | hist->SetMinimum(0.00001); // minimum y-value - to be improved |
| 517 | hist->SetMaximum(0.34499); // maximum y-value - to be improved |
| 518 | //hist->SetLineStyle(7); |
| 519 | return hist; |
| 520 | |
| 521 | } // end of TH1D* StyleHistogramVsPtUp() |
| 522 | |
| 523 | // ===================================================================================== |
| 524 | |
| 525 | TH1D* StyleHistogramVsPtDown() |
| 526 | { |
| 527 | // Style histogram: |
| 528 | |
| 529 | Int_t nPtBins = 100; // to be improved - hardwired 100 |
| 530 | Double_t ptMin = 0.; // in [GeV/c] // to be improved - hardwired 0 |
| 531 | Double_t ptMax = 5.044; // in [GeV/c] // to be improved - hardwired 5.75 |
| 532 | TString xTitle = "p_{t} (GeV/#font[72]{c})"; |
| 533 | TString yTitle = "v_{2}{4}"; |
| 534 | |
| 535 | TH1D* hist = new TH1D("","",nPtBins,ptMin,ptMax); |
| 536 | hist->SetXTitle(xTitle.Data()); |
| 537 | hist->SetYTitle(yTitle.Data()); |
| 538 | hist->SetMinimum(0.00001); // minimum y-value - to be improved |
| 539 | hist->SetMaximum(0.29999); // maximum y-value - to be improved |
| 540 | //hist->SetLineStyle(7); |
| 541 | return hist; |
| 542 | |
| 543 | } // end of TH1D* StyleHistogramVsPtDown() |
| 544 | |
| 545 | // ===================================================================================== |
| 546 | |
| 547 | |
| 548 | void RemoveZeroPoints(TGraphErrors *ge) |
| 549 | { |
| 550 | // Remove zero points from TGraphErrors. |
| 551 | |
| 552 | if(!ge){return;} |
| 553 | |
| 554 | Int_t nPoints = ge->GetN(); |
| 555 | Double_t x = 0.; |
| 556 | Double_t y = 0.; |
| 557 | for(Int_t p=0;p<nPoints;p++) |
| 558 | { |
| 559 | ge->GetPoint(p,x,y); |
| 560 | if(TMath::Abs(y)<1.e-10) |
| 561 | { |
| 562 | ge->RemovePoint(p); |
| 563 | cout<<Form(" WARNING (%s): point %d is < 1.e-10 and it was removed from the plot !!!!",ge->GetName(),p+1)<<endl; |
| 564 | } |
| 565 | } // end of for(Int_t p=0;p<nPoints;p++) |
| 566 | |
| 567 | cout<<endl; |
| 568 | return; |
| 569 | |
| 570 | } // end of void RemoveZeroPoints(TGraphErrors *ge) |
| 571 | |
| 572 | // ===================================================================================== |
| 573 | |
| 574 | void ShiftAlongXaxis(TGraphErrors *ge, Double_t shift) |
| 575 | { |
| 576 | // Shift original TGraphErrors along x-axis by amount determined by 'shift'. |
| 577 | |
| 578 | if(!ge) |
| 579 | { |
| 580 | printf("\n WARNING: ge is NULL in ShiftAlongXaxis() !!!! \n\n"); |
| 581 | return; |
| 582 | } |
| 583 | |
| 584 | Int_t nPoints = ge->GetN(); |
| 585 | Double_t x = 0.; |
| 586 | Double_t y = 0.; |
| 587 | for(Int_t p=0;p<nPoints;p++) |
| 588 | { |
| 589 | ge->GetPoint(p,x,y); |
| 590 | x+=shift; |
| 591 | ge->SetPoint(p,x,y); |
| 592 | } // end of for(Int_t p=0;p<nPoints;p++) |
| 593 | |
| 594 | } // end of void ShiftAlongXaxis(TGraphErrors *ge, Double_t shift) |
| 595 | |
| 596 | // ===================================================================================== |
| 597 | |
| 598 | void SetFlowStyle() |
| 599 | { |
| 600 | // Set style which will affect all plots. |
| 601 | |
| 602 | gStyle->Reset(); |
| 603 | // gStyle->SetOptitle(0); |
| 604 | // gStyle->SetOptStat(0); |
| 605 | //gStyle->SetOptDate(1); |
| 606 | // gStyle->SetPalette(8,0); // (1,0) |
| 607 | gStyle->SetPalette(1); // (1,0) |
| 608 | gStyle->SetDrawBorder(0); |
| 609 | // gStyle->SetFillColor(0); // kills palete ??? |
| 610 | gStyle->SetCanvasColor(0); |
| 611 | gStyle->SetPadColor(0); |
| 612 | // gStyle->SetFillColor(0); // otherwize it affects Fill colors later |
| 613 | gStyle->SetFrameFillColor(0); |
| 614 | gStyle->SetCanvasBorderMode(0); |
| 615 | gStyle->SetFrameLineWidth(2); |
| 616 | // gStyle->SetFrameFillStyle(4000); |
| 617 | gStyle->SetPadBorderMode(0); |
| 618 | gStyle->SetPadTickX(1); |
| 619 | gStyle->SetPadTickY(1); |
| 620 | gStyle->SetPadBottomMargin(0.15); |
| 621 | gStyle->SetPadLeftMargin(0.15); |
| 622 | gStyle->SetHistLineWidth(2); |
| 623 | gStyle->SetFuncWidth(2); |
| 624 | gStyle->SetLineWidth(2); |
| 625 | gStyle->SetLabelSize(0.05,"xyz"); |
| 626 | gStyle->SetLabelOffset(0.01,"y"); |
| 627 | gStyle->SetLabelColor(kBlack,"xyz"); |
| 628 | gStyle->SetTitleSize(0.06,"xyz"); |
| 629 | gStyle->SetTitleOffset(1.3,"y"); |
| 630 | gStyle->SetTitleFillColor(0); |
| 631 | gStyle->SetLineWidth(2); |
| 632 | gStyle->SetHistLineColor(1); |
| 633 | gStyle->SetTextColor(1); |
| 634 | gStyle->SetTitleTextColor(1); |
| 635 | TGaxis::SetMaxDigits(4); |
| 636 | gStyle->SetOptStat(0); // removes stat. box from all histos |
| 637 | gROOT->ForceStyle(); |
| 638 | |
| 639 | } // end of void SetFlowStyle() |
| 640 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame