diff --git a/glauber/Fitter.cpp b/glauber/Fitter.cpp
index 99fed417f9258f4fbd42f40e748fcc0790687ff7..b0c930652446a5086dad7e9d5c93909cb2992af9 100644
--- a/glauber/Fitter.cpp
+++ b/glauber/Fitter.cpp
@@ -28,7 +28,8 @@ Glauber::Fitter::Fitter(std::unique_ptr<TTree> tree)
 void Glauber::Fitter::Init(int nEntries)
 {
     if ( nEntries < 0 || nEntries > fSimTree->GetEntries() ){
-        std::cout << "Init: *** ERROR - Number of Entries < 0 or less that number of entries in input tree" << std::endl;
+        std::cout << "Init: *** ERROR - number of entries < 0 or less that number of entries in input tree" << std::endl;
+        std::cout << "Init: *** number of entries in input tree = " << fSimTree->GetEntries() << std::endl;        
         exit(EXIT_FAILURE);
     }
     
diff --git a/glauber/main.cpp b/glauber/main.cpp
index e48a1d2f61de81f86a3103b89a310535ae4a1b53..5b3dbba9ba216ee8356f8e75f041f2f4478fab1c 100644
--- a/glauber/main.cpp
+++ b/glauber/main.cpp
@@ -27,23 +27,28 @@ int main(int argc, char *argv[])
     // Modify this part according to your needs
     // *****************************************
 
-    const Int_t min_bin = 50;
+    ///  |   mode    |   function for Na     |
+    ///  |   Default |  Npart + (1-f)*Ncoll  |
+    ///  |   PSD     |     f - Npart         |
+    ///  |   Npart   |     Npart^f           |
+    ///  |   Ncoll   |     Ncoll^f           |
     const TString mode = "Default";
     
-    const TString glauber_filename = "/home/vklochkov/Data/glauber/pbpb/30agev/glau_pbpb_ntuple_signn_31.0_7.6AGeV_CM_30AGeV_LC.root";
+    const TString glauber_filename = "../input/glauber_auau_sigma_30_100k.root";   // input files
+    const TString glauber_treename = "nt_Au_Au";
     const TString in_filename = "../input/test_input.root";
     const TString histoname = "hMreco";
 
-    const Int_t max_bin = 400;
-    const Int_t nevents = 199000;
+    const Int_t min_bin = 50;      // not fitting low multiplicity region due to trigger bias, etc
+    const Int_t max_bin = 10000;   // very large number to fit the whole histo
+    const Int_t nevents = 100000;
 
-    const TString treename = "nt_Pb_Pb";
     const TString outdir = ".";
     // *****************************************
     // *****************************************
 
     std::unique_ptr<TFile> glauber_file{ TFile::Open(glauber_filename, "read") };
-    std::unique_ptr<TTree> glauber_tree{ (TTree*) glauber_file->Get(treename) };
+    std::unique_ptr<TTree> glauber_tree{ (TTree*) glauber_file->Get(glauber_treename) };
 
     std::unique_ptr<TFile> f{TFile::Open(in_filename)};    
     TH1F *hdata = (TH1F*)f->Get(histoname);
diff --git a/input/glauber_auau_sigma_30_100k.root b/input/glauber_auau_sigma_30_100k.root
new file mode 100644
index 0000000000000000000000000000000000000000..7260febba6d50456daebd78ae93b227ccc49ca2b
Binary files /dev/null and b/input/glauber_auau_sigma_30_100k.root differ