/**
This Macro   
1. sums hits of different tornado plots corresponding to different pixels 
2. calculate the ratio of this sum over the number of triggers * number of tornado plots
   this ratio is 1 for values of CalDel/VcThr that give a fully efficienct response
The aim is to find those values of CalDel/VcThr

Need to specify
1. See "Declare Constants"

Notes
*/
/////
//   To run: root -l -q -b goodCalDelVcThr.cc+ > goodValues.txt
/////
/////
//   Prepare Root and Roofit
/////
#include "TFile.h"
#include "TH1F.h"
#include "TH2F.h"
#include "THStack.h"
#include "TLegend.h"
#include "TTree.h"
#include "TTreePlayer.h"
#include "TStyle.h"
#include "TGaxis.h"
#include "TPaveLabel.h"
#include "TCanvas.h"
#include "TGraphErrors.h"
#include "TGraphAsymmErrors.h"
#include "TEfficiency.h"
#include "TBranch.h"
#include "TLatex.h"
#include "TROOT.h"
#include "Math/DistFunc.h"
#include <fstream>
#include <iostream>
#include <iomanip>
#include <string>
using namespace std;
/////
//   Declare constants
/////
const string channel = "mFar21";
const string pathFile = "/cmsnfsbrilpro/brilpro/plt/cmsplt/interface/trimCalib/"+channel+"/tornado/";
const string histName_ = "DACDAC_FEDChannel_10_ROC_0";
const char *filesName[] = {
"fullDetectorDACDAC_col13row20.root",
"fullDetectorDACDAC_col26row40.root",
"fullDetectorDACDAC_col39row60.root",
"fullDetectorDACDAC_col39row20.root",
"fullDetectorDACDAC_col13row60.root"
};
//PLT specific
const int rocIdIni = 0;
const int rocIdFin = 3; 
const int maxValCalDel = 151;
const int maxValVcThr = 161;
const int numTriggers = 5;
/////
//   Declare functions 
/////
TFile* Call_TFile(string fileName);
void setTDRStyle();
/////
//   Main function
/////
void goodCalDelVcThr(){
 setTDRStyle();
 vector<string> fileNames(filesName, filesName + sizeof(filesName)/sizeof(filesName[0]));
 int numFiles = fileNames.size();
 string nameCan = "goodCalDelVcThr_roc";
 //TCanvas* can = new TCanvas(nameCan.c_str(),nameCan.c_str(),100,100,1300,4900);
 //can->Divide(rocIdFin-rocIdIni,1);
 //TPaveLabel *title = new TPaveLabel(0.1,0.98,0.9,0.99,("Percentage of good fit in "+to_string(numFiles)+" scans").c_str(),"brndc");
 //title->Draw(); 
 int canPos = 0;
 //Loop over different rocs
 for(int r=rocIdIni; r<rocIdFin; r++){
  //can->cd(r);
  string histName = histName_+to_string(r);
  TCanvas* can = new TCanvas((nameCan+to_string(r)).c_str(),(nameCan+to_string(r)).c_str(),100,100,2100,1600);
  TH2F* htot = new TH2F(histName.c_str(),histName.c_str(),maxValCalDel,0,maxValCalDel,maxValVcThr,0,maxValVcThr); 
  htot->GetYaxis()->SetTitle("VcThr");
  htot->GetXaxis()->SetTitle("CalDel");
  //Loop over different files
  for(int fn=0; fn<numFiles; fn++){
   TFile* file = Call_TFile(fileNames[fn]);
   if(file->GetListOfKeys()->Contains(histName.c_str())){
    TH2I* h = (TH2I*)file->Get(histName.c_str()); 
    htot->Add(h,100/(numTriggers*float(numFiles)));
   }
   delete file;
  }//Loop over different files
  htot->Draw("COLZ");
  //htot->Draw("textsame");
  htot->SetMinimum(0);
  htot->SetMaximum(100);
  //print info
  cout<<"Good CalDel VcThr values for "<<histName_+to_string(r)<<endl;
  for(int col=0; col<maxValCalDel; col++) for(int row=0; row<maxValVcThr; row++) if(htot->GetBinContent(col,row)==100) cout<<col<<" "<<row<<endl;
  cout<<endl;
  can->SaveAs(((nameCan+to_string(r))+".png").c_str());
 }//Loop over different rocs
 //can->SaveAs((nameCan+".png").c_str());
}
/////
//   Declare functions 
/////
/////
//   Call TFile to be read
/////
TFile* Call_TFile(string fileName) {
 string fN = pathFile+fileName;
 TFile* file = new TFile(fN.c_str());
 return file;
}
/////
//   Set setTDRStyle_modified (from link https://twiki.cern.ch/twiki/pub/CMS/TRK10001/setTDRStyle_modified.C)
/////
void setTDRStyle(){
  TStyle *tdrStyle = new TStyle("tdrStyle","Style for P-TDR");

  // For the canvas:
  tdrStyle->SetCanvasBorderMode(0);
  tdrStyle->SetCanvasColor(kWhite);
  tdrStyle->SetCanvasDefH(600); //Height of canvas
  tdrStyle->SetCanvasDefW(600); //Width of canvas
  tdrStyle->SetCanvasDefX(0);   //POsition on screen
  tdrStyle->SetCanvasDefY(0);
  // For the Pad:
  tdrStyle->SetPadBorderMode(0);
  // tdrStyle->SetPadBorderSize(Width_t size = 1);
  tdrStyle->SetPadColor(kWhite);
  tdrStyle->SetPadGridX(false);
  tdrStyle->SetPadGridY(false);
  tdrStyle->SetGridColor(0);
  tdrStyle->SetGridStyle(3);
  tdrStyle->SetGridWidth(1);
  // For the frame:
  tdrStyle->SetFrameBorderMode(0);
  tdrStyle->SetFrameBorderSize(1);
  tdrStyle->SetFrameFillColor(0);
  tdrStyle->SetFrameFillStyle(0);
  tdrStyle->SetFrameLineColor(1);
  tdrStyle->SetFrameLineStyle(1);
  tdrStyle->SetFrameLineWidth(1);

  // For the histo:
  tdrStyle->SetHistFillColor(0);
  // tdrStyle->SetHistFillStyle(0);
  tdrStyle->SetHistLineColor(1);
  tdrStyle->SetHistLineStyle(0);
  tdrStyle->SetHistLineWidth(1);
  // tdrStyle->SetLegoInnerR(Float_t rad = 0.5);
  // tdrStyle->SetNumberContours(Int_t number = 20);
//  tdrStyle->SetEndErrorSize(0);
  tdrStyle->SetErrorX(0.);
//  tdrStyle->SetErrorMarker(20);

  tdrStyle->SetMarkerStyle(20);

  //For the fit/function:
  tdrStyle->SetOptFit(1);
  tdrStyle->SetFitFormat("5.4g");
  //tdrStyle->SetFuncColor(1);
  tdrStyle->SetFuncStyle(1);
  tdrStyle->SetFuncWidth(1);

  //For the date:
  tdrStyle->SetOptDate(0);
  // tdrStyle->SetDateX(Float_t x = 0.01);
  // tdrStyle->SetDateY(Float_t y = 0.01);

  // For the statistics box:
  tdrStyle->SetOptFile(0);
  tdrStyle->SetOptStat(""); // To display the mean and RMS:   SetOptStat("mr");
  //tdrStyle->SetStatColor(kWhite);
  tdrStyle->SetStatColor(kGray);
  tdrStyle->SetStatFont(42);

  tdrStyle->SetStatTextColor(1);
  tdrStyle->SetStatFormat("6.4g");
  tdrStyle->SetStatBorderSize(0);
  tdrStyle->SetStatX(0.9); //Starting position on X axis
  tdrStyle->SetStatY(1.); //Starting position on Y axis
  tdrStyle->SetStatFontSize(0.025); //Vertical Size
  tdrStyle->SetStatW(0.15); //Horizontal size 
  // tdrStyle->SetStatStyle(Style_t style = 1001);

  // Margins:
  tdrStyle->SetPadTopMargin(0.1);
  tdrStyle->SetPadBottomMargin(0.125);
  tdrStyle->SetPadLeftMargin(0.115);
  tdrStyle->SetPadRightMargin(0.135);

  // For the Global title:

  //  tdrStyle->SetOptTitle(0);
  tdrStyle->SetTitleFont(50);
  tdrStyle->SetTitleColor(1);
  tdrStyle->SetTitleTextColor(1);
  tdrStyle->SetTitleFillColor(10);
  tdrStyle->SetTitleFontSize(0.06);
  //center position
  tdrStyle->SetTitleX(0.1f);
  tdrStyle->SetTitleW(0.8f);
  // tdrStyle->SetTitleH(0); // Set the height of the title box
  // tdrStyle->SetTitleW(0); // Set the width of the title box
  // tdrStyle->SetTitleX(0); // Set the position of the title box
  // tdrStyle->SetTitleY(0.985); // Set the position of the title box
  // tdrStyle->SetTitleStyle(Style_t style = 1001);
  tdrStyle->SetTitleBorderSize(0);

  // For the axis titles:
  tdrStyle->SetTitleColor(1, "XYZ");
  tdrStyle->SetTitleFont(42, "XYZ");
  tdrStyle->SetTitleSize(0.06, "XYZ");
  // tdrStyle->SetTitleXSize(Float_t size = 0.02); // Another way to set the size?
  // tdrStyle->SetTitleYSize(Float_t size = 0.02);
  tdrStyle->SetTitleXOffset(0.9);
  tdrStyle->SetTitleYOffset(0.9);
  // tdrStyle->SetTitleOffset(1.1, "Y"); // Another way to set the Offset

  // For the axis labels:

  tdrStyle->SetLabelColor(1, "XYZ");
  tdrStyle->SetLabelFont(42, "XYZ");
  tdrStyle->SetLabelOffset(0.007, "XYZ");
  tdrStyle->SetLabelSize(0.05, "XYZ");

  // For the axis:

  tdrStyle->SetAxisColor(1, "XYZ");
  tdrStyle->SetStripDecimals(kTRUE);
  tdrStyle->SetTickLength(0.03, "XYZ");
  tdrStyle->SetNdivisions(510, "XYZ");
  tdrStyle->SetPadTickX(1);  // To get tick marks on the opposite side of the frame
  tdrStyle->SetPadTickY(1);

  // Change for log plots:
  tdrStyle->SetOptLogx(0);
  tdrStyle->SetOptLogy(0);
  tdrStyle->SetOptLogz(0);
  // Postscript options:
  // tdrStyle->SetPaperSize(15.,15.);
  // tdrStyle->SetLineScalePS(Float_t scale = 3);
  // tdrStyle->SetLineStyleString(Int_t i, const char* text);
  // tdrStyle->SetHeaderPS(const char* header);
  // tdrStyle->SetTitlePS(const char* pstitle);
  // tdrStyle->SetBarOffset(Float_t baroff = 0.5);
  // tdrStyle->SetBarWidth(Float_t barwidth = 0.5);
  //tdrStyle->SetPaintTextFormat(const char* format = "g");
  //tdrStyle->SetPalette(0,0);
  // tdrStyle->SetTimeOffset(Double_t toffset);
  // tdrStyle->SetHistMinimumZero(kTRUE);
  tdrStyle->SetPaintTextFormat("3.1f");

  tdrStyle->cd();
}