Changeset 19904 for trunk/FACT++/src

Timestamp:

12/15/19 16:23:56 (5 years ago)

Author:

tbretz

Message:

Implemented the Flux Sigma and control plots to check if a point should be replaced by an UL.

File:

• trunk/FACT++/src/spectrum.cc (modified) (4 diffs)

trunk/FACT++/src/spectrum.cc

@@ -2060 +2060 @@
             map<size_t, double> rolke_int;
 
+            if (verbose>0)
+                cout << "Bin     Excess     Significance          Flux ErrFlux" << endl;
+
             for (auto ir=res13.cbegin(); ir!=res13.cend(); ir++)
             {
@@ -2094 +2097 @@
                 if (verbose>0)
                 {
-                    cout << setw(5)  << row["center"]  << ": ";
-                    cout << setw(10) << row["Excess"]  << " ";
-                    cout << setw(10) << row["Flux"]    << " ";
-                    cout << setw(10) << row["ErrFlux"] << " ";
-                    cout << setw(10) << row["Significance"] << '\n';
+                    cout << setw(5)  << row["center"] << ":";
+                    cout << " " << setw(10) << row["Excess"];
+                    cout << " " << setw(10) << row["Significance"];
+                    cout << " " << setw(10) << row["Flux"];
+                    cout << " " << setw(10) << row["ErrFlux"];
+                    cout << endl;
                 }
             }
@@ -2169 +2173 @@
             WriteHistogram(connection, hist);
 
+            hist.name = "SigmaFlux";
+            hist.v    = "SigmaFlux";
+            hist.err  = "";
+            WriteHistogram(connection, hist);
+
             if (*ib=="Energy")
             {
@@ -2177 +2186 @@
                 hist.v     = "FluxI";
                 hist.err   = "ErrFluxI";
+                WriteHistogram(connection, hist);
+
+                hist.dir   = "Data/Energy/Integral";
+                hist.name  = "SigmaFlux";
+                hist.v     = "SigmaFluxI";
+                hist.err   = "";
                 WriteHistogram(connection, hist);