In [ ]:

    

import scipy.io
import numpy
import os

percorsoIN = "/home/protoss/Documenti/TESI/DATI/candCoinc.mat"
candy = scipy.io.loadmat(percorsoIN)

candy1 = candy['candy'][0,0]['cand1']
candy2 = candy['candy'][0,0]['cand2']

frequenze = candy1[0]

critRat1 = candy1[5]
critRat2 = candy2[5]


    

In [ ]:

    

from matplotlib import pyplot
%matplotlib notebook


figura = pyplot.figure()
scat = figura.add_subplot(111)
figura, scat = pyplot.subplots(figsize=(10, 7))
scat.scatter(frequenze,critRat1, s = 5, label='H')
scat.scatter(frequenze,critRat2, s = 5, label = 'L')
pyplot.semilogx()
pyplot.legend(loc='upper left');

pyplot.show()


    

In [ ]:

    

filtro1 = numpy.where(critRat1 > 4)
filtro2 = numpy.where(critRat2 > 4)

critRatMedio = (critRat1+critRat2)/2
filtromedio = numpy.where(critRatMedio > 4)

critRat1filtro1 = critRat1[filtro1]
critRat2filtro1 = critRat2[filtro1]
critRat1filtro2 = critRat1[filtro2]
critRat2filtro2 = critRat2[filtro2]

critRat1filtromedio = critRat1[filtromedio]
critRat2filtromedio = critRat2[filtromedio]

annullatore1 = numpy.where(critRat1 <= 4)
annullatore2 = numpy.where(critRat2 <= 4)
critRat1[annullatore1] = 0
critRat1[annullatore2] = 0
critRat2[annullatore1] = 0
critRat2[annullatore2] = 0



critRat1filtroboth = critRat1#[filtroboth]
critRat2filtroboth = critRat2#[filtroboth]


print(critRat1filtro1.size, critRat1filtro2.size, critRat1filtroboth.size, critRat2filtro1.size, critRat2filtro2.size, critRat2filtroboth.size)


    

In [ ]:

    

print(critRat1filtromedio.size, critRat2filtromedio.size)

from matplotlib import pyplot
%matplotlib notebook

pyplot.figure(figsize=(10, 8))
a = pyplot.scatter(numpy.arange(critRat1filtromedio.size),critRat1filtromedio, s = 10)
b = pyplot.scatter(numpy.arange(critRat2filtromedio.size),critRat2filtromedio, s = 10)
#pyplot.semilogx()

pyplot.show()


    

In [ ]:

    

from matplotlib import pyplot
%matplotlib notebook

pyplot.figure(figsize=(10, 8))
a = pyplot.scatter(numpy.arange(critRat1filtro1.size),critRat1filtro1, s = 10)
b = pyplot.scatter(numpy.arange(critRat2filtro1.size),critRat2filtro1, s = 10)
#pyplot.semilogx()

pyplot.show()


    

In [ ]:

    

from matplotlib import pyplot
%matplotlib notebook

pyplot.figure(figsize=(10, 8))
a = pyplot.scatter(numpy.arange(critRat1filtro2.size),critRat1filtro2, s = 10)
b = pyplot.scatter(numpy.arange(critRat2filtro2.size),critRat2filtro2, s = 10)
#pyplot.semilogx()

pyplot.show()


    

In [ ]:

    

from matplotlib import pyplot
%matplotlib notebook

pyplot.figure(figsize=(10, 8))
a = pyplot.scatter(numpy.arange(critRat1filtroboth.size),critRat1filtroboth, s = 10)
b = pyplot.scatter(numpy.arange(critRat2filtroboth.size),critRat2filtroboth, s = 10)
#pyplot.semilogx()

pyplot.show()


    

In [ ]: