Chronicle/TimeLine/0709FirstPlasma/Shot17withPhotodiode/rsrc/basicdiagn/PlasmaDetect.pl

#!/bin/perl

do "parameters.pl";


# Open the file for reading.
open DATA, "Iplasma_dp.csv" or die "can't open file $!";
my @data = <DATA>;
close (DATA);
system("echo 0 > Plasma;rm -f PlasmaStart PlasmaEnd PlasmaStartAdvanced PlasmaEndDelayed");

@Ipl1=0;@Ipl2=0;@Ipl3=0;@Ipl4=0;



for ($i = 1; $i < $LineCount-500; ++$i)  
{
    @row = split(';',$data[$i]);$Ipl1=$row[1];
    @row = split(';',$data[$i+100]);$TimeStart=$row[0];$Ipl2=$row[1];
    @row = split(';',$data[$i+150]);$Ipl3=$row[1];
    @row = split(';',$data[$i+200]);$Ipl4=$row[1];
    if ( $Ipl2 > $Ipl1+$PlCurrentRise && $Ipl3 > $Ipl1+$PlCurrentRise && $Ipl4 > $Ipl1+$PlCurrentRise ){ $StartRow=$i+100;system ("echo $TimeStart > PlasmaStart;echo $TimeStart*9/10 > PlasmaStartAdvanced;echo $StartRow > PlasmaStartRow");$i=$LineCount};
}

for ($i = $LineCount-500; $i > 1 ; --$i)  
{
    @row = split(';',$data[$i]);$Ipl1=$row[1];
    @row = split(';',$data[$i-100]);$TimeEnd=$row[0];$Ipl2=$row[1];
    @row = split(';',$data[$i-150]);$Ipl3=$row[1];
    @row = split(';',$data[$i-200]);$Ipl4=$row[1];
    if ( $Ipl2 > $Ipl1+$PlCurrentRise && $Ipl3 > $Ipl1+$PlCurrentRise && $Ipl4 > $Ipl1+$PlCurrentRise ){$EndRow=$i-100;$dT=$TimeEnd-$TimeStart; system ("echo $TimeEnd > PlasmaEnd;echo $TimeEnd*11/10 > PlasmaEndDelayed;echo $dT > PlasmaTimeLength;echo $EndRow > PlasmaEndRow;echo 1 > Plasma");$i=1};
}