development/doxygen/CDCSimpleSimulation_8cc_source.html

/**************************************************************************

 * 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.                  *

 **************************************************************************/


#include <tracking/trackFindingCDC/sim/CDCSimpleSimulation.h>


#include <tracking/trackingUtilities/eventdata/tracks/CDCTrack.h>


#include <tracking/trackingUtilities/eventdata/hits/CDCRecoHit3D.h>

#include <tracking/trackingUtilities/eventdata/hits/CDCRLWireHit.h>

#include <tracking/trackingUtilities/eventdata/hits/CDCWireHit.h>


#include <cdc/topology/CDCWire.h>

#include <cdc/topology/CDCWireLayer.h>

#include <cdc/topology/CDCWireTopology.h>


#include <tracking/trackingUtilities/utilities/VectorRange.h>


#include <framework/gearbox/Const.h>

#include <framework/logging/Logger.h>


#include <Math/Vector3D.h>

#include <Math/Vector2D.h>

#include <TRandom.h>


#include <algorithm>


using namespace Belle2;

using namespace CDC;

using namespace TrackFindingCDC;

using namespace TrackingUtilities;


ConstVectorRange<CDCWireHit> CDCSimpleSimulation::getWireHits() const

{

  if (m_sharedWireHits) {

    return {m_sharedWireHits->begin(), m_sharedWireHits->end()};

  } else {

    return ConstVectorRange<CDCWireHit>();

  }

}


CDCTrack CDCSimpleSimulation::simulate(const CDCTrajectory3D& trajectory3D)

{

  return std::move(simulate(std::vector<CDCTrajectory3D>(1, trajectory3D)).front());

}


std::vector<CDCTrack> CDCSimpleSimulation::simulate(const std::vector<CDCTrajectory3D>& trajectories3D)

{

  std::vector<SimpleSimHit> simpleSimHits;

  const size_t nMCTracks = trajectories3D.size();


  for (size_t iMCTrack = 0;  iMCTrack < nMCTracks; ++iMCTrack) {

    const CDCTrajectory3D& trajectory3D = trajectories3D[iMCTrack];


    const UncertainHelix& localHelix = trajectory3D.getLocalHelix();

    const ROOT::Math::XYZVector& localOrigin = trajectory3D.getLocalOrigin();


    Helix globalHelix = localHelix;

    const double arcLength2DOffset = globalHelix.passiveMoveBy(-localOrigin);

    std::vector<SimpleSimHit> simpleSimHitsForTrajectory = createHits(globalHelix, arcLength2DOffset);


    for (SimpleSimHit& simpleSimHit : simpleSimHitsForTrajectory) {

      simpleSimHit.m_iMCTrack = iMCTrack;

      simpleSimHits.push_back(simpleSimHit);

    }

  }


  std::vector<CDCTrack> mcTracks = constructMCTracks(nMCTracks, std::move(simpleSimHits));


  for (size_t iMCTrack = 0; iMCTrack < nMCTracks; ++iMCTrack) {

    CDCTrack& mcTrack = mcTracks[iMCTrack];

    CDCTrajectory3D mcTrajectory = trajectories3D[iMCTrack];

    if (not mcTrack.empty()) {

      mcTrajectory.setLocalOrigin(mcTrack.front().getRecoPos3D());

      mcTrack.setStartTrajectory3D(mcTrajectory);

      mcTrajectory.setLocalOrigin(mcTrack.back().getRecoPos3D());

      mcTrack.setEndTrajectory3D(mcTrajectory);

    } else {

      mcTrack.setStartTrajectory3D(mcTrajectory);

      mcTrack.setEndTrajectory3D(mcTrajectory);

    }

  }

  return mcTracks;

}


std::vector<CDCTrack>


CDCSimpleSimulation::constructMCTracks(int nMCTracks, std::vector<SimpleSimHit> simpleSimHits)

