Image datacube
Image temp
number summed, bg1, bg2
number width, height, depth

GetOneLabeledImageWithPrompt("Please select the rotated datacube.","Spectrum Grabber","Datacube",datacube)

get3dsize(datacube,width,height,depth)

for (number n=255; n<=274; n++)
{

temp = slice2( datacube, 0, 0, n, 0, width, 1, 1, height, 1 )

number tempsum = sum(temp)

summed = summed + tempsum
}

for (number n=238; n<=247; n++)
{
temp = slice2(datacube,0,0,n,0,width,1,1,height,1)
number tempsum = sum(temp)
bg1= bg1+tempsum
}

for (number n=288; n<=317; n++)
{
temp = slice2(datacube,0,0,n,0,width,1,1,height,1)
number tempsum = sum(temp)
bg2= bg2+tempsum
}

number output = summed - 20*((bg1+bg2)/40)

result("Non-subtracted Au Ma counts:" + summed + "\n")
result("Background subtracted Au Ma counts:" + output + "\n")