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=1006; n<=1030; n++)
{

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

number tempsum = sum(temp)

summed = summed + tempsum
}

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

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

number output = summed - 25*(((bg2/30)+(bg1/30))/2)

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