release-05-02-19/doxygen/TOPRingPlotterModule_8h_source.html

/**************************************************************************

* BASF2 (Belle Analysis Framework 2)                                     *

* Copyright(C) 2020 - Belle II Collaboration                             *

*                                                                        *

* Author: The Belle II Collaboration                                     *

* Contributors: Umberto Tamponi                                          *

*                                                                        *

* This software is provided "as is" without any warranty.                *

**************************************************************************/


#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>

#include <top/reconstruction/TOPreco.h>


using namespace Belle2::Variable;


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;

    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;

    TH2F* m_hitMapMCK = nullptr;

    TH2F* m_hitMapMCPI =  nullptr;

    TH2F* m_hitMapMCP =  nullptr;

    TH2F* m_hitMapMCE =  nullptr;

    TH2F* m_hitMapMCMU =  nullptr;

    TH2F* m_pdfAsHisto =  nullptr;

  };

}