GaussianStateConversions Namespace Reference

Functions
MultiGaussianState< 5 >	multiGaussianStateFromTSOS (const TrajectoryStateOnSurface &tsos)
MultiGaussianState< 3 >	multiGaussianStateFromVertex (const VertexState aState)
TrajectoryStateOnSurface	tsosFromMultiGaussianState (const MultiGaussianState< 5 > &multiState, const TrajectoryStateOnSurface &refTsos)
VertexState	vertexFromMultiGaussianState (const MultiGaussianState< 3 > &multiState)

Function Documentation

multiGaussianStateFromTSOS()

MultiGaussianState< 5 > GaussianStateConversions::multiGaussianStateFromTSOS

(

const TrajectoryStateOnSurface &

tsos

)

Definition at line 10 of file TsosGaussianStateConversions.cc.

References makeMuonMisalignmentScenario::components, TrajectoryStateOnSurface::components(), and TrajectoryStateOnSurface::isValid().

Referenced by MultiTrajectoryStateMerger::merge().

                                                                                         {
    if (!tsos.isValid())
      return MultiGaussianState<5>();

    using SingleStatePtr = std::shared_ptr<SingleGaussianState<5>>;
    auto const& components = tsos.components();
    MultiGaussianState<5>::SingleStateContainer singleStates;
    singleStates.reserve(components.size());
    for (auto const& ic : components) {
      if (ic.isValid()) {
        auto sgs = std::make_shared<SingleGaussianState<5>>(
            ic.localParameters().vector(), ic.localError().matrix(), ic.weight());
        singleStates.push_back(sgs);
      }
    }
    return MultiGaussianState<5>(singleStates);
  }

multiGaussianStateFromVertex()

MultiGaussianState< 3 > GaussianStateConversions::multiGaussianStateFromVertex

(

const VertexState

aState

)

Definition at line 8 of file VertexGaussianStateConversions.cc.

References makeMuonMisalignmentScenario::components, and VertexState::components().

Referenced by GsfVertexMerger::merge().

                                                                               {
    typedef std::shared_ptr<SingleGaussianState<3> > SingleStatePtr;
    const std::vector<VertexState> components = aState.components();
    MultiGaussianState<3>::SingleStateContainer singleStates;
    singleStates.reserve(components.size());
    for (std::vector<VertexState>::const_iterator ic = components.begin(); ic != components.end(); ic++) {
      if (ic->isValid()) {
        GlobalPoint pos(ic->position());
        AlgebraicVector3 parameters;
        parameters(0) = pos.x();
        parameters(1) = pos.y();
        parameters(2) = pos.z();
        SingleStatePtr sgs(new SingleGaussianState<3>(parameters, ic->error().matrix(), ic->weightInMixture()));
        singleStates.push_back(sgs);
      }
    }
    return MultiGaussianState<3>(singleStates);
  }

tsosFromMultiGaussianState()

TrajectoryStateOnSurface GaussianStateConversions::tsosFromMultiGaussianState	(	const MultiGaussianState< 5 > &	multiState,
		const TrajectoryStateOnSurface &	refTsos
	)

Definition at line 28 of file TsosGaussianStateConversions.cc.

References PostProcessorHGCAL_cfi::charged, makeMuonMisalignmentScenario::components, TrajectoryStateOnSurface::components(), TrajectoryStateOnSurface::magneticField(), MultiGaussianStateTransform::multiState(), ALPAKA_ACCELERATOR_NAMESPACE::ecal::reconstruction::internal::barrel::side(), TrajectoryStateOnSurface::surface(), and TrajectoryStateOnSurface::surfaceSide().

Referenced by MultiTrajectoryStateMerger::merge().

                                                                                               {
    if (multiState.components().empty())
      return TrajectoryStateOnSurface();
    const Surface& surface = refTsos.surface();
    SurfaceSide side = refTsos.surfaceSide();
    const MagneticField* field = refTsos.magneticField();
    auto const& refTsos1 = refTsos.components().front();
    auto pzSign = refTsos1.localParameters().pzSign();
    bool charged = refTsos1.charge() != 0;

    auto const& singleStates = multiState.components();
    std::vector<TrajectoryStateOnSurface> components;
    components.reserve(singleStates.size());
    for (auto const& ic : singleStates) {
      //require states to be positive-definite
      if (double det = 0; (*ic).covariance().Det2(det) && det > 0) {
        components.emplace_back((*ic).weight(),
                                LocalTrajectoryParameters((*ic).mean(), pzSign, charged),
                                LocalTrajectoryError((*ic).covariance()),
                                surface,
                                field,
                                side);
      }
    }
    return components.empty()
               ? TrajectoryStateOnSurface()
               : TrajectoryStateOnSurface((BasicTrajectoryState*)new BasicMultiTrajectoryState(components));
  }

vertexFromMultiGaussianState()

VertexState GaussianStateConversions::vertexFromMultiGaussianState

(

const MultiGaussianState< 3 > &

multiState

)

Definition at line 27 of file VertexGaussianStateConversions.cc.

References makeMuonMisalignmentScenario::components, relativeConstraints::error, MultiGaussianStateTransform::multiState(), and position.

Referenced by GsfVertexMerger::merge().

                                                                                    {
    if (multiState.components().empty())
      return VertexState();

    const MultiGaussianState<3>::SingleStateContainer& singleStates = multiState.components();
    std::vector<VertexState> components;
    components.reserve(singleStates.size());
    for (MultiGaussianState<3>::SingleStateContainer::const_iterator ic = singleStates.begin();
         ic != singleStates.end();
         ic++) {
      const AlgebraicVector3& par = (**ic).mean();
      GlobalPoint position(par(0), par(1), par(2));
      const AlgebraicSymMatrix33& cov = (**ic).covariance();
      GlobalError error(cov(0, 0), cov(1, 0), cov(2, 0), cov(1, 1), cov(2, 1), cov(2, 2));
      components.push_back(VertexState(position, error, (**ic).weight()));
    }
    return VertexState(new BasicMultiVertexState(components));
  }