← Back to branch summary

~mbogomilov/maus/devel3

« back to all changes in this revision

Viewing changes to src/legacy/Optics/TransferMapCalculator.cc

  • Committer: Durga Rajaram
  • Date: 2014-01-14 07:07:02 UTC
  • mfrom: (659.1.80 relcand)
  • Revision ID: durga@fnal.gov-20140114070702-2l1fuj1w6rraw7xe
Tags: MAUS-v0.7.6
MAUS-v0.7.6

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/legacy/Optics/TransferMapCalculator.cc
75
75
 
76
76
TransferMap TransferMapCalculator::GetPolynomialTransferMap(std::vector<PhaseSpaceVector> in, std::vector<PhaseSpaceVector> out, int order)
77
77
{
78
 
  if(in.size() < 1 || out.size() != in.size()) throw(Squeal(Squeal::recoverable, "Badly conditioned input", "TransferMapCalculator::GetPolynomialTransferMap"));
 
78
  if(in.size() < 1 || out.size() != in.size()) throw(MAUS::Exception(MAUS::Exception::recoverable, "Badly conditioned input", "TransferMapCalculator::GetPolynomialTransferMap"));
79
79
  std::vector< std::vector<double> > in_v (in.size(), std::vector<double>(6,0.));
80
80
  std::vector< std::vector<double> > out_v(in.size(), std::vector<double>(6,0.));
81
81
  for(size_t i=0; i<in.size(); i++) 
105
105
  MAUS::PolynomialMap* pvec = MAUS::PolynomialMap::Chi2SweepingLeastSquaresFitVariableWalls(
106
106
    *trackingOutput, order, std::vector< MAUS::PolynomialMap::PolynomialCoefficient >(), 
107
107
                               chi2Max, deltaV, deltaFactor, maxNumberOfSteps, deltaMaxV);
108
 
  if(!pvec) throw(Squeal(Squeal::recoverable, "Failed to make any polynomial fit at all - try tweaking control parameters", "TransferMapCalculator::GetSweepingPolynomialTransferMap"));
 
108
  if(!pvec) throw(MAUS::Exception(MAUS::Exception::recoverable, "Failed to make any polynomial fit at all - try tweaking control parameters", "TransferMapCalculator::GetSweepingPolynomialTransferMap"));
109
109
  for(size_t i=0; i<deltaV.size(); i++) delta[i] = deltaV[i];
110
110
  TransferMap       map;
111
111
  map.SetPolynomialMap(pvec);