def METHOD_NAME(self):
    """
    Measures n-shots keeping all the probabilities fixed
    """
    n = self.N_shots()
    if self.T1_sigma() != 0:
        T1 = np.random.normal(self.T1(), self.T1_sigma(), 1)[0]
    else:
        T1 = self.T1()
    p_relax = (1-np.exp(-self.tau_d()/T1))
    state = self.state()
    P_RB = self.P_RB()
    F_discr = self.F_discr()
    measured_shots = np.empty(n)
    for i in range(n):
        # Relaxation
        if state == -1:
            if np.random.rand() < p_relax:
                state *= -1
        if np.random.rand() < P_RB:
            state *= -1
        # Readout
        if np.random.rand() > F_discr:
            measured_shots[i] = -1*state
        else:
            measured_shots[i] = state
    return measured_shots