ATLASCuts.hpp

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file ATLASCuts.hpp

// This file is part of the Acts project.
//
// Copyright (C) 2023 CERN for the benefit of the Acts project
//
// This Source Code Form is subject to the terms of the Mozilla Public
// License, v. 2.0. If a copy of the MPL was not distributed with this
// file, You can obtain one at http://mozilla.org/MPL/2.0/.

#pragma once

#include "Acts/Seeding/IExperimentCuts.hpp"

#include <algorithm>

namespace Acts {
template <typename SpacePoint>
class ATLASCuts : public IExperimentCuts<SpacePoint> {
 public:
  float seedWeight(const InternalSpacePoint<SpacePoint>& bottom,
                   const InternalSpacePoint<SpacePoint>& middle,
                   const InternalSpacePoint<SpacePoint>& top) const override;
  bool singleSeedCut(
      float weight, const InternalSpacePoint<SpacePoint>& bottom,
      const InternalSpacePoint<SpacePoint>& /*middle*/,
      const InternalSpacePoint<SpacePoint>& /*top*/) const override;

  std::vector<typename CandidatesForMiddleSp<
      const InternalSpacePoint<SpacePoint>>::value_type>
  cutPerMiddleSP(std::vector<typename CandidatesForMiddleSp<
                     const InternalSpacePoint<SpacePoint>>::value_type>
                     seedCandidates) const override;
};

template <typename SpacePoint>
float ATLASCuts<SpacePoint>::seedWeight(
    const InternalSpacePoint<SpacePoint>& bottom,
    const InternalSpacePoint<SpacePoint>& /*middle*/,
    const InternalSpacePoint<SpacePoint>& top) const {
  float weight = 0;
  if (bottom.radius() > 150) {
    weight = 400;
  }
  if (top.radius() < 150) {
    weight = 200;
  }
  return weight;
}

template <typename SpacePoint>
bool ATLASCuts<SpacePoint>::singleSeedCut(
    float weight, const InternalSpacePoint<SpacePoint>& b,
    const InternalSpacePoint<SpacePoint>& /*m*/,
    const InternalSpacePoint<SpacePoint>& /*t*/) const {
  return !(b.radius() > 150. && weight < 380.);
}

template <typename SpacePoint>
std::vector<typename CandidatesForMiddleSp<
    const InternalSpacePoint<SpacePoint>>::value_type>
ATLASCuts<SpacePoint>::cutPerMiddleSP(
    std::vector<typename CandidatesForMiddleSp<
        const InternalSpacePoint<SpacePoint>>::value_type>
        seedCandidates) const {
  std::vector<typename CandidatesForMiddleSp<
      const InternalSpacePoint<SpacePoint>>::value_type>
      newSeedsVector;
  if (seedCandidates.size() <= 1) {
    return seedCandidates;
  }

  newSeedsVector.push_back(std::move(seedCandidates[0]));
  std::size_t itLength = std::min(seedCandidates.size(), std::size_t(5));
  // don't cut first element
  for (std::size_t i(1); i < itLength; i++) {
    float weight = seedCandidates[i].weight;
    const auto& bottom = seedCandidates[i].bottom;
    if (weight > 200. or bottom->radius() > 43.) {
      newSeedsVector.push_back(std::move(seedCandidates[i]));
    }
  }
  return newSeedsVector;
}
}  // namespace Acts