← Back to branch summary

~christopher-hunt08/maus/maus_radiation_det

« back to all changes in this revision

Viewing changes to src/common_cpp/Optics/PolynomialTransferMap.cc

  • Committer: Christopher Hunt
  • Date: 2014-01-28 11:09:46 UTC
  • mfrom: (1019.5.38 tracker_devel)
  • Revision ID: christopher.hunt08@imperial.ac.uk-20140128110946-vy5z1vax8gshxtxm
Merged Again

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/common_cpp/Optics/PolynomialTransferMap.cc
24
24
#include "Maths/PolynomialMap.hh"
25
25
#include "src/common_cpp/Optics/CovarianceMatrix.hh"
26
26
#include "src/common_cpp/Optics/PhaseSpaceVector.hh"
 
27
#include "Recon/Global/ParticleOpticalVector.hh"
27
28
 
28
29
namespace MAUS {
 
30
 
 
31
using MAUS::recon::global::ParticleOpticalVector;
 
32
 
29
33
// ##############################
30
34
//  PolynomialTransferMap public
31
35
// ##############################
101
105
            << "Input Vector = " << std::endl << input_vector << std::endl;
102
106
  // operate on the input vector with the polynomial map
103
107
  MAUS::PhaseSpaceVector output_vector;
 
108
  #if 1
 
109
  polynomial_map_->F(input_vector-(*reference_trajectory_in_), output_vector);
 
110
  output_vector += (*reference_trajectory_out_);
 
111
  #else
104
112
  polynomial_map_->F(input_vector, output_vector);
 
113
  #endif
105
114
  std::cout << "DEBUG PolynomialTransferMap::Transport() "
106
115
            << "Transported Vector = " << std::endl << output_vector
107
116
            << std::endl;