portelli/LatAnalyze
1
0
Fork 0
mirror of https://github.com/aportelli/LatAnalyze.git synced 2025-04-11 03:20:46 +01:00
Code Releases Activity

Last commit before sync with original repo

Browse Source
This commit is contained in:
Andrew Zhen Ning Yong 2019-07-29 19:39:14 +08:00
parent 4cab49cd92
commit 2096238ff6
3 changed files with 42 additions and 26 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
physics/2pt-fit.cpp
View File

@ -30,7 +30,7 @@ int main(int argc, char *argv[])
opt.addOption("s", "shift" , OptParser::OptType::value , true, opt.addOption("s", "shift" , OptParser::OptType::value , true,
"time variable shift", "0"); "time variable shift", "0");
opt.addOption("m", "model" , OptParser::OptType::value , true, opt.addOption("m", "model" , OptParser::OptType::value , true,
"fit model (exp|exp2|exp3|sinh|cosh|cosh2|cosh3|cotanh|explin|const|<interpreter code>)", "cosh"); "fit model (exp|exp2|exp3|sinh|cosh|cosh2|cosh3|coth|explin|const|<interpreter code>)", "cosh");
opt.addOption("" , "nPar" , OptParser::OptType::value , true, opt.addOption("" , "nPar" , OptParser::OptType::value , true,
"number of model parameters for custom models " "number of model parameters for custom models "
"(-1 if irrelevant)", "-1"); "(-1 if irrelevant)", "-1");
@ -186,13 +186,16 @@ int main(int argc, char *argv[])
+ p[5]*(exp(-p[2]*x[0])+exp(-p[4]*(nt-x[0]))); + p[5]*(exp(-p[2]*x[0])+exp(-p[4]*(nt-x[0])));
}, 1, nPar); }, 1, nPar);
} }
else if (model == "cotanh") else if (model == "coth")
{ {
cotanhModel = true; cotanhModel = true;
nPar = 2; nPar = 2;
mod.setFunction([nt](const double *x, const double *p) mod.setFunction([nt](const double *x, const double *p)
{ {
return p[0] + p[1]/tanh(0.1892*(nt/2-x[0])); // 0.1892 is the meson eff mass in lattice unit from 2pt fit // cout << "Using effective pion mass, am = 0.1892" << endl;
// kaon physical mass:0.28856
// pion physical mass: 0.08046
return p[0] + p[1]/tanh(0.28856*(nt/2-x[0])); // 0.1892 is the meson eff mass in lattice unit from 2pt fit; 0.08046 is the meson mass at physical point
}, 1, nPar); }, 1, nPar);
} }
else if (model == "explin") else if (model == "explin")
@ -270,8 +273,13 @@ int main(int argc, char *argv[])
} }
else if(cotanhModel) else if(cotanhModel)
{ {
init(0) = 0.5*(data.y(0, 0)[central]+ data.y(nt-1, 0)[central]); // init(0) = 0.5*(data.y(0, 0)[central]+ data.y(nt-1, 0)[central]);
init(1) = 0.5*(data.y(0, 0)[central]- data.y(nt-1, 0)[central]); // init(1) = 0.5*(data.y(0, 0)[central]- data.y(nt-1, 0)[central]);
init(0) = -0.1677;
init(1) = 0.02526;
cout << init(0) << endl;
} }
else else
{ {
@ -282,6 +290,7 @@ int main(int argc, char *argv[])
{ {
init(p) = 2*init(p - 2); init(p) = 2*init(p - 2);
init(p + 1) = init(p - 1)/2.; init(p + 1) = init(p - 1)/2.;
} }
// set limits for minimiser ////////////// // set limits for minimiser //////////////
for (Index p = 0; p < nPar; p += 2) for (Index p = 0; p < nPar; p += 2)
@ -299,9 +308,16 @@ int main(int argc, char *argv[])
} }
else if(cotanhModel) else if(cotanhModel)
{ {
globMin.setLowLimit(p,-0.5); // globMin.setLowLimit(p,-0.5);
locMin.setLowLimit(p, -0.3); // locMin.setLowLimit(p, -0.3);
globMin.setHighLimit(p,1); // globMin.setHighLimit(p,1);
globMin.setLowLimit(p,-1);
locMin.setLowLimit(p, -0.3);
globMin.setHighLimit(p,0);
globMin.setLowLimit(p+1,-0.12);
locMin.setLowLimit(p+1, -0.1);
globMin.setHighLimit(p+1,0.9);
} }
else else
{ {

18
physics/eff-mass.cpp
View File

@ -1,12 +1,12 @@
#include <LatCore/OptParser.hpp> #include <LatAnalyze/Core/OptParser.hpp>
#include <LatAnalyze/CompiledModel.hpp> #include <LatAnalyze/Functional/CompiledModel.hpp>
#include <LatAnalyze/Io.hpp> #include <LatAnalyze/Io/Io.hpp>
#include <LatAnalyze/MatSample.hpp> #include <LatAnalyze/Statistics/MatSample.hpp>
#include <LatAnalyze/Math.hpp> #include <LatAnalyze/Core/Math.hpp>
#include <LatAnalyze/MinuitMinimizer.hpp> #include <LatAnalyze/Numerical/MinuitMinimizer.hpp>
#include <LatAnalyze/NloptMinimizer.hpp> #include <LatAnalyze/Numerical/NloptMinimizer.hpp>
#include <LatAnalyze/Plot.hpp> #include <LatAnalyze/Core/Plot.hpp>
#include <LatAnalyze/XYSampleData.hpp> #include <LatAnalyze/Statistics/XYSampleData.hpp>
using namespace std; using namespace std;
using namespace Latan; using namespace Latan;

18
utils/plot.cpp
View File

@ -1,12 +1,12 @@
#include <LatCore/OptParser.hpp> #include <LatAnalyze/Core/OptParser.hpp>
#include <LatAnalyze/CompiledModel.hpp> #include <LatAnalyze/Functional/CompiledModel.hpp>
#include <LatAnalyze/Io.hpp> #include <LatAnalyze/Io/Io.hpp>
#include <LatAnalyze/MatSample.hpp> #include <LatAnalyze/Statistics/MatSample.hpp>
#include <LatAnalyze/Math.hpp> #include <LatAnalyze/Core/Math.hpp>
#include <LatAnalyze/MinuitMinimizer.hpp> #include <LatAnalyze/Numerical/MinuitMinimizer.hpp>
#include <LatAnalyze/NloptMinimizer.hpp> #include <LatAnalyze/Numerical/NloptMinimizer.hpp>
#include <LatAnalyze/Plot.hpp> #include <LatAnalyze/Core/Plot.hpp>
#include <LatAnalyze/XYSampleData.hpp> #include <LatAnalyze/Statistics/XYSampleData.hpp>
using namespace std; using namespace std;
using namespace Latan; using namespace Latan;