DATCONSVDClusterLoaderAndPreparer.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 <string>
#include <vector>
 
namespace Belle2 {
  class ModuleParamList;
  class SVDCluster;
  class VxdID;
 
  class DATCONSVDClusterLoaderAndPreparer : public
    TrackFindingCDC::Findlet<const SVDCluster, const SVDCluster, std::pair<VxdID, std::pair<long, long>>,
        std::pair<VxdID, std::pair<long, long>>> {
    using Super = TrackFindingCDC::Findlet<const SVDCluster, const SVDCluster, std::pair<VxdID, std::pair<long, long>>,
          std::pair<VxdID, std::pair<long, long>>>;
 
  public:
    DATCONSVDClusterLoaderAndPreparer();
 
    void exposeParameters(ModuleParamList* moduleParamList, const std::string& prefix) override;
 
    void initialize() override;
 
    void apply(const std::vector<SVDCluster>& uClusters, const std::vector<SVDCluster>& vClusters,
               std::vector<std::pair<VxdID, std::pair<long, long>>>& uHits,
               std::vector<std::pair<VxdID, std::pair<long, long>>>& vHits) override;
 
  private:
    // Parameters
    uint m_param_maxClustersPerLayer = 200;
 
    std::array<uint, 8> m_nClusterPerLayer = {0};
 
    // ATTENTION: all the values below are hardcoded and taken from svd/data/SVD-Components.xml
    const std::array<int, 4> m_const_SVDRadii = {38990, 80000, 104000, 135150};
    const std::array<int, 4> m_const_RPhiShiftsOfLayers = { -4680, -10780, -11036, -19076};
    const std::array<double, 4> m_const_InitialAngle = {(342.195 - 360.) / 180. * M_PI, 8. / 180. * M_PI, -8. / 180. * M_PI, -4. / 180. * M_PI};
    const std::array<double, 4> m_const_AngleStep = {2. * M_PI / 7., M_PI / 5., M_PI / 6., M_PI / 8.};
    const std::array<double, 4> m_const_CosSlantedAngles = {1, cos(0.207694180987), cos(0.279252680319), cos(0.368264472171)};
 
    const std::array<int, 2> m_const_ZShiftL3 = {92350, -32650};
    const std::array<int, 3> m_const_ZShiftL4 = {149042, 24760, -100240};
    const std::array<int, 4> m_const_ZShiftL5 = {233754, 110560, -14440, -139440};
    const std::array<int, 5> m_const_ZShiftL6 = {303471, 182060, 57060, -67940, -192940};
 
  };
}