I believe the issue is in your feature collection. I can't access it as you don't have it shared, so can't investigate what is going on there.
But I have tested the script with a random geometry placed on the map, and had no issues.
Check that your feature collection doesn't have too many features to it.
Try using raster calculator to calculate a single summarizing statistic based on all 15 bands.
You will have to have each band separated as a one-band raster dataset or as a raster layer in the map.
I have never used this this, but the Band Collection Statistics tools sounds useful.
MODIS assumes the Earth is a sphere with a R = 6.3781e6 m radius. The world is subdivided into 18 vertical tiles and 36 horizontal tiles, each tile has 2400 cells of 500m x 500m.
We can code this as:
CELLS = 2400
VERTICAL_TILES = 18
HORIZONTAL_TILES = 36
EARTH_RADIUS = 6371007.181
EARTH_WIDTH = 2 * math.pi * EARTH_RADIUS
TILE_WIDTH = ...
I recommend resampling your raster data to match the higher resolution dataset (i.e. smaller pixel dimensions) using gdal.
First run the gdalinfo utility on both files to determine the resolution
then choose the smaller set of values (call these xres and yres) of the two and use gdalwarp to ...