#!/usr/bin/perl

use lib 'd:/Programs/Perl/lib';

use strict;

use JFile;
use Sci::Ellipsometry;

my $InFile = "SI111025.spe";
my $el = new Ellipsometry;
if(!$el->Read($InFile)) {
	print "Error: Can not read [$InFile].\n";
	exit;
}

my $FileType    = $el->FileType();
my $Title       = $el->Title();
my $Angle       = $el->IncidentAngle();
my $FileType    = $el->FileType();
my $pLabelArray = $el->LabelArray();
my $pDataArray  = $el->DataArray();
my $nDataArray  = $el->nDataArray();

print "Title: $Title\n";
print "Angle: $Angle\n";
print "nDataArray: $nDataArray\n";
my $nData = $el->nData();
print "nData: $nData\n";
for(my $i = 0 ; $i < $nDataArray ; $i++) {
	print "$pLabelArray->[$i]\t";
}
print "\n";

for(my $i = 0 ; $i < $nData ; $i++) {
	for(my $j = 0 ; $j < $nDataArray ; $j++) {
		print "$pDataArray->[$j]->[$i]\t";
	}
	print "\n";
}

print "\n";
#$el->SetIncidentAngle(70);
print "Angle: $Angle\n";
$nData = 10;
my $E     = $el->GetData('eV');
my $Psi   = $el->GetData('Psi');
my $Delta = $el->GetData('Delta');
my $e1    = $el->GetData('e1');
my $e2    = $el->GetData('e2');
print "E\tPsi\tDelta\te1\te2\n";
for(my $i = 0 ; $i < $nData ; $i++) {
#	my ($n, $k)     = $el->PsiDeltaToNK($Psi->[$i], $Delta->[$i]);
	my ($ee1, $ee2) = $el->PsiDeltaToEps($Psi->[$i], $Delta->[$i]);
	my ($n, $k)     = $el->EpsToNK($ee1, $ee2);
	print "$E->[$i]\t$Psi->[$i]\t$Delta->[$i]\t$e1->[$i]\t$e2->[$i]\t$n\t$k\t$ee1\t$ee2\n";
#	($ee1, $ee2) = $el->NKToEps($n, $k);
#	print "e=$ee1\t$ee2\n";
}