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 "common/Benchmark.h"
#include "dataobjects/FitGeoParams.h"
#include "dataobjects/Point.h"

// C++ headers
#include <cmath>
#include <vector>

// this in c++17 will become
// namespace Herd::FitDigitizationTools {
namespace Herd {
namespace 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;

  std::vector<std::pair<Point, double>> effectivePathFracInCores(const Point &enPoint, const Point &exPoint,
                                                                 const FitGeoParams &geoParams) const;

  void setMatWidth(float matWidth);

private:
  inline double getFiberX(unsigned int iLayer, int iFiber) const {
    return m_fibrePitch * (iFiber + 0.5 * (iLayer % 2));
  }
  inline double getFiberZ(unsigned int iLayer) const { return (iLayer - 0.5 * (m_nlayers - 1)) * m_deltaZ; }

  inline double firstFiber(unsigned int iLayer) const { return ((iLayer % 2) == 0) ? m_firstFiber0 : m_firstFiber1; }

  inline double lastFiber(unsigned int iLayer) const { return ((iLayer % 2) == 0) ? -m_firstFiber0 : m_lastFiber1; }

  void initNFibers(unsigned int nFibers);

  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 FitDigitizationTools
} // namespace Herd

#endif