SegmentNetworkAnalyzerModule.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 <tracking/trackFindingVXD/segmentNetwork/DirectedNodeNetworkContainer.h>
#include <tracking/trackFindingVXD/segmentNetwork/Segment.h>
#include <tracking/trackFindingVXD/segmentNetwork/TrackNode.h>
 
#include <framework/datastore/StoreObjPtr.h>
#include <framework/datastore/StoreArray.h>
#include <mdst/dataobjects/MCParticle.h>
#include <framework/core/Module.h>
 
#include <TFile.h>
#include <TTree.h>
 
#include <string>
 
namespace Belle2 {
  class SegmentNetworkAnalyzerModule : public Module {
 
  public:

    struct RootVariables {
      std::vector<double> phi{}; 
      std::vector<double> theta{}; 
      std::vector<double> pT{}; 
      std::vector<int> signal{}; 
      std::vector<int> passed{}; 
      std::vector<int> pdg{}; 
      std::vector<int> virtualIP{}; 
      unsigned networkSize{}; 
      unsigned networkConnections{}; 
    };

    SegmentNetworkAnalyzerModule();

    void initialize() override;

    void event() override;

    void terminate() override;
 
  private:

    std::string m_PARAMnetworkName;

    std::string m_PARAMrootFileName;

    Belle2::StoreObjPtr<Belle2::DirectedNodeNetworkContainer> m_network;

    Belle2::StoreArray<Belle2::MCParticle> m_mcParticles;

    TFile* m_rFilePtr;

    TTree* m_treePtr;

    RootVariables m_rootVariables;

    void makeBranches();

    void analyzeCombination(const Belle2::Segment<Belle2::TrackNode>& outer, const Belle2::Segment<Belle2::TrackNode>& inner,
                            bool passed);

    template<typename EntryType, typename MetaInfoType>
    static size_t getNConnections(Belle2::DirectedNodeNetwork<EntryType, MetaInfoType>& network);
 
  };

}