{


  // Sort the hits along side their wire hits

  std::stable_sort(simpleSimHits.begin(), simpleSimHits.end(),

  [](const SimpleSimHit & lhs, const SimpleSimHit & rhs) -> bool {

    return lhs.m_wireHit < rhs.m_wireHit;

  });


  // Discard multiple hits on the same wire up to the maximal exceeding the maximal desired number

  if (m_maxNHitOnWire > 0) {

    const CDCWire* lastWire = nullptr;

    size_t nSameWire = 0;

    const size_t maxNHitOnWire = m_maxNHitOnWire;


    auto exceedsMaxNHitOnWire =

    [&lastWire, &nSameWire, maxNHitOnWire](const SimpleSimHit & simpleSimHit) -> bool {


      if (&(simpleSimHit.m_wireHit.getWire()) == lastWire)

      {

        ++nSameWire;

      } else

      {

        nSameWire = 1;

        lastWire = &(simpleSimHit.m_wireHit.getWire());

      }

      return nSameWire > maxNHitOnWire ? true : false;

    };


    auto itLast = std::remove_if(simpleSimHits.begin(), simpleSimHits.end(), exceedsMaxNHitOnWire);

    simpleSimHits.erase(itLast, simpleSimHits.end());

  }


  // Write the created hits and move them to the their storage place.

  {

    std::vector<CDCWireHit> wireHits;

    wireHits.reserve(simpleSimHits.size());

    for (SimpleSimHit& simpleSimHit : simpleSimHits) {

      wireHits.push_back(simpleSimHit.m_wireHit);

    }


    B2ASSERT("WireHits should be sorted as a result from sorting the SimpleSimHits. "

             "Algorithms may relay on the sorting o the WireHits",

             std::is_sorted(wireHits.begin(), wireHits.end()));


    m_sharedWireHits.reset(new const std::vector<CDCWireHit>(std::move(wireHits)));

  }


  // TODO: Decide if the EventMeta should be incremented after write.


  // Now construct the tracks.

  std::vector<CDCTrack> mcTracks;

  mcTracks.resize(nMCTracks);

  ConstVectorRange<CDCWireHit> wireHits = getWireHits();

  const size_t nWireHits = wireHits.size();


  for (size_t iWireHit = 0; iWireHit < nWireHits; ++iWireHit) {

    const CDCWireHit& wireHit = wireHits[iWireHit];

    const SimpleSimHit& simpleSimHit = simpleSimHits[iWireHit];


    CDCTrack& mcTrack = mcTracks[simpleSimHit.m_iMCTrack];


    CDCRLWireHit rlWireHit(&wireHit, simpleSimHit.m_rlInfo);

    CDCRecoHit3D recoHit3D(rlWireHit, simpleSimHit.m_pos3D, simpleSimHit.m_arcLength2D);

    mcTrack.push_back(recoHit3D);

  }


  for (CDCTrack& mcTrack : mcTracks) {

    mcTrack.sortByArcLength2D();

  }


  return mcTracks;

}


std::vector<CDCSimpleSimulation::SimpleSimHit>


CDCSimpleSimulation::createHits(const Helix& globalHelix,

                                double arcLength2DOffset) const

