void RunGlauberEstimator(TString flag="")
{
    gROOT->LoadMacro("CentralityEstimator.cxx+g");
    
    CentralityEstimator *centEst = new CentralityEstimator();
    
    TString infilename ;
    TString treename ;
    TString outfilename ;    


    if(flag == "AuAu"){
        infilename = "./root/GlauberMC_Au3Au3_ntuple_23.6mb_r6.5541_a0.523_m0.9.root";
        treename = "nt_Au3_Au3";
        outfilename = "./root/GlauberMC_Au3Au3_ntuple_23.6mb_r6.5541_a0.523_m0.9_withEstimator.root";    
        
        centEst->AddEstimator("TOF"                  , -1.64069e-06 , 0.202237 , 6.36083  );
        //f = -1.64069e-06    mu = 0.202237    k = 6.36083    Chi2Min = 50.1303    
        centEst->AddEstimator("RPC"               ,   -1.64069e-06  , 0.503523  , 29.063 ); 
        // f = -1.64069e-06    mu = 0.503523    k = 29.063    Chi2Min = 39.0985
        centEst->AddEstimator("PrimaryParticleCand"  , -1.3825e-07 , 0.185581  , 7.60447 );
        //f = -1.3825e-07    mu = 0.185581    k = 7.60447    Chi2Min = 89.24
        centEst->AddEstimator("SelectedParticleCand" ,  -1.3825e-07  ,0.215636 , 5.18754  );
        //f = -1.3825e-07    mu = 0.215636    k = 5.18754    Chi2Min = 45.4244
        
    }
    else if(flag == "AgAg"){
        infilename = "./root/GlauberMC_AgAg_ntuple_25.0mb_m0.9AgAg-snn25.0-md0.9-nd0.0-rc1-smax99.0.root";
        treename = "nt_Ag_Ag";
        outfilename = "./root/GlauberMC_AgAg_ntuple_25.0mb_m0.9AgAg-snn25.0-md0.9-nd0.0-rc1-smax99.0_withEstimator.root"; 
        
        centEst->AddEstimator("TOFRPC"               ,  1           ,0.798707 , 120.746  );
        //f = 1    mu = 0.798707    k = 120.746    Chi2Min = 76.0805
        centEst->AddEstimator("TOF"                  , -7.07424e-06 , 0.274859 , 9.55442  );
        //f = -7.07424e-06    mu = 0.274859    k = 9.55442    Chi2Min = 23.0534
        centEst->AddEstimator("RPC"               ,   -1.96366e-06  , 0.604601  , 55.284 ); 
        //f = -1.96366e-06    mu = 0.604601    k = 55.284    Chi2Min = 33.0516
        centEst->AddEstimator("PrimaryParticleCand"  , -1.27035e-06 , 0.168996 , 6.23886 );
        //f = -1.27035e-06    mu = 0.168996    k = 6.23886    Chi2Min = 49.7326
        centEst->AddEstimator("SelectedParticleCand" , -1.78391e-06 , 0.221611 , 8.55049  );
        //f = -1.78391e-06    mu = 0.221611    k = 8.55049    Chi2Min = 108.999
    }
    else  if(flag == "AuC"){
        infilename = "./root/GlauberMC_AgC_ntuple_25mb_ntot.root";
        treename = "nt_Ag_C";
        outfilename = "./root/GlauberMC_AgC_ntuple_25mb_withEstimator.root"; 
    }
    else  if(flag == "AgC"){
        infilename = "./root/GlauberMC_Au3C_ntuple_23.6mb_ntot.root";
        treename = "nt_Au3_C";
        outfilename = "./root/GlauberMC_Au3C_ntuple_25mb_withEstimator.root"; 
    }
    
    //centEst->SetMode("Default");
    centEst->SetInputTree ( infilename, treename);
    centEst->SetOutputTree ( outfilename );
    centEst->LoadInputTree();
    //centEst->SetNpartMax (400);   //not used in Default-Setting
    //centEst->SetNcollMax (1100);
    centEst->Run();

    cout << "Estimator calulaterd " << "\r" << endl;
    return;
}