--- branches/development/src/mdParser/MDTreeParser.cpp	2012/08/08 18:45:52	1775
+++ branches/development/src/mdParser/MDTreeParser.cpp	2012/08/09 15:52:59	1776
@@ -1,4 +1,4 @@
-/* $ANTLR 2.7.7 (20110725): "MDTreeParser.g" -> "MDTreeParser.cpp"$ */
+/* $ANTLR 2.7.7 (20120725): "MDTreeParser.g" -> "MDTreeParser.cpp"$ */
 #include "MDTreeParser.hpp"
 #include <antlr/Token.hpp>
 #include <antlr/AST.hpp>
@@ -141,11 +141,11 @@ void MDTreeParser::componentblock(ANTLR_USE_NAMESPACE(
 	ANTLR_USE_NAMESPACE(antlr)RefAST componentblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t9 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t10 = _t;
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp2_AST_in = _t;
 		match(_t,COMPONENT);
 		_t = _t->getFirstChild();
-#line 65 "MDTreeParser.g"
+#line 72 "MDTreeParser.g"
 		Component* currComponet = new Component(); blockStack.push(currComponet);
 #line 151 "MDTreeParser.cpp"
 		{ // ( ... )*
@@ -157,18 +157,18 @@ void MDTreeParser::componentblock(ANTLR_USE_NAMESPACE(
 				_t = _retTree;
 			}
 			else {
-				goto _loop11;
+				goto _loop12;
 			}
 			
 		}
-		_loop11:;
+		_loop12:;
 		} // ( ... )*
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp3_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t9;
+		_t = __t10;
 		_t = _t->getNextSibling();
-#line 67 "MDTreeParser.g"
+#line 74 "MDTreeParser.g"
 		blockStack.top()->validate();blockStack.pop(); currConf->addComponent(currComponet);
 #line 174 "MDTreeParser.cpp"
 	}
@@ -184,11 +184,11 @@ void MDTreeParser::moleculeblock(ANTLR_USE_NAMESPACE(a
 	ANTLR_USE_NAMESPACE(antlr)RefAST moleculeblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t33 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t34 = _t;
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp4_AST_in = _t;
 		match(_t,MOLECULE);
 		_t = _t->getFirstChild();
-#line 96 "MDTreeParser.g"
+#line 103 "MDTreeParser.g"
 		MoleculeStamp* currMoleculeStamp = new MoleculeStamp(); blockStack.push(currMoleculeStamp);
 #line 194 "MDTreeParser.cpp"
 		{ // ( ... )*
@@ -200,18 +200,18 @@ void MDTreeParser::moleculeblock(ANTLR_USE_NAMESPACE(a
 				_t = _retTree;
 			}
 			else {
-				goto _loop35;
+				goto _loop36;
 			}
 			
 		}
-		_loop35:;
+		_loop36:;
 		} // ( ... )*
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp5_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t33;
+		_t = __t34;
 		_t = _t->getNextSibling();
-#line 98 "MDTreeParser.g"
+#line 105 "MDTreeParser.g"
 		blockStack.top()->validate(); blockStack.pop(); currConf->addMoleculeStamp(currMoleculeStamp);
 #line 217 "MDTreeParser.cpp"
 	}
@@ -227,11 +227,11 @@ void MDTreeParser::zconstraintblock(ANTLR_USE_NAMESPAC
 	ANTLR_USE_NAMESPACE(antlr)RefAST zconstraintblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t13 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t14 = _t;
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp6_AST_in = _t;
 		match(_t,ZCONSTRAINT);
 		_t = _t->getFirstChild();
-#line 70 "MDTreeParser.g"
+#line 77 "MDTreeParser.g"
 		ZConsStamp* currZConsStamp = new ZConsStamp(); blockStack.push(currZConsStamp);
 #line 237 "MDTreeParser.cpp"
 		{ // ( ... )*
@@ -243,18 +243,18 @@ void MDTreeParser::zconstraintblock(ANTLR_USE_NAMESPAC
 				_t = _retTree;
 			}
 			else {
-				goto _loop15;
+				goto _loop16;
 			}
 			
 		}
-		_loop15:;
+		_loop16:;
 		} // ( ... )*
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp7_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t13;
+		_t = __t14;
 		_t = _t->getNextSibling();
-#line 72 "MDTreeParser.g"
+#line 79 "MDTreeParser.g"
 		blockStack.top()->validate();blockStack.pop(); currConf->addZConsStamp(currZConsStamp);
 #line 260 "MDTreeParser.cpp"
 	}
@@ -270,11 +270,11 @@ void MDTreeParser::restraintblock(ANTLR_USE_NAMESPACE(
 	ANTLR_USE_NAMESPACE(antlr)RefAST restraintblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t17 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t18 = _t;
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp8_AST_in = _t;
 		match(_t,RESTRAINT);
 		_t = _t->getFirstChild();
-#line 75 "MDTreeParser.g"
+#line 82 "MDTreeParser.g"
 		RestraintStamp* currRestraintStamp = new RestraintStamp(); blockStack.push(currRestraintStamp);
 #line 280 "MDTreeParser.cpp"
 		{ // ( ... )*
@@ -286,18 +286,18 @@ void MDTreeParser::restraintblock(ANTLR_USE_NAMESPACE(
 				_t = _retTree;
 			}
 			else {
-				goto _loop19;
+				goto _loop20;
 			}
 			
 		}
-		_loop19:;
+		_loop20:;
 		} // ( ... )*
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp9_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t17;
+		_t = __t18;
 		_t = _t->getNextSibling();
-#line 77 "MDTreeParser.g"
+#line 84 "MDTreeParser.g"
 		blockStack.top()->validate();blockStack.pop(); currConf->addRestraintStamp(currRestraintStamp);
 #line 303 "MDTreeParser.cpp"
 	}
@@ -313,11 +313,11 @@ void MDTreeParser::flucqblock(ANTLR_USE_NAMESPACE(antl
 	ANTLR_USE_NAMESPACE(antlr)RefAST flucqblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t21 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t22 = _t;
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp10_AST_in = _t;
 		match(_t,FLUCQ);
 		_t = _t->getFirstChild();
-#line 80 "MDTreeParser.g"
+#line 87 "MDTreeParser.g"
 		FluctuatingChargeParameters* flucQpars = new FluctuatingChargeParameters(); blockStack.push(flucQpars);
 #line 323 "MDTreeParser.cpp"
 		{ // ( ... )*
@@ -329,18 +329,18 @@ void MDTreeParser::flucqblock(ANTLR_USE_NAMESPACE(antl
 				_t = _retTree;
 			}
 			else {
-				goto _loop23;
+				goto _loop24;
 			}
 			
 		}
-		_loop23:;
+		_loop24:;
 		} // ( ... )*
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp11_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t21;
+		_t = __t22;
 		_t = _t->getNextSibling();
-#line 82 "MDTreeParser.g"
+#line 89 "MDTreeParser.g"
 		blockStack.top()->validate();blockStack.pop(); currConf->addFluctuatingChargeParameters(flucQpars);
 #line 346 "MDTreeParser.cpp"
 	}
@@ -356,11 +356,11 @@ void MDTreeParser::rnemdblock(ANTLR_USE_NAMESPACE(antl
 	ANTLR_USE_NAMESPACE(antlr)RefAST rnemdblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t25 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t26 = _t;
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp12_AST_in = _t;
 		match(_t,RNEMD);
 		_t = _t->getFirstChild();
-#line 85 "MDTreeParser.g"
+#line 92 "MDTreeParser.g"
 		RNEMDParameters* rnemdPars = new RNEMDParameters(); blockStack.push(rnemdPars);
 #line 366 "MDTreeParser.cpp"
 		{ // ( ... )*
@@ -372,18 +372,18 @@ void MDTreeParser::rnemdblock(ANTLR_USE_NAMESPACE(antl
 				_t = _retTree;
 			}
 			else {
-				goto _loop27;
+				goto _loop28;
 			}
 			
 		}
-		_loop27:;
+		_loop28:;
 		} // ( ... )*
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp13_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t25;
+		_t = __t26;
 		_t = _t->getNextSibling();
-#line 87 "MDTreeParser.g"
+#line 94 "MDTreeParser.g"
 		blockStack.top()->validate();blockStack.pop(); currConf->addRNEMDParameters(rnemdPars);
 #line 389 "MDTreeParser.cpp"
 	}
@@ -399,11 +399,11 @@ void MDTreeParser::minimizerblock(ANTLR_USE_NAMESPACE(
 	ANTLR_USE_NAMESPACE(antlr)RefAST minimizerblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t29 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t30 = _t;
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp14_AST_in = _t;
 		match(_t,MINIMIZER);
 		_t = _t->getFirstChild();
-#line 90 "MDTreeParser.g"
+#line 97 "MDTreeParser.g"
 		MinimizerParameters* minimizerPars = new MinimizerParameters(); blockStack.push(minimizerPars);
 #line 409 "MDTreeParser.cpp"
 		{ // ( ... )*
@@ -415,18 +415,18 @@ void MDTreeParser::minimizerblock(ANTLR_USE_NAMESPACE(
 				_t = _retTree;
 			}
 			else {
-				goto _loop31;
+				goto _loop32;
 			}
 			
 		}
-		_loop31:;
+		_loop32:;
 		} // ( ... )*
 		ANTLR_USE_NAMESPACE(antlr)RefAST tmp15_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t29;
+		_t = __t30;
 		_t = _t->getNextSibling();
-#line 92 "MDTreeParser.g"
+#line 99 "MDTreeParser.g"
 		blockStack.top()->validate();blockStack.pop(); currConf->addMinimizerParameters(minimizerPars);
 #line 432 "MDTreeParser.cpp"
 	}
@@ -447,9 +447,10 @@ void MDTreeParser::constant(ANTLR_USE_NAMESPACE(antlr)
 #line 50 "MDTreeParser.g"
 	
 	int ival;
-	RealType dval;
+	RealType dval, x, y, z;
+	Vector3d dvec;
 	
-#line 453 "MDTreeParser.cpp"
+#line 454 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -460,9 +461,9 @@ void MDTreeParser::constant(ANTLR_USE_NAMESPACE(antlr)
 		{
 			ival=intConst(_t);
 			_t = _retTree;
-#line 55 "MDTreeParser.g"
+#line 56 "MDTreeParser.g"
 			blockStack.top()->assign(id->getText(), ival);
-#line 466 "MDTreeParser.cpp"
+#line 467 "MDTreeParser.cpp"
 			break;
 		}
 		case NUM_FLOAT:
@@ -470,19 +471,51 @@ void MDTreeParser::constant(ANTLR_USE_NAMESPACE(antlr)
 		{
 			dval=floatConst(_t);
 			_t = _retTree;
-#line 56 "MDTreeParser.g"
+#line 57 "MDTreeParser.g"
 			blockStack.top()->assign(id->getText(), dval);
-#line 476 "MDTreeParser.cpp"
+#line 477 "MDTreeParser.cpp"
 			break;
 		}
+		case LPAREN:
+		{
+			ANTLR_USE_NAMESPACE(antlr)RefAST __t8 = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp16_AST_in = _t;
+			match(_t,LPAREN);
+			_t = _t->getFirstChild();
+			x=doubleNumber(_t);
+			_t = _retTree;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp17_AST_in = _t;
+			match(_t,COMMA);
+			_t = _t->getNextSibling();
+			y=doubleNumber(_t);
+			_t = _retTree;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp18_AST_in = _t;
+			match(_t,COMMA);
+			_t = _t->getNextSibling();
+			z=doubleNumber(_t);
+			_t = _retTree;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp19_AST_in = _t;
+			match(_t,RPAREN);
+			_t = _t->getNextSibling();
+			_t = __t8;
+			_t = _t->getNextSibling();
+#line 59 "MDTreeParser.g"
+			dvec.x() = x; 
+			dvec.y() = y; 
+			dvec.z() = z; 
+			blockStack.top()->assign(id->getText(), dvec);
+			
+#line 509 "MDTreeParser.cpp"
+			break;
+		}
 		case ID:
 		{
 			str1 = _t;
 			match(_t,ID);
 			_t = _t->getNextSibling();
-#line 57 "MDTreeParser.g"
+#line 64 "MDTreeParser.g"
 			blockStack.top()->assign(id->getText(), str1->getText());
-#line 486 "MDTreeParser.cpp"
+#line 519 "MDTreeParser.cpp"
 			break;
 		}
 		case StringLiteral:
@@ -490,12 +523,12 @@ void MDTreeParser::constant(ANTLR_USE_NAMESPACE(antlr)
 			str2 = _t;
 			match(_t,StringLiteral);
 			_t = _t->getNextSibling();
-#line 58 "MDTreeParser.g"
+#line 65 "MDTreeParser.g"
 			std::string s =  str2->getText();
 			s = s.substr(1, s.length()-2);
 			blockStack.top()->assign(id->getText(),s);
 			
-#line 499 "MDTreeParser.cpp"
+#line 532 "MDTreeParser.cpp"
 			break;
 		}
 		default:
@@ -513,9 +546,9 @@ int  MDTreeParser::intConst(ANTLR_USE_NAMESPACE(antlr)
 }
 
 int  MDTreeParser::intConst(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
-#line 286 "MDTreeParser.g"
+#line 294 "MDTreeParser.g"
 	int ival;
-#line 519 "MDTreeParser.cpp"
+#line 552 "MDTreeParser.cpp"
 	ANTLR_USE_NAMESPACE(antlr)RefAST intConst_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	ANTLR_USE_NAMESPACE(antlr)RefAST i1 = ANTLR_USE_NAMESPACE(antlr)nullAST;
 	ANTLR_USE_NAMESPACE(antlr)RefAST i2 = ANTLR_USE_NAMESPACE(antlr)nullAST;
@@ -529,9 +562,9 @@ int  MDTreeParser::intConst(ANTLR_USE_NAMESPACE(antlr)
 			i1 = _t;
 			match(_t,NUM_INT);
 			_t = _t->getNextSibling();
-#line 287 "MDTreeParser.g"
+#line 295 "MDTreeParser.g"
 			ival = lexi_cast<int>(i1->getText());
-#line 535 "MDTreeParser.cpp"
+#line 568 "MDTreeParser.cpp"
 			break;
 		}
 		case NUM_LONG:
@@ -539,9 +572,9 @@ int  MDTreeParser::intConst(ANTLR_USE_NAMESPACE(antlr)
 			i2 = _t;
 			match(_t,NUM_LONG);
 			_t = _t->getNextSibling();
-#line 288 "MDTreeParser.g"
+#line 296 "MDTreeParser.g"
 			ival = lexi_cast<int>(i2->getText());
-#line 545 "MDTreeParser.cpp"
+#line 578 "MDTreeParser.cpp"
 			break;
 		}
 		default:
@@ -560,9 +593,9 @@ RealType  MDTreeParser::floatConst(ANTLR_USE_NAMESPACE
 }
 
 RealType  MDTreeParser::floatConst(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
-#line 300 "MDTreeParser.g"
+#line 306 "MDTreeParser.g"
 	RealType dval;
-#line 566 "MDTreeParser.cpp"
+#line 599 "MDTreeParser.cpp"
 	ANTLR_USE_NAMESPACE(antlr)RefAST floatConst_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	ANTLR_USE_NAMESPACE(antlr)RefAST d1 = ANTLR_USE_NAMESPACE(antlr)nullAST;
 	ANTLR_USE_NAMESPACE(antlr)RefAST d2 = ANTLR_USE_NAMESPACE(antlr)nullAST;
@@ -576,9 +609,9 @@ RealType  MDTreeParser::floatConst(ANTLR_USE_NAMESPACE
 			d1 = _t;
 			match(_t,NUM_FLOAT);
 			_t = _t->getNextSibling();
-#line 301 "MDTreeParser.g"
+#line 307 "MDTreeParser.g"
 			dval = lexi_cast<RealType>(d1->getText());
-#line 582 "MDTreeParser.cpp"
+#line 615 "MDTreeParser.cpp"
 			break;
 		}
 		case NUM_DOUBLE:
@@ -586,11 +619,60 @@ RealType  MDTreeParser::floatConst(ANTLR_USE_NAMESPACE
 			d2 = _t;
 			match(_t,NUM_DOUBLE);
 			_t = _t->getNextSibling();
-#line 302 "MDTreeParser.g"
+#line 308 "MDTreeParser.g"
 			dval = lexi_cast<RealType>(d2->getText());
-#line 592 "MDTreeParser.cpp"
+#line 625 "MDTreeParser.cpp"
+			break;
+		}
+		default:
+		{
+			throw ANTLR_USE_NAMESPACE(antlr)NoViableAltException(_t);
+		}
+		}
+	}
+	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
+		reportError(ex);
+		if ( _t != ANTLR_USE_NAMESPACE(antlr)nullAST )
+			_t = _t->getNextSibling();
+	}
+	_retTree = _t;
+	return dval;
+}
+
+RealType  MDTreeParser::doubleNumber(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
+#line 300 "MDTreeParser.g"
+	RealType dval;
+#line 646 "MDTreeParser.cpp"
+	ANTLR_USE_NAMESPACE(antlr)RefAST doubleNumber_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
+	ANTLR_USE_NAMESPACE(antlr)RefAST ic = ANTLR_USE_NAMESPACE(antlr)nullAST;
+	ANTLR_USE_NAMESPACE(antlr)RefAST fc = ANTLR_USE_NAMESPACE(antlr)nullAST;
+	
+	try {      // for error handling
+		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
+			_t = ASTNULL;
+		switch ( _t->getType()) {
+		case NUM_INT:
+		case NUM_LONG:
+		{
+			ic = (_t == ASTNULL) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
+			intConst(_t);
+			_t = _retTree;
+#line 301 "MDTreeParser.g"
+			dval = lexi_cast<RealType>(ic->getText());
+#line 663 "MDTreeParser.cpp"
 			break;
 		}
+		case NUM_FLOAT:
+		case NUM_DOUBLE:
+		{
+			fc = (_t == ASTNULL) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
+			floatConst(_t);
+			_t = _retTree;
+#line 302 "MDTreeParser.g"
+			dval = lexi_cast<RealType>(fc->getText());
+#line 674 "MDTreeParser.cpp"
+			break;
+		}
 		default:
 		{
 			throw ANTLR_USE_NAMESPACE(antlr)NoViableAltException(_t);
@@ -683,22 +765,22 @@ void MDTreeParser::atomblock(ANTLR_USE_NAMESPACE(antlr
 
 void MDTreeParser::atomblock(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST atomblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 112 "MDTreeParser.g"
+#line 119 "MDTreeParser.g"
 	
 	int index;
 	
-#line 691 "MDTreeParser.cpp"
+#line 773 "MDTreeParser.cpp"
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t38 = _t;
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp16_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t39 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp20_AST_in = _t;
 		match(_t,ATOM);
 		_t = _t->getFirstChild();
 		index=intConst(_t);
 		_t = _retTree;
-#line 116 "MDTreeParser.g"
+#line 123 "MDTreeParser.g"
 		AtomStamp* currAtomStamp = new AtomStamp(index); blockStack.push(currAtomStamp);
-#line 702 "MDTreeParser.cpp"
+#line 784 "MDTreeParser.cpp"
 		{ // ( ... )*
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -708,25 +790,25 @@ void MDTreeParser::atomblock(ANTLR_USE_NAMESPACE(antlr
 				_t = _retTree;
 			}
 			else {
-				goto _loop40;
+				goto _loop41;
 			}
 			
 		}
-		_loop40:;
+		_loop41:;
 		} // ( ... )*
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp17_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp21_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t38;
+		_t = __t39;
 		_t = _t->getNextSibling();
-#line 118 "MDTreeParser.g"
+#line 125 "MDTreeParser.g"
 		
 		blockStack.top()->validate();
 		blockStack.pop(); 
 		MoleculeStamp* currMoleculeStamp = static_cast<MoleculeStamp*>(blockStack.top());
 		currMoleculeStamp->addAtomStamp(currAtomStamp); 
 		
-#line 730 "MDTreeParser.cpp"
+#line 812 "MDTreeParser.cpp"
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
 		reportError(ex);
@@ -740,13 +822,13 @@ void MDTreeParser::bondblock(ANTLR_USE_NAMESPACE(antlr
 	ANTLR_USE_NAMESPACE(antlr)RefAST bondblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t45 = _t;
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp18_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t46 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp22_AST_in = _t;
 		match(_t,BOND);
 		_t = _t->getFirstChild();
-#line 138 "MDTreeParser.g"
+#line 145 "MDTreeParser.g"
 		BondStamp* currBondStamp = new BondStamp(); blockStack.push(currBondStamp);
-#line 750 "MDTreeParser.cpp"
+#line 832 "MDTreeParser.cpp"
 		{ // ( ... )*
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -756,24 +838,24 @@ void MDTreeParser::bondblock(ANTLR_USE_NAMESPACE(antlr
 				_t = _retTree;
 			}
 			else {
-				goto _loop47;
+				goto _loop48;
 			}
 			
 		}
-		_loop47:;
+		_loop48:;
 		} // ( ... )*
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp19_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp23_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t45;
+		_t = __t46;
 		_t = _t->getNextSibling();
-#line 140 "MDTreeParser.g"
+#line 147 "MDTreeParser.g"
 		
 		blockStack.pop(); 
 		MoleculeStamp* currMoleculeStamp = static_cast<MoleculeStamp*>(blockStack.top());
 		currMoleculeStamp->addBondStamp(currBondStamp); 
 		
-#line 777 "MDTreeParser.cpp"
+#line 859 "MDTreeParser.cpp"
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
 		reportError(ex);
@@ -787,13 +869,13 @@ void MDTreeParser::bendblock(ANTLR_USE_NAMESPACE(antlr
 	ANTLR_USE_NAMESPACE(antlr)RefAST bendblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t51 = _t;
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp20_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t52 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp24_AST_in = _t;
 		match(_t,BEND);
 		_t = _t->getFirstChild();
-#line 156 "MDTreeParser.g"
+#line 163 "MDTreeParser.g"
 		BendStamp* currBendStamp = new BendStamp(); blockStack.push(currBendStamp);
-#line 797 "MDTreeParser.cpp"
+#line 879 "MDTreeParser.cpp"
 		{ // ( ... )*
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -803,25 +885,25 @@ void MDTreeParser::bendblock(ANTLR_USE_NAMESPACE(antlr
 				_t = _retTree;
 			}
 			else {
-				goto _loop53;
+				goto _loop54;
 			}
 			
 		}
-		_loop53:;
+		_loop54:;
 		} // ( ... )*
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp21_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp25_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t51;
+		_t = __t52;
 		_t = _t->getNextSibling();
-#line 158 "MDTreeParser.g"
+#line 165 "MDTreeParser.g"
 		
 		blockStack.top()->validate();
 		blockStack.pop(); 
 		MoleculeStamp* currMoleculeStamp = static_cast<MoleculeStamp*>(blockStack.top());
 		currMoleculeStamp->addBendStamp(currBendStamp); 
 		
-#line 825 "MDTreeParser.cpp"
+#line 907 "MDTreeParser.cpp"
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
 		reportError(ex);
@@ -835,13 +917,13 @@ void MDTreeParser::torsionblock(ANTLR_USE_NAMESPACE(an
 	ANTLR_USE_NAMESPACE(antlr)RefAST torsionblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t57 = _t;
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp22_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t58 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp26_AST_in = _t;
 		match(_t,TORSION);
 		_t = _t->getFirstChild();
-#line 175 "MDTreeParser.g"
+#line 182 "MDTreeParser.g"
 		TorsionStamp* currTorsionStamp = new TorsionStamp(); blockStack.push(currTorsionStamp);
-#line 845 "MDTreeParser.cpp"
+#line 927 "MDTreeParser.cpp"
 		{ // ( ... )*
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -851,25 +933,25 @@ void MDTreeParser::torsionblock(ANTLR_USE_NAMESPACE(an
 				_t = _retTree;
 			}
 			else {
-				goto _loop59;
+				goto _loop60;
 			}
 			
 		}
-		_loop59:;
+		_loop60:;
 		} // ( ... )*
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp23_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp27_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t57;
+		_t = __t58;
 		_t = _t->getNextSibling();
-#line 177 "MDTreeParser.g"
+#line 184 "MDTreeParser.g"
 		
 		blockStack.top()->validate();
 		blockStack.pop(); 
 		MoleculeStamp* currMoleculeStamp = static_cast<MoleculeStamp*>(blockStack.top());
 		currMoleculeStamp->addTorsionStamp(currTorsionStamp); 
 		
-#line 873 "MDTreeParser.cpp"
+#line 955 "MDTreeParser.cpp"
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
 		reportError(ex);
@@ -883,13 +965,13 @@ void MDTreeParser::inversionblock(ANTLR_USE_NAMESPACE(
 	ANTLR_USE_NAMESPACE(antlr)RefAST inversionblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t63 = _t;
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp24_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t64 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp28_AST_in = _t;
 		match(_t,INVERSION);
 		_t = _t->getFirstChild();
-#line 194 "MDTreeParser.g"
+#line 201 "MDTreeParser.g"
 		InversionStamp* currInversionStamp = new InversionStamp(); blockStack.push(currInversionStamp);
-#line 893 "MDTreeParser.cpp"
+#line 975 "MDTreeParser.cpp"
 		{ // ( ... )*
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -899,25 +981,25 @@ void MDTreeParser::inversionblock(ANTLR_USE_NAMESPACE(
 				_t = _retTree;
 			}
 			else {
-				goto _loop65;
+				goto _loop66;
 			}
 			
 		}
-		_loop65:;
+		_loop66:;
 		} // ( ... )*
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp25_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp29_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t63;
+		_t = __t64;
 		_t = _t->getNextSibling();
-#line 196 "MDTreeParser.g"
+#line 203 "MDTreeParser.g"
 		
 		blockStack.top()->validate();
 		blockStack.pop(); 
 		MoleculeStamp* currMoleculeStamp = static_cast<MoleculeStamp*>(blockStack.top());
 		currMoleculeStamp->addInversionStamp(currInversionStamp); 
 		
-#line 921 "MDTreeParser.cpp"
+#line 1003 "MDTreeParser.cpp"
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
 		reportError(ex);
@@ -929,22 +1011,22 @@ void MDTreeParser::rigidbodyblock(ANTLR_USE_NAMESPACE(
 
 void MDTreeParser::rigidbodyblock(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST rigidbodyblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 213 "MDTreeParser.g"
+#line 220 "MDTreeParser.g"
 	
 	int index;
 	
-#line 937 "MDTreeParser.cpp"
+#line 1019 "MDTreeParser.cpp"
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t69 = _t;
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp26_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t70 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp30_AST_in = _t;
 		match(_t,RIGIDBODY);
 		_t = _t->getFirstChild();
 		index=intConst(_t);
 		_t = _retTree;
-#line 217 "MDTreeParser.g"
+#line 224 "MDTreeParser.g"
 		RigidBodyStamp* currRigidBodyStamp = new RigidBodyStamp(index); blockStack.push(currRigidBodyStamp);
-#line 948 "MDTreeParser.cpp"
+#line 1030 "MDTreeParser.cpp"
 		{ // ( ... )*
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -954,25 +1036,25 @@ void MDTreeParser::rigidbodyblock(ANTLR_USE_NAMESPACE(
 				_t = _retTree;
 			}
 			else {
-				goto _loop71;
+				goto _loop72;
 			}
 			
 		}
-		_loop71:;
+		_loop72:;
 		} // ( ... )*
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp27_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp31_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t69;
+		_t = __t70;
 		_t = _t->getNextSibling();
-#line 219 "MDTreeParser.g"
+#line 226 "MDTreeParser.g"
 		
 		blockStack.top()->validate();
 		blockStack.pop(); 
 		MoleculeStamp* currMoleculeStamp = static_cast<MoleculeStamp*>(blockStack.top());
 		currMoleculeStamp->addRigidBodyStamp(currRigidBodyStamp); 
 		
-#line 976 "MDTreeParser.cpp"
+#line 1058 "MDTreeParser.cpp"
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
 		reportError(ex);
@@ -986,13 +1068,13 @@ void MDTreeParser::cutoffgroupblock(ANTLR_USE_NAMESPAC
 	ANTLR_USE_NAMESPACE(antlr)RefAST cutoffgroupblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t75 = _t;
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp28_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t76 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp32_AST_in = _t;
 		match(_t,CUTOFFGROUP);
 		_t = _t->getFirstChild();
-#line 236 "MDTreeParser.g"
+#line 243 "MDTreeParser.g"
 		CutoffGroupStamp* currCutoffGroupStamp = new CutoffGroupStamp(); blockStack.push(currCutoffGroupStamp);
-#line 996 "MDTreeParser.cpp"
+#line 1078 "MDTreeParser.cpp"
 		{ // ( ... )*
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -1002,25 +1084,25 @@ void MDTreeParser::cutoffgroupblock(ANTLR_USE_NAMESPAC
 				_t = _retTree;
 			}
 			else {
-				goto _loop77;
+				goto _loop78;
 			}
 			
 		}
-		_loop77:;
+		_loop78:;
 		} // ( ... )*
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp29_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp33_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t75;
+		_t = __t76;
 		_t = _t->getNextSibling();
-#line 238 "MDTreeParser.g"
+#line 245 "MDTreeParser.g"
 		
 		blockStack.top()->validate();
 		blockStack.pop(); 
 		MoleculeStamp* currMoleculeStamp = static_cast<MoleculeStamp*>(blockStack.top());
 		currMoleculeStamp->addCutoffGroupStamp(currCutoffGroupStamp); 
 		
-#line 1024 "MDTreeParser.cpp"
+#line 1106 "MDTreeParser.cpp"
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
 		reportError(ex);
@@ -1032,20 +1114,20 @@ void MDTreeParser::fragmentblock(ANTLR_USE_NAMESPACE(a
 
 void MDTreeParser::fragmentblock(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST fragmentblock_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 255 "MDTreeParser.g"
+#line 262 "MDTreeParser.g"
 	int ival;
-#line 1038 "MDTreeParser.cpp"
+#line 1120 "MDTreeParser.cpp"
 	
 	try {      // for error handling
-		ANTLR_USE_NAMESPACE(antlr)RefAST __t81 = _t;
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp30_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST __t82 = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp34_AST_in = _t;
 		match(_t,FRAGMENT);
 		_t = _t->getFirstChild();
 		ival=intConst(_t);
 		_t = _retTree;
-#line 256 "MDTreeParser.g"
+#line 263 "MDTreeParser.g"
 		FragmentStamp* currFragmentStamp = new FragmentStamp(ival); blockStack.push(currFragmentStamp);
-#line 1049 "MDTreeParser.cpp"
+#line 1131 "MDTreeParser.cpp"
 		{ // ( ... )*
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -1055,25 +1137,25 @@ void MDTreeParser::fragmentblock(ANTLR_USE_NAMESPACE(a
 				_t = _retTree;
 			}
 			else {
-				goto _loop83;
+				goto _loop84;
 			}
 			
 		}
-		_loop83:;
+		_loop84:;
 		} // ( ... )*
-		ANTLR_USE_NAMESPACE(antlr)RefAST tmp31_AST_in = _t;
+		ANTLR_USE_NAMESPACE(antlr)RefAST tmp35_AST_in = _t;
 		match(_t,ENDBLOCK);
 		_t = _t->getNextSibling();
-		_t = __t81;
+		_t = __t82;
 		_t = _t->getNextSibling();
-#line 258 "MDTreeParser.g"
+#line 265 "MDTreeParser.g"
 		
 		blockStack.top()->validate();
 		blockStack.pop(); 
 		MoleculeStamp* currMoleculeStamp = static_cast<MoleculeStamp*>(blockStack.top());
 		currMoleculeStamp->addFragmentStamp(currFragmentStamp); 
 		
-#line 1077 "MDTreeParser.cpp"
+#line 1159 "MDTreeParser.cpp"
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
 		reportError(ex);
@@ -1085,13 +1167,13 @@ void MDTreeParser::atomstatement(ANTLR_USE_NAMESPACE(a
 
 void MDTreeParser::atomstatement(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST atomstatement_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 126 "MDTreeParser.g"
+#line 133 "MDTreeParser.g"
 	
 	vector<RealType> dvec;
 	AtomStamp* currAtomStamp =  static_cast<AtomStamp*>(blockStack.top());
 	
 	
-#line 1095 "MDTreeParser.cpp"
+#line 1177 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -1105,32 +1187,32 @@ void MDTreeParser::atomstatement(ANTLR_USE_NAMESPACE(a
 		}
 		case POSITION:
 		{
-			ANTLR_USE_NAMESPACE(antlr)RefAST __t42 = _t;
-			ANTLR_USE_NAMESPACE(antlr)RefAST tmp32_AST_in = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST __t43 = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp36_AST_in = _t;
 			match(_t,POSITION);
 			_t = _t->getFirstChild();
 			dvec=doubleNumberTuple(_t);
 			_t = _retTree;
-			_t = __t42;
+			_t = __t43;
 			_t = _t->getNextSibling();
-#line 133 "MDTreeParser.g"
+#line 140 "MDTreeParser.g"
 			currAtomStamp->setPosition(dvec);
-#line 1119 "MDTreeParser.cpp"
+#line 1201 "MDTreeParser.cpp"
 			break;
 		}
 		case ORIENTATION:
 		{
-			ANTLR_USE_NAMESPACE(antlr)RefAST __t43 = _t;
-			ANTLR_USE_NAMESPACE(antlr)RefAST tmp33_AST_in = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST __t44 = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp37_AST_in = _t;
 			match(_t,ORIENTATION);
 			_t = _t->getFirstChild();
 			dvec=doubleNumberTuple(_t);
 			_t = _retTree;
-			_t = __t43;
-			_t = _t->getNextSibling();
-#line 134 "MDTreeParser.g"
+			_t = __t44;
+			_t = _t->getNextSibling();
+#line 141 "MDTreeParser.g"
 			currAtomStamp->setOrientation(dvec);
-#line 1134 "MDTreeParser.cpp"
+#line 1216 "MDTreeParser.cpp"
 			break;
 		}
 		default:
@@ -1148,36 +1230,36 @@ vector<RealType>  MDTreeParser::doubleNumberTuple(ANTL
 }
 
 vector<RealType>  MDTreeParser::doubleNumberTuple(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
-#line 271 "MDTreeParser.g"
+#line 278 "MDTreeParser.g"
 	vector<RealType> dvec;
-#line 1154 "MDTreeParser.cpp"
+#line 1236 "MDTreeParser.cpp"
 	ANTLR_USE_NAMESPACE(antlr)RefAST doubleNumberTuple_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 271 "MDTreeParser.g"
+#line 278 "MDTreeParser.g"
 	
 	RealType dval;
 	
-#line 1160 "MDTreeParser.cpp"
+#line 1242 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		{ // ( ... )+
-		int _cnt87=0;
+		int _cnt88=0;
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
 				_t = ASTNULL;
 			if (((_t->getType() >= NUM_INT && _t->getType() <= NUM_DOUBLE))) {
 				dval=doubleNumber(_t);
 				_t = _retTree;
-#line 275 "MDTreeParser.g"
+#line 282 "MDTreeParser.g"
 				dvec.push_back(dval);
-#line 1173 "MDTreeParser.cpp"
+#line 1255 "MDTreeParser.cpp"
 			}
 			else {
-				if ( _cnt87>=1 ) { goto _loop87; } else {throw ANTLR_USE_NAMESPACE(antlr)NoViableAltException(_t);}
+				if ( _cnt88>=1 ) { goto _loop88; } else {throw ANTLR_USE_NAMESPACE(antlr)NoViableAltException(_t);}
 			}
 			
-			_cnt87++;
+			_cnt88++;
 		}
-		_loop87:;
+		_loop88:;
 		}  // ( ... )+
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
@@ -1191,12 +1273,12 @@ void MDTreeParser::bondstatement(ANTLR_USE_NAMESPACE(a
 
 void MDTreeParser::bondstatement(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST bondstatement_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 147 "MDTreeParser.g"
+#line 154 "MDTreeParser.g"
 	
 	vector<int> ivec;
 	BondStamp* currBondStamp = static_cast<BondStamp*>(blockStack.top());
 	
-#line 1200 "MDTreeParser.cpp"
+#line 1282 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -1210,17 +1292,17 @@ void MDTreeParser::bondstatement(ANTLR_USE_NAMESPACE(a
 		}
 		case MEMBERS:
 		{
-			ANTLR_USE_NAMESPACE(antlr)RefAST __t49 = _t;
-			ANTLR_USE_NAMESPACE(antlr)RefAST tmp34_AST_in = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST __t50 = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp38_AST_in = _t;
 			match(_t,MEMBERS);
 			_t = _t->getFirstChild();
 			ivec=inttuple(_t);
 			_t = _retTree;
-			_t = __t49;
+			_t = __t50;
 			_t = _t->getNextSibling();
-#line 153 "MDTreeParser.g"
+#line 160 "MDTreeParser.g"
 			currBondStamp->setMembers(ivec);
-#line 1224 "MDTreeParser.cpp"
+#line 1306 "MDTreeParser.cpp"
 			break;
 		}
 		default:
@@ -1238,36 +1320,36 @@ vector<int>  MDTreeParser::inttuple(ANTLR_USE_NAMESPAC
 }
 
 vector<int>  MDTreeParser::inttuple(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
-#line 278 "MDTreeParser.g"
+#line 286 "MDTreeParser.g"
 	vector<int> ivec;
-#line 1244 "MDTreeParser.cpp"
+#line 1326 "MDTreeParser.cpp"
 	ANTLR_USE_NAMESPACE(antlr)RefAST inttuple_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 278 "MDTreeParser.g"
+#line 286 "MDTreeParser.g"
 	
 	int ival;
 	
-#line 1250 "MDTreeParser.cpp"
+#line 1332 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		{ // ( ... )+
-		int _cnt90=0;
+		int _cnt91=0;
 		for (;;) {
 			if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
 				_t = ASTNULL;
 			if ((_t->getType() == NUM_INT || _t->getType() == NUM_LONG)) {
 				ival=intConst(_t);
 				_t = _retTree;
-#line 282 "MDTreeParser.g"
+#line 290 "MDTreeParser.g"
 				ivec.push_back(ival);
-#line 1263 "MDTreeParser.cpp"
+#line 1345 "MDTreeParser.cpp"
 			}
 			else {
-				if ( _cnt90>=1 ) { goto _loop90; } else {throw ANTLR_USE_NAMESPACE(antlr)NoViableAltException(_t);}
+				if ( _cnt91>=1 ) { goto _loop91; } else {throw ANTLR_USE_NAMESPACE(antlr)NoViableAltException(_t);}
 			}
 			
-			_cnt90++;
+			_cnt91++;
 		}
-		_loop90:;
+		_loop91:;
 		}  // ( ... )+
 	}
 	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
@@ -1281,12 +1363,12 @@ void MDTreeParser::bendstatement(ANTLR_USE_NAMESPACE(a
 
 void MDTreeParser::bendstatement(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST bendstatement_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 166 "MDTreeParser.g"
+#line 173 "MDTreeParser.g"
 	
 	vector<int> ivec;
 	BendStamp* currBendStamp = static_cast<BendStamp*>(blockStack.top());
 	
-#line 1290 "MDTreeParser.cpp"
+#line 1372 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -1300,17 +1382,17 @@ void MDTreeParser::bendstatement(ANTLR_USE_NAMESPACE(a
 		}
 		case MEMBERS:
 		{
-			ANTLR_USE_NAMESPACE(antlr)RefAST __t55 = _t;
-			ANTLR_USE_NAMESPACE(antlr)RefAST tmp35_AST_in = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST __t56 = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp39_AST_in = _t;
 			match(_t,MEMBERS);
 			_t = _t->getFirstChild();
 			ivec=inttuple(_t);
 			_t = _retTree;
-			_t = __t55;
+			_t = __t56;
 			_t = _t->getNextSibling();
-#line 172 "MDTreeParser.g"
+#line 179 "MDTreeParser.g"
 			currBendStamp->setMembers(ivec);
-#line 1314 "MDTreeParser.cpp"
+#line 1396 "MDTreeParser.cpp"
 			break;
 		}
 		default:
@@ -1329,12 +1411,12 @@ void MDTreeParser::torsionstatement(ANTLR_USE_NAMESPAC
 
 void MDTreeParser::torsionstatement(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST torsionstatement_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 185 "MDTreeParser.g"
+#line 192 "MDTreeParser.g"
 	
 	vector<int> ivec;
 	TorsionStamp* currTorsionStamp = static_cast<TorsionStamp*>(blockStack.top());
 	
-#line 1338 "MDTreeParser.cpp"
+#line 1420 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -1348,17 +1430,17 @@ void MDTreeParser::torsionstatement(ANTLR_USE_NAMESPAC
 		}
 		case MEMBERS:
 		{
-			ANTLR_USE_NAMESPACE(antlr)RefAST __t61 = _t;
-			ANTLR_USE_NAMESPACE(antlr)RefAST tmp36_AST_in = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST __t62 = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp40_AST_in = _t;
 			match(_t,MEMBERS);
 			_t = _t->getFirstChild();
 			ivec=inttuple(_t);
 			_t = _retTree;
-			_t = __t61;
+			_t = __t62;
 			_t = _t->getNextSibling();
-#line 191 "MDTreeParser.g"
+#line 198 "MDTreeParser.g"
 			currTorsionStamp->setMembers(ivec);
-#line 1362 "MDTreeParser.cpp"
+#line 1444 "MDTreeParser.cpp"
 			break;
 		}
 		default:
@@ -1377,12 +1459,12 @@ void MDTreeParser::inversionstatement(ANTLR_USE_NAMESP
 
 void MDTreeParser::inversionstatement(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST inversionstatement_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 204 "MDTreeParser.g"
+#line 211 "MDTreeParser.g"
 	
 	int icent;
 	InversionStamp* currInversionStamp = static_cast<InversionStamp*>(blockStack.top());
 	
-#line 1386 "MDTreeParser.cpp"
+#line 1468 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -1396,17 +1478,17 @@ void MDTreeParser::inversionstatement(ANTLR_USE_NAMESP
 		}
 		case CENTER:
 		{
-			ANTLR_USE_NAMESPACE(antlr)RefAST __t67 = _t;
-			ANTLR_USE_NAMESPACE(antlr)RefAST tmp37_AST_in = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST __t68 = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp41_AST_in = _t;
 			match(_t,CENTER);
 			_t = _t->getFirstChild();
 			icent=intConst(_t);
 			_t = _retTree;
-			_t = __t67;
+			_t = __t68;
 			_t = _t->getNextSibling();
-#line 210 "MDTreeParser.g"
+#line 217 "MDTreeParser.g"
 			currInversionStamp->setCenter(icent);
-#line 1410 "MDTreeParser.cpp"
+#line 1492 "MDTreeParser.cpp"
 			break;
 		}
 		default:
@@ -1425,12 +1507,12 @@ void MDTreeParser::rigidbodystatement(ANTLR_USE_NAMESP
 
 void MDTreeParser::rigidbodystatement(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST rigidbodystatement_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 227 "MDTreeParser.g"
+#line 234 "MDTreeParser.g"
 	
 	vector<int> ivec;
 	RigidBodyStamp* currRigidBodyStamp = static_cast<RigidBodyStamp*>(blockStack.top());
 	
-#line 1434 "MDTreeParser.cpp"
+#line 1516 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -1444,17 +1526,17 @@ void MDTreeParser::rigidbodystatement(ANTLR_USE_NAMESP
 		}
 		case MEMBERS:
 		{
-			ANTLR_USE_NAMESPACE(antlr)RefAST __t73 = _t;
-			ANTLR_USE_NAMESPACE(antlr)RefAST tmp38_AST_in = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST __t74 = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp42_AST_in = _t;
 			match(_t,MEMBERS);
 			_t = _t->getFirstChild();
 			ivec=inttuple(_t);
 			_t = _retTree;
-			_t = __t73;
+			_t = __t74;
 			_t = _t->getNextSibling();
-#line 233 "MDTreeParser.g"
+#line 240 "MDTreeParser.g"
 			currRigidBodyStamp->setMembers(ivec);
-#line 1458 "MDTreeParser.cpp"
+#line 1540 "MDTreeParser.cpp"
 			break;
 		}
 		default:
@@ -1473,12 +1555,12 @@ void MDTreeParser::cutoffgroupstatement(ANTLR_USE_NAME
 
 void MDTreeParser::cutoffgroupstatement(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
 	ANTLR_USE_NAMESPACE(antlr)RefAST cutoffgroupstatement_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-#line 246 "MDTreeParser.g"
+#line 253 "MDTreeParser.g"
 	
 	vector<int> ivec;
 	CutoffGroupStamp* currCutoffGroupStamp = static_cast<CutoffGroupStamp*>(blockStack.top());
 	
-#line 1482 "MDTreeParser.cpp"
+#line 1564 "MDTreeParser.cpp"
 	
 	try {      // for error handling
 		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
@@ -1492,17 +1574,17 @@ void MDTreeParser::cutoffgroupstatement(ANTLR_USE_NAME
 		}
 		case MEMBERS:
 		{
-			ANTLR_USE_NAMESPACE(antlr)RefAST __t79 = _t;
-			ANTLR_USE_NAMESPACE(antlr)RefAST tmp39_AST_in = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST __t80 = _t;
+			ANTLR_USE_NAMESPACE(antlr)RefAST tmp43_AST_in = _t;
 			match(_t,MEMBERS);
 			_t = _t->getFirstChild();
 			ivec=inttuple(_t);
 			_t = _retTree;
-			_t = __t79;
+			_t = __t80;
 			_t = _t->getNextSibling();
-#line 252 "MDTreeParser.g"
+#line 259 "MDTreeParser.g"
 			currCutoffGroupStamp->setMembers(ivec);
-#line 1506 "MDTreeParser.cpp"
+#line 1588 "MDTreeParser.cpp"
 			break;
 		}
 		default:
@@ -1534,55 +1616,6 @@ void MDTreeParser::fragmentstatement(ANTLR_USE_NAMESPA
 	_retTree = _t;
 }
 
-RealType  MDTreeParser::doubleNumber(ANTLR_USE_NAMESPACE(antlr)RefAST _t) {
-#line 292 "MDTreeParser.g"
-	RealType dval;
-#line 1541 "MDTreeParser.cpp"
-	ANTLR_USE_NAMESPACE(antlr)RefAST doubleNumber_AST_in = (_t == ANTLR_USE_NAMESPACE(antlr)RefAST(ASTNULL)) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-	ANTLR_USE_NAMESPACE(antlr)RefAST ic = ANTLR_USE_NAMESPACE(antlr)nullAST;
-	ANTLR_USE_NAMESPACE(antlr)RefAST fc = ANTLR_USE_NAMESPACE(antlr)nullAST;
-	
-	try {      // for error handling
-		if (_t == ANTLR_USE_NAMESPACE(antlr)nullAST )
-			_t = ASTNULL;
-		switch ( _t->getType()) {
-		case NUM_INT:
-		case NUM_LONG:
-		{
-			ic = (_t == ASTNULL) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-			intConst(_t);
-			_t = _retTree;
-#line 294 "MDTreeParser.g"
-			dval = lexi_cast<RealType>(ic->getText());
-#line 1558 "MDTreeParser.cpp"
-			break;
-		}
-		case NUM_FLOAT:
-		case NUM_DOUBLE:
-		{
-			fc = (_t == ASTNULL) ? ANTLR_USE_NAMESPACE(antlr)nullAST : _t;
-			floatConst(_t);
-			_t = _retTree;
-#line 295 "MDTreeParser.g"
-			dval = lexi_cast<RealType>(fc->getText());
-#line 1569 "MDTreeParser.cpp"
-			break;
-		}
-		default:
-		{
-			throw ANTLR_USE_NAMESPACE(antlr)NoViableAltException(_t);
-		}
-		}
-	}
-	catch (ANTLR_USE_NAMESPACE(antlr)RecognitionException& ex) {
-		reportError(ex);
-		if ( _t != ANTLR_USE_NAMESPACE(antlr)nullAST )
-			_t = _t->getNextSibling();
-	}
-	_retTree = _t;
-	return dval;
-}
-
 void MDTreeParser::initializeASTFactory( ANTLR_USE_NAMESPACE(antlr)ASTFactory& )
 {
 }