{


  std::vector<SimpleSimHit> simpleSimHits;


  CDCWireTopology& wireTopology = CDCWireTopology::getInstance();

  const double outerWallCylinderR = wireTopology.getOuterCylindricalR();


  const double minR = globalHelix.minimalCylindricalR();

  const double maxR = globalHelix.maximalCylindricalR();


  const double globalArcLength2DToOuterWall = globalHelix.arcLength2DToCylindricalR(outerWallCylinderR);

  const double localArcLength2DToOuterWall = arcLength2DOffset + globalArcLength2DToOuterWall;


  if (localArcLength2DToOuterWall < 0) {

    // Trajectory starts outside the CDC and initially flys away from it

    // Do not try to createHit hits for it

    B2WARNING("Simple simulation got trajectory outside CDC that moves away from the detector.");

    return simpleSimHits;

  }


  // Two dimensional arc length where the trajectory

  // * leaves the outer wall of the CDC or

  // * made a full circle (cut off for curlers)

  const bool isCurler = std::isnan(localArcLength2DToOuterWall);

  const double perimeterXY = globalHelix.perimeterXY();

  const double maxArcLength2D = isCurler ? fabs(perimeterXY) : localArcLength2DToOuterWall;


  if (isCurler) {

    B2INFO("Simulating curler");

  }


  for (const CDCWireLayer& wireLayer : wireTopology.getWireLayers()) {

    double outerR = wireLayer.getOuterCylindricalR();

    double innerR = wireLayer.getInnerCylindricalR();


    if ((maxR < innerR) or (outerR < minR)) {

      // Trajectory does not reaching the layer

      continue;

    }


    double centerR = (std::min(outerR, maxR) + std::max(innerR, minR)) / 2;


    double globalArcLength2D = globalHelix.arcLength2DToCylindricalR(centerR);

    double localArcLength2D = arcLength2DOffset + globalArcLength2D;


    std::vector<SimpleSimHit> simpleSimHitsInLayer;

    if (localArcLength2D > 0 and localArcLength2D < maxArcLength2D) {


      ROOT::Math::XYZVector pos3DAtLayer = globalHelix.atArcLength2D(globalArcLength2D);

      const CDCWire& closestWire =  wireLayer.getClosestWire(pos3DAtLayer);


      simpleSimHitsInLayer = createHitsForLayer(closestWire, globalHelix, arcLength2DOffset);


      for (SimpleSimHit& simpleSimHit : simpleSimHitsInLayer) {

        if (simpleSimHit.m_arcLength2D < maxArcLength2D) {

          simpleSimHits.push_back(simpleSimHit);

        }

      }

    } else {

      B2INFO("Arc length to long");

    }


    bool oneSegment = outerR > maxR or innerR < minR;

    if (not oneSegment) {


      // Check the second branch for more hits

      double secondGlobalArcLength2D = -globalArcLength2D;

      double secondArcLength2DOffset = arcLength2DOffset;

      double secondLocalArcLength2D = secondArcLength2DOffset + secondGlobalArcLength2D;


      if (isCurler and secondLocalArcLength2D < 0) {

        secondLocalArcLength2D += perimeterXY;

        secondArcLength2DOffset += perimeterXY;

        secondGlobalArcLength2D += perimeterXY;

      }


      if (secondLocalArcLength2D > 0 and secondLocalArcLength2D < maxArcLength2D) {

        ROOT::Math::XYZVector pos3DAtLayer = globalHelix.atArcLength2D(secondGlobalArcLength2D);

        const CDCWire& closestWire =  wireLayer.getClosestWire(pos3DAtLayer);


        // Check again if the wire has been hit before

        bool wireAlreadyHit = false;

        for (const SimpleSimHit& simpleSimHit : simpleSimHits) {

          if (simpleSimHit.m_wireHit.isOnWire(closestWire)) {

            wireAlreadyHit = true;

          }

        }

        if (not wireAlreadyHit) {

          std::vector<SimpleSimHit> secondSimpleSimHitsInLayer =

            createHitsForLayer(closestWire, globalHelix, secondArcLength2DOffset);


          for (SimpleSimHit& simpleSimHit : secondSimpleSimHitsInLayer) {

            if (simpleSimHit.m_arcLength2D < maxArcLength2D) {

              simpleSimHits.push_back(simpleSimHit);

            }

          }

        }

      }

    }

  }


  return simpleSimHits;

}


std::vector<CDCSimpleSimulation::SimpleSimHit>


CDCSimpleSimulation::createHitsForLayer(const CDCWire& nearWire,

                                        const Helix& globalHelix,

                                        double arcLength2DOffset) const

{

  std::vector<SimpleSimHit> result;


  SimpleSimHit simpleSimHit = createHitForCell(nearWire, globalHelix, arcLength2DOffset);

  if (not std::isnan(simpleSimHit.m_wireHit.getRefDriftLength())) {

    result.push_back(simpleSimHit);

  }


  const CDCWire* ccwWire = nearWire.getNeighborCCW();

  while (true) {

    SimpleSimHit simpleSimHitForWire = createHitForCell(*ccwWire, globalHelix, arcLength2DOffset);

    if (std::isnan(simpleSimHitForWire.m_arcLength2D) or

        std::isnan(simpleSimHitForWire.m_trueDriftLength)) {

      break;

    }

    result.push_back(simpleSimHitForWire);

    ccwWire = ccwWire->getNeighborCCW();

  }


  const CDCWire* cwWire = nearWire.getNeighborCW();

  while (true) {

    SimpleSimHit simpleSimHitForWire = createHitForCell(*cwWire, globalHelix, arcLength2DOffset);

    if (std::isnan(simpleSimHitForWire.m_arcLength2D) or

        std::isnan(simpleSimHitForWire.m_trueDriftLength)) {

      break;

    }

    result.push_back(simpleSimHitForWire);

    cwWire = cwWire->getNeighborCW();

  }


  return result;

}


CDCSimpleSimulation::SimpleSimHit


CDCSimpleSimulation::createHitForCell(const CDCWire& wire,

                                      const Helix& globalHelix,

                                      double arcLength2DOffset) const

