Suppose that ω 1 , . . . , ω N are positive real numbers and x is a complex number with positive real part. The multiple Barnes-Euler zeta function ζ E,N (s, x; ω) with parameter vector ω = (ω 1 , . . . , ω N ) is defined as a deformation of the Barnes multiple zeta function as followsIn this paper, based on the fermionic p-adic integral, we define the padic analogue of multiple Barnes-Euler zeta function ζ E,N (s, x; ω) which we denote by ζ p,E,N (s, x; ω). We prove several properties of ζ p,E,N (s, x; ω), including the convergent Laurent series expansion, the distribution formula, the difference equation, the reflection functional equation and the derivative formula. By computing the values of this kind of p-adic zeta function at nonpositive integers, we show that it interpolates the higher order Euler polynomials E N,n (x; ω) p-adically.Furthermore, we define the corresponding multiple p-adic Diamond-Euler Log Gamma function. We also show that the multiple p-adic Diamond-Euler Log Gamma function Log Γ D,E,N (x; ω) has an integral representation by the multiple fermionic p-adic integral, and it satisfies the distribution formula, the difference equation, the reflection functional equation, the derivative formula and also the Stirling's series expansions.