Dear Material Studio users,

I have a system that includes several 6-rings. Following MD simulations, I aimed to conduct post-processing analyses, one of which requires obtaining the normal vectors of each 6-ring plane versus time.

So, I made BestFitPlane out of all 6-rings (please see the attached figure) and then wrote the following code:

#!perl

use strict;
use Getopt::Long;
use MaterialsScript qw(:all);

my \$doc = \$Documents{"new0.xtd"};
my \$std = Documents->New("NormaleVectors.std");
\$std->ColumnHeading(0) = "Nx";
\$std->ColumnHeading(1) = "Ny";
\$std->ColumnHeading(2) = "Nz"
my \$bestFitPlane = \$doc->DisplayRange->BestFitPlane;
my \$planes = \$bestFitPlane->BestFitPlane;
for (my \$frame=1; \$frame<=\$doc->Trajectory->NumFrames; ++\$frame){
\$doc->Trajectory->CurrentFrame = \$frame;
my \$BFPCounter = 0;
my \$i=0
foreach my \$plane (@{\$doc->DisplayRange->BestFitPlane}){{++\$BFPCounter}  
\$std->Cell(\$i,0+(\$frame-1)*3) = \$plane->UnitNormalVector->X;
\$std->Cell(\$i,1+(\$frame-1)*3) = \$plane->UnitNormalVector->Y;
\$std->Cell(\$i,2+(\$frame-1)*3) = \$plane->UnitNormalVector->Z;
\$i++;} 
}

However, upon running the code, I got this error:

"There is more than one object in the "BestFitPlane"s collection (function/property "BestFitPlane") at -e line 13."

Would you be able to kindly advise me on what I have missed in the code? Thank you, SJ