double Fs = 25.e6; // sampling frequency
double MemoryDepth = 10000.; // memory depth
double NBins = 10000.;

//const char *fname = "Apr05-02-24-08_16bit.root";  // Fs = 200 MHz. no LPF.
//const char *fname = "May19-12-52-14_16bit.root";  // Fs = 50 MHz. LPF.
//const char *fname = "data/May19-12-59-28_16bit.root";  // Fs = 50 MHz. no LPF.
//const char *fname = "data/Jun02-02-39-16_16bit.root";  // Fs = 25 MHz
//const char *fname = "data/Jun02-03-03-17_16bit.root";  // Fs = 25 MHz, door closed.
//const char *fname = "data/Jun03-03-08-41_16bit.root";  // Fs = 25 MHz.
const char *fname = "data/Jun10-12-45-25_16bit.root";  // Fs = 25 MHz.

int run()
{


// Check if (Oscillo)Scope is already loaded
if (!TClass::GetDict("Scope")) {
gROOT->ProcessLine(".L ScopeClass.C");
//gROOT->ProcessLine(".L CryoClass.C");
//gROOT->ProcessLine(".L ABCClass.C");
}

//Use ABC Class
//ABC *ABC1 = new ABC;
//ABC1->singlehisto();
//ABC1->PlotSomething();
//delete ABC1;




// Use Scope Class
Scope *Scope1 = new Scope;
//Scope1->Print();
//Scope1->PlotSomething();
Scope1->FFTScopeTraces();
//Scope1->PlotOneScopeTrace();
//Scope1->PlotJustNoiseTemperature();
//Scope1->PlotNoiseTemperature();
//Scope1->PlotNoiseTemperatureFull();
//Scope1->PlotPowerOverlap();
delete Scope1;


// Use Cryo Class
//Cryo *Cryo2 = new Cryo;
//Cryo2->Print();
//Cryo2->GetTemperature();
//Cryo2->GetRoomTemperature();
//delete Cryo2;


printf("Back in main.\n");
return 0;

}