{

  double arcLength2D = globalHelix.arcLength2DToXY(wire.getRefPos2D());

  if ((arcLength2D + arcLength2DOffset) < 0) {

    arcLength2D += globalHelix.perimeterXY();

  }


  ROOT::Math::XYZVector pos3D = globalHelix.atArcLength2D(arcLength2D);


  ROOT::Math::XYZVector correctedPos3D = pos3D;

  ROOT::Math::XYVector correctedWirePos(wire.getWirePos2DAtZ(correctedPos3D.z()));

  double correctedArcLength2D = arcLength2D;


  for (int c_Iter = 0; c_Iter < 2; c_Iter++) {

    // Iterate the extrapolation to the stereo shifted position.

    correctedWirePos = wire.getWirePos2DAtZ(correctedPos3D.z());

    correctedArcLength2D = globalHelix.arcLength2DToXY(correctedWirePos);


    if ((correctedArcLength2D + arcLength2DOffset) < 0) {

      correctedArcLength2D += globalHelix.perimeterXY();

    }

    correctedPos3D = globalHelix.atArcLength2D(correctedArcLength2D);

  }


  const double trueDriftLength = wire.getDriftLength(correctedPos3D);

  const double smearedDriftLength = trueDriftLength + gRandom->Gaus(0, m_driftLengthSigma);


  double delayTime = getEventTime();

  if (m_addTOFDelay) {

    double arcLength3D = hypot2(1, globalHelix.tanLambda()) * (correctedArcLength2D + arcLength2DOffset);

    delayTime += arcLength3D / Const::speedOfLight;

  }


  if (m_addInWireSignalDelay) {

    double backwardZ = wire.getBackwardZ();

    // Position where wire has been hit

    ROOT::Math::XYZVector wirePos = wire.getClosest(correctedPos3D);

    double distanceToBack = (wirePos.z() - backwardZ) * hypot2(1, wire.getTanStereoAngle());


    delayTime += distanceToBack / m_propSpeed;

  }


  double measuredDriftLength = smearedDriftLength + delayTime * m_driftSpeed;


  ERightLeft rlInfo = globalHelix.circleXY().isRightOrLeft(correctedWirePos);


  // if (not std::isnan(trueDriftLength)){

  //   B2INFO("Delay time " << delayTime);

  //   B2INFO("True dirft length " << trueDriftLength);

  //   B2INFO("Measured drift length " << measuredDriftLength);

  //   B2INFO("Absolute deviation " << measuredDriftLength - trueDriftLength);

  //   B2INFO("Relative deviation " << (measuredDriftLength / trueDriftLength - 1) * 100 << "%");

  // }


  return SimpleSimHit{

    CDCWireHit(wire.getWireID(), measuredDriftLength, m_driftLengthVariance),

    0,

    rlInfo,

    correctedPos3D,

    correctedArcLength2D,

    trueDriftLength

  };

}


std::vector<CDCTrack>


CDCSimpleSimulation::loadPreparedEvent()

