Documentation / Manuel utilisateur en C++ :

The constructor prototype used with an analytic function is:

// Create a TMorris object with an analytic function
TMorris(TDataServer *tds,const char *fcn,Int_t nreplica, Int_t level, Double_t delta=0, TString sout="")
TMorris(TDataServer *tds,void *fcn,TString sinp,TString sout,Int_t nreplica, Int_t level, Double_t delta=0)

This constructor takes six arguments:

Here is an example of how to use the constructor with an analytic function:

Int_t nreplica = 3;
Int_t nlevel = 10;
TMorris * scmo = new TMorris(tds, "flowrateModel", nreplica, nlevel); 

The constructor prototype used with a code is:

// Create a TMorris object with a code
TMorris(TDataServer *tds,TCode *fcode, Int_t nreplica, Int_t level,  Double_t delta=0)

This constructor takes five arguments:

Here is an example using this constructor for the flowrate case:

The constructor prototype used with a runner object is:

// Create a TMorris object with a runner
TMorris(TDataServer *tds,TRun *run, Int_t nreplica, Int_t level,  Double_t delta=0)

This constructor takes five arguments:

Here is an example using this constructor for the flowrate code in sequential mode:

// The input file
TKeyScript infile("flowrate_input_with_keys.in");
// provide the input and their key
infile.setInputs(8, tds->getAttribute("rw"), "Rw", tds->getAttribute("r"), "R",
tds->getAttribute("tu"), "Tu", tds->getAttribute("tl"), "Tl", tds->getAttribute("hu"), "Hu",
tds->getAttribute("hl"), "Hl", tds->getAttribute("l"), "L", tds->getAttribute("kw"), "Kw");

TAttribute yhat("yhat");
// The output file of the code
TKeyResult outfile("_output_flowrate_withKey_.dat");
// The attribute in the output file
outfile.addOutput(&yhat, "yhat");

// Instanciation de mon code
TCodeEval code("flowrate -s -k");
// Adding the intput/output file to the code
code.addInputFile(&infile);
code.addOutputFile(&outfile);

TSequentialRun run(&code);
run.startSlave();
if(run.onMaster())
{
    Int_t nreplicas = 3;
    Int_t nlevels = 10;
    TMorris * scmoR = new TMorris(tds, &run, nreplicas, nlevels);
    // ...
}

The two constructors before are used in the cases where the computation of the model outputs are run from Uranie. However it is possible to compute the outputs of the model outside of Uranie then load them in a TDataServer object (via a file) and use that TDataServer object to compute the finite differences indices.

The constructor prototype used with a TDataServer object already containing the simulations is:

//  Create a TMorris object with already filled TDS
TMorris(TDataServer *tds,  const char *inp,  const char *out,  Option_t *option="")

This constructor takes three arguments:

Below is an example of the constructor with a TDataServer object filled:

// Define the DataServer
TDataServer *tds = new TDataServer();
tds->fileDataRead("_morris_launching_.dat");

TMorris * scmo = new TMorris(tds, "rw:r:tu:tl:hu:hl:l:kw", "flowrateModel"); 

The first method use the method getValue of the TMorris object specifying the order of the extract value ("mustar", "mu" ou "sigma"), the related input and possibly more selected options.

double hl_mustar_index = scmo->getValue("mustar","hl");
double hl_mu_index = scmo->getValue("mu","hl");
double hl_sigma_index = scmo->getValue("sigma","hl");

The second method uses 3 steps to extract an index:

Below is an example of extraction of the index for hl in our flowrate case:

ntresu->Scan("*","Input==\"hl\"");
ntresu->GetEntry(5);
Double_t Hl_mu_Indexe =  ntresu->GetLeaf("mu")->GetValue(); 

The third method uses 2 steps to extract an index:

Below is another example of extraction of the index for "rw" in our flowrate case:

ntresu->Draw("mustar","Input==\"rw\"","goff");
Double_t Rw_mustar_Indexe = ntresu->GetV1()[0];