#pragma rtGlobals=1		// Use modern global access method.
#include ":IFDL Procedures:IFDL"
#include ":IFDL Procedures:Apply filter"

Function Master()

     variable SingleDouble
     SingleDouble=1
     
    Prompt SingleDouble,"1: for single A-tag towed, 2: for double A-tag towed"
    DoPrompt "FOR TOWED SURVEY TO COUNT ANIMALS",SingleDouble
    if(V_Flag)
    Abort "The user pressed Cancel"
    endif

   if (SingleDouble==1)
       Data_loadingCSV()
       PulseInterval()
       DrawSingle()
       
   elseif  (SingleDouble==2)
       Data_loadingCSVB()
       PulseIntervalB()
       Data_loadingCSV()
       PulseInterval()
       DrawDouble()
    
    else
       print "unexpected towed type. Let me assume single type"
       Data_loadingCSV()
       PulseInterval()
       DrawSingle()
    endif
    
    Bunpu()
    
end

//-----------------------------------------------------------
//Data loading
Function Data_loadingCSV()
   string /G InputFileName
   
   print "*****************************"
   print "******                             ******"
   print "****** select FRONT data  ******"
   print "******                             ******"
   print "*****************************"
   LoadWave/G/D/L={0, 0, 0, 0, 0}

   rename wave0, Igtime
   rename wave1, SPL1
   rename wave2, SPL2
   rename wave3, td

   InputFileName=S_fileName
   //print InputFileName

   print "FRONT data loading completed"
   print ""
   
end

//-----------------------------------------------------------
//Data loading
Function Data_loadingCSVB()
   string /G InputFileName

   print "****************************"
   print "******                            ******"
   print "****** select REAR data  ******"
   print "******                            ******"
   print "****************************"
   
   LoadWave/G/D/L={0, 0, 0, 0, 0}

   rename wave0, IgtimeB
   rename wave1, SPL1B
   rename wave2, SPL2B
   rename wave3, tdB

   InputFileName=S_fileName
   //print InputFileName

   print "REAR data loading completed"
   print ""
   
end


//-----------------------------------------------------------
//create inter-pulse interval
Function PulseInterval()
wave Igtime
variable p, m, N, j, dlt, scale
string w

m=numpnts(Igtime) //total number of pulses

Make/D/O/N=(m) PlsInt
  
    j=0
    do
		PlsInt[j]=(Igtime[j+1]*100000-Igtime[j]*100000)/100
		j+=1
	while (j<m-1)
	

End

//-----------------------------------------------------------
//create inter-pulse interval
Function PulseIntervalB()
wave IgtimeB
variable p, m, N, j, dlt, scale
string w

m=numpnts(IgtimeB) //total number of pulses

Make/D/O/N=(m) PlsIntB
  
    j=0
    do
		PlsIntB[j]=(IgtimeB[j+1]*100000-IgtimeB[j]*100000)/100
		j+=1
	while (j<m-1)
	

End

//----------------------------------------------------------

Function DrawDouble()
    wave SPL1, SPL2, SPL1B, SPL2B
    //Sound pressure
    //SPL1=SPL1*0.0771
    //SPL2=SPL2*0.0771
    //SPL1B=SPL1B*0.0771
    //SPL2B=SPL2B*0.0771

    
    Display/T SPL1 vs Igtime
    ModifyGraph axisEnab(left)={0.8, 1}
    Label left "Sound Pressure (relative to Pa)\\Z12"
    ModifyGraph mode(SPL1)=2,lsize(SPL1)=1.5,rgb(SPL1)=(0,0,0)
    AppendToGraph/T SPL1B vs IgtimeB
    ModifyGraph mode(SPL1B)=2,lsize(SPL1B)=1.5
    //Time difference
    AppendToGraph/L=td/T td vs Igtime
    ModifyGraph axisEnab(td)={0.32,0.65}
    ModifyGraph mode(td)=2,rgb(td)=(0,0,0)
    ModifyGraph lsize(td)=1.5
    SetAxis td -450,450 
    Label td "time difference"
    AppendToGraph/L=td/T tdB vs IgtimeB
    ModifyGraph mode(tdB)=2
    ModifyGraph lsize(tdB)=1.5
    //Inter pulse interval
    AppendToGraph/L=PlsInt/T PlsInt vs Igtime
    ModifyGraph axisEnab(PlsInt)={0,0.3};DelayUpdate
    Label PlsInt "PI (ms)";DelayUpdate
    SetAxis PlsInt 1,300 
    ModifyGraph log(PlsInt)=1;DelayUpdate
    ModifyGraph mode(PlsInt)=2,lsize(PlsInt)=2,rgb(PlsInt)=(0,15872,65280)
    AppendToGraph/L=PlsInt/T PlsIntB vs IgtimeB
    ModifyGraph mode(PlsIntB)=2,lsize(PlsIntB)=1.5
 
    ShowInfo

end

//----------------------------------------------------------

Function DrawSingle()
    wave SPL1, SPL2, SPL1B, SPL2B
    //Sound pressure in Pa
    //SPL1=SPL1*0.0771
    //SPL2=SPL2*0.0771

    //Sound pressure
    Display/T SPL1 vs Igtime
    ModifyGraph axisEnab(left)={0.8, 1}
    Label left "Sound Pressure (relative to Pa)\\Z12"
    ModifyGraph mode(SPL1)=1
    //Time difference
    AppendToGraph/L=td/T td vs Igtime
    ModifyGraph axisEnab(td)={0.32,0.75}
    ModifyGraph mode(td)=2,rgb(td)=(0,0,0)
    ModifyGraph lsize(td)=2
    SetAxis td -450,450 
    Label td "time difference"
    //Inter pulse interval
    AppendToGraph/L=PlsInt/T PlsInt vs Igtime
    ModifyGraph axisEnab(PlsInt)={0,0.3};DelayUpdate
    Label PlsInt "PI (ms)";DelayUpdate
    SetAxis PlsInt 1,300 
    ModifyGraph log(PlsInt)=1;DelayUpdate
    ModifyGraph mode(PlsInt)=2,lsize(PlsInt)=2,rgb(PlsInt)=(0,15872,65280)
 
end

//-----------------------------------------------------------
//eliminate low intensity signals, surface reflections and unusual change of inter-pulse intervals

Function Bunpu()
wave SPL1, SPL2, td, PlsInt, NofP

Make/O/N=1 SPL1Hist
histogram /B={0,1,1000}  SPL1,  SPL1Hist
display SPL1Hist
Label left "\\Z12number of events"
ModifyGraph fSize(left)=12
ModifyGraph fSize=12
Label bottom "\\Z12sound pressure 1"

Make/O/N=1 SPL2Hist
histogram /B={0,1,1000}  SPL2,  SPL2Hist
display SPL2Hist
Label left "\\Z12number of events"
ModifyGraph fSize(left)=12
ModifyGraph fSize=12
Label bottom "\\Z12sound pressure 2"

Make/O/N=1 tdHist
histogram /B={-511,1,1025}  td,  tdHist
display tdHist
Label left "\\Z12number of events"
ModifyGraph fSize(left)=12
ModifyGraph fSize=12
ModifyGraph log(left)=1;DelayUpdate
Label bottom "\\Z12time difference(count)"

Make/O/N=1 PlsIntHist
histogram /B={0,1,200}  PlsInt,  PlsIntHist
display PlsIntHist
Label left "\\Z12number of events"
ModifyGraph fSize(left)=12
ModifyGraph fSize=12
Label bottom "\\Z12pulse interval (ms)"


end