#ifndef drawAsymSysErrOme_h
#define drawAsymSysErrOme_h
const int nsysSc=9;  //The number of the source of systematic errors, such as accidentals, coherent beam energy, cut conditions
double Asym[5][9][3][2][3]; //i("standard", "accidentals", "coherent BE", "MMsq", "Delta_phi", "ME","MMXp", "IM_2gamma"), j("cuta", "cutb", "cutc"), k(pi0, eta), l(TOT, Ran, TSubR)
double Err[5][9][3][2][3];
double Asym_m2g_bin[5][9][3][6]; //i("standard", "accidentals", "coherent BE", "MMsq", "Delta_phi", "ME","MMXp", "IM_2gamma"), j("cuta", "cutb", "cutc"), l(TOT, Ran, TSubR)
double Err_m2g_bin[5][9][3][6];
double Errtmp;
TString hN[5][9][3][2][3]; //The different Name of ijkl
double NPERP_a[5][9][3][2][3];
double NPARA_a[5][9][3][2][3];
double Norm_a[5][9][3][2][3];
//double NPERP_btmp;
double NPARA_btmp;
double Asymtmp;

double frac[5][9][3][2][3];
double Asymsig[5][9][3][2][3]; 
double Asymbkg[5][9][3][2][3];
 
double Sigma_pi0[5][9][3];// m(rebin = 1,2,4,5,10), i("standard", "accidentals", "coherent BE", "MMsq", "Delta_phi", "ME","MMXp", "IM_2gamma"), j("cuta", "cutb", "cutc"),
double Sigma_eta[5][9][3];
double Err_Sigma_pi0[5][9][3];
double Err_Sigma_eta[5][9][3];

#endif