{

  const size_t nMCTracks = 2;

  std::vector<SimpleSimHit> simpleSimHits;

  simpleSimHits.reserve(128 + 64); // First plus second mc track


  // First MC track

  size_t iMCTrack = 0;


  // SL 6

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 4, 251), 0.104), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 4, 250), 0.272), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 4, 249), 0.488), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 4, 248), 0.764), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 3, 247), 0.9), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 4, 247), 1.024), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 3, 246), 0.64), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 3, 245), 0.304), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 3, 244), 0.012), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 3, 243), 0.352), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 3, 242), 0.74), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 2, 241), 0.46), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 2, 240), 0.02), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 2, 239), 0.46), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 2, 238), 0.884), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 1, 238), 1.104), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 1, 237), 0.612), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 1, 236), 0.12), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 1, 235), 0.356), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 1, 234), 0.884), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 0, 235), 1.032), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 0, 234), 0.52), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 0, 233), 0.06), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 0, 232), 0.62), iMCTrack, ERightLeft::c_Left});


  // SL 5

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 5, 206), 1.116), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 5, 205), 0.464), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 5, 204), 0.168), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 4, 204), 1.08), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 4, 203), 0.392), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 4, 202), 0.304), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 3, 201), 0.968), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 3, 200), 0.252), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 3, 199), 0.476), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 2, 199), 0.736), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 2, 198), 0.008), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 2, 197), 0.788), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 1, 197), 1.188), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 1, 196), 0.404), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 1, 195), 0.356), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 0, 195), 0.74), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 0, 194), 0.04), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 0, 193), 0.832), iMCTrack, ERightLeft::c_Left});


  // SL 4

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 5, 173), 0.692), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 5, 172), 0.22), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 5, 171), 1.132), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 4, 172), 0.816), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 4, 171), 0.136), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 4, 170), 1.048), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 3, 170), 0.884), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 3, 169), 0.032), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 3, 168), 0.96), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 2, 169), 0.972), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 2, 168), 0.044), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 2, 167), 0.872), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 1, 167), 1.004), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 1, 166), 0.1), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 1, 165), 0.828), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 0, 166), 1.004), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 0, 165), 0.084), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 0, 164), 0.82), iMCTrack, ERightLeft::c_Left});


  // SL 3

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 5, 145), 0.508), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 5, 144), 0.5), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 4, 145), 1.348), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 4, 144), 0.292), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 4, 143), 0.68), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 3, 143), 1.136), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 3, 142), 0.12), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 3, 141), 0.872), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 2, 142), 0.96), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 2, 141), 0.036), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 1, 140), 0.756), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 1, 139), 0.204), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 0, 139), 0.588), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 0, 138), 0.332), iMCTrack, ERightLeft::c_Left});


  // SL 2

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 5, 116), 1.1), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 5, 115), 0.008), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 5, 114), 1.048), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 4, 115), 0.712), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 4, 114), 0.316), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 3, 113), 0.428), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 3, 112), 0.572), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 2, 112), 0.188), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 2, 111), 0.776), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 1, 111), 0.92), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 1, 110), 0.024), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 1, 109), 0.928), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 0, 110), 0.776), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 0, 109), 0.116), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 0, 108), 0.992), iMCTrack, ERightLeft::c_Left});


  // SL 1

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 5, 87), 0.664), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 5, 86), 0.3), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 4, 86), 0.504), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 4, 85), 0.424), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 3, 85), 1.256), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 3, 84), 0.388), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 3, 83), 0.5), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 2, 84), 1.128), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 2, 83), 0.28), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 2, 82), 0.532), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 1, 82), 1.084), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 1, 81), 0.3), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 1, 80), 0.472), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 0, 81), 1.124), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 0, 80), 0.428), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 0, 79), 0.296), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 0, 78), 0.972), iMCTrack, ERightLeft::c_Left});


  // SL 0

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 7, 81), 0.192), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 7, 80), 0.452), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 6, 80), 0.596), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 6, 79), 0.024), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 6, 78), 0.66), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 5, 79), 0.388), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 5, 78), 0.184), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 4, 77), 0.296), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 4, 76), 0.244), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 3, 76), 0.268), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 3, 75), 0.212), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 2, 74), 0.316), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 2, 73), 0.112), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 2, 72), 0.588), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 1, 73), 0.464), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 1, 72), 0.028), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 0, 70), 0.284), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 0, 69), 0.088), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 0, 68), 0.416), iMCTrack, ERightLeft::c_Left});


  // Second MC track

  iMCTrack = 1;


  // SL 0

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 1, 140), 0.308), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 2, 139), 0.08), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 3, 139), 0.16), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 4, 139), 0.404), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 4, 138), 0.38), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 5, 139), 0.132), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 6, 138), 0.108), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 7, 139), 0.48), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(0, 7, 138), 0.424), iMCTrack, ERightLeft::c_Right});


  // SL 1

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 0, 136), 0.532), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 0, 135), 0.452), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 1, 135), 0.396), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 2, 135), 0.26), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 1, 134), 0.64), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 3, 134), 0.092), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 4, 134), 0.16), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(1, 5, 133), 0.524), iMCTrack, ERightLeft::c_Right});


  // SL 2

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 0, 163), 0.228), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 1, 162), 0.356), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 2, 163), 0.776), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 2, 162), 0.46), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 3, 162), 0.744), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 3, 161), 0.58), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 4, 162), 0.656), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 4, 161), 0.68), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 5, 161), 0.568), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(2, 5, 160), 0.812), iMCTrack, ERightLeft::c_Right});


  // SL 3

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 0, 190), 0.54), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 1, 188), 0.688), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 2, 188), 0.656), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 3, 188), 0.664), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 3, 187), 0.68), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 4, 188), 0.724), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 4, 187), 0.656), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(3, 5, 186), 0.676), iMCTrack, ERightLeft::c_Right});


  // SL 4

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 0, 211), 0.42), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 0, 210), 0.872), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 1, 210), 0.548), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 1, 209), 0.796), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 2, 210), 0.716), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 2, 209), 0.656), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 3, 209), 0.856), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 4, 209), 1.056), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 4, 208), 0.36), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(4, 5, 207), 0.232), iMCTrack, ERightLeft::c_Right});


  // SL 5

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 0, 231), 0.224), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 0, 230), 1.088), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 1, 230), 0.452), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 1, 229), 0.912), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 2, 230), 0.72), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 2, 229), 0.632), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 3, 229), 1.016), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 3, 228), 0.34), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 4, 228), 0.04), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 5, 227), 0.22), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(5, 5, 226), 1.196), iMCTrack, ERightLeft::c_Right});


  // SL 6

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 0, 254), 0.104), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 1, 253), 0.504), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 1, 252), 0.78), iMCTrack, ERightLeft::c_Right});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 2, 253), 0.968), iMCTrack, ERightLeft::c_Left});

  simpleSimHits.push_back(SimpleSimHit{CDCWireHit(WireID(6, 2, 252), 0.332), iMCTrack, ERightLeft::c_Right});


  std::vector<CDCTrack> mcTracks = constructMCTracks(nMCTracks, std::move(simpleSimHits));

  return mcTracks;

}


