Abstract. The aim of this paper is to study the boundedness of different conical square functions that arise naturally from second order divergence form degenerate elliptic operators. More precisely, let L w = −w −1 div(w A ∇) where w ∈ A 2 and A is an n × n bounded, complex-valued, uniformly elliptic matrix. D. Cruz-Uribe and C. Rios solved the L 2 (w)-Kato square root problem obtaining that √ L w is equivalent to the gradient on L 2 (w). The same authors in collaboration with the second named author of this paper studied the L p (w)-boundedness of operators that are naturally associated with L w , such as the functional calculus, Riesz transforms, or vertical square functions. The theory developed admitted also weighted estimates (i.e., estimates in L p (vdw) for v ∈ A ∞ (w)), and in particular a class of "degeneracy" weights w was found in such a way that the classical L 2 -Kato problem can be solved. In this paper, continuing this line of research, and also that originated in some recent results by the second and third named authors of the current paper, we study the boundedness on L p (w) and on L p (vdw), with v ∈ A ∞ (w), of the conical square functions that one can construct using the heat or Poisson semigroup associated with L w . As a consequence of our methods, we find a class of degeneracy weights w for which L 2 -estimates for these conical square functions hold. This opens the door to the study of weighted and unweighted Hardy spaces and of boundary value problems associated with L w .