#ifndef functions_hadcalib_h
#define functions_hadcalib_h

#include "TProfile2D.h"
#include "TString.h"
#include <TH2.h>
#include <TStyle.h>
#include <TCanvas.h>
#include <iostream>
#include <cmath>
#include <iomanip>
#include <string>
#include "functions_hadcalib.h"
#include <TTree.h>
#include <fstream>
#include <sstream>
#include <stdlib.h>
#include <math.h>
#include <TH3F.h>

void printvector(double N_vector[], TString ttext, const int N); 

void printvector(int N_vector[], TString ttext, const int N);

void logaxisbins(const Int_t nbins, Double_t lowbin, Double_t upbin, Double_t* logbins);

Int_t find_seedcluster(Int_t Nclusters, Float_t* cluster_energy);

Double_t make_cone(Int_t Nclusters, Float_t* cluster_energy, Float_t* eta_cl, Float_t* phi_cl, Double_t coneR, Int_t* coneclusters, Int_t* cellcluster_topo, Int_t* clusters_in_cone);

Double_t make_cones(Int_t Nclusters, Float_t* cluster_energy, Float_t* eta_cl, Float_t* phi_cl, Double_t coneR, Int_t* coneclusters, Int_t* cellcluster_topo, Int_t* Nclusters_in_cone, Int_t* excluded_clusters);

Double_t get_weight(TProfile2D *weights_table, Double_t Econe, Double_t E_dens);

Double_t get_weight_old(TH2F *weights_table, Double_t Econe, Double_t E_dens);

void estimate_EMfraction(Float_t* cl_EM_frac, Float_t* cl_e_topo, Float_t* cl_et_topo, Float_t* cl_m2_lambda_topo, Float_t* cl_m2_r_topo, Float_t* cl_center_lambda_topo, Int_t Nclusters, Double_t Ecl_cut, Double_t* c);

void CalibHit_cells_to_cluster(Float_t* cell_var, Float_t* cellweight_topo, Double_t* cluster_var, Int_t* cellcluster_topo,Int_t Ncells, Int_t Nclusters);

Double_t EMfrac_from_histo(TProfile2D *EM_histo, Double_t log_densmom, Double_t log_lambdac, char opt, Double_t EMnoise);

void cellcluster_array(Int_t* cellcluster, Int_t Nclusters, Int_t* ncells_topoC);

#endif // #ifdef hadcalib_mc_v1204_SinglePions_all_cxx