#ifndef source_h
#define source_h


//C++ PRE DEFINED LIBRARY
#include <iostream>
#include <vector>
#include <fstream>
#include <string>
#include <iomanip>
#include <math.h>
//ROOT libraries
//#include "TApplication.h"
//#include "TLine.h"
//root-config --cflags --glibs
#include "TAttMarker.h"
#include "TPolyMarker.h"
#include "TLegend.h"
#include "TLegendEntry.h"
#include "TGraph.h"
#include "TMultiGraph.h"
#include "TCanvas.h"
#include "TLatex.h"
#include "TDirectory.h"
#include "TAxis.h"
#include "TGaxis.h"
#include "TFile.h"
#include "TPaveText.h"
#include "TObject.h"
#include "TStyle.h"
#include "TAttLine.h"
#include "TAttFill.h"
#include "TPave.h"
#include "TPad.h"
#include "TMath.h"
#include "TString.h"
#include "TCollection.h"
#include "RConfigure.h"
#include "TSystem.h"
#include "Math/Minimizer.h"
#include "Minuit2/MnUserParameterState.h"

// GLOBAL VARIABLES
const double PI = 3.141592653589793;
const double Avogadro   =  6.022E23; //atoms/mole
const double c_light    =  2.99793E8; //m/s
const double c_square   =  931.5 ;  //MeV/u
const double FourPiRsquare= 4.0* PI * (2.818*1.E-15) *  (2.818*1.E-15); //m^2

//================================
//charge of the particles
const double z_electron = 1.0 ; // = z in the formula
//////////////////////////////////////////////////////////
//mass of the particles
const double m_electron =  0.00054858 ; //in akb, u;  9.11E-31;  //kg
//===================Medium Material=============
 double A_Si       = 28.0  ; 
 double rho_Si     = 2.33*1.E3 ;
 double N_Si       = (rho_Si*Avogadro)/ (A_Si*1.0E-3);
 double Z_Si       = 14.0 ;
 double MEP_Si     = 172.*1.0E-6 ;

#endif /* source_h */