eclLOMModule.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
/* Basf2 headers. */
#include <framework/datastore/StoreArray.h>
#include <framework/core/Module.h>
#include <mdst/dataobjects/MCParticle.h>
#include <trg/ecl/dataobjects/TRGECLWaveform.h>
 
/* ROOT headers. */
#include <TTree.h>
#include <TFile.h>
#include <TH1D.h>
#include <TH2D.h>
 
namespace Belle2 {
  namespace ECL {

    class ECLLOMModule : public Module {
    public:

      ECLLOMModule();

      virtual ~ECLLOMModule() override;

      virtual void initialize() override;

      virtual void event() override;

      virtual void terminate() override;
 
    protected:
 
    private:

      void get_MCparticles();

      void get_waveforms();

      void calculate_discr_output();

      void calculate_amplitudes();

      bool calculate_BE_quality(int iSample);

      bool calculate_FE_quality(int iSample);

      void calculate_coincidence(int iSample);

      void clear_lom_data();
 
      //module parameters
      std::string m_testFileName; 
      double m_thresholdFE = 0.0; 
      double m_thresholdBE = 0.0; 
      double m_thresholdBkg = 0.0; 
      double m_discrTime = 0.0; 
      bool m_includeInnerFE = false; 
      bool m_saveSignal = false; 
 
      //input data
      int m_evtNum; 
      double m_BE_Waveform_100ns[16][64] = {}; 
      double m_FE_Waveform_100ns[16][64] = {}; 
      double m_mcen[2] = {}; 
      double m_mcth[2] = {}; 
      double m_mcph[2] = {}; 
      double m_com_en[2] = {}; 
      double m_com_th[2] = {}; 
      double m_com_ph[2] = {}; 

      StoreArray<MCParticle> m_MCParticles;

      StoreArray<TRGECLWaveform> m_TrgEclWaveforms;
 
      //important output
      bool m_isBhabha = false; 
 
      int m_BhNum = 0; 
      double m_BE_Amplitude[16] = {}; 
      double m_FE_Amplitude[16] = {}; 
      double m_BESum_Amplitude[16] = {}; 
      double m_FESum_Amplitude[16] = {}; 
      double m_BE_Pedal[16] = {}; 
      double m_FE_Pedal[16] = {}; 
      double m_FESum_MaxAmp = 0.0; 
      double m_BESum_MaxAmp = 0.0; 
      int m_FESum_MaxId = 0; 
      int m_BESum_MaxId = 0; 
 
      //validation output
      TH2D* m_h2Coin = nullptr; 
      TH2D* m_h2SumCoin =
        nullptr; 
      TH2D* m_h2FEAmp = nullptr; 
      TH2D* m_h2BEAmp = nullptr; 
      TH1D* m_h1BEHits =
        nullptr; 
      TH1D* m_h1FEHits =
        nullptr; 
 
      //internal variables
      TFile* m_testfile = nullptr; 
      TTree* m_testtree = nullptr; 
      int m_NSamples = 0; 
      double m_BE_Waveform_10ns[16][631] = {}; 
      double m_FE_Waveform_10ns[16][631] = {}; 
      double m_BESum_Waveform_10ns[16][631] = {}; 
      double m_FESum_Waveform_10ns[16][631] = {}; 
      bool m_FESum_Discr[16][631] = {}; 
      bool m_BESum_Discr[16][631] = {}; 
      bool m_FEQual_Discr[16][631] = {}; 
      bool m_BEQual_Discr[16][631] = {}; 
      bool m_isBhabhaPatternFE = false; 
      bool m_isBhabhaPatternBE = false; 
      int m_CoincidenceMatrix[16][16] = {}; 
      int m_SumCoincidenceMatrix[16][16] = {}; 
      int m_CoincidenceCounterMatrix[16][16] = {}; 
      int m_SumCoincidenceCounterMatrix[16][16] = {}; 
    };
  }//namespace ECL
}//namespace Belle2