Belle2::CDC::CDCWireLayer
Class representing a sense wire layer in the central drift chamber.
Definition CDCWireLayer.h:42

Belle2::CDC::CDCWireTopology
Class representing the sense wire arrangement in the whole of the central drift chamber.
Definition CDCWireTopology.h:44

Belle2::CDC::CDCWireTopology::getWireLayers
const std::vector< CDCWireLayer > & getWireLayers() const
Getter for the underlying storing layer vector.
Definition CDCWireTopology.h:149

Belle2::CDC::CDCWireTopology::getInstance
static CDCWireTopology & getInstance()
Getter for the singleton instance of the wire topology.
Definition CDCWireTopology.cc:21

Belle2::CDC::CDCWireTopology::getOuterCylindricalR
double getOuterCylindricalR() const
Getter for the outer radius of the outer most wire layer.
Definition CDCWireTopology.h:471

Belle2::CDC::CDCWire
Class representing a sense wire in the central drift chamber.
Definition CDCWire.h:50

Belle2::CDC::CDCWire::getRefPos2D
const ROOT::Math::XYVector & getRefPos2D() const
Getter for the wire reference position for 2D tracking Gives the wire's reference position projected ...
Definition CDCWire.h:221

Belle2::CDC::CDCWire::getWireID
const WireID & getWireID() const
Getter for the wire id.
Definition CDCWire.h:114

Belle2::CDC::CDCWire::getNeighborCCW
TrackingUtilities::MayBePtr< const CDCWire > getNeighborCCW() const
Gives the closest neighbor in the counterclockwise direction - always exists.
Definition CDCWire.cc:159

Belle2::CDC::CDCWire::getWirePos2DAtZ
ROOT::Math::XYVector getWirePos2DAtZ(const double z) const
Gives the xy projected position of the wire at the given z coordinate.
Definition CDCWire.h:184

Belle2::CDC::CDCWire::getDriftLength
double getDriftLength(const ROOT::Math::XYZVector &pos3D) const
Calculates the straight drift length from the position to the wire This is essentially the same as th...
Definition CDCWire.h:206

Belle2::CDC::CDCWire::getNeighborCW
TrackingUtilities::MayBePtr< const CDCWire > getNeighborCW() const
Gives the closest neighbor in the clockwise direction - always exists.
Definition CDCWire.cc:164

Belle2::CDC::CDCWire::getTanStereoAngle
double getTanStereoAngle() const
Getter for the tangents of the stereo angle of the wire.
Definition CDCWire.h:232

Belle2::CDC::CDCWire::getBackwardZ
double getBackwardZ() const
Getter for the z coordinate at the backward joint points of the wires.
Definition CDCWire.h:272

Belle2::CDC::CDCWire::getClosest
ROOT::Math::XYZVector getClosest(const ROOT::Math::XYZVector &pos3D) const
Calculates the closest approach in the wire to the position.
Definition CDCWire.h:196

Belle2::CDC::Helix
Helix parameter class.
Definition Helix.h:48

