summaryrefslogtreecommitdiff
path: root/src/main/Quadrupole.cc
diff options
context:
space:
mode:
authorJakob Odersky <jodersky@gmail.com>2011-04-24 13:31:11 +0000
committerJakob Odersky <jodersky@gmail.com>2011-04-24 13:31:11 +0000
commita91ecfe471c4659e65be573e0217c8b72d5fa958 (patch)
treeddb16859fdd206b5978902f2d965bfc7bc900a94 /src/main/Quadrupole.cc
parenta80f5cb7b3d8e4cf3e48c55895b3ea23a940b8ec (diff)
downloadvhc-a91ecfe471c4659e65be573e0217c8b72d5fa958.tar.gz
vhc-a91ecfe471c4659e65be573e0217c8b72d5fa958.tar.bz2
vhc-a91ecfe471c4659e65be573e0217c8b72d5fa958.zip
*Ajout d'un pointeur sur les element precedents.
*Implementation de la methode 'close()' dans l'accelerateur. *Restructuration des methodes geometriques des elements. -Ajout d'une methode isBefore(): determine si un point est contenu dans l'espace avant le plan de la face d'entree. -Ajout d'une methode isBeside(): determine si un point est plus proche de la trajectoire ideal que le rayon de section. -Ajout d'une methode isAfter(): determine si un point est contenu dans l'espace apres le plan de la face de sortie. -Ajout d'une methode contains(): determine si un point est contenu dans l'element (contains = !(isBefore || isBeside || isAfter)) -La methode isPast() (passe_au_suivant selon projet) a ete supprimee, elle est a remplacer par isAfter(). -La methode hitWall() (heurte_bord selon projet) a ete supprimee, elle est a remplacer par isBeside(). *Resolution d'un bug dans la methode isBeside() de StraightElement, avant une particule ne pouvait jamais heuter le bord! *Ajout des cibles gui et gui-build dans eclipse.
Diffstat (limited to 'src/main/Quadrupole.cc')
-rw-r--r--src/main/Quadrupole.cc6
1 files changed, 2 insertions, 4 deletions
diff --git a/src/main/Quadrupole.cc b/src/main/Quadrupole.cc
index 911ae18..ff0e952 100644
--- a/src/main/Quadrupole.cc
+++ b/src/main/Quadrupole.cc
@@ -19,14 +19,12 @@ Quadrupole::Quadrupole(const Vector3D& entry, const Vector3D& exit, double secti
Quadrupole::~Quadrupole() {};
Vector3D Quadrupole::magneticFieldAt(const Vector3D& position) const {
+ if (!Quadrupole::contains(position)) return Vector3D::Null;
Vector3D x = position - getEntryPosition();
Vector3D d = getDiagonal().unit();
Vector3D y = x - x.dot(d) * d;
Vector3D u = Vector3D::k.cross(d);
- if ((x.dot(d) >= 0) && ((position - getExitPosition()).dot(-d) >= 0))
- return focalizingCoefficient * (y.dot(u) * Vector3D::k + position.getZ() * u);
- else
- return Vector3D::Null;
+ return focalizingCoefficient * (y.dot(u) * Vector3D::k + position.getZ() * u);
}
double Quadrupole::getFocalizingCoefficient() const {