#include "TSystem.h"
#include "TROOT.h"
#include "TFile.h"
#include "TH1.h"
#include "TH1D.h"
#include "TCanvas.h"
#include "TPad.h"
#include "TProfile.h"
#include "THistPainter.h"
#include "TPaveStats.h"
#include "THStack.h"
#include "TList.h"
#include "TDirectory.h"
#include "THashTable.h"
#include "TStyle.h"
#include "TLegend.h"

#include <math.h>
#include <iostream>
#include <vector>
#include <stdio.h>
#include <fstream>
#include <stdlib.h>

   
void read(){

        static const int nBCID = 3564; 

              
	//histograms of <f>_bcid - f_ped
        TH1D* h1_firCorrection_0f1 = new TH1D("h1_firCorrection_0f1","firCorrection_0FIR_{bcid}-FIR_{ped}",2548,-500,2048); 
	TH1D* h1_firCorrection_0f2 = new TH1D("h1_firCorrection_0f2","firCorrection_0FIR_{bcid}-FIR_{ped}",2548,-500,2048);


   h1_firCorrection_0f1->SetBinContent(151,1);
   h1_firCorrection_0f1->SetBinContent(152,1);
   h1_firCorrection_0f1->SetBinContent(154,1);
   h1_firCorrection_0f1->SetBinContent(155,1);
   h1_firCorrection_0f1->SetBinContent(156,1);
   h1_firCorrection_0f1->SetBinContent(158,2);
   h1_firCorrection_0f1->SetBinContent(159,3);
   h1_firCorrection_0f1->SetBinContent(160,9);
   h1_firCorrection_0f1->SetBinContent(161,5);
   h1_firCorrection_0f1->SetBinContent(162,2);
   h1_firCorrection_0f1->SetBinContent(163,2);
   h1_firCorrection_0f1->SetBinContent(164,3);
   h1_firCorrection_0f1->SetBinContent(165,2);
   h1_firCorrection_0f1->SetBinContent(166,2);
   h1_firCorrection_0f1->SetBinContent(167,2);
   h1_firCorrection_0f1->SetBinContent(168,1);
   h1_firCorrection_0f1->SetBinContent(372,1);
   h1_firCorrection_0f1->SetBinContent(376,1);
   h1_firCorrection_0f1->SetBinContent(378,1);
   h1_firCorrection_0f1->SetBinContent(379,2);
   h1_firCorrection_0f1->SetBinContent(380,5);
   h1_firCorrection_0f1->SetBinContent(381,2);
   h1_firCorrection_0f1->SetBinContent(382,1);
   h1_firCorrection_0f1->SetBinContent(383,1);
   h1_firCorrection_0f1->SetBinContent(384,4);
   h1_firCorrection_0f1->SetBinContent(385,3);
   h1_firCorrection_0f1->SetBinContent(386,3);
   h1_firCorrection_0f1->SetBinContent(387,2);
   h1_firCorrection_0f1->SetBinContent(388,2);
   h1_firCorrection_0f1->SetBinContent(389,3);
   h1_firCorrection_0f1->SetBinContent(390,3);
   h1_firCorrection_0f1->SetBinContent(392,1);
   h1_firCorrection_0f1->SetBinContent(393,2);
   h1_firCorrection_0f1->SetBinContent(397,1);
   h1_firCorrection_0f1->SetBinContent(425,2);
   h1_firCorrection_0f1->SetBinContent(429,2);
   h1_firCorrection_0f1->SetBinContent(430,3);
   h1_firCorrection_0f1->SetBinContent(431,2);
   h1_firCorrection_0f1->SetBinContent(432,6);
   h1_firCorrection_0f1->SetBinContent(433,5);
   h1_firCorrection_0f1->SetBinContent(434,11);
   h1_firCorrection_0f1->SetBinContent(435,18);
   h1_firCorrection_0f1->SetBinContent(436,16);
   h1_firCorrection_0f1->SetBinContent(437,26);
   h1_firCorrection_0f1->SetBinContent(438,33);
   h1_firCorrection_0f1->SetBinContent(439,32);
   h1_firCorrection_0f1->SetBinContent(440,52);
   h1_firCorrection_0f1->SetBinContent(441,54);
   h1_firCorrection_0f1->SetBinContent(442,80);
   h1_firCorrection_0f1->SetBinContent(443,82);
   h1_firCorrection_0f1->SetBinContent(444,114);
   h1_firCorrection_0f1->SetBinContent(445,117);
   h1_firCorrection_0f1->SetBinContent(446,145);
   h1_firCorrection_0f1->SetBinContent(447,151);
   h1_firCorrection_0f1->SetBinContent(448,133);
   h1_firCorrection_0f1->SetBinContent(449,149);
   h1_firCorrection_0f1->SetBinContent(450,144);
   h1_firCorrection_0f1->SetBinContent(451,120);
   h1_firCorrection_0f1->SetBinContent(452,128);
   h1_firCorrection_0f1->SetBinContent(453,117);
   h1_firCorrection_0f1->SetBinContent(454,94);
   h1_firCorrection_0f1->SetBinContent(455,94);
   h1_firCorrection_0f1->SetBinContent(456,74);
   h1_firCorrection_0f1->SetBinContent(457,72);
   h1_firCorrection_0f1->SetBinContent(458,52);
   h1_firCorrection_0f1->SetBinContent(459,51);
   h1_firCorrection_0f1->SetBinContent(460,26);
   h1_firCorrection_0f1->SetBinContent(461,23);
   h1_firCorrection_0f1->SetBinContent(462,15);
   h1_firCorrection_0f1->SetBinContent(463,15);
   h1_firCorrection_0f1->SetBinContent(464,8);
   h1_firCorrection_0f1->SetBinContent(465,6);
   h1_firCorrection_0f1->SetBinContent(466,8);
   h1_firCorrection_0f1->SetBinContent(467,5);
   h1_firCorrection_0f1->SetBinContent(468,1);
   h1_firCorrection_0f1->SetBinContent(469,5);
   h1_firCorrection_0f1->SetBinContent(470,4);
   h1_firCorrection_0f1->SetBinContent(471,3);
   h1_firCorrection_0f1->SetBinContent(472,2);
   h1_firCorrection_0f1->SetBinContent(473,2);
   h1_firCorrection_0f1->SetBinContent(474,2);
   h1_firCorrection_0f1->SetBinContent(475,3);
   h1_firCorrection_0f1->SetBinContent(477,3);
   h1_firCorrection_0f1->SetBinContent(479,3);
   h1_firCorrection_0f1->SetBinContent(480,2);
   h1_firCorrection_0f1->SetBinContent(481,1);
   h1_firCorrection_0f1->SetBinContent(482,1);
   h1_firCorrection_0f1->SetBinContent(503,1);
   h1_firCorrection_0f1->SetBinContent(504,2);
   h1_firCorrection_0f1->SetBinContent(505,1);
   h1_firCorrection_0f1->SetBinContent(506,1);
   h1_firCorrection_0f1->SetBinContent(509,3);
   h1_firCorrection_0f1->SetBinContent(510,1);
   h1_firCorrection_0f1->SetBinContent(511,4);
   h1_firCorrection_0f1->SetBinContent(512,4);
   h1_firCorrection_0f1->SetBinContent(513,3);
   h1_firCorrection_0f1->SetBinContent(514,3);
   h1_firCorrection_0f1->SetBinContent(515,2);
   h1_firCorrection_0f1->SetBinContent(516,2);
   h1_firCorrection_0f1->SetBinContent(518,4);
   h1_firCorrection_0f1->SetBinContent(519,4);
   h1_firCorrection_0f1->SetBinContent(521,1);
   h1_firCorrection_0f1->SetBinContent(525,2);
   h1_firCorrection_0f1->SetBinContent(574,1);
   h1_firCorrection_0f1->SetBinContent(575,3);
   h1_firCorrection_0f1->SetBinContent(577,2);
   h1_firCorrection_0f1->SetBinContent(578,3);
   h1_firCorrection_0f1->SetBinContent(579,2);
   h1_firCorrection_0f1->SetBinContent(580,3);
   h1_firCorrection_0f1->SetBinContent(581,1);
   h1_firCorrection_0f1->SetBinContent(582,1);
   h1_firCorrection_0f1->SetBinContent(583,2);
   h1_firCorrection_0f1->SetBinContent(584,1);
   h1_firCorrection_0f1->SetBinContent(585,2);
   h1_firCorrection_0f1->SetBinContent(586,1);
   h1_firCorrection_0f1->SetBinContent(588,4);
   h1_firCorrection_0f1->SetBinContent(589,3);
   h1_firCorrection_0f1->SetBinContent(590,2);
   h1_firCorrection_0f1->SetBinContent(592,1);
   h1_firCorrection_0f1->SetBinContent(593,2);
   h1_firCorrection_0f1->SetBinContent(594,2);
   h1_firCorrection_0f1->SetBinContent(595,1);
   h1_firCorrection_0f1->SetBinContent(598,1);
   h1_firCorrection_0f1->SetBinContent(660,1);
   h1_firCorrection_0f1->SetBinContent(661,3);
   h1_firCorrection_0f1->SetBinContent(662,1);
   h1_firCorrection_0f1->SetBinContent(663,1);
   h1_firCorrection_0f1->SetBinContent(664,1);
   h1_firCorrection_0f1->SetBinContent(665,3);
   h1_firCorrection_0f1->SetBinContent(666,1);
   h1_firCorrection_0f1->SetBinContent(667,2);
   h1_firCorrection_0f1->SetBinContent(668,2);
   h1_firCorrection_0f1->SetBinContent(669,2);
   h1_firCorrection_0f1->SetBinContent(670,2);
   h1_firCorrection_0f1->SetBinContent(671,1);
   h1_firCorrection_0f1->SetBinContent(672,4);
   h1_firCorrection_0f1->SetBinContent(673,1);
   h1_firCorrection_0f1->SetBinContent(674,1);
   h1_firCorrection_0f1->SetBinContent(675,2);
   h1_firCorrection_0f1->SetBinContent(676,2);
   h1_firCorrection_0f1->SetBinContent(679,2);
   h1_firCorrection_0f1->SetBinContent(680,1);
   h1_firCorrection_0f1->SetBinContent(681,3);
   h1_firCorrection_0f1->SetBinContent(684,1);
   h1_firCorrection_0f1->SetBinContent(687,1);
   h1_firCorrection_0f1->SetBinContent(757,1);
   h1_firCorrection_0f1->SetBinContent(761,3);
   h1_firCorrection_0f1->SetBinContent(762,6);
   h1_firCorrection_0f1->SetBinContent(763,1);
   h1_firCorrection_0f1->SetBinContent(764,1);
   h1_firCorrection_0f1->SetBinContent(765,1);
   h1_firCorrection_0f1->SetBinContent(766,4);
   h1_firCorrection_0f1->SetBinContent(767,3);
   h1_firCorrection_0f1->SetBinContent(768,1);
   h1_firCorrection_0f1->SetBinContent(769,2);
   h1_firCorrection_0f1->SetBinContent(770,3);
   h1_firCorrection_0f1->SetBinContent(771,2);
   h1_firCorrection_0f1->SetBinContent(773,1);
   h1_firCorrection_0f1->SetBinContent(774,2);
   h1_firCorrection_0f1->SetBinContent(775,3);
   h1_firCorrection_0f1->SetBinContent(776,1);
   h1_firCorrection_0f1->SetBinContent(777,1);
   h1_firCorrection_0f1->SetBinContent(780,2);
   h1_firCorrection_0f1->SetBinContent(896,1);
   h1_firCorrection_0f1->SetBinContent(899,3);
   h1_firCorrection_0f1->SetBinContent(900,2);
   h1_firCorrection_0f1->SetBinContent(901,3);
   h1_firCorrection_0f1->SetBinContent(902,1);
   h1_firCorrection_0f1->SetBinContent(903,2);
   h1_firCorrection_0f1->SetBinContent(904,1);
   h1_firCorrection_0f1->SetBinContent(906,3);
   h1_firCorrection_0f1->SetBinContent(907,3);
   h1_firCorrection_0f1->SetBinContent(908,5);
   h1_firCorrection_0f1->SetBinContent(909,8);
   h1_firCorrection_0f1->SetBinContent(910,1);
   h1_firCorrection_0f1->SetBinContent(911,1);
   h1_firCorrection_0f1->SetBinContent(913,1);
   h1_firCorrection_0f1->SetBinContent(914,1);
   h1_firCorrection_0f1->SetBinContent(916,1);
   h1_firCorrection_0f1->SetBinContent(917,1);
   h1_firCorrection_0f1->SetBinContent(1119,1);
   h1_firCorrection_0f1->SetBinContent(1123,1);
   h1_firCorrection_0f1->SetBinContent(1124,1);
   h1_firCorrection_0f1->SetBinContent(1125,2);
   h1_firCorrection_0f1->SetBinContent(1126,1);
   h1_firCorrection_0f1->SetBinContent(1127,3);
   h1_firCorrection_0f1->SetBinContent(1128,5);
   h1_firCorrection_0f1->SetBinContent(1129,6);
   h1_firCorrection_0f1->SetBinContent(1130,1);
   h1_firCorrection_0f1->SetBinContent(1131,2);
   h1_firCorrection_0f1->SetBinContent(1132,3);
   h1_firCorrection_0f1->SetBinContent(1133,2);
   h1_firCorrection_0f1->SetBinContent(1134,1);
   h1_firCorrection_0f1->SetBinContent(1135,3);
   h1_firCorrection_0f1->SetBinContent(1136,1);
   h1_firCorrection_0f1->SetBinContent(1137,2);
   h1_firCorrection_0f1->SetBinContent(1138,2);
   h1_firCorrection_0f1->SetBinContent(1139,1);
   h1_firCorrection_0f1->SetBinContent(1171,1);
   h1_firCorrection_0f1->SetBinContent(1172,1);
   h1_firCorrection_0f1->SetBinContent(1173,1);
   h1_firCorrection_0f1->SetBinContent(1174,2);
   h1_firCorrection_0f1->SetBinContent(1175,2);
   h1_firCorrection_0f1->SetBinContent(1177,3);
   h1_firCorrection_0f1->SetBinContent(1178,1);
   h1_firCorrection_0f1->SetBinContent(1179,1);
   h1_firCorrection_0f1->SetBinContent(1180,7);
   h1_firCorrection_0f1->SetBinContent(1181,1);
   h1_firCorrection_0f1->SetBinContent(1182,2);
   h1_firCorrection_0f1->SetBinContent(1184,1);
   h1_firCorrection_0f1->SetBinContent(1185,1);
   h1_firCorrection_0f1->SetBinContent(1188,2);
   h1_firCorrection_0f1->SetBinContent(1189,1);
   h1_firCorrection_0f1->SetBinContent(1237,2);
   h1_firCorrection_0f1->SetBinContent(1238,1);
   h1_firCorrection_0f1->SetBinContent(1241,2);
   h1_firCorrection_0f1->SetBinContent(1242,2);
   h1_firCorrection_0f1->SetBinContent(1243,2);
   h1_firCorrection_0f1->SetBinContent(1246,1);
   h1_firCorrection_0f1->SetBinContent(1250,1);
   h1_firCorrection_0f1->SetBinContent(1370,2);
   h1_firCorrection_0f1->SetBinContent(1373,3);
   h1_firCorrection_0f1->SetBinContent(1374,3);
   h1_firCorrection_0f1->SetBinContent(1375,2);
   h1_firCorrection_0f1->SetBinContent(1376,2);
   h1_firCorrection_0f1->SetBinContent(1377,4);
   h1_firCorrection_0f1->SetBinContent(1378,4);
   h1_firCorrection_0f1->SetBinContent(1379,2);
   h1_firCorrection_0f1->SetBinContent(1380,1);
   h1_firCorrection_0f1->SetBinContent(1381,3);
   h1_firCorrection_0f1->SetBinContent(1382,1);
   h1_firCorrection_0f1->SetBinContent(1383,2);
   h1_firCorrection_0f1->SetBinContent(1384,4);
   h1_firCorrection_0f1->SetBinContent(1386,1);
   h1_firCorrection_0f1->SetBinContent(1387,2);
   h1_firCorrection_0f1->SetBinContent(1388,2);
   h1_firCorrection_0f1->SetBinContent(1436,1);
   h1_firCorrection_0f1->SetBinContent(1438,2);
   h1_firCorrection_0f1->SetBinContent(1439,1);
   h1_firCorrection_0f1->SetBinContent(1440,1);
   h1_firCorrection_0f1->SetBinContent(1441,2);
   h1_firCorrection_0f1->SetBinContent(1443,2);
   h1_firCorrection_0f1->SetBinContent(1444,4);
   h1_firCorrection_0f1->SetBinContent(1446,4);
   h1_firCorrection_0f1->SetBinContent(1447,2);
   h1_firCorrection_0f1->SetBinContent(1448,1);
   h1_firCorrection_0f1->SetBinContent(1449,1);
   h1_firCorrection_0f1->SetBinContent(1450,1);
   h1_firCorrection_0f1->SetBinContent(1451,2);
   h1_firCorrection_0f1->SetBinContent(1453,2);
   h1_firCorrection_0f1->SetBinContent(1456,1);
   h1_firCorrection_0f1->SetBinContent(1458,1);
   h1_firCorrection_0f1->SetBinContent(1461,1);
   h1_firCorrection_0f1->SetBinContent(1462,2);
   h1_firCorrection_0f1->SetBinContent(1463,1);
   h1_firCorrection_0f1->SetBinContent(1466,3);
   h1_firCorrection_0f1->SetBinContent(1468,1);
   h1_firCorrection_0f1->SetBinContent(1469,1);
   h1_firCorrection_0f1->SetBinContent(1474,1);
   h1_firCorrection_0f1->SetEntries(2736);

   h1_firCorrection_0f2->SetBinContent(504,1);
   h1_firCorrection_0f2->SetBinContent(505,1);
   h1_firCorrection_0f2->SetBinContent(506,2);
   h1_firCorrection_0f2->SetBinContent(507,18);
   h1_firCorrection_0f2->SetBinContent(508,84);
   h1_firCorrection_0f2->SetBinContent(509,229);
   h1_firCorrection_0f2->SetBinContent(510,444);
   h1_firCorrection_0f2->SetBinContent(511,559);
   h1_firCorrection_0f2->SetBinContent(512,485);
   h1_firCorrection_0f2->SetBinContent(513,269);
   h1_firCorrection_0f2->SetBinContent(514,138);
   h1_firCorrection_0f2->SetBinContent(515,61);
   h1_firCorrection_0f2->SetBinContent(516,42);
   h1_firCorrection_0f2->SetBinContent(517,12);
   h1_firCorrection_0f2->SetBinContent(518,6);
   h1_firCorrection_0f2->SetBinContent(519,5);
   h1_firCorrection_0f2->SetBinContent(524,2);
   h1_firCorrection_0f2->SetBinContent(525,4);
   h1_firCorrection_0f2->SetBinContent(526,10);
   h1_firCorrection_0f2->SetBinContent(527,6);
   h1_firCorrection_0f2->SetBinContent(528,9);
   h1_firCorrection_0f2->SetBinContent(529,11);
   h1_firCorrection_0f2->SetBinContent(530,5);
   h1_firCorrection_0f2->SetBinContent(531,10);
   h1_firCorrection_0f2->SetBinContent(532,6);
   h1_firCorrection_0f2->SetBinContent(533,10);
   h1_firCorrection_0f2->SetBinContent(534,11);
   h1_firCorrection_0f2->SetBinContent(535,5);
   h1_firCorrection_0f2->SetBinContent(536,9);
   h1_firCorrection_0f2->SetBinContent(537,6);
   h1_firCorrection_0f2->SetBinContent(538,8);
   h1_firCorrection_0f2->SetBinContent(539,8);
   h1_firCorrection_0f2->SetBinContent(540,11);
   h1_firCorrection_0f2->SetBinContent(541,4);
   h1_firCorrection_0f2->SetBinContent(542,7);
   h1_firCorrection_0f2->SetBinContent(543,6);
   h1_firCorrection_0f2->SetBinContent(545,4);
   h1_firCorrection_0f2->SetBinContent(548,3);
   h1_firCorrection_0f2->SetBinContent(549,3);
   h1_firCorrection_0f2->SetBinContent(550,6);
   h1_firCorrection_0f2->SetBinContent(551,9);
   h1_firCorrection_0f2->SetBinContent(552,13);
   h1_firCorrection_0f2->SetBinContent(553,16);
   h1_firCorrection_0f2->SetBinContent(554,14);
   h1_firCorrection_0f2->SetBinContent(555,8);
   h1_firCorrection_0f2->SetBinContent(556,1);
   h1_firCorrection_0f2->SetBinContent(557,1);
   h1_firCorrection_0f2->SetBinContent(558,1);
   h1_firCorrection_0f2->SetBinContent(559,2);
   h1_firCorrection_0f2->SetBinContent(561,4);
   h1_firCorrection_0f2->SetBinContent(562,4);
   h1_firCorrection_0f2->SetBinContent(563,7);
   h1_firCorrection_0f2->SetBinContent(564,8);
   h1_firCorrection_0f2->SetBinContent(565,13);
   h1_firCorrection_0f2->SetBinContent(566,4);
   h1_firCorrection_0f2->SetBinContent(567,15);
   h1_firCorrection_0f2->SetBinContent(568,13);
   h1_firCorrection_0f2->SetBinContent(569,14);
   h1_firCorrection_0f2->SetBinContent(570,13);
   h1_firCorrection_0f2->SetBinContent(571,9);
   h1_firCorrection_0f2->SetBinContent(572,4);
   h1_firCorrection_0f2->SetBinContent(573,6);
   h1_firCorrection_0f2->SetBinContent(574,3);
   h1_firCorrection_0f2->SetBinContent(575,4);
   h1_firCorrection_0f2->SetBinContent(576,4);
   h1_firCorrection_0f2->SetBinContent(577,7);
   h1_firCorrection_0f2->SetBinContent(578,4);
   h1_firCorrection_0f2->SetBinContent(579,3);
   h1_firCorrection_0f2->SetBinContent(581,2);
   h1_firCorrection_0f2->SetBinContent(583,2);
   h1_firCorrection_0f2->SetBinContent(584,3);
   h1_firCorrection_0f2->SetBinContent(585,3);
   h1_firCorrection_0f2->SetBinContent(586,1);
   h1_firCorrection_0f2->SetBinContent(587,1);
   h1_firCorrection_0f2->SetEntries(2736);
	

	THStack *hs = new THStack("hs", "<f>_bcid - f_ped Histograms for matched and autocorrelation modes");

	h1_firCorrection_0f1->SetLineColor(1);
	hs->Add(h1_firCorrection_0f1,"sames");
	h1_firCorrection_0f2->SetLineColor(2);
	hs->Add(h1_firCorrection_0f2,"sames");

	//Get bin position of peak (maximum value in histogram)
	Int_t binmaxf1 = h1_firCorrection_0f1->GetMaximumBin();//bin position
	TString xf1 = TString::Itoa(h1_firCorrection_0f1->GetXaxis()->GetBinCenter(binmaxf1),10);//value position
	Int_t binmaxf2 = h1_firCorrection_0f2->GetMaximumBin();
	TString xf2 = TString::Itoa(h1_firCorrection_0f2->GetXaxis()->GetBinCenter(binmaxf2),10);//value position

	//Plot of firCorrection for f1 and f2, both in a same canvas

	TCanvas *c1 = new TCanvas("c1", "FirCorrection <f>_bcid-f_ped Histograms", 900, 700);
	c1->cd();

	//c1->SetLogy(1);
	hs->Draw("ehist");
	c1->Modified();
	c1->Update();

	TPaveStats *st1 = (TPaveStats*)hs->GetStack()->FindObject("h1_firCorrection_0f1")->FindObject("stats");
	TPaveStats *st2 = (TPaveStats*)hs->GetStack()->FindObject("h1_firCorrection_0f2")->FindObject("stats");
	h1_firCorrection_0f1->SetBit(TH1::kNoStats);
	h1_firCorrection_0f2->SetBit(TH1::kNoStats);
	st1->AddText("Xpeakf1	"+xf1); 	
	st2->AddText("Xpeakf2	"+xf2);
	st1->Print();	
	st2->Print();
	st1->SetX1NDC(.77); st1->SetX2NDC(.97);
	st1->SetY1NDC(.75); st1->SetY2NDC(.95);
	st2->SetX1NDC(.77); st2->SetX2NDC(.97);
	st2->SetY1NDC(.55); st2->SetY2NDC(.75);
	st1->SetOptStat(1111);
	st2->SetOptStat(1111);
	st1->Draw("same");
	st2->Draw("same");
	c1->Modified();
	c1->Update();

}