OptMatrix

 
+	inline void GetTotalInteractionAsMatrixOpt(Interaction::SubType st,
+					    Matrix3** mat,
+					    energy* xx,energy* yy, energy* zz) const {
+	    *mat=mOptList[st].matrix;
+	    *xx=mOptList[st].xx;
+	    *yy=mOptList[st].yy;
+	    *zz=mOptList[st].zz;
+	}
+	inline void GetTotalInteractionAsTraceOpt(Interaction::SubType st, energy* trace) const {
+	    *trace=mOptList[st].trace;
+	}
+
     private:
         std::vector<Interaction*> mInter;
         Vector3l mPosition;
         std::string mLabel;
         long mElement;
         std::vector<long> mIsotopes;
+
+	struct OptRecord {
+	    Matrix3 matrix;
+	    energy xx;
+	    energy yy;
+	    energy zz;
+	    energy trace;
+	};
+	std::vector<OptRecord> mOptList;
     };