I have a finite reflection (or Coxeter) group defined abstractly through the standard presentation
$$(s_i s_j)^{c_{ij}}=1$$
For each of its elements I want to find the number of reduced words equal to it. Is this possible in GAP or one of its packages?

Here is some GAP code using KBMAG to do this example. I won't try and explain in detail what it is doing, but please let me know if anything is unclear.
As I said, this approach will become too slow if you want to do larger examples with longer words. There is a better method for testing words in Coxeter groups for irreducibility but I am not aware of any GAP implementation.