The problems of characterizing quadratic stability and computing an estimation of the domain of attraction for saturated hybrid systems are addressed. Hybrid systems presenting saturations and nested saturations on signals involved in both the continuous-time and the discrete-time dynamics are considered. Geometrical characterizations of local and global quadratic stability are provided. Computation oriented conditions for quadratic stability are given in form of convex constraints.