Changeset 19951 for trunk/Mars/hawc/plot_trace.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/plot_trace.C (modified) (8 diffs)

trunk/Mars/hawc/plot_trace.C

@@ +1 @@
+#include <TStyle.h>
+#include <TCanvas.h>
+#include <TSystem.h>
+#include <TF1.h>
+#include <TProfile.h>
+#include <TProfile2D.h>
+#include <TMath.h>
+#include <TGraph.h>
+#include <TLine.h>
+#include <TFitResultPtr.h>
+#include <TFitResult.h>
+#include <TFile.h>
+#include <TLine.h>
+
+#include <cstdio>
+#include <stdio.h>
+#include <stdint.h>
+
+#include "Getline.h"
+#include "MH.h"
+#include "MArrayI.h"
+#include "MDirIter.h"
+#include "MFillH.h"
+#include "MEvtLoop.h"
+#include "MCamEvent.h"
+#include "MGeomApply.h"
+#include "MTaskList.h"
+#include "MParList.h"
+#include "MContinue.h"
+#include "MBinning.h"
+#include "MDrsCalibApply.h"
+#include "MDrsCalibration.h"
+#include "MRawFitsRead.h"
+#include "MStatusDisplay.h"
+#include "MTaskInteractive.h"
+#include "MPedestalSubtractedEvt.h"
+#include "MGeomCamFAMOUS.h"
+#include "MRawRunHeader.h"
+#include "MPedestalCam.h"
+#include "MPedestalPix.h"
+#include "MParameters.h"
+
 // ==========================================================================
 // ============ see plot_callisto function at the end of the file ===========
@@ -21 +63 @@
 
 // Create a histogram for the trace
-TH1D hist("Trace", "Waveform", 1024, -0.25, 511.75);
+//TH1D hist("Trace", "Waveform", 1024, -0.25, 511.75); //range = [-0.25,511.75]ns
+TH1D hist("Trace", "Waveform (Pixel 5, Event 1)", 1024, -0.5, 1024-0.5);
 
 // All 'data members' that are required globally
@@ -34 +77 @@
 {
     fEvt = (MPedestalSubtractedEvt*)plist->FindObject("MPedestalSubtractedEvt");
-    if (!fEvt)
-    {
-        gfLog << err << "MPedestalSubtractedEvt not found... aborting." << endl;
-        return kALSE;
-    }
+
+    //if (!fEvt)
+    //{
+    //    gfLog << err << "MPedestalSubtractedEvt not found... aborting." << endl;
+    //    return kFALSE;
+    //}
 
     // Create a canvas and plot the histogram into it
-    c = new TCanvas;
+    c = new TCanvas();
     hist.SetStats(kFALSE);
     hist.Draw();
-
-    TF1 constant("zero", "[0]", -0.25, 511.75);
+    hist.SetXTitle("Samples");
+    hist.SetYTitle("Amplitude [a.u.]");
+
+    //TF1 constant("zero", "[0]", -0.25, 511.75);
+    TF1 constant("zero", "[0]", -0.5, 1024);
     constant.SetParameter(0, 0);
     constant.SetLineColor(kBlack);
@@ -59 +106 @@
     constant.DrawCopy("same");
 
+
     return kTRUE;
 }
@@ -64 +112 @@
 Int_t Process()
 {
-    int pixel_index = 0;
+    int pixel_index = 6; //plot trace only of this pixel
+    int skipstart = 5;
+    int skipend = 10;
 
     // Number of samples per pixel (usually 1024)
@@ -73 +123 @@
     hist.Reset();
     for (int i=0; i<numsamples; i++)
-        hist.Fill(i*0.5, cal[i]); // Convert sample to nano-second
+        hist.Fill(i, cal[i]); // 0.5ns = 1 sample
 
     // Reset min/max range
-    hist.SetMinimum();
-    hist.SetMaximum();
-
-    // Set a reasonable range
-    if (hist.GetMinimum()>-50)
-        hist.SetMinimum(-50);
-    if (hist.GetMaximum()<250)
-        hist.SetMaximum(250);
+    hist.SetMinimum(-50);
+    hist.SetMaximum(50);
+
+    //// Set a reasonable range
+    //if (hist.GetMinimum()>-30)
+    //    hist.SetMinimum(-20);
+    //if (hist.GetMaximum()<60)
+    //    hist.SetMaximum(50);
+
+
+
+    TLine *line = new TLine(skipstart, 50, skipstart, -50);
+    line->SetLineStyle(kSolid);
+    line->SetLineColor(kRed);
+    line->Draw();
+
+    TLine *line2 = new TLine(1024-skipend, 50, 1024-skipend, -50);
+    line2->SetLineStyle(kSolid);
+    line2->SetLineColor(kRed);
+    line2->Draw();
 
     // Signal root to update the canvas/pad
     c->Modified();
     c->Update();
+    //c->SaveAs("/home/giangdo/Documents/Master/MA/pulse/data_single_pe/traces/pixel3/event_0.png");
+
 
     bool rc = HandleInput();
@@ -131 +195 @@
 
     // mapping file (found in Mars/hawc)
-    const char *mmap = usemap ? "../hawc/FAMOUSmap171218.txt" : NULL;
+    const char *mmap = usemap ? "/home/giangdo/Documents/Master/MA/Software/Mars/hawc/HAWCsEyemap181214.txt" : NULL;
 
     // ======================================================
@@ -193 +257 @@
     MEvtLoop loop(gSystem->BaseName(datafile));
     loop.SetParList(&plist);
- 
-    if (!loop.Eventloop())
+
+    if (!loop.Eventloop(100))
         return 4;