source: fact/tools/pyscripts/sandbox/dneise/rdietlic/periods2.py@ 13408

# Programm zur Bestimmung der Genauigkeit des konvergierten Algorithmus
#
# Remo Dietlicher
#
#

from optparse import OptionParser
import pyfact
from myhisto import *
from hist import *
import numpy as np
import numpy.random as rnd
from ROOT import *
from time import time

def periods2(Data, CellTime, h):

        NROI  = 1024                            # length of the DRS pipeline
        fsampling = 2.                          # sampling frequency
        freq = 250.                             # testfrequency
        P_nom = 1000./freq                      # nominal Period due to testfrequency


        nomCellTime = np.linspace(0., 1024., 1025)/fsampling


        # FUNCTION TO DETERMINE CROSSINGS

        def Crossing(Mean, rCellTime):
                TimeXing = np.zeros(NROI)
                MeanXing = np.zeros(NROI)
                NumXing = 0

                for i in range(NROI-1):
                        if ((Data[i] > Mean) & (Data[i+1] < Mean)):
                                MeanXing[NumXing] = i

                                FirstCell = rCellTime[i]
                                SecondCell = rCellTime[i+1]

                                TimeXing[NumXing] = FirstCell+(SecondCell-FirstCell)/(1.-Data[i+1]/(Data[i]))*(1.-Mean/(Data[i]))
                                NumXing += 1

                return MeanXing, TimeXing

        def CalculatePeriods(rCellTime, name):

                Period = np.zeros(126)

                MeanXing, TimeXing = Crossing(np.average(Data), rCellTime)

                for i in range(int(126)):
                        Period[i] = TimeXing[i+1] - TimeXing[i]

                for val in Period:
                        h.dict[str(name)].Fill(val)

                return

        CalculatePeriods(CellTime, "avperiods")
        CalculatePeriods(nomCellTime, "avperiods0")
        
        return