Changeset 32 in TRACY3 for trunk/tracy/tracy/src/nsls-ii_lib.cc

Timestamp:

Apr 9, 2014, 3:50:11 PM (10 years ago)

Author:

zhangj

Message:

active the transport of the twiss functions and orbits of the transfer line.

File:

• trunk/tracy/tracy/src/nsls-ii_lib.cc (modified) (9 diffs)

trunk/tracy/tracy/src/nsls-ii_lib.cc

@@ -763 +763 @@
                       1.5  sextupole
                       2.0  bpm 
-      Fnum        family index
+      Fnum        family index of the BPM.
       all         true, print the cod at all elements
-                  false, print the cod at the family elements
+                  false, print the cod at the BPM elements, the name of the BPM must be defined in the "*.prm" file.
    Output:
        none
@@ -784 +784 @@
 void prt_cod(const char *file_name, const int Fnum, const bool all)
 {
-  long      i;
+  long      i=0L;
   double    code = 0.0;
   FILE      *outf;
-  long      FORLIM;
+  long      FORLIM=0L;
   struct    tm *newtime;
 
@@ -828 +828 @@
       }
       /* COD is in local coordinates */
-      fprintf(outf, "%4ld %.*s %6.2f %4.1f %6.3f %6.3f %6.3f %6.3f"
-              " %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f\n",
+      fprintf(outf, "%4ld %.*s %6.2f %4.1f %6.4f %6.4f %6.4f %6.4f"
+              " %6.4f %6.4f %6.4f %6.4f %6.4f %6.4f\n",
               i, SymbolLength, Cell[i].Elem.PName, Cell[i].S, code,
               Cell[i].Beta[X_], Cell[i].Nu[X_],
@@ -842 +842 @@
 }
 
+/****************************************************************
+ * void prt_beampos(const char *file_name)
+ * 
+ *  print the orbits at all lattice elements.
+ * 
+ * Comments:   
+ *     To use this command, the BPM name should be defined in the
+ *  "*.prm" file with the command 
+ *                bpm_name   BPM_names
+ *
+ * 
+ ****************************************************************/
 
 void prt_beampos(const char *file_name)
 {
-  int       k;
+  int       i=0;
   ofstream  outf;
 
   file_wr(outf, file_name);
 
-  outf << "# k  s  name x   px    y   py   delta ct" << endl;
+  outf << "#   i     name              s[m]     x[mm]   px    y[mm]   py   delta ct" << endl;
   outf << "#" << endl;
 
-  for (k = 0; k <= globval.Cell_nLoc; k++)
-    outf << scientific << setprecision(5)
-         << setw(5) << k << setw(11) << Cell[k].Elem.PName
-         << setw(13) << Cell[k].BeamPos << endl;
-
+  
+  for (i = 0; i <= globval.Cell_nLoc; i++){
+   if (Cell[i].Fnum == globval.bpm)
+        outf << scientific << setprecision(4)
+          << setw(5) << i << setw(20) << Cell[i].Elem.PName 
+          << setw(5) << Cell[i].S 
+          << setw(13) << Cell[i].BeamPos*1e3 << endl;
+  }     
   outf.close();
 }
@@ -874 +889 @@
   // the dT and roll angle are all printed out
 {
-  int i,j;
-  int  n_girders;
-  int  gs_Fnum, ge_Fnum;
-  int  gs_nKid, ge_nKid;
-  int  dip_Fnum,dip_nKid;
-  int  loc, loc_gs, loc_ge;
+  int i=0,j=0;
+  int  n_girders=0;
+  int  gs_Fnum=0, ge_Fnum=0;
+  int  gs_nKid=0, ge_nKid=0;
+  int  dip_Fnum=0,dip_nKid=0;
+  int  loc=0, loc_gs=0, loc_ge=0;
   char * name;
-  double s;
+  double s=0.0;
   double PdSsys[2], PdSrms[2], PdSrnd[2], dS[2], dT[2];
   fstream fout;
@@ -5285 +5300 @@
   Purpose:
            Calculate momentum compact factor up to third order.
-           Method:
-           track the orbit offset c*tau at the first element,
-           at 11 different energy offset(-10-3 to 10-3), then 
-           use polynomal to fit the momentum compact faction factor
-           up to 3rd order.
-           The initial coorinates are (x_co,px_co,y_co,py_co,delta,0).
-           
+           Method:
+           track the orbit offset c*tau at the first element,
+           at 11 different energy offset(-10-3 to 10-3), then 
+           use polynomal to fit the momentum compact faction factor
+           up to 3rd order.
+          The initial coorinates are (x_co,px_co,y_co,py_co,delta,0).
+          
  Input:
        none
@@ -5307 +5322 @@
        change the initial coorinates from (0 0 0 0 delta 0) to 
               (x_co,px_co,y_co,py_co,delta,0).  i.e., track
-              around the off-momentum close orbit but not zero orbit.      
+       around the off-momentum close orbit but not zero orbit.      
 
 *******************************************************************/
@@ -5320 +5335 @@
  //const double  d_delta  = 1e-4;
 
-  int       i, n;
-  long int  lastpos;  // last tracking position when the particle is stable 
-  double    delta[2*n_points+1], alphac[2*n_points+1], sigma;
+  int       i=0, n=0;
+  long int  lastpos=0L;  // last tracking position when the particle is stable 
+  double    delta[2*n_points+1], alphac[2*n_points+1], sigma=0.0;
   Vector    x, b;
   CellType  Cell;
-  bool      cod;
+  bool      cod=false;
   
   globval.pathlength = false;
   getelem(globval.Cell_nLoc, &Cell); 
-  n = 0;
   
   for (i = -n_points; i <= n_points; i++) {
@@ -5353 +5367 @@
 
 
+/**********************************************************
+ * float f_bend(float b0L[])
+ * 
+ * 
+ * 
+ **********************************************************/
 float f_bend(float b0L[])
 {