double precision onshellcutoff
      parameter(onshellcutoff=1d-8)