Fix Drift

src/Algorithms/ThickTracker.cpp

@@ -462,6 +462,10 @@ void ThickTracker::execute() {
                         - ( py*py )
                         - 1./( beta0 * beta0 * gamma0 * gamma0 ),order+1
                 );
+                
+                Drift* drift_p = dynamic_cast<Drift*> (element.get());
+                
+                msg << "drift NSLICES: " << drift_p->getNSlices() << endl;
 
                 break;
             }

src/Classic/AbsBeamline/Drift.cpp

@@ -28,17 +28,20 @@ extern Inform *gmsg;
 // ------------------------------------------------------------------------
 
 Drift::Drift():
-    Component()
+    Component(),
+    nSlices_m(1)
 { }
 
 
 Drift::Drift(const Drift &right):
-    Component(right)
+    Component(right),
+    nSlices_m(right.nSlices_m)
 { }
 
 
 Drift::Drift(const std::string &name):
-    Component(name) {
+    Component(name),
+    nSlices_m(1) {
 
 }
 

tests/Maps/map.in

@@ -14,12 +14,12 @@ REAL rf     = 50.6328e6;  //need to be confirmed
 
 VALUE,{gamma,brho,Edes,beta,gambet};
 
-QP1: QUADRUPOLE, L=1.0, K1= 2.00, ELEMEDGE=0.000;
+QP1: QUADRUPOLE, L=1.0, K1= 2.00, ELEMEDGE=0.000, NSLICES=3;
 QP2: QUADRUPOLE, L=1.0, K1=-3.00, ELEMEDGE=1.000;
-QP3: QUADRUPOLE, L=1.0, K1= 4.00, ELEMEDGE=2.000;
+QP3: QUADRUPOLE, L=1.0, K1= 4.00, ELEMEDGE=2.000, NSLICES=4;
 
 
-D: DRIFT, L=2.0, ELEMEDGE=3.000 ;
+D: DRIFT, L=2.0, ELEMEDGE=3.000, NSLICES=4;
 
 //BENDS: RBend, K0 = 1.000, FMAPFN = "1DPROFILE1-DEFAULT", ELEMEDGE = 5.0, DESIGNENERGY = Edes*1e3, L = 1.0, GAP = 0.02;
 
@@ -30,7 +30,7 @@ BENDS: SBend,   ANGLE = 30.0*Pi / 180.0, FMAPFN = "1DPROFILE1-DEFAULT", ELEMEDGE
 S11: SEXTUPOLE, L=0.4, K2=0.00134, ELEMEDGE=4.000;
 
 //QUADTEST: LINE=( QP1, QP2, D, QP3);
-QUADTEST: LINE=( QP1 );
+QUADTEST: LINE=( QP1, QP2, D, QP3 );
 
 D1: DISTRIBUTION, TYPE=GAUSS,
 SIGMAX=  5*1.0e-03,   SIGMAPX= 0.696*1.0e06,  CORRX= 0.0,