KLMTimeCollectorModule.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
 
/* KLM headers. */
#include <klm/calibration/KLMTimeAlgorithm.h>
#include <klm/bklm/geometry/GeometryPar.h>
#include <klm/dataobjects/KLMHit2d.h>
#include <klm/dbobjects/KLMChannelStatus.h>
#include <klm/eklm/geometry/GeometryData.h>
#include <klm/eklm/geometry/TransformData.h>
 
/* Basf2 headers */
#include <calibration/CalibrationCollectorModule.h>
#include <framework/dataobjects/EventT0.h>
#include <framework/datastore/StoreArray.h>
#include <framework/datastore/StoreObjPtr.h>
#include <mdst/dataobjects/Track.h>
#include <tracking/dataobjects/ExtHit.h>
 
/* ROOT headers */
#include <TH1D.h>
#include <TH1I.h>
#include <TH2D.h>
#include <TTree.h>
 
namespace Belle2 {
  class KLMTimeCollectorModule : public CalibrationCollectorModule {
 
  public:
 
    KLMTimeCollectorModule();
 
    virtual ~KLMTimeCollectorModule();
 
    void prepare() override;
 
    void collect() override;
 
    void finish() override;
 
  private:
 
    void collectRPC(RelationVector<KLMHit2d>&);
 
    void collectScint(RelationVector<KLMHit2d>&);
 
    void collectScintEnd(const RelationVector<KLMHit2d>&);
 
    std::pair<ExtHit*, ExtHit*> matchExt(
      KLMChannelNumber channelID, std::multimap<unsigned int, ExtHit>&);
 
    void storeDistDiff(ROOT::Math::XYZVector&);
 
    bool m_useEvtT0;
 
    bool m_IgnoreBackwardPropagation;
 
    std::string m_inputListName;
 
    std::multimap<unsigned int, ExtHit> m_vExtHits_RPC;
 
    std::multimap<unsigned int, ExtHit> m_vExtHits;
 
    StoreArray<Track> m_tracks;
 
    StoreObjPtr<EventT0> m_eventT0;
 
    KLMTimeAlgorithm::Event m_Event;
 
    const bklm::GeometryPar* m_geoParB;
 
    const EKLM::GeometryData* m_geoParE;
 
    EKLM::TransformData* m_TransformData;
 
    const KLMElementNumbers* m_elementNum;
 
    DBObjPtr<KLMChannelStatus> m_channelStatus;
 
    /* Monitor histograms. */
 
    TH1D* m_HeventT0_0;
 
    TH1D* m_HeventT0_1;
 
    TH1I* m_HnumTrack;
 
    TH1I* m_HnKLMHit2dOfTrack;
 
    TH1D* m_HpositionDiff;
 
    TH1D* m_HpositionXDiff;
 
    TH1D* m_HpositionYDiff;
 
    TH1D* m_HpositionZDiff;
 
    TH2D* m_HflyTimeB;
 
    TH2D* m_HflyTimeE;
 
    /* Sum number of digits collected. */
 
    TH1I* m_HnumDigit_scint_end;
 
    TH1I* m_HnumDigit_scint;
 
    TH1I* m_HnumDigit_rpc;
 
    TTree* m_outTree;
 
  };
}