TOPRingPlotterModule.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
 
#include <framework/core/Module.h>
#include <TTree.h>
#include <TH2F.h>
#include <TFile.h>
#include <analysis/VariableManager/Manager.h>
#include <mdst/dataobjects/Track.h>
 
namespace Belle2 {
  class TOPRingPlotterModule : public Module {
 
  public:
 
    TOPRingPlotterModule();
 
    void initialize() override;
 
 
    void event() override;
 
 
    void terminate() override;
 
 
  private:
 
    void fillPDF(Belle2::Const::ChargedStable, const Track*, TH2F*, short*, float*,  int&, int&);
 
    void resetTree();
 
 
    static const int m_MaxPhotons = 5000; 
    static const int m_MaxPDFPhotons = 5000; 
    // steering parameters
    std::string m_particleList = "pi+:all";  
    std::string m_outputName = "TOPRings.root";  
    int m_toyNumber = 1; 
    bool m_saveHistograms = false;  
    bool m_saveLLScan = false; 
    TFile* m_outputFile = nullptr;   
    TTree* m_tree = nullptr;   
    std::vector<double> m_branchAddresses = {};
    std::vector<Variable::Manager::FunctionPtr> m_functions = {};
    std::vector<std::string> m_variables = {};
 
    std::vector<short> m_pdgHyp = {11, 13, 211, 321, 2212}; 
    float m_digitTime[m_MaxPhotons] = {0}; 
    short m_digitChannel[m_MaxPhotons] = {0}; 
    short m_digitPixel[m_MaxPhotons] = {0}; 
    short m_digitPixelCol[m_MaxPhotons] = {0}; 
    short m_digitPixelRow[m_MaxPhotons] = {0}; 
    float m_digitAmplitude[m_MaxPhotons] = {0}; 
    float m_digitWidth[m_MaxPhotons] = {0}; 
    short m_digitSlot[m_MaxPhotons] = {0}; 
    short m_digitBoardstack[m_MaxPhotons] = {0}; 
    short m_digitCarrier[m_MaxPhotons] = {0}; 
    short m_digitAsic[m_MaxPhotons] = {0}; 
    short m_digitQuality[m_MaxPhotons] = {0}; 
    short m_digitASICChannel[m_MaxPhotons] = {0}; 
    short m_digitPMTNumber[m_MaxPhotons] = {0}; 
    short m_nDigits = 0; 
    short m_pdfPixelE[m_MaxPDFPhotons] = {0}; 
    float m_pdfTimeE[m_MaxPDFPhotons] = {0}; 
    short m_pdfPixelMU[m_MaxPDFPhotons] = {0}; 
    float m_pdfTimeMU[m_MaxPDFPhotons] = {0}; 
    short m_pdfPixelPI[m_MaxPDFPhotons] = {0}; 
    float m_pdfTimePI[m_MaxPDFPhotons] = {0}; 
    short m_pdfPixelK[m_MaxPDFPhotons] = {0}; 
    float m_pdfTimeK[m_MaxPDFPhotons] = {0}; 
    short m_pdfPixelP[m_MaxPDFPhotons] = {0}; 
    float m_pdfTimeP[m_MaxPDFPhotons] = {0}; 
    int m_pdfSamplesP = 0; 
    int m_pdfToysP = 0; 
    int m_pdfSamplesK = 0; 
    int m_pdfToysK = 0; 
    int m_pdfSamplesPI = 0; 
    int m_pdfToysPI = 0; 
    int m_pdfSamplesMU = 0; 
    int m_pdfToysMU = 0; 
    int m_pdfSamplesE = 0; 
    int m_pdfToysE = 0; 
    int m_nScanPoints = 0; 
    float m_scanMass[10000] = {0}; 
    float m_scanLL[10000] = {0}; 
    TH2F* m_hitMapMCK = nullptr; 
    TH2F* m_hitMapMCPI =  nullptr; 
    TH2F* m_hitMapMCP =  nullptr; 
    TH2F* m_hitMapMCE =  nullptr; 
    TH2F* m_hitMapMCMU =  nullptr; 
    TH2F* m_pdfAsHisto =  nullptr; 
  };
}