The eikonal reaction theory (ERT) proposed lately is a method of calculating one-neutron removal reactions at intermediate incident energies in which Coulomb breakup is treated accurately with the continuum discretized coupled-channels method. ERT is extended to two-neutron removal reactions. ERT reproduces measured oneand two-neutron removal cross sections for 6 He scattering on 12 C and 208 Pb targets at 240 MeV/nucleon and also on a 28 Si target at 52 MeV/nucleon. For the heavier target in which Coulomb breakup is important, ERT yields much better agreement with the measured cross sections than the Glauber model. PACS numbers: 24.10.Eq, 25.60.Gc, 25.70.De Introduction. Removal reactions are a quite useful tool for investigating structure of valence nucleons in weakly-bound nuclei such as one-and two-neutron halo nuclei. Spectroscopic factors and orbital angular momenta of valence nucleons in incident nuclei can be deduced from the removal reactions; see for example Ref. [1]. In particular, two-neutron removal reactions are crucial for analyzing two-neutron correlations between valence neutrons. The two-neutron removal reactions were investigated for light targets [2][3][4] with the Glauber model [5].The Glauber model is based on the eikonal and adiabatic approximations. The theoretical foundation of the model is shown in Ref. [6]. Once Coulomb breakup is taken into account in the Glauber model, the calculated removal cross sections diverge because of the adiabatic approximation. For this reason, the model has been applied to light targets. Lately a way of making Coulomb corrections to Glauber-model calculations was proposed [7,8]; the divergent dipole component of the eikonal Coulomb phase is replaced by that estimated by the first-order perturbation.The Coulomb problem is solved by the eikonal reaction theory (ERT) [17,18] in which Coulomb breakup is treated accurately with the continuum discretized coupledchannels method (CDCC) [19][20][21]. For one-proton andneutron removal cross sections of deuteron scattering at 200 MeV/nucleon, the Glauber-model results are found to be largely deviated from the ERT results for heavier targets [18]. In Ref. [17], ERT was applied to recently measured one-neutron removal cross sections for 31 Ne scattering from 12 C and 208 Pb targets at 230 MeV/nucleon [22]. Spectroscopic factors and asymptotic normalization coefficients for the 30 Ne + n bound system are consistently deduced from the measured cross sections for both the light and heavy targets. The analysis for both light and heavy thus makes it possible to determine spectroscopic factors and asymptotic normalization coefficients of valence neutrons definitely.Scattering of three-body projectiles such as 6 He can be described by the four-body model composed of three constituents of projectile and a target. Four-body CDCC [9-13] is a method of treating projectile breakup in the four-body scattering. Four-body CDCC was applied so far to elastic scattering and exclusive breakup reactions of 6 He from 12 C and 208 Pb targ...