Changeset 94 in svn for trunk/src/BFieldProp.cc

Timestamp:

Dec 12, 2008, 5:32:29 PM (16 years ago)

Author:

severine ovyn

Message:

Add frog plus cleaning + bugs remove

File:

• trunk/src/BFieldProp.cc (modified) (3 diffs)

trunk/src/BFieldProp.cc

@@ -43 +43 @@
   
   //out of trackibg coverage?
-  if(sqrt(Xvertex1*Xvertex1+Yvertex1*Yvertex1) > TRACKING_RADIUS){return;}
-  if(fabs(Zvertex1) > TRACKING_LENGTH){return;}
+  if(sqrt(Xvertex1*Xvertex1+Yvertex1*Yvertex1) > TRACK_radius){return;}
+  if(fabs(Zvertex1) > TRACK_length){return;}
   
   float Px = Part->Px;
@@ -62 +62 @@
      double vz = pz/M;
 
-     double Bx = BFIELD_X;
-     double By = BFIELD_Y;
-     double Bz = BFIELD_Z;
+     double Bx = TRACK_bfield_x;
+     double By = TRACK_bfield_y;
+     double Bz = TRACK_bfield_z;
 
      double ax =  (q/M)*(Bz*vy - By*vz);
@@ -99 +99 @@
         z  += vz*dt;
 
-       if( (x*x+y*y) > TRACKING_RADIUS*TRACKING_RADIUS ){ x /= (x*x+y*y)/(TRACKING_RADIUS*TRACKING_RADIUS); y /= (x*x+y*y)/(TRACKING_RADIUS*TRACKING_RADIUS); break;}
-       if( fabs(z)>TRACKING_LENGTH)break;
+       if( (x*x+y*y) > TRACK_radius*TRACK_radius ){ x /= (x*x+y*y)/(TRACK_radius*TRACK_radius); y /= (x*x+y*y)/(TRACK_radius*TRACK_radius); break;}
+       if( fabs(z)>TRACK_length)break;
 
        xold = x;