Skip to content
Snippets Groups Projects
Commit e76662df authored by Marc Pelizaeus's avatar Marc Pelizaeus
Browse files

Add AICc to App

parent eb23491e
No related branches found
No related tags found
No related merge requests found
......@@ -208,16 +208,19 @@ int main(int __argc,char *__argv[]){
double evtWeightSumData = theJpsiGamXEventListPtr->NoOfWeightedDataEvts();
double BICcriterion=2.*theLh+noOfFreeFitParams*log(evtWeightSumData);
double AICcriterion=2.*theLh+2.*noOfFreeFitParams;
double AICccriterion=AICcriterion+2.*noOfFreeFitParams*(noOfFreeFitParams+1)/(evtWeightSumData-noOfFreeFitParams-1);
Info << "noOfFreeFitParams:\t" <<noOfFreeFitParams;
Info << "evtWeightSumData:\t" <<evtWeightSumData;
Info << "BIC:\t" << BICcriterion << endmsg;
Info << "AIC:\t" << AICcriterion << endmsg;
Info << "AICc:\t" << AICccriterion << endmsg;
std::string qaSummaryFileName = "qaSummary" + jobOption + ".dat";
std::ofstream theQaStream ( qaSummaryFileName.c_str() );
theQaStream << "BIC\t" << BICcriterion << "\n";
theQaStream << "AIC\t" << AICcriterion << "\n";
theQaStream << "AICa\t" << AICcriterion << "\n";
theQaStream << "AICc\t" << AICccriterion << "\n";
theQaStream << "logLh\t" << theLh << "\n";
theQaStream << "free parameter\t" << noOfFreeFitParams << "\n";
theQaStream.close();
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment