///////////////////////////////////////////////////////////////////
// ROOT macro to study the behaviour of signal detection related //
// to (astrophysical) transient phenomena.                       //
// This macro makes use of the generic NcAstrolab facility       //
// of the ROOT based NCFS-Pack framework.                        //
// See https://nick-ve.github.io/ncfs/docs for the documentation.//
//                                                               //
// Note that the corresponding ROOT and NCFS packages have to be //
// correctly installed before one can use this macro.            //
// For information about the installation of these packages,     //
// see https://nick-ve.github.io/ncfs                            //
//                                                               //
// To execute the macro, just invoke ($ = prompt)                //
// $ root -b -q burst-model.cc                                   //
//                                                               //
// Several standard histograms, depending on the user selections,//
// are automatically generated and may be written into a ROOT    //
// output file for later inspection, as indicated below.         //
//                                                               //
//--- Nick van Eijndhoven, IIHE-VUB, Brussel, March 2, 2020      //
///////////////////////////////////////////////////////////////////
{
 NcVersion version;
 version.Data();

 // The virtual lab for a cosmic transient phenomena study
 NcAstrolab2 lab;

 lab.SetExperiment("IceCube");
 lab.SetUT("11-04-2020","12:00:00.0",0); // Fixed fictative analysis date
 lab.SetRandomiser(-1); // Use the UT timestamp to generate a seed

 Float_t pi=acos(-1.);

 // User defined parameter settings
 Float_t bkgrate=-0.003/(2.*pi);
 lab.SetBurstParameter("Bkgrate",bkgrate);
 lab.SetBurstParameter("Nbkg",2);
 lab.SetBurstParameter("Grbnu",-0.02);
 lab.SetBurstParameter("Ezcor",1);
 lab.SetBurstParameter("Declmin",5);
 lab.SetBurstParameter("Declmax",85);
 lab.SetBurstParameter("Dtnu",500);
 lab.SetBurstParameter("Dtnus",-0.5);
 lab.SetBurstParameter("Datype",0);
 lab.SetBurstParameter("Dawin",5);
 lab.SetBurstParameter("Tunits",2);
 lab.SetBurstParameter("Tmin",-3600);
 lab.SetBurstParameter("Tmax",3600);

 // Obtain burst locations, durations etc.
 Bool_t scale=kFALSE;
 lab.GenBurstGCNdata(500,"GRB",scale); // Generate fictative burst GCN data

 // Generate burst signals according to the user defined parameters
 lab.GenBurstSignals();

 // Provide a listing of the first 10 stored burst positions
 cout << endl;
 lab.ListSignals("equ","J",1,"T",10);
 cout << endl;

 lab.ListBurstHistograms();
}