OPAL (Object Oriented Parallel Accelerator Library)  2.2.0
OPAL
CoordinateTransform.h
Go to the documentation of this file.
1 /*
2  * Copyright (c) 2018, Martin Duy Tat
3  * All rights reserved.
4  * Redistribution and use in source and binary forms, with or without
5  * modification, are permitted provided that the following conditions are met:
6  * 1. Redistributions of source code must retain the above copyright notice,
7  * this list of conditions and the following disclaimer.
8  * 2. Redistributions in binary form must reproduce the above copyright notice,
9  * this list of conditions and the following disclaimer in the documentation
10  * and/or other materials provided with the distribution.
11  * 3. Neither the name of STFC nor the names of its contributors may be used to
12  * endorse or promote products derived from this software without specific
13  * prior written permission.
14  *
15  * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
16  * AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17  * IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18  * ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
19  * LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
20  * CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
21  * SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
22  * INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
23  * CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
24  * ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
25  * POSSIBILITY OF SUCH DAMAGE.
26  */
27 
28 #ifndef COORDINATE_TRANSFORM_H
29 #define COORDINATE_TRANSFORM_H
30 
48 #include <vector>
49 
50 namespace coordinatetransform {
51 
52 struct myParams {
53  double s_0;
54  double lambdaleft;
55  double lambdaright;
56  double rho;
57 };
58 
60 public:
62  CoordinateTransform() = delete;
73  CoordinateTransform(const double &xlab,
74  const double &ylab,
75  const double &zlab,
76  const double &s_0,
77  const double &lambdaleft,
78  const double &lambdaright,
79  const double &rho);
81  CoordinateTransform(const CoordinateTransform &transform);
87  std::vector<double> getTransformation() const;
91  std::vector<double> calcReferenceTrajectory(const double &s) const;
95  std::vector<double> getUnitTangentVector(const double &s) const;
96 private:
101  void calcSCoordinate(const double &xlab, const double &ylab);
106  void calcXCoordinate(const double &xlab, const double &ylab);
114  void transformFromEntranceCoordinates(std::vector<double> &coordinates,
115  const double &boundingBoxLength);
116  double s_0_m;
117  double lambdaleft_m;
119  double rho_m;
120  double x_m;
121  double z_m;
122  double s_m;
123  static const double error;
124  static const int workspaceSize;
125  static const int algorithm;
126 };
127 
129 double getUnitTangentVectorX(double s, void *p);
131 double getUnitTangentVectorY(double s, void *p);
132 
133 }
134 #endif
double getUnitTangentVectorY(double s, void *p)
void calcSCoordinate(const double &xlab, const double &ylab)
std::vector< double > calcReferenceTrajectory(const double &s) const
CoordinateTransform & operator=(const CoordinateTransform &transform)
void transformFromEntranceCoordinates(std::vector< double > &coordinates, const double &boundingBoxLength)
std::vector< double > getTransformation() const
std::vector< double > getUnitTangentVector(const double &s) const
void calcXCoordinate(const double &xlab, const double &ylab)
double getUnitTangentVectorX(double s, void *p)