  The study of the anomalous non-resonant $HH$ production is presented from the
point of view of effective field theory (EFT) parameterization. The searches
for $pp \rightarrow HH \rightarrow \gamma \gamma WW$, $pp \rightarrow HH
\rightarrow \gamma \gamma bb$ and $pp \rightarrow HH \rightarrow bbbb$
processes by ATLAS experiment at the Large Hadron Collider (LHC) are
reproduced. The selection efficiencies and exclusion limits on the $HH$
production cross sections for different Beyond Standard Model (BSM) EFT
benchmark models are obtained using the fast Monte-Carlo simulation of the
ATLAS detector.