Belle2::Const::speedOfLight
static const double speedOfLight
[cm/ns]
Definition Const.h:695

Belle2::TrackFindingCDC::CDCSimpleSimulation::getWireHits
TrackingUtilities::ConstVectorRange< TrackingUtilities::CDCWireHit > getWireHits() const
Getter for the wire hits created in the simulation.
Definition CDCSimpleSimulation.cc:37

Belle2::TrackFindingCDC::CDCSimpleSimulation::m_sharedWireHits
std::shared_ptr< const std::vector< TrackingUtilities::CDCWireHit > > m_sharedWireHits
Space for the memory of the generated wire hits.
Definition CDCSimpleSimulation.h:167

Belle2::TrackFindingCDC::CDCSimpleSimulation::constructMCTracks
std::vector< TrackingUtilities::CDCTrack > constructMCTracks(int nMCTracks, std::vector< SimpleSimHit > simpleSimHits)
Creates TrackingUtilities::CDCWireHits and uses them to construct the true TrackingUtilities::CDCTrac...
Definition CDCSimpleSimulation.cc:93

Belle2::TrackFindingCDC::CDCSimpleSimulation::m_addInWireSignalDelay
bool m_addInWireSignalDelay
Switch to activate the in wire signal delay.
Definition CDCSimpleSimulation.h:188

Belle2::TrackFindingCDC::CDCSimpleSimulation::m_propSpeed
double m_propSpeed
Electrical current propagation speed in the wires.
Definition CDCSimpleSimulation.h:198

Belle2::TrackFindingCDC::CDCSimpleSimulation::createHits
std::vector< SimpleSimHit > createHits(const TrackingUtilities::Helix &globalHelix, double arcLength2DOffset) const
Generate hits for the given helix in starting from the two dimensional arc length.
Definition CDCSimpleSimulation.cc:171

Belle2::TrackFindingCDC::CDCSimpleSimulation::simulate
std::vector< TrackingUtilities::CDCTrack > simulate(const std::vector< TrackingUtilities::CDCTrajectory3D > &trajectories3D)
Propagates the trajectories through the CDC as without energy loss until they first leave the CDC.
Definition CDCSimpleSimulation.cc:52

Belle2::TrackFindingCDC::CDCSimpleSimulation::getEventTime
double getEventTime() const
Getter for a global event time offset.
Definition CDCSimpleSimulation.h:142

Belle2::TrackFindingCDC::CDCSimpleSimulation::m_driftLengthVariance
double m_driftLengthVariance
Variance by which the drift length should be smeared.
Definition CDCSimpleSimulation.h:192

Belle2::TrackFindingCDC::CDCSimpleSimulation::m_addTOFDelay
bool m_addTOFDelay
Switch to activate the addition of the time of flight.
Definition CDCSimpleSimulation.h:185

Belle2::TrackFindingCDC::CDCSimpleSimulation::m_maxNHitOnWire
int m_maxNHitOnWire
Maximal number of hits allowed on each wire (0 means all).
Definition CDCSimpleSimulation.h:179

Belle2::TrackFindingCDC::CDCSimpleSimulation::m_driftLengthSigma
double m_driftLengthSigma
Standard deviation by which the drift length should be smeared.
Definition CDCSimpleSimulation.h:195

Belle2::TrackFindingCDC::CDCSimpleSimulation::loadPreparedEvent
std::vector< TrackingUtilities::CDCTrack > loadPreparedEvent()
Fills the wire hits with a hard coded event from the real simulation.
Definition CDCSimpleSimulation.cc:387

Belle2::TrackFindingCDC::CDCSimpleSimulation::createHitForCell
SimpleSimHit createHitForCell(const CDC::CDCWire &wire, const TrackingUtilities::Helix &globalHelix, double arcLength2DOffset) const
Generate a hit for the given wire.
Definition CDCSimpleSimulation.cc:319

Belle2::TrackFindingCDC::CDCSimpleSimulation::createHitsForLayer
std::vector< SimpleSimHit > createHitsForLayer(const CDC::CDCWire &nearWire, const TrackingUtilities::Helix &globalHelix, double arcLength2DOffset) const
Generate connected hits for wires in the same layer close to the given wire.
Definition CDCSimpleSimulation.cc:279

Belle2::TrackFindingCDC::CDCSimpleSimulation::m_driftSpeed
double m_driftSpeed
Electron drift speed in the cdc gas.
Definition CDCSimpleSimulation.h:201

