DATCONSVDClusterizer.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/trackFindingCDC/findlets/base/Findlet.h>
#include <framework/datastore/StoreArray.h>
#include <vxd/geometry/GeoCache.h>
#include <svd/calibration/SVDNoiseCalibrations.h>
 
#include <string>
#include <vector>
#include <fstream>
#include <numeric>
 
namespace Belle2 {
  class ModuleParamList;
  class DATCONSVDDigit;
  class SVDCluster;
  struct DATCONSVDClusterCandidate;
 
  class DATCONSVDClusterizer : public TrackFindingCDC::Findlet<const DATCONSVDDigit, SVDCluster> {
    using Super = TrackFindingCDC::Findlet<const DATCONSVDDigit, SVDCluster>;
 
    typedef std::map<int, float> SimpleSVDNoiseMap;
 
  public:
    DATCONSVDClusterizer();
 
    void exposeParameters(ModuleParamList* moduleParamList, const std::string& prefix) override;
 
    void initialize() override;
 
    void beginRun() override;
 
    void apply(const std::vector<DATCONSVDDigit>& digits, std::vector<SVDCluster>& clusters) override;
 
  private:
    void fillDATCONSVDNoiseMap();
 
    float calculateSNR(DATCONSVDDigit digit);
 
    void saveCluster(DATCONSVDClusterCandidate& clusterCand, std::vector<SVDCluster>& clusters);
 
    SVDNoiseCalibrations m_NoiseCal;
 
    // Parameters
    SimpleSVDNoiseMap m_svdNoiseMap;
 
    std::string m_param_noiseMapfileName = "noiseMap.txt";
    bool m_param_writeNoiseMapsToFile = false;
    bool m_param_isU = true;
 
    bool m_param_saveClusterToDataStore = false;
    std::string m_param_storeSVDClustersName = "DATCONSVDClusters";
    StoreArray<SVDCluster> m_storeSVDClusters;
 
    unsigned short m_param_maxiClusterSize = 20;
    float m_param_noiseCut = 4;
    float m_param_requiredSNRstrip = 5;
    float m_param_requiredSNRcluster = 5;
 
    const VXD::GeoCache& m_geoCache = VXD::GeoCache::getInstance();
 
  };
}