Changeset 19951 for trunk/Mars/hawc/fit_singles.C

Timestamp:

04/29/20 21:48:05 (5 years ago)

Author:

giangdo

Message:

script to synchronise the data from HAWC's Eye and HAWC

File:

• trunk/Mars/hawc/fit_singles.C (modified) (7 diffs)

trunk/Mars/hawc/fit_singles.C

@@ -93 +93 @@
 
 // Print function parameters
-void PrintFunc(TF1 &f, float integration_window=30)
+void PrintFunc(TF1 &f, float integration_window=40)
 {
     //cout << "--------------------" << endl;
@@ -113 +113 @@
 // The parameters for the function are the filenam from the gain extraction
 // and the output filename
-int fit_singles(const char *filename = "20191027_006_006.root",
-                const char *outfile  = "20190227_006_006-fit.root",
+int fit_singles(const char *filename = "/home/giangdo/Documents/Master/MA/pulse/data_single_pe/extracted/20191027_006_006.root",
+                const char *outfile  = "/home/giangdo/Documents/Master/MA/pulse/data_single_pe/extracted/20191027_006_006_fit.root",
                 bool fixednoise=false)
 {
     // Read the mapping between bias channels and hardware channels
     PixelMap pmap;
-    if (!pmap.Read("../hawc/FAMOUSmap171218.txt"))
-    {
-        cout << "FAMOUSmap171218.txt not found." << endl;
+    if (!pmap.Read("/home/giangdo/Documents/Master/MA/Software/Mars/hawc/HAWCsEyemap181214.txt"))
+    {
+        cout << "HAWCsEyemap181214.txt not found." << endl;
         return 1;
     }
 
     // It is more correct to fit the integral, but this is super
-    // slow, especially for 1440 pixel. To get that, one would have
+    // slow. To get that, one would have
     // to analytically integrate and fit this function.
-    // (Currently the integral is switched off for the 1440 individual
-    // spectra in all cases)
     bool fast = false; // Switch off using integral
 
     // Values which should be read from the file but not available atm
-    Int_t integration_window = 30;
+    Int_t integration_window = 40;
 
     // Map for which pixel shall be plotted and which not
@@ -266 +264 @@
     TH1F hCoeffR2    ("CoeffR2",    "Coefficient R",          90, -1,     2);
 
-    // Histigram for the sum of all spectrums
-    TH1F hSum("Sum1", "Signal spectrum of all pixels",
+    // Histogram for the sum of all spectrums
+    TH1F hSum("Sum1", "Signal spectrum of all pixels (excl. blind pixels)",
               hsignal->GetNbinsY(), hsignal->GetYaxis()->GetXmin(),  hsignal->GetYaxis()->GetXmax());
     hSum.SetXTitle("pulse integral [mV sample]");
@@ -275 +273 @@
 
     // Histogram for the sum of all pixels excluding the ones with faulty fits
-    TH1F hSumClear1("SumC1", "Signal spectrum of all pixels (incl TM)",
+    TH1F hSumClear1("SumC1", "Signal spectrum of all pixels (excl. faulty and blind pixels)",
                     hsignal->GetNbinsY(), hsignal->GetYaxis()->GetXmin(),  hsignal->GetYaxis()->GetXmax());
     hSumClear1.SetXTitle("pulse integral [mV sample]");
@@ -303 +301 @@
     hPulse.SetStats(false);
 
-    // Spektrum for the normalized individual spectra
-    TH1F hSumScale1("SumScale1", "Signal spectrum of all pixels (incl TM)",  120, 0.05, 12.05);
+    // Spectrum for the normalized individual spectra
+    TH1F hSumScale1("SumScale1", "Signal spectrum of all pixels (excl. blind pixels)",  120, 0.05, 12.05);
     hSumScale1.SetXTitle("pulse integral [pe]");
     hSumScale1.SetYTitle("Rate");
@@ -332 +330 @@
         if (usePixel[pixel]==0)
             continue;
+
+        //exclude the pixels without cones
+        if (pixel == 61)
+           continue;
+        if (pixel == 62)
+           continue;
+        if (pixel == 63)
+           continue;
+
+
+
 
         TH1D *hist = hsignal->ProjectionY("proj", pixel+1, pixel+1);
@@ -434 +443 @@
         if (usePixel[pixel]==0)
             continue;
+
+        // Exclude pixels without cones
+        if (pixel == 61)
+           continue;
+        if (pixel == 62)
+           continue;
+        if (pixel == 63)
+           continue;
 
         //Projectipon of a certain Pixel out of Ampl.Spectrum