def test():
data = FitxerRugRealitzacio("rug.dat")
mitjana_ = data.mitjana()
bd = BurstDetector(data)
bd.set_threshold(mitjana_)
bursts = bd.bursts()
bi = BurstIntegrator(bursts)
print bi
def test():
data = FitxerEneRealitzacio("./ene.dat")
print data
mitjana_ = data.mitjana()
bd = BurstDetector(data)
bd.set_threshold(mitjana_)
bursts = bd.bursts()
print "mitjana:", mitjana_
bi = BurstIntegrator(bursts)
print bi
tl = TimeLapses(bursts)
print tl
def main(directori):
stressmf = os.path.join(directori, "stressm.dat")
data = RealitzacioDStressm(stressmf)
threshold = data.threshold()
bd = BurstDetector(data)
bd.set_threshold(threshold)
bursts = bd.bursts()
bi = BurstIntegrator(bursts)
bi.integrate_with_x_from(column=3)
outfile = os.path.join(directori, "burst_dstressm.dat")
outfile = open(outfile, 'w')
print >> outfile, bi
print >> sys.stderr, "printing", outfile.name
outfile.close()
def main(directori, threshold=None, outfile="burst.dat"):
fitxer = os.path.join(directori, "rug.dat")
data = FitxerRugRealitzacio(fitxer)
# amb aquest donava bé!!!
#if threshold is None:
#threshold = data.mitjana() / 10. # el / 10. el copio de energy
threshold = data.threshold()
bd = BurstDetector(data)
bd.set_threshold(threshold)
bursts = bd.bursts()
bi = BurstIntegrator(bursts)
bi.integrate_with_x_from(column=3)
outfile = os.path.join(directori, outfile)
assert_file_is_empty (outfile)
outfile = open(outfile,"w")
print >> outfile, bi
print >> sys.stderr, "printing", outfile.name
outfile.close()
def main(directori, threshold=None, outfile="burst_dW.dat", outfile_time_lapses="time_lapses.dat"):
fitxer = os.path.join(directori, "eneNow.dat")
data = FitxerEneRealitzacio(fitxer)
# if threshold is None:
# threshold = data.mitjana()/10.
threshold = data.threshold()
bd = BurstDetector(data)
bd.set_threshold(threshold)
bursts = bd.bursts()
bi = BurstIntegrator(bursts)
bi.integrate_with_x_from(column=3)
bm = BurstMaxims(bursts)
outfile = os.path.join(directori, outfile)
assert_file_is_empty(outfile)
outfile = open(outfile, "w")
print >> outfile, bi
print >> sys.stderr, "printing", outfile.name
outfile.close()
outfile = os.path.join(directori, "burst_maxims.dat")
assert_file_is_empty(outfile)
outfile = open(outfile, "w")
print >> outfile, bm
print >> sys.stderr, "printing", outfile.name
outfile.close()
tl = TimeLapses(bursts)
outfile = os.path.join(directori, outfile_time_lapses)
assert_file_is_empty(outfile)
outfile = open(outfile, "w")
print >> outfile, tl
print >> sys.stderr, "printing", outfile.name
outfile.close()