Belle2::TrackingUtilities::CDCRLWireHit
Class representing an oriented hit wire including a hypotheses whether the causing track passes left ...
Definition CDCRLWireHit.h:46

Belle2::TrackingUtilities::CDCRecoHit3D
Class representing a three dimensional reconstructed hit.
Definition CDCRecoHit3D.h:59

Belle2::TrackingUtilities::CDCTrack
Class representing a sequence of three dimensional reconstructed hits.
Definition CDCTrack.h:37

Belle2::TrackingUtilities::CDCTrack::setStartTrajectory3D
void setStartTrajectory3D(const CDCTrajectory3D &startTrajectory3D)
Setter for the two dimensional trajectory.
Definition CDCTrack.h:94

Belle2::TrackingUtilities::CDCTrack::setEndTrajectory3D
void setEndTrajectory3D(const CDCTrajectory3D &endTrajectory3D)
Setter for the three dimensional trajectory.
Definition CDCTrack.h:101

Belle2::TrackingUtilities::CDCTrajectory3D
Particle full three dimensional trajectory.
Definition CDCTrajectory3D.h:46

Belle2::TrackingUtilities::CDCTrajectory3D::setLocalOrigin
double setLocalOrigin(const ROOT::Math::XYZVector &localOrigin)
Setter for the origin of the local coordinate system.
Definition CDCTrajectory3D.cc:368

Belle2::TrackingUtilities::CDCTrajectory3D::getLocalHelix
const UncertainHelix & getLocalHelix() const
Getter for the helix in local coordinates.
Definition CDCTrajectory3D.h:340

Belle2::TrackingUtilities::CDCTrajectory3D::getLocalOrigin
const ROOT::Math::XYZVector & getLocalOrigin() const
Getter for the origin of the local coordinate system.
Definition CDCTrajectory3D.h:352

Belle2::TrackingUtilities::CDCWireHit
Class representing a hit wire in the central drift chamber.
Definition CDCWireHit.h:56

Belle2::TrackingUtilities::CDCWireHit::getRefDriftLength
double getRefDriftLength() const
Getter for the drift length at the reference position of the wire.
Definition CDCWireHit.h:225

Belle2::TrackingUtilities::Range::size
std::size_t size() const
Returns the total number of objects in this range.
Definition Range.h:76

Belle2::UncertainHelix
This class represents an ideal helix in perigee parameterization including the covariance matrix of t...
Definition UncertainHelix.h:36

Belle2::WireID
Class to identify a wire inside the CDC.
Definition WireID.h:34

Belle2::CDC
Definition CrudeT0CalibrationAlgorithm.h:20

Belle2
Abstract base class for different kinds of events.
Definition MillepedeAlgorithm.h:17

Belle2::TrackFindingCDC::CDCSimpleSimulation::SimpleSimHit
Structure to accommodate information about the individual hits during the simulation.
Definition CDCSimpleSimulation.h:53

Belle2::TrackFindingCDC::CDCSimpleSimulation::SimpleSimHit::m_arcLength2D
double m_arcLength2D
Memory for the true two dimensional arc length on the helix to this hit.
Definition CDCSimpleSimulation.h:98

Belle2::TrackFindingCDC::CDCSimpleSimulation::SimpleSimHit::m_trueDriftLength
double m_trueDriftLength
Memory for the true drift length from the true position to the wire.
Definition CDCSimpleSimulation.h:101

Belle2::TrackFindingCDC::CDCSimpleSimulation::SimpleSimHit::m_rlInfo
TrackingUtilities::ERightLeft m_rlInfo
Memory for the true right left passage information.
Definition CDCSimpleSimulation.h:92

Belle2::TrackFindingCDC::CDCSimpleSimulation::SimpleSimHit::m_iMCTrack
size_t m_iMCTrack
Memory for the true index of the track this hit is contained in.
Definition CDCSimpleSimulation.h:89

Belle2::TrackFindingCDC::CDCSimpleSimulation::SimpleSimHit::m_pos3D
ROOT::Math::XYZVector m_pos3D
Memory for the true position on the track closest to the wire.
Definition CDCSimpleSimulation.h:95

Belle2::TrackFindingCDC::CDCSimpleSimulation::SimpleSimHit::m_wireHit
TrackingUtilities::CDCWireHit m_wireHit
Memory for the wire hit instance that will be given to the reconstruction.
Definition CDCSimpleSimulation.h:86