The experimental study of relativistic heavy-ion collisions provides a strong evidence of the formation of hot and dense Quark-Gluon Plasma (QGP) state of matter at the center of the collision zone [1]}, [2]}, [3]}, [4]}, [5]}. The relativistic hydrodynamic models present a satisfactory description of fireball evolution as the copious production and large harmonic flow of hadrons have been successfully explained by these frameworks [6]}, [7]}, [8]}, [9]}, [10]}, [11]}, [12]}, [13]}, [14]}, [15]}, [16]}, [17]}, [18]}, [19]}. The past studies have shown that nuclear deformation can be investigated from the anisotropic flow of produced hadrons at intermediate and high collision energies  [20]}, [21]}, [22]}, [23]}. The nuclear collisions at such energies provide us a unique way to investigate the intrinsic substructure of nuclei as the time during which the participating nuclei pass through each other is too small for any slower nuclear excitation to occur and thus one can obtain a projection of overlying nuclear distributions on the fireball.
