KinkFinderModule.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
//Object with performing the actual algorithm:
#include <tracking/kinkFinding/fitter/KinkFitter.h>
 
#include <tracking/dataobjects/RecoTrack.h>
#include <mdst/dataobjects/Track.h>
 
#include <framework/datastore/StoreArray.h>
#include <framework/core/Module.h>
 
#include <tracking/dbobjects/KinkFinderParameters.h>
 
#include <string>
#include <memory>
 
namespace Belle2 {
  class KinkFinderModule : public Module {
 
  public:
 
    KinkFinderModule();
 
    ~KinkFinderModule() override = default;
 
    void initialize() override;
 
    void beginRun() override;
 
    void event() override;
 
    void terminate() override;
 
  private:
 
    void fitAndStore(const Track* trackMother, const Track* trackDaughter, const short filterFlag);
 
    bool kinkFitterModeSplitTrack();
 
    bool ifInCDC(const ROOT::Math::XYZVector& pos);
 
    bool preFilterMotherTracks(Track const* const track);
 
    bool preFilterDaughterTracks(Track const* const track);
 
    bool preFilterTracksToSplit(Track const* const track);
 
    short isTrackPairSelected(const Track* motherTrack, const Track* daughterTrack);
 
    struct CDCForwardBackwardWallLine {
 
      void setLine(const double tangent, const double offset)
      {
        m_tangent = tangent;
        m_offset = offset;
      }
 
      double getLine(const double x) const
      {
        return m_tangent * x + m_offset;
      }
 
    private:
      double m_tangent; 
      double m_offset; 
    };
 
    std::string m_arrayNameTrack;  
    StoreArray<Track> m_tracks;  
 
    std::unique_ptr<KinkFitter> m_kinkFitter;  
    std::string m_arrayNameRecoTrack;  
    std::string m_arrayNameCopiedRecoTrack;  
    std::string m_arrayNameTFResult;  
    std::string m_arrayNameKink;  
 
    DBObjPtr<KinkFinderParameters> m_kinkFinderParameters; 
 
    // CDC and SVD geometry variables
    CDCForwardBackwardWallLine m_cdcForwardBottomWall; 
    CDCForwardBackwardWallLine m_cdcForwardTopWall; 
    CDCForwardBackwardWallLine m_cdcBackwardBottomWall; 
    CDCForwardBackwardWallLine m_cdcBackwardTopWall; 
    static constexpr double m_cdcInnerWithFirstLayerWall = 15; 
    static constexpr double m_cdcInnerWallWithoutFirstLayer = 17; 
    static constexpr double m_cdcInnerWallWithoutFirstTwoLayers = 18; 
    static constexpr double m_cdcOuterWall = 112; 
    static constexpr double m_svdBeforeOuterLayer = 12; 
 
    // counter for KinkFinder statistics
    unsigned int m_allStored = 0;    
    unsigned int m_f1Stored = 0;    
    unsigned int m_f2Stored = 0;    
    unsigned int m_f3Stored = 0;    
    unsigned int m_f4Stored = 0;    
    unsigned int m_f5Stored = 0;    
  };
}