SpikeStats.cpp

Go to the documentation of this file.

#include "CMtuFormat.h"
#include "CDespike.h"
#include "CSpikeStats.h"
#include "CTsSpectrum.h"
#include <iostream>
#include <fstream>
#include <string>

using namespace std;
int main(void)
{
        CMtuFormat TestData;
        CSpikeStats SpikeStatsX;
        CSpikeStats SpikeStatsY;
        ofstream spikestats,statsum;
        titspikeindex spikeposout;
        string infilename;
        string temp;
        double width;
        
        do
        {
        cout << "Inputfile: ";
        cin >> infilename;
        }
        while (TestData.GetData(infilename.c_str()));
        temp = infilename + "_spikestats.dat";
        spikestats.open(temp.c_str());
        temp = infilename + "_statsum.dat";
        statsum.open(temp.c_str());
        SpikeStatsX.nValueBins = 1000;
        SpikeStatsX.nDerivBins = 1000;
        SpikeStatsX.samplerate = TestData.samplerate;
        
        
        SpikeStatsX.Data = &(TestData.Ex.data);
        SpikeStatsX.Time = &(TestData.t);
        SpikeStatsX.AnalyseSpikes();    
        
        vector<int>::iterator occ = SpikeStatsX.AbsoluteValueOcc.begin();
        vector<double>::iterator bin = SpikeStatsX.AbsoluteValueBin.begin();
        double sum = 0;
        width = (*(bin+1) - *bin);
        while (occ != SpikeStatsX.AbsoluteValueOcc.end() )
        {
                spikestats.precision(15);
                sum += *occ/SpikeStatsX.AbsoluteTotalOcc;
                //statsum << *bin/(SpikeStatsX.absolutemax - SpikeStatsX.absolutemin) << " " << sum << endl;
                statsum << *bin << " " << sum << endl;
                //spikestats << *bin/(SpikeStatsX.absolutemax - SpikeStatsX.absolutemin) << " " << *occ/(SpikeStatsX.AbsoluteTotalOcc*
                //      width/(SpikeStatsX.absolutemax - SpikeStatsX.absolutemin)) << endl;
                spikestats << *bin << " " << *occ/(SpikeStatsX.AbsoluteTotalOcc*
                        width) << endl;
                occ++;
                bin++;
        }
        spikestats.close();
        statsum.close();
        cout << " done" << endl << flush;
}