{
#include <fstream.h>
gROOT->LoadMacro("AtlasUtils.C");
ATLAS->SetPadTickX(1);
ATLAS->SetPadTickY(0);
ATLAS->SetPadRightMargin(0.12);
ATLAS->SetPadLeftMargin(0.12);
c1 = new TCanvas("c1","A Simple Graph",600,10,700,500);
//c1->SetFillStyle(4100);
//TDatime db(2010,01,01,00,00,00);
TDatime db(2010,01,01,00,00,00);
int t1 = db.Convert();
//TDatime db(2011,04,01,00,00,00);
TDatime db(2021,12,31,00,00,00);
int t2 = db.Convert();
gStyle->SetOptStat(0);
//c1->SetGrid();
// create the arrays for the points
Double_t ymax= 350;
double Tscale=ymax*0.1/20., Toffset = ymax*0.8;
//double Tscale=0.01/10., Toffset = 0.02;
Int_t n = 3000;
double intL[n],Temp[4][n],Temp_plot[4][n], Tbias=3.6;
double time[n], Ttime[n];
double weekL[n], energy[n], weekLL[n];
double leak[4][n],dleak[4][n],leakop[4][n],dleakop[4][n],
leak3[5][n],dleak3[5][n];
int week,year1,month1,day1,year2,month2,day2; 
//--------readin luminosity data-------------------
int iweek=1, iT=0;
ifstream fin("lhcall_7_intL.txt");
while(fin>>year1>>month1>>day1>>year2>>month2>>day2>>
weekL[iweek]>>intL[iweek] )
{
//cout<<year1<<" "<<month1<<" "<<day1<<" "<<year2<<" "
//<<month2<<" "<<day2<<" "<<weekL[iweek]<<" "<<intL[iweek]<<endl;
  intL[iweek]= intL[iweek]/1000000.; 
  TDatime d1(year1,month1,day1,12,0,0);
  int time1=d1.Convert();
  TDatime d2(year2,month2,day2,12,0,0);
  time[iweek]=d2.Convert();
  Ttime[iT]=time1;
  weekLL[iT] = weekL[iweek];
  iT++;
  Ttime[iT]=time[iweek];
  weekLL[iT] = weekL[iweek];
  iT++;
  iweek++;
}
fin.close();
cout<<"number of time intervals= "<< iweek << endl;
cout<<"last week: from "<<year1<<" "<<month1<<" "<<day1<<" to "<<
year2<<" "<<month2<<" "<<day2<< endl;



//-----------------------------------------------------
//
leg0 = new TLegend(0.20 ,0.55 ,0.55,0.75);
for (int iw=0; iw<iT; iw++) weekLL[iw]/=50000.;

  hu = new TGraph(iweek,time,intL);
  hu->SetMaximum(ymax);

  hu->GetYaxis()->SetTitle("Integrated luminosity (fb^{-1})");
  //hu->GetXaxis()->SetTitle("DATE");
  hu->GetXaxis()->SetTitleOffset(1.4);
  hu->GetXaxis()->SetLabelOffset(0.045);
  hu->GetXaxis()->SetLabelSize(0.04);
  hu->GetYaxis()->SetTitleSize(0.05);
  hu->GetYaxis()->SetTitleOffset(1.2);
  hu->GetYaxis()->SetLabelOffset(0.01);
//  hu->GetXaxis()->SetRangeUser(13,45);
  hu->GetXaxis()->SetTimeDisplay(1);
  hu->GetXaxis()->SetTimeFormat("#splitline{  %Y}{%b %d}");
  hu->GetXaxis()->SetLimits(t1,t2);
//  hu->GetXaxis()->SetNdivisions(206);
  hu->GetXaxis()->SetNdivisions(206);
  hu->Draw("AL");
leg0->AddEntry(hu,"Lum/week","l");
gL = new TGraph(iT,Ttime,weekLL);
gL->SetLineStyle(1); gL->SetLineColor(1);gL->Draw("SAME LF");
leg0->AddEntry(gL,"int-Lum (fb^{-1})","l");
//
leg0->SetFillColor(10); leg0->SetBorderSize(0); leg0->SetTextSize(0.04);
//leg0->Draw();
//draw an axis on the right side
TGaxis *axis = new TGaxis(t2,0.,t2, 120,0.,6.,205,"+L");
axis->SetLabelOffset(0.02);
axis->SetTitleOffset(1.0);
axis->SetTitleSize(0.05);
axis->SetLabelSize(0.05);
axis->SetLabelFont(42);
axis->SetTitleFont(42);
axis->Draw();
TDatime db(2022,07,15,00,00,00);
int t4 = db.Convert();
legA = new TLatex(t4 ,140 ,"Luminosity (fb^{-1}/week)");
legA->SetTextSize(0.045);
legA->SetTextAngle(90);
legA->Draw();
//-- write last day
TDatime db(2021,09,01,00,00,00);int t4 = db.Convert();
legC = new TText(t4 ,-24. ,"2022");
legC->SetTextSize(0.040);legC->Draw();
TDatime db(2021,08,01,00,00,00);int t4 = db.Convert();
legB = new TText(t4 ,-42 ,"Jan 01");
legB->SetTextSize(0.040);legB->Draw();
//-- write start day
TDatime db(2009,09,01,00,00,00);int t4 = db.Convert();
legB = new TText(t4 ,-24 ,"2010");
legB->SetTextSize(0.040);legB->Draw();
TDatime db(2009,08,01,00,00,00);int t4 = db.Convert();
legC = new TText(t4 ,-42 ,"Jan 01");
legC->SetTextSize(0.040);legC->Draw();
//-- write energy
TDatime db(2011,01,01,00,00,00);int t4 = db.Convert();
legB = new TText(t4 ,320 ,"7 TeV");
legB->SetTextSize(0.040);legB->Draw();
TDatime db(2018,01,01,00,00,00);int t4 = db.Convert();
legB = new TText(t4 ,320 ,"14 TeV");
legB->SetTextSize(0.040);legB->Draw();
// arrows
TDatime db(2010,01,01,00,00,00);int t4 = db.Convert();
TDatime db(2013,01,01,00,00,00);int t5 = db.Convert();
ar = new TArrow(t4,315.,t5,315.,0.02,"<>");
ar->SetArrowSize(0.015);ar->SetLineWidth(1);ar->Draw(); 
TDatime db(2015,01,01,00,00,00);int t4 = db.Convert();
TDatime db(2021,12,31,00,00,00);int t5 = db.Convert();
ar = new TArrow(t4,315.,t5,315.,0.02,"<>");
ar->SetArrowSize(0.015);ar->SetLineWidth(1);ar->Draw();

c1->Print("Luminsity_profile.png");
c1->Print("Luminsity_profile.pdf");
}