flaredetector: allow bigger delta to mean

This commit is contained in:
2024-04-25 15:50:28 +02:00
parent 31dfaebc3e
commit e835964857
+3 -2
View File
@@ -57,6 +57,7 @@ def getFlareRange(flux, peakIndex):
maxInd = np.inf maxInd = np.inf
maxDelta = 0.005 maxDelta = 0.005
maxDeltaToMean = 2*maxDelta
meanVal = 1 meanVal = 1
maxBacklog = 20 maxBacklog = 20
maxForwardlog = 100 maxForwardlog = 100
@@ -66,7 +67,7 @@ def getFlareRange(flux, peakIndex):
prevFlux = flux[peakIndex] prevFlux = flux[peakIndex]
for i in range(peakIndex-1, peakIndex - maxBacklog, -1): for i in range(peakIndex-1, peakIndex - maxBacklog, -1):
delta = prevFlux - flux[i] delta = prevFlux - flux[i]
if(delta < maxDelta and abs(flux[i] - meanVal) < maxDelta): if(delta < maxDelta and abs(flux[i] - meanVal) < maxDeltaToMean):
deltaCount += 1 deltaCount += 1
if(deltaCount > 3): if(deltaCount > 3):
minInd = i minInd = i
@@ -89,7 +90,7 @@ def getFlareRange(flux, peakIndex):
prevFlux = flux[peakIndex] prevFlux = flux[peakIndex]
for i in range(peakIndex+1, peakIndex + maxForwardlog, 1): for i in range(peakIndex+1, peakIndex + maxForwardlog, 1):
delta = prevFlux - flux[i] delta = prevFlux - flux[i]
if(delta < maxDelta and abs(flux[i] - meanVal) < maxDelta): if(delta < maxDelta and abs(flux[i] - meanVal) < maxDeltaToMean):
deltaCount += 1 deltaCount += 1
if(deltaCount > 3): if(deltaCount > 3):
maxInd = i maxInd = i