void mass_T1(){ TCanvas *c1 = new TCanvas("median_T1_FullFit","median_T1_FullFit", 200, 100, 600, 600); double number_of_energy_bounds = 50; double Emean[5][50] = { {0.}, {2.80471,3.92638,4.31826,4.59535,4.82018,5.01898,5.19626,5.35873,5.51289,5.65855,5.79894,5.9355,6.06924,6.20289,6.33568,6.46857,6.60092,6.73338,6.86745,7.00427,7.14206,7.28318,7.4269,7.5717,7.72009,7.87265,8.03008,8.19437,8.3653,8.5421,8.72564,8.91631,9.11745,9.33035,9.55376,9.78953,10.0411,10.3101,10.5982,10.9105,11.2523,11.6317,12.0547,12.5332,13.0802,13.7161,14.4876,15.4535,16.7678,18.7008}, {2.29742,2.91467,3.28068,3.57719,3.83316,4.06328,4.27425,4.47018,4.65492,4.83087,5.0016,5.16585,5.32485,5.48131,5.63527,5.78768,5.93924,6.09064,6.24192,6.39178,6.54267,6.69682,6.85372,7.01285,7.17595,7.34312,7.51377,7.68817,7.86925,8.0559,8.24968,8.45046,8.66086,8.88004,9.10943,9.35304,9.61309,9.88706,10.1835,10.5033,10.8493,11.2296,11.6467,12.1153,12.6539,13.2865,14.0471,15.0206,16.3894,18.494}, {1.62784,1.74471,1.82771,1.89934,1.96523,2.0271,2.08661,2.14502,2.20199,2.25748,2.3125,2.36744,2.42246,2.4774,2.53242,2.58788,2.64384,2.69986,2.75686,2.81532,2.87483,2.93527,2.99726,3.06077,3.12578,3.19277,3.26174,3.33267,3.40571,3.4807,3.55863,3.63959,3.72359,3.81155,3.9034,4.00044,4.10291,4.21069,4.325,4.44845,4.58229,4.72768,4.8887,5.07122,5.2809,5.53105,5.84709,6.2813,7.01908,9.79533}, {1.67516,1.83167,1.94067,2.03325,2.11606,2.19229,2.26522,2.33514,2.40265,2.46861,2.53306,2.59656,2.65902,2.72101,2.78346,2.84647,2.90944,2.97297,3.03743,3.10236,3.16836,3.2354,3.30288,3.37187,3.44237,3.51386,3.58683,3.66186,3.73937,3.81882,3.9008,3.98639,4.0748,4.16665,4.26275,4.3632,4.4691,4.582,4.70195,4.83035,4.97135,5.12514,5.2956,5.48988,5.71659,5.99072,6.33829,6.82753,7.67935,10.6342} }; double median[5][50] = { {0.}, {126.934,137.321,135.299,134.912,135.133,135.318,135.214,135.365,135.192,135.364,135.533,135.201,135.13,135.387,135.408,135.242,135.483,135.305,135.202,135.074,134.828,134.888,134.811,135.104,134.906,135.067,134.787,134.651,134.791,134.792,134.966,134.763,134.834,134.927,134.631,134.974,134.93,134.734,135.19,135.202,135.55,135.591,135.558,135.572,135.709,135.499,135.266,134.765,134.164,135.914}, {134.002,140.522,142.547,142.173,141.864,140.217,137.906,138.092,135.889,136.019,134.924,134.939,135.207,134.887,135.436,135.591,135.422,135.727,135.228,135.831,136.093,135.496,135.611,135.988,135.681,135.904,135.352,135.243,135.588,135.112,135.095,135.326,135.294,135.391,135.208,135.411,134.856,134.98,135.342,135.284,134.943,135.262,135.201,135.247,135.54,135.383,135.315,135.466,134.239,135.731}, {134.553,134.651,135.094,135.612,135.243,135.578,135.144,135.531,135.597,135.522,135.664,135.535,135.94,135.566,135.512,135.64,135.65,135.378,135.2,135.143,135.19,135.019,135.092,135.094,135.17,134.861,135.039,134.873,135.003,134.877,134.829,134.764,134.643,134.937,134.804,134.947,134.946,134.968,134.936,134.963,135.333,135.107,135.182,135.291,135.216,135.338,134.709,134.396,133.818,135.446}, {129.911,132.407,133.022,134.691,134.624,136.584,134.159,137.087,135.174,135.913,135.816,137.133,136.412,136.749,135.233,137.609,135.781,135.906,136.327,135.875,135.174,135.405,134.816,134.487,134.435,135.156,134.64,134.16,135.273,133.618,134.479,134.342,134.713,134.666,134.451,134.052,134.643,134.326,134.424,134.185,134.606,134.863,134.388,134.607,134.523,133.621,132.749,130.152,126.136,132.773} }; double median_Error[5][50] = { {0.}, {1.13159,0.287092,0.223377,0.189021,0.176366,0.17302,0.162633,0.15422,0.147725,0.171788,0.142454,0.14272,0.139349,0.141459,0.13487,0.135045,0.140518,0.132211,0.132111,0.128276,0.133715,0.127375,0.129223,0.125835,0.128363,0.125795,0.124352,0.123496,0.122798,0.120372,0.11695,0.11857,0.120818,0.119503,0.117381,0.115714,0.114017,0.115024,0.113313,0.115268,0.113117,0.113415,0.114484,0.116915,0.116473,0.120115,0.12354,0.123824,0.134119,0.148153}, {5,5,5,5,5,5,5,5,0.403427,0.372505,0.338748,0.345899,0.312533,0.293239,0.277104,0.280365,0.271901,0.272051,0.258336,0.243804,0.292383,0.243236,0.234446,0.226362,0.302695,0.222779,0.370137,0.213141,0.210636,0.21708,0.206494,0.203113,0.202106,0.196262,0.188554,0.187267,0.321769,0.251485,0.186146,0.216013,0.209829,0.305763,0.163389,0.167265,0.166272,0.163377,0.179215,0.181766,0.192578,0.250175}, {0.429163,0.386277,0.350372,0.321243,0.331839,0.31292,0.296106,0.269761,0.272137,0.501366,0.260827,0.246684,0.240125,0.230346,0.23401,0.229521,0.252999,0.207044,0.271316,0.201515,0.189008,0.184576,0.239848,0.176357,0.178422,0.169122,0.214328,0.160277,0.193249,0.151126,0.185365,0.1808,0.183472,0.171398,0.129372,0.127931,0.125759,0.119639,0.114681,0.116116,0.111249,0.110316,0.107245,0.105608,0.10236,0.102514,0.0987777,0.100013,0.100617,0.142933}, {0.906258,0.796369,0.829269,0.75501,0.667252,0.734101,0.664262,0.703425,0.658122,0.615963,0.668561,0.604411,0.646601,0.545502,0.560123,0.723382,0.599128,0.966132,1.69011,0.967526,1.02043,0.504325,0.735808,0.780305,0.502215,0.687845,0.454895,0.43745,0.461373,0.410648,0.397913,0.385691,0.406209,0.358387,0.350169,0.354825,0.33356,0.330613,0.30821,0.311457,0.301826,0.329717,0.296845,0.300918,0.322605,0.333141,0.343361,0.336338,0.383666,10} }; double median_unscaled[5][50] = { {0.}, {118.223,132.498,132.784,133.796,135.146,136.208,136.818,137.613,137.969,138.664,139.276,139.376,139.578,140.211,140.508,140.698,141.17,141.284,141.441,141.527,141.45,141.713,141.873,142.371,142.303,142.625,142.493,142.512,142.778,142.93,143.225,143.152,143.363,143.614,143.485,144.014,144.2,144.21,144.973,145.305,146.1,146.674,147.294,148.193,149.523,150.839,152.684,154.997,158.097,163.062}, {129.628,141.253,144.161,144.42,144.73,143.482,141.35,141.995,140.028,140.36,139.552,139.761,140.215,140.177,140.868,141.154,141.17,141.614,141.313,142.023,142.379,142.01,142.164,142.739,142.555,142.957,142.48,142.545,143.072,142.722,142.868,143.346,143.524,143.887,143.954,144.428,144.174,144.71,145.501,145.901,146.161,147.194,147.914,149.045,150.646,152.187,154.261,157.47,160.52,167.337}, {157.035,156.452,155.768,156.114,155.056,155.098,153.924,153.923,153.733,153.245,153.138,152.672,152.749,152.132,151.791,151.531,151.245,150.821,150.26,150.096,149.915,149.5,149.339,149.214,149.165,148.674,148.676,148.368,148.44,148.225,148.066,147.96,147.826,148.053,147.901,148.085,148.175,148.288,148.404,148.556,149.216,149.093,149.552,150.044,150.478,151.292,151.503,152.461,153.592,156.496}, {141.216,143.684,144.141,146.272,146.522,148.149,146.286,148.47,146.918,147.998,147.378,149.501,148.758,149.104,147.286,149.86,148.075,148.504,148.627,149.024,147.849,148.255,147.986,147.871,147.583,148.323,148.099,147.454,148.694,147.214,148.688,148.353,148.863,149.436,149.458,149.028,149.901,149.97,150.642,150.583,151.303,152.353,152.446,153.309,154.564,155.204,156.4,156.31,159.031,164.985} }; double Median_Error_unscaled[5][50] = { {0.}, {2.71501,0.276311,0.216149,0.191191,0.169004,0.165964,0.164761,0.165564,0.160481,0.152979,0.14651,0.144791,0.145385,0.144037,0.14033,0.139973,0.137314,0.139705,0.135944,0.134739,0.133848,0.135697,0.13537,0.132176,0.134152,0.131921,0.13249,0.130208,0.129008,0.12801,0.127073,0.126556,0.129594,0.124281,0.125538,0.123585,0.123858,0.122498,0.120145,0.122991,0.121109,0.123108,0.123647,0.127681,0.127698,0.128225,0.134186,0.13904,0.155404,0.182334}, {5,5,5,5,5,5,5,5,0.421548,0.377007,0.353346,0.358905,0.326309,0.30437,0.290729,0.286574,0.280249,0.281927,0.259749,0.254051,0.254026,0.254928,0.244184,0.237975,0.239696,0.236895,0.236108,0.224987,0.222911,0.230129,0.219001,0.218014,0.213831,0.207827,0.201406,0.198781,0.203197,0.198976,0.199104,0.191375,0.181498,0.188715,0.181357,0.185066,0.184806,0.191665,0.218492,0.222782,0.262959,0.437238}, {0.45784,0.438902,0.397275,0.37411,0.384381,0.347289,0.328722,0.31155,0.303213,0.303583,0.308329,0.278569,0.271793,0.26268,0.259974,0.255611,0.236465,0.244853,0.234295,0.2303,0.224685,0.211654,0.209112,0.200843,0.207574,0.193339,0.19161,0.180786,0.176594,0.171317,0.16634,0.160743,0.157105,0.150549,0.147738,0.148265,0.142739,0.13481,0.132337,0.130672,0.126542,0.124088,0.120039,0.119858,0.119363,0.11692,0.11691,0.118954,0.122145,0.167564}, {0.907181,0.813054,0.817571,0.774266,0.711823,0.697419,0.76054,0.737793,0.690691,0.652779,0.694117,0.668644,0.689009,0.597388,0.624854,0.816489,0.690194,0.620593,0.680406,0.631305,0.613342,0.60539,0.58603,0.660159,0.623975,0.557819,0.569879,0.522793,0.564598,0.526309,0.517058,0.456479,0.506818,0.465329,0.444835,0.420356,0.399001,0.39098,0.371708,0.370448,0.369952,0.391234,0.366304,0.366291,0.391154,0.399652,0.438672,0.429727,0.433934,10} }; double e[50] = { 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0. }; gStyle->SetOptFit(1111); gStyle->SetStatY(0); gStyle->SetStatX(0); gStyle->SetStatW(0); gStyle->SetStatH(0); TGraphErrors *g1 = new TGraphErrors(number_of_energy_bounds, Emean[1], median[1], e, median_Error[1]); g1->SetName("g1"); g1->SetMarkerStyle(8); g1->SetMarkerColor(9); g1->Fit("pol5","","",4.0,17.0); g1->GetFunction("pol5")->SetLineColor(9); TGraphErrors *g5 = new TGraphErrors(number_of_energy_bounds, Emean[1], median_unscaled[1], e, Median_Error_unscaled[1]); g5->SetName("g5"); g5->SetMarkerStyle(8); g5->Fit("pol4","","",4.0,17.0); TMultiGraph *g = new TMultiGraph(); g->Add(g1); g->Add(g5); g->Draw("AP"); g->GetXaxis()->SetLimits(0., 22.); g->SetMinimum(100.); g->SetMaximum(180.); g->SetTitle("median_FullFit"); g->GetXaxis()->SetTitle("#pi^{0} Energy [GeV]"); g->GetYaxis()->SetTitle("Median mass [MeV]"); auto legend = new TLegend(0.55, 0.8, 0.9, 0.9); legend->AddEntry("g1", "Scaled full mass spectra fit median T1", "p,E"); legend->AddEntry("g5", "Unscaled full mass spectra fit median T1", "p,E"); legend->Draw(); TLine *line = new TLine(0, 135, 100, 135); line->SetLineStyle(2); line->Draw(); c1->SetGrid(1, 1); c1->Update(); return c1; }