/*
      //
      // ---------------- values for thick optics (lowb.coll) ----------------- 
      //
      // --- TCTH.4R5.B2 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000, 
			   betx =  1587.445970, // value at start of collimator
			   bety =   598.119830,
			   alfx =    -1.759930,
			   alfy =    -9.800047,
			   betax_mkd = 534.46,
			   betay_mkd = 0.0,
			   collname = "TCTH.4R5.B2",
			   icoll = 52,
			   angle = 0.0);
      
      // --- TCP.C6L7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx =  150.939,
			   bety =   82.5329,
			   alfx =    2.05448,
			   alfy =   -1.1503,
			   betax_mkd = 542.112,
			   betay_mkd = 0.0,
			   collname = "TCP.C6L7.B1",
			   icoll = 11,
			   angle = 0.0);
      
      // --- TCSG.B4L7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx =  141.193,
			   bety =   129.752,
			   alfx =    1.44029,
			   alfy =   -1.23866,
			   betax_mkd = 542.112,
			   betay_mkd = 0.0,
			   collname = "TCSG.B4L7.B1",
			   icoll = 26,
			   angle = 0.0);

      // --- TCSM.B4L7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx =  132.747,
			   bety =   137.36,
			   alfx =    1.37496,
			   alfy =   -1.29725,
			   betax_mkd = 542.112,
			   betay_mkd = 0.0,
			   collname = "TCSM.B4L7.B1",
			   icoll = 27,
			   angle = 0.0);

      // --- TCSG.6R7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx =  338.612,
			   bety =   47.3792,
			   alfx =    2.84627,
			   alfy =   0.028661,
			   betax_mkd = 542.112,
			   betay_mkd = 0.0,
			   collname = "TCSG.6R7.B1",
			   icoll = 44,
			   angle = 0.9E-02);

      // --- TCSM.6R7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx =  321.776,
			   bety =   47.3973,
			   alfx =    2.76564,
			   alfy =   -0.03471,
			   betax_mkd = 542.112,
			   betay_mkd = 0.0,
			   collname = "TCSM.6R7.B1",
			   icoll = 45,
			   angle = 0.9E-02);
      */      

      //
      // ---------------- values for thin optics (lowb.coll) ----------------- 
      //
      // --- TCTH.4R5.B2 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000, 
			   betx =  1646.65,
			   bety =   624.456,
			   alfx =    -1.51735,
			   alfy =    -9.93876,
			   betax_mkd = 513.834834,
			   betay_mkd = 0.0,
			   collname = "TCTH.4R5.B2",
			   icoll = 52,
			   angle = 0.0);
      
      // --- TCP.C6L7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx =  153.6886405004573,
			   bety =   81.09542777148614,
			   alfx =    2.07677,
			   alfy =   -1.13737,
			   betax_mkd = 561.419,
			   betay_mkd = 0.0,
			   collname = "TCP.C6L7.B1",
			   icoll = 11,
			   angle = 0.0);

      // --- TCSG.B4L7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx = 136.982 ,
			   bety = 133.156,
			   alfx =  1.38302,
			   alfy = -1.25464 ,
			   betax_mkd = 561.419,
			   betay_mkd = 0.0,
			   collname = "TCSG.B4L7.B1",
			   icoll = 26,
			   angle = 0.0);

      // --- TCSM.B4L7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx = 131.535,
			   bety = 138.252,
			   alfx =  1.34049,
			   alfy = -1.29331,
			   betax_mkd = 561.419,
			   betay_mkd = 0.0,
			   collname = "TCSM.B4L7.B1",
			   icoll = 27,
			   angle = 0.0);

      // --- TCSG.6R7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx =  340.822,
			   bety =   46.419,
			   alfx =    2.89808,
			   alfy =   0.026019,
			   betax_mkd = 561.419,
			   betay_mkd = 0.0,
			   collname = "TCSG.6R7.B1",
			   icoll = 44,
			   angle = 0.9E-02);

      // --- TCSM.6R7.B1 --- input
      ParticleDistribution(xk,0,seed+nd,cut,
			   energy = 7000,
			   betx =  329.34,
			   bety =   46.4012,
			   alfx =    2.84292,
			   alfy =   -0.017096,
			   betax_mkd = 561.419,
			   betay_mkd = 0.0,
			   collname = "TCSM.6R7.B1",
			   icoll = 45,
			   angle = 0.9E-02);