#!/bin/perl

do "parameters.pl";

$Mirnov_S_Eff=0.0038;#[m]
$MirnovPosition=93;#[mm]
#$MirnovPosition=9,3;#[cm]


# Open the data-file for reading.
open DATA, "Nidatap.lvm" or die "cant open Nidatap.lvm";
@data=<DATA>;
close (DATA);



for ($i=0;$i<$LineCount;$i++)
{
	@row=split('	',$data[$i]);
	$bt[$i]=$row[2];
	$mc1[$i]=$row[5];
	$mc5[$i]=$row[6];
	$mc9[$i]=$row[7];
	$mc13[$i]=$row[8];
}


for ($i=495;$i<505;$i++)
{
	$mc1[$i]=0;
	$mc5[$i]=0;
	$mc9[$i]=0;
	$mc13[$i]=0;
}

$aver1=0;
$aver5=0;
$aver9=0;
$aver13=0;
for ($i=0;$i<$TriggerLine;$i++)
{
	$aver1+=$mc1[$i];
	$aver5+=$mc5[$i];
	$aver9+=$mc9[$i];
	$aver13+=$mc13[$i];
}
$aver1=$aver1/$TriggerLine;
$aver5=$aver5/$TriggerLine;
$aver9=$aver9/$TriggerLine;
$aver13=$aver13/$TriggerLine;


#Pro� je tu tady tohle??? Je to �patn�, �e�eno v n�sleduj�c� "for i"
#for ($i=0;$i<$LineCount;$i++)
#{
#	$mc1[$i]+=$Bt[$i]/(0.0038*166.736);
#}


$integr1=0;
$integr5=0;
$integr9=0;
$integr13=0;
for ($i=0;$i<$LineCount;$i++)
{
	$integr1+=($mc1[$i]-$aver1)+$bt[$i]/(205);#166.736);#upravena konstanta!! potom presneji zmerit
	$mc1[$i]=$integr1*$DeltaTs/$Mirnov_S_Eff;
	$integr5+=($mc5[$i]-$aver5)+$bt[$i]/76.68;
	$mc5[$i]=$integr5*$DeltaTs/$Mirnov_S_Eff;
	$integr9+=($mc9[$i]-$aver9)-$bt[$i]/(90);#105.598);#to same jako u mc1
	$mc9[$i]=$integr9*$DeltaTs/$Mirnov_S_Eff;
	$integr13+=($mc13[$i]-$aver13);
	$mc13[$i]=$integr13*$DeltaTs/$Mirnov_S_Eff;
}

$aver1=0;
$aver5=0;
$aver9=0;
$aver13=0;
for ($i=2949;$i<$LineCount;$i++)
{
	$aver1+=$mc1[$i];
	$aver5+=$mc5[$i];
	$aver9+=$mc9[$i];
	$aver13+=$mc13[$i];
}
$aver1=$aver1/50;
$aver5=$aver5/50;
$aver9=$aver9/50;
$aver13=$aver13/50;

for ($i=499;$i<2950;$i++)
{
	$mc1[$i]-=$aver1*($i-499)/2450;
	$mc5[$i]-=$aver5*($i-499)/2450;
	$mc9[$i]-=$aver9*($i-499)/2450;
	$mc13[$i]-=$aver13*($i-499)/2450;
}

for ($i=2950;$i<3000;$i++)
{
	$mc1[$i]-=$aver1;
	$mc5[$i]-=$aver5;
	$mc9[$i]-=$aver9;
	$mc3[$i]-=$aver13;
}




open CIL1, ">mc1.txt";
open CIL5, ">mc5.txt";
open CIL9, ">mc9.txt";
open CIL13, ">mc13.txt";
for ($i=0;$i<$LineCount;$i++)
{
	print CIL1 $mc1[$i]."\n";
	print CIL5 $mc5[$i]."\n";
	print CIL9 $mc9[$i]."\n";
	print CIL13 $mc13[$i]."\n";
}
close (CIL1);
close (CIL5);
close (CIL9);
close (CIL13);

open CIL, ">displaceR.txt";
open CIL1, ">displaceZ.txt";
for ($i=0;$i<$LineCount;$i++)
{
	$displaceR=($mc1[$i]-$mc9[$i])*$MirnovPosition/($mc1[$i]+$mc9[$i]);
	$displaceZ=($mc5[$i]-$mc13[$i])*$MirnovPosition/($mc5[$i]+$mc13[$i]);
	print CIL $displaceR."\n";
	print CIL1 $displaceZ."\n";
}
close (CIL);
close (CIL1);

