EstimateCredibility.cpp File Reference

Functions
int	main (int argc, char **argv)

Variables
namespace	filesystem = std::experimental::filesystem

Function Documentation

main()

int main	(	int	argc,
		char **	argv
	)

                               {
  if(argc != 3) {
    std::cerr << "Usage: 1.This 2.Directory 3.Nsample\n";
    std::exit(EX_USAGE);
  }
  constexpr int NMAX = 20;
  Integer_t const Nsample = std::atoll(argv[2]);
  std::string const baseDir(argv[1]);
  std::string const OutDIR(baseDir+"/Output");
  filesystem::create_directory(OutDIR);
  std::cout << "Directory: " << baseDir << "\n";
 
  std::string   InFname(OutDIR+"/Stat_MaxFluc.txt");
  std::ifstream InFs(InFname); checkIsFileOpen(InFs, InFname);
  std::string line;
  skip_header(InFs);
  debug_print("# Input: " << InFname);
 
  int n, Nmax=0, Nmin=NMAX;
  double dtemp;
  std::vector<double> Mean(NMAX+1, NAN), Stddev(NMAX+1, NAN);
  std::vector<Integer_t> Ndata(NMAX+1, NAN);
  try {
    for(getline(InFs, line); !InFs.eof(); getline(InFs, line)) {
      std::istringstream inStream(line);
      inStream >> n;
      inStream >> Mean[n] >> Stddev[n] >> dtemp >> Ndata[n];
      Stddev[n] /= sqrt((double)Ndata[n]);
      Nmax = (n > Nmax? n: Nmax);
      Nmin = (n < Nmin? n: Nmin);
    }
  } catch(std::exception& e) {
    if( !InFs.eof() )
    std::cerr << "Error(" << InFs.eof() << InFs.fail() << InFs.bad() << "):" << e.what() << std::endl;
  }
  InFs.close();
  for(n = Nmax;n >= Nmin; --n) {
    std::cout << n << "\t" << Mean[n] << "\t" << Stddev[n] << std::endl;
  }
 
  bool flag;
  std::vector<double> Rand(NMAX+1);
  std::vector<double> ProbDec(NMAX,0.0), ProbInc(NMAX,0.0);
  std::vector<double> ProbDecMax(NMAX,0.0), ProbIncMax(NMAX,0.0);
  for(size_t sample = 0;sample < Nsample; ++sample) {
    for(n = Nmin;n <= Nmax; ++n) Rand[n] = Mean[n] + Stddev[n]*Gaussian();
    for(int mod=0; mod <= 1; ++mod) {
      flag = true;
      for(n = Nmax-2;n >= Nmin; --n) {
        if(n%2 != mod) continue;
        if(Rand[n]<Rand[n+2] || flag==false) flag = false;
        else ProbDec[n] += 1.0/(double)Nsample;
        if(Rand[n] > Rand[(Nmax/2)*2-mod]) ProbDecMax[n] += 1.0/(double)Nsample;
        debug_print("# sample(" << sample << ") Rand[" << n+2 << "] = " << Rand[n+2]);
      }
      debug_print("# sample(" << sample << ") Rand[" << n << "] = " << Rand[n]);
      debug_print("");
 
      flag = true;
      for(n = Nmax-2;n >= Nmin; --n) {
        if(n%2 != mod) continue;
        if(Rand[n]>Rand[n+2] || flag==false) flag = false;
        else ProbInc[n] += 1.0/(double)Nsample;
        if(Rand[n] < Rand[(Nmax/2)*2-mod]) ProbIncMax[n] += 1.0/(double)Nsample;
      }
    }
    debug_print("");
  }
 
  std::string   OutFname(OutDIR+"/Stat_MaxFluc_Credibitity.txt");
  std::ofstream OutFs(OutFname); checkIsFileOpen(OutFs, OutFname);
  OutFs << "# Nsample = " << Nsample << "\n";
  OutFs << "# 1.Start 2.End 3.ProbDec 4.ProbInc 5.ProbDecMax 6.ProbIncMax\n";
  for(int mod=0; mod <= 1; ++mod) {
    for(n = Nmax-2;n >= Nmin; --n) {
      if(n%2 != mod) continue;
      OutFs << n       << " "
            << (Nmax/2)*2-mod << " "
            << ProbDec[n] << " "
            << ProbInc[n] << " "
            << ProbDecMax[n] << " "
            << ProbIncMax[n] << "\n";
    }
    OutFs << "\n";
  }
  OutFs.close();
 
 
  return 0;
}

Variable Documentation

filesystem

namespace filesystem = std::experimental::filesystem