Changeset 32 in TRACY3 for trunk/tracy/tracy/src/soleilcommon.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/soleilcommon.cc (modified) (5 diffs)

trunk/tracy/tracy/src/soleilcommon.cc

@@ -153 +153 @@
     /*check whether the RF cavity hormonic number is defined or not*/
     if(Cell[Elem_GetPos(globval.cav, 1)].Elem.C->Ph == 0){
-      printf("soleilcommon():    Error!!! Please define RF harmonic number in the lattice! \n"
-              "Because this parameter will be used to calculate stable phase for positive / negative momentum compaction factor\n");
+      printf("soleilcommon():    Error!!! Please define the RF harmonic number in the lattice! \n"
+              "Because this parameter will be used to calculate the stable phase for the lattice with the "
+              "positive / negative momentum compaction factor\n");
       exit_(1);
     }
@@ -169 +170 @@
     globval.emittance  = false; /* emittance  on/off */
     globval.pathlength = false; /* Path lengthening computation */
-    globval.CODimax    = 40L;   /* maximum number of iterations for COD algo */
+    globval.CODimax    = 40L;   /* maximum number of iterations for COD algorithm */
     globval.dPcommon   = 1e-10; /* Common energy step for energy differentiation */
     globval.delta_RF  = RFacceptance;/* energy acceptance for SOLEIL */
 
-   /* define x/y physical aperture, use the default values: +- 1 meter  */
+   /* define x/y physical aperture (vacuum chamber), use the default values: +- 1 meter  */
     globval.Aperture_on = false;
     
@@ -183 +184 @@
   else 
   { // for transfer lines  
-    /* Initial settings : */
-    beta[0]  = 8.1;
-    alpha[0] = 0.0;
-    beta[1]  = 8.1;
-    alpha[1] = 0.0;
-    eta[0]   = 0.0;
-    etap[0]  = 0.0;
-    eta[1]   = 0.0;
-    etap[1]  = 0.0;
-
+    /* Initial optical parameters of the transfer line: */
+    beta[0]  = 34.46;   // beta_x
+    beta[1]  = 33.94;   // beta_y
+    alpha[0] = -4.24;   // alpha_x
+    alpha[1] = -4.34;   // alpha_y
+    eta[0]   = 0.0;   // eta_x
+    etap[0]  = 0.0;   // etap_x
+    eta[1]   = 0.0;   // eta_y
+    etap[1]  = 0.0;   // etap_y
+
+    // set the initial orbit of the transfer line is zero
      for (i = 0; i < ss_dim; i++) {
     {
@@ -200 +202 @@
     dP = codvect[4];
 
+    codvect[0]=1e-4;
+    
     /* Defines global variables for Tracy code */
     globval.MatMeth = false;  /* matrix method */
@@ -214 +218 @@
     globval.Aperture_on = false;
     
-    
+    // transfer the optical functions and the orbit
     TransTwiss(alpha, beta, eta, etap, codvect);
+  
+    // print location, twiss parameters and close orbit/orbit at all elements position to a file
+      prt_cod("cod.out", globval.bpm, true);
     }
   }//finish read the transfer line