SANCinterface.cxx

#include <iostream>
#include "SANCtable.h"
using std::cout;

int main()
{
        SANCtable f1("table1-1.txt"),f2("table2-2.txt");
        f1.setFlavor(1);
        f2.setFlavor(2);

        //Uncomment for born-level computation
        //f1.setBornLevel(true);
        //f2.setBornLevel(true);

        //Write with fixed precision
        f1.setFixedLength(8);
        f2.setFixedLength(8);

        //Set the basic variables. Only needed once before all computation
        SANCtable::setFlags();
        //The dimensions must match the tauola interface
        SANCtable::setDimensions(100,100,100,21);
        //Define the sqrt mass range for all three data sets
        SANCtable::setRanges(6,17000,85,110,160,220);

        //Add header and additional information to files
        if(!f1.addHeader())
        {
                cout<<"Ranges or dimensions not set.\n";
                return -1;
        }
        if(!f1.addFile("SancLib_v1_02/lib.txt"))
        {
                cout<<"No info file on electrowek library SANC (file SancLib_v1_02/lib.txt missing).\n";
                return -1;
        }

        if(!f1.addFile("var.dump"))
        {
                cout<<"No initialization variables of SANC (file var.dump missing).\n";
                return -1;
        }
        //Same for the second file
        if(!f2.addHeader() || !f2.addFile("SancLib_v1_02/lib.txt") || !f2.addFile("var.dump"))
        {
                cout<<"Second file I/O operations failed.\n";
                return -1;
        }

        //Start computation
        f1.addRange(1,true); //The first range is a logarithmic range
        f1.addRange(2);
        f1.addRange(3);
        f1.close();

        f2.addRange(1,true); //The first range is a logarithmic range
        f2.addRange(2);
        f2.addRange(3);
        f2.close();
        return 0;
}