pdl> $photon_list = new PDL::RandVar->sample(500)->reshape(2,250)*10
pdl> $histogram = zeroes(10,10)
pdl> $histogram->range($photon_list,1)++;  #not what you wanted