TOPTBCComparatorModule.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Umberto Tamponi (tamponi@to.infn.it)                     *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <framework/core/HistoModule.h>
#include "TH1F.h"
#include "TH2F.h"
#include <string>
#include <vector>
 
namespace Belle2 {
  class TOPTBCComparatorModule : public HistoModule {
  public:
    TOPTBCComparatorModule();
 
    ~TOPTBCComparatorModule() {};
 
    void defineHisto() override;
 
    void initialize() override;
 
    void beginRun() override;
 
    void event() override;
 
    void endRun() override;
 
    void terminate() override;
 
    int analyzeCalFile();
 
    int makeComparisons();
 
 
 
    int parseSlotAndScrodIDfromFileName(std::string);
 
    int parseInputDirectoryLine(std::string);
 
    TH1F* calculateHistoRatio(TH1F*, TH1F*, TH1F*);
 
    TH2F* calculateHistoRatio(TH2F*, TH2F*, TH2F*);
 
 
 
 
  private:
    // steering parameters
    std::string m_inputDirectoryList =
      ""; 
    bool m_compareToPreviousSet =
      true; 
    std::string m_outputFile = ""; 
    short m_minCalPulses = 200; 
    short m_numSamples = 256; 
    // utilities
    std::string m_calSetDirectory; 
    TFile*  m_calSetFile =
      nullptr; 
    std::string m_calSetLabel; 
    short m_slotID = -1; 
    short m_boardstackID = -1; 
    short m_scrodID = -1; 
    short m_calSetID = 0; 
    short m_totCalSets = 0; 
    // Delta T plots, slot-by-slot
    std::vector<TH1F*>
    m_slotAverageDeltaT[16]; 
    std::vector<TH1F*>
    m_slotSigmaDeltaT[16]; 
    std::vector<TH2F*>
    m_slotAverageDeltaTMap[16]; 
    std::vector<TH2F*>
    m_slotSigmaDeltaTMap[16]; 
    // Average timing plots, all the detector at once
    std::vector<TH1F*>
    m_topAverageDeltaT; 
    std::vector<TH1F*>
    m_topSigmaDeltaT; 
    // Occupancy plots, slot-by-slot
    std::vector<TH1F*>
    m_slotSampleOccupancy[16]; 
    std::vector<TH1F*>
    m_slotEmptySamples[16]; 
    std::vector<TH2F*>
    m_slotSampleOccupancyMap[16]; 
    std::vector<TH2F*>
    m_slotEmptySamplesMap[16]; 
    // Occupancy plot, all the detector at once
    std::vector<TH1F*>
    m_topSampleOccupancy; 
    // Delta T ratio plots, slot-by-slot
    std::vector<TH1F*>
    m_slotAverageDeltaTComparison[16]; 
    std::vector<TH1F*>
    m_slotSigmaDeltaTComparison[16]; 
    std::vector<TH2F*>
    m_slotAverageDeltaTMapComparison[16]; 
    std::vector<TH2F*>
    m_slotSigmaDeltaTMapComparison[16]; 
    // Delta T ratio plots, whole detector
    std::vector<TH1F*>
    m_topAverageDeltaTComparison; 
    std::vector<TH1F*>
    m_topSigmaDeltaTComparison; 
    std::vector<TH1F*>
    m_topSampleOccupancyComparison; 
  };
 
} //namespace Belle2