#ifndef detector_offset_h
#define detector_offset_h


#include <vector>

#include "TROOT.h"
#include "TNtuple.h"
#include "TString.h"

#include "TFile.h"
#include "TCanvas.h"
#include "TH1F.h"
#include "TGraph.h"
#include "TGraphErrors.h"

#include "TChain.h"
#include "TLine.h"

#include "TH2F.h"
#include "TCut.h"
#include "TStyle.h"

/*--------------------------------------------------*/
/// List of functions


void Init();

void Add_Root_File(Int_t* grp);

void Cer_plot(TChain*chain, TString plot_title, TCut data_cut, Float_t threshhold, Float_t bin_min, Float_t bin_max, Int_t bin_num);

void Beta_plot(TChain*chain, TString plot_title, TCut data_cut, Float_t threshhold, Float_t bin_min, Float_t bin_max, Int_t bin_num);

void Aero_plot(TChain*chain, TString plot_title, TCut data_cut, Float_t threshhold, Float_t bin_min, Float_t bin_max, Int_t bin_num); 

void check_plots(TChain*chain, TString parameter, TString log);

void plot_cointime_finebin(TChain*chain, TString plot_title, Float_t bin_min, Float_t bin_max, Int_t bin_num);

void plot_beta_finebin(TChain*chain, TString plot_title, Float_t bin_min, Float_t bin_max, Int_t bin_num);

void plot_beta_vs_cointime_finebin(TChain*chain, TString plot_title, Float_t x_bin_min, Float_t x_bin_max, Int_t x_bin_num, Float_t y_bin_min, Float_t y_bin_max, Int_t y_bin_num);

Int_t plot_beta_vs_cointime_finebin_cut(TChain*chain, TCut cut, TString plot_title, Float_t x_bin_min, Float_t x_bin_max, Int_t x_bin_num, Float_t y_bin_min, Float_t y_bin_max, Int_t y_bin_num);

void plot_beta_vs_cointime_finebin_pretty(TChain*chain, TString plot_title, Float_t x_bin_min, Float_t x_bin_max, Int_t x_bin_num, Float_t y_bin_min, Float_t y_bin_max, Int_t y_bin_num);


//void Coin_time_loop(Int_t* grp);
void Coin_time_loop(std::vector<int> grp);



void Run_by_Run(Int_t run_num, Float_t x_bin_min, Float_t x_bin_max, Int_t x_bin_num, Float_t y_bin_min, Float_t y_bin_max, Int_t y_bin_num); 



/// hsbeta float
Int_t Calculate_efficiency(TH2F* myhist_2D);

/// hsbeta with fixed value
Int_t Calculate_efficiency(TH2F* myhist_2D, Float_t hsbeta_limit);

void DrawGoodBox(float x_center, float y_center);
void DrawRandomBox(float x_center, float y_center);

void Histogram_output(TH2F *hist);

bool Check_exists(std::vector<int> vectora, int number);

void Write_overall();

Float_t Calculate_tail_x (Float_t xx, Float_t yy);

Float_t Calculate_miss_x (Float_t x_boundary, Float_t y_boudary, Float_t y_pos);

TChain* globchain1; 		
TChain* globchain2;		
TChain* globchain3;	
TChain* globchain4;	
TChain* globchain_total; 



TFile* file_out;

TDirectory* test_dir;              
TDirectory* plot_dir_Q2_2p4_EB_4p2;
TDirectory* plot_dir_Q2_6p5_EB_5p2;
TDirectory* plot_dir_Q2_4p4_EB_3p7;
TDirectory* plot_dir_Q2_5p4_EB_4p7;



TH2F* Q2_2p4_EB_4p2_combine; 
TH2F* Q2_6p5_EB_5p2_combine;
TH2F* Q2_4p4_EB_3p7_combine;
TH2F* Q2_5p4_EB_4p7_combine;
TH2F* total_hbeta_cointime;

float beta_cut;             
float global_centered_hsbeta_cut;
float gradi;                 

#endif