Page MenuHomeHEPForge
Files F11221952

BasicMultiTrajectoryState.cc
No OneTemporary
Actions

BasicMultiTrajectoryState.cc
View Options

#include "TrackingTools/GsfTools/interface/BasicMultiTrajectoryState.h"
#include "DataFormats/GeometrySurface/interface/BoundPlane.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
using namespace SurfaceSideDefinition;
BasicMultiTrajectoryState::BasicMultiTrajectoryState( const std::vector<TSOS>& tsvec) :
theCombinedStateUp2Date( false)
{
theStates.reserve(tsvec.size());
for (std::vector<TSOS>::const_iterator i=tsvec.begin(); i!=tsvec.end(); i++) {
if (!i->isValid()) {
throw cms::Exception("LogicError") << "MultiTrajectoryState constructed with invalid state";
}
if (i->hasError() != tsvec.front().hasError()) {
throw cms::Exception("LogicError") << "MultiTrajectoryState mixes states with and without errors";
}
if ( &i->surface() != &tsvec.front().surface()) {
throw cms::Exception("LogicError") << "MultiTrajectoryState mixes states with different surfaces";
}
if ( i->surfaceSide() != tsvec.front().surfaceSide()) {
throw cms::Exception("LogicError")
<< "MultiTrajectoryState mixes states defined before and after material";
}
if ( i->localParameters().pzSign()*tsvec.front().localParameters().pzSign()<0. ) {
throw cms::Exception("LogicError")
<< "MultiTrajectoryState mixes states with different signs of local p_z";
}
if ( i==tsvec.begin() ) {
// only accept planes!!
const BoundPlane* bp = dynamic_cast<const BoundPlane*>(&i->surface());
if ( bp==0 )
throw cms::Exception("LogicError") << "MultiTrajectoryState constructed on cylinder";
}
theStates.push_back( *i);
}
}
void BasicMultiTrajectoryState::checkCombinedState() const
{
if (theCombinedStateUp2Date) return;
theCombinedState = theCombiner.combine( theStates);
theCombinedStateUp2Date = true;
}
double BasicMultiTrajectoryState::weight() const {
if (theStates.empty()) {
edm::LogError("BasicMultiTrajectoryState")
<< "Asking for weight of empty MultiTrajectoryState, returning zero!";
return 0.;
}
double sumw = 0.;
for (std::vector<TSOS>::const_iterator it = theStates.begin(); it != theStates.end(); it++) {
sumw += it->weight();
}
return sumw;
}
void BasicMultiTrajectoryState::rescaleError(double factor) {
if (theStates.empty()) {
edm::LogError("BasicMultiTrajectoryState") << "Trying to rescale errors of empty MultiTrajectoryState!";
return;
}
for (std::vector<TSOS>::iterator it = theStates.begin(); it != theStates.end(); it++) {
it->rescaleError(factor);
}
theCombinedStateUp2Date = false;
}
const MagneticField*
BasicMultiTrajectoryState::magneticField () const
{
//
// Magnetic field should be identical in all components:
// avoid forcing the combination of states and take value from 1st component!
//
if (theStates.empty()) {
edm::LogError("BasicMultiTrajectoryState")
<< "Asking for magneticField of empty MultiTrajectoryState, returning null pointer!";
return 0;
}
return theStates.front().magneticField();
}
SurfaceSide
BasicMultiTrajectoryState::surfaceSide () const
{
//
// SurfaceSide should be identical in all components:
// avoid forcing the combination of states and take value from 1st component!
//
if (theStates.empty()) {
edm::LogError("BasicMultiTrajectoryState")
<< "Asking for magneticField of empty MultiTrajectoryState, returning atCenterOfSurface!";
return atCenterOfSurface;
}
return theStates.front().surfaceSide();
}
void
BasicMultiTrajectoryState::
update( const LocalTrajectoryParameters& p,
const Surface& aSurface,
const MagneticField* field,
const SurfaceSide side)
{
throw cms::Exception("LogicError",
"BasicMultiTrajectoryState::update(LocalTrajectoryParameters, Surface, ...) called even if canUpdateLocalParameters() is false");
}
void
BasicMultiTrajectoryState::
update( const LocalTrajectoryParameters& p,
const LocalTrajectoryError& err,
const Surface& aSurface,
const MagneticField* field,
const SurfaceSide side,
double weight)
{
throw cms::Exception("LogicError",
"BasicMultiTrajectoryState::update(LocalTrajectoryParameters, LocalTrajectoryError, ...) called even if canUpdateLocalParameters() is false");
}

File Metadata

Mime Type
text/x-c
Expires
Wed, May 14, 11:09 AM (20 h, 15 m)
Storage Engine
blob
Storage Format
Raw Data
Storage Handle
5111361
Default Alt Text
BasicMultiTrajectoryState.cc (4 KB)

Event Timeline