ToPXDExtrapolator.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>
#include <cmath>
 
namespace Belle2 {
  class ModuleParamList;
  class VxdID;
  class PXDIntercept;

  class ToPXDExtrapolator : public TrackFindingCDC::Findlet<const std::pair<double, double>, const std::pair<double, double>,
    std::pair<VxdID, long>, std::pair<VxdID, long>> {
    using Super = TrackFindingCDC::Findlet<const std::pair<double, double>, const std::pair<double, double>,
          std::pair<VxdID, long>, std::pair<VxdID, long>>;
 
  public:
    ToPXDExtrapolator();

    void exposeParameters(ModuleParamList* moduleParamList, const std::string& prefix) override;

    void initialize() override;

    void apply(const std::vector<std::pair<double, double>>& uTracks, const std::vector<std::pair<double, double>>& vTracks,
               std::vector<std::pair<VxdID, long>>& uExtrapolations, std::vector<std::pair<VxdID, long>>& vExtrapolations) override;
 
  private:
    void extrapolateUTrack(const double trackPhi, const double trackRadius, const uint layer,
                           std::vector<std::pair<VxdID, long>>& uExtrapolations);

    void extrapolateVTrack(const long tanLambda, const uint layer,
                           std::vector<std::pair<VxdID, long>>& vExtrapolations);

    double m_param_phiCutL1 = M_PI / 3.;
    double m_param_phiCutL2 = M_PI / 4.;

    bool m_param_createPXDIntercepts = false;
    std::string m_param_PXDInterceptStoreArrayName = "DATCONFPGAPXDIntercepts";
    StoreArray<PXDIntercept> m_pxdIntercepts;
 
    // ATTENTION: hard coded values taken and derived from pxd/data/PXD-Components.xml
    const long centerZShiftLayer1[2] = {36825500, -8825500};
    const long centerZShiftLayer2[2] = {50145500, -12145500};
    const long layerRadius[2] = {14285, 22121};
    const long sensorLength[2] = {44800000, 61440000};
    const long sensorMinY  = -3600000;
    const long sensorMaxY  = 8900000;
    const long shiftY      = (sensorMaxY + sensorMinY) / 2;

    const uint laddersPerLayer[2] = {8, 12};
 
  };

}