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

double GainCheck = -89.14;  // dBm/Hz
double RoomTemperature = 304.5; // K, temperature inside shielded room.

string filestem =
"Jun09-11-36-28_342776GHz";  // Fs = 25 MHz., -88.86, 304.1, 34.2776 GHz  (34.3024)
//"Jun08-01-07-54_342800GHz";  // Fs = 25 MHz., -88.975, 304.55, 34.280 GHz (34.3)
//"Jun08-11-47-14_342824GHz";  // Fs = 25 MHz., -88.94, 304.35, 34.2824 GHz (34.2976)
//"Jun14-01-31-01_342850GHz";  // Fs = 25 MHz., -88.5, 303.45, 34.2850 GHz (34.295)
//"Jun09-02-47-29_342877GHz";  // Fs = 25 MHz., -89.13, 304.6, 34.2877 GHz (34.2923)
//"Jun14-12-23-12_342900GHz";  // Fs = 25 MHz., -88.465, 303.35, 34.2900 GHz, 10.2k triggers.
//"Jun14-11-12-39_342925GHz";  // Fs = 25 MHz., -88.465, 302.85, 34.2925 GHz
//"Jun10-02-00-43_342950GHz";  // Fs = 25 MHz., -89.14, 304.5, 34.2950 GHz
//"Jun10-12-45-25_342975GHz";  // Fs = 25 MHz., -88.90, 304.5, 34.2975 GHz
//"Jun10-11-38-35_343003GHz";  // Fs = 25 MHz., -88.91, 304.1, 34.3003 GHz
//"Jun09-04-09-24_343025GHz";  // Fs = 25 MHz., -88.86, 304.6, 34.3025 GHz (? - check logbook)
//"Jun14-02-35-25_background"; // Fs = 25 MHz., -88.55, 303.6

//"Jul12-12-03-20"; // Fs = 50 MHz, MemDepth=10000.  Room temp. waveguide at input to room temp. receiver.

string histoname = filestem;
string filename="data/"+filestem+"_16bit.root"; 
string filename_background="data/Jun14-02-35-25_background_16bit.root";
const char *fname = filename.c_str();
const char *fname_background = filename_background.c_str();
const char *hname = histoname.c_str();



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->fitexample();
//ABC1->singlehisto();
//ABC1->PlotSomething();
//delete ABC1;




// Use Scope Class
//Scope *Scope1 = new Scope;
//Scope1->Print();
//Scope1->FittingTest();
//Scope1->Read(fname);
//Scope1->PlotGains();
//Scope1->PlotKnee();
//Scope1->PlotSomething();
//Scope1->PlotOneScopeTrace();
//Scope1->PlotNoiseTemperature(450, fname_background, 0, 100);  // 90 or 360 or 450 for Q, I, or Q+I, filename, scope trace, scope trace.
//Scope1->PlotIvsQ(fname_background, 0, 2);  // filename, scope trace, scope trace.
//Scope1->PlotIvsQfft(fname_background, 0, 10);  // filename, scope trace, scope trace.
//delete Scope1;

//}



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




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

}