package etomica.modules.dcvgcmd;

import etomica.atom.Atom;
import etomica.atom.AtomList;
import etomica.atom.iterator.AtomIteratorListSimple;
import etomica.config.Conformation;
import etomica.space.Space;

/**
 * @author nsives, msellers
 * 
 * Sets carbon atom coordinates in tubular form.
 *  
 */

public class ConformationTube extends Conformation {

	int atomsPerRing;

	int ctr;

	double tubeRadius;

	double A, B, Q, dz;

	double dtheta, theta, theta0;

	public ConformationTube(Space space, int atomsPerRing) {
		
		super(space);	
		this.atomsPerRing = atomsPerRing;

		dtheta = Math.PI * 2.0 / atomsPerRing;
		
		
		A= 1.42;

		B= 0.705;
		
		//Distance between ORTHO carbons in hexagonal ring, x-z plane
		Q = 2.45;

		//Trig operation on right triangle formed, x-y plane
		tubeRadius = Q / 2.0/Math.sin(dtheta / 2);
		
	}	

	public void initializePositions(AtomList atomList) {

		int size = atomList.size();
		if (size == 0)
			return;

		AtomIteratorListSimple atomIterator = new AtomIteratorListSimple();
		atomIterator.setList(atomList);
		atomIterator.reset();

		int N = 0;
		double theta, theta0, x, y, z, dz;
		int ctr;
		
//		tubeRadius = atomsPerRing*.5;
		x = tubeRadius;
		y = 0;
		z = 0;
		dz = 0;
		theta = 0;
		theta0 = 0;
		ctr = 0;

		while (atomIterator.hasNext()) {
			Atom a = atomIterator.nextAtom();

			a.coord.position().setX(0, x);
			a.coord.position().setX(1, y);
			a.coord.position().setX(2, z);

			N = N + 1;

			if (N == atomsPerRing) {
				N = 0;
				switch (ctr) {
				case (0):
					dz = A;
					theta0 = 0;
					ctr = 1;
					break;

				case (1):
					dz = B;
					theta0 = dtheta / 2;
					ctr = 2;
					break;

				case (2):
					dz = A;
					theta0 = dtheta / 2;
					ctr = 3;
					break;

				case (3):
					dz = B;
					theta0 = 0;
					ctr = 0;
					break;
				}
				z = z + dz;
			}

			theta = theta0 + N * dtheta;

			x = tubeRadius * Math.cos(theta);
			y = tubeRadius * Math.sin(theta);
		}
	}

}
@@ -374,19 +374,8 @@
 g;%0A%0A
-%09int ctr;%0A%0A
 %09dou
@@ -406,20 +406,16 @@
  A, B, Q
-, dz
 ;%0A%0A%09doub
@@ -427,23 +427,8 @@
 heta
-, theta, theta0
 ;%0A%0A%09

remove unused variables