DistributionFiberTool.h

Go to the documentation of this file.

/*
 * DistributionFiberTool.h
 *
 *  This tool calculates the fraction of the total path length
 *  of a track in the core of each fibre that was hit along with
 *  the position of the fibre.
 *
 *  Created on 03 Mar 2020
 *      Author: Valerio Formato (based on code by Junjing Wang)
 */
 
#ifndef HERD_DISTRIBUTIONFIBERTOOL_H_
#define HERD_DISTRIBUTIONFIBERTOOL_H_
 
// HerdSoftware headers
#include "herd/analysis/common/Benchmark.h"
#include "herd/analysis/dataobjects/FitGeoParams.h"
#include "herd/analysis/dataobjects/Point.h"
 
// C++ headers
#include <cmath>
#include <vector>
 
// this in c++17 will become
// namespace Herd::FitDigitizationTools {
namespace Herd::FitDigitizationTools {
 
namespace MCFT {
 
class Vec2 : public Vec3D {
public:
  Vec2() : Vec3D(std::numeric_limits<float>::max(), std::numeric_limits<float>::max(), 0) {}
  Vec2(float u, float v) : Vec3D(u, v, 0) {}
 
  Vec2(const Vec3D &vec) : Vec2(vec[RefFrame::Coo::X], vec[RefFrame::Coo::Y]) {}
 
  // float u() { return (*this)[RefFrame::Coo::X]; }
  // float v() { return (*this)[RefFrame::Coo::Y]; }
};
 
double distanceToSquared(const Vec2 &v, const Vec2 &p);
double distanceTo(const Vec2 &v, const Vec2 &p);
double minDistanceFromLineSegmentToPoint(const Vec2 &v, const Vec2 &w, const Vec2 &p);
double minDistanceFromLineToPoint(const Vec2 &v, const Vec2 &w, const Vec2 &p);
 
inline double getXOnPath(const Vec2 &entry, const Vec2 &exit, double z) {
  return entry[RefFrame::Coo::X] + (exit[RefFrame::Coo::X] - entry[RefFrame::Coo::X]) /
                                       (exit[RefFrame::Coo::Y] - entry[RefFrame::Coo::Y]) *
                                       (z - entry[RefFrame::Coo::Y]);
}
 
inline double getZOnPath(const Vec2 &entry, const Vec2 &exit, double x) {
  return entry[RefFrame::Coo::Y] + (exit[RefFrame::Coo::Y] - entry[RefFrame::Coo::Y]) /
                                       (exit[RefFrame::Coo::X] - entry[RefFrame::Coo::X]) *
                                       (x - entry[RefFrame::Coo::X]);
}
 
} // namespace MCFT
 
class DistributionFiberTool {
public:
  DistributionFiberTool();
  ~DistributionFiberTool() = default;
 
  [[nodiscard]] std::vector<std::pair<Point, double>>
  effectivePathFracInCores(const Point &enPoint, const Point &exPoint, const FitGeoParams &geoParams) const;
 
  void setMatWidth(float matWidth);
 
private:
  [[nodiscard]] inline double getFiberX(unsigned int iLayer, int iFiber) const {
    return m_fibrePitch * (iFiber + 0.5 * (iLayer % 2));
  }
  [[nodiscard]] inline double getFiberZ(unsigned int iLayer) const {
    return (iLayer - 0.5 * (m_nlayers - 1)) * m_deltaZ;
  }
 
  [[nodiscard]] inline double firstFiber(unsigned int iLayer) const {
    return ((iLayer % 2) == 0) ? m_firstFiber0 : m_firstFiber1;
  }
 
  [[nodiscard]] inline double lastFiber(unsigned int iLayer) const {
    return ((iLayer % 2) == 0) ? -m_firstFiber0 : m_lastFiber1;
  }
 
  void initNFibers(unsigned int nFibers);
 
  [[nodiscard]] std::pair<int, int> getAdjacentFiberIds(double x, double dx_over_dz, unsigned int iLayer) const;
 
  // const double m_fibreDiameter = 0.025;   // Diameter of the fiber in the fiber mats
  const double m_fibreCoreRadius = 0.011; // Radius of the fiber core
  const double m_fibrePitch = 0.0275;     // Distance between two neighbour fibre centers in x-direction
  const double m_deltaZ = 0.021;          // Pitch in z-direction
  // junjing 0809
  int m_nlayers = 6; // Number of layers
 
  //     ()()()
  //    ()()()()
  //       ()()()
  //      ()()()()
  //       ()()()
  //      ()()()()
 
  const double m_crossTalkProb = 0.193;   // Probability of cross talk between fibers
  const double m_CTmodelTouching = 0.163; // Fraction of cross talk going to the each of the 4 neighbors touching fibers
  const double m_CTmodelSide = 0.147;     // Fraction of cross talk going to each  of the 2 neighbors on the sides
  const double m_CTmodel2layers = 0.028;  // Fraction of cross talk going to each of the 2 neighbors on the same column
  //*/
 
  double m_minZ;
  unsigned int m_nFibers = 355;
  int m_firstFiber0;
  int m_firstFiber1;
  int m_lastFiber1;
};
} // namespace Herd::FitDigitizationTools
 
#endif