Recursive algorithm for the generation of relativistic kinematics for collisions and decays with regularizations of sharp peaks

We describe an algorithm for the generation of relativistic kinematics for collision and decay processes with multiparticle final states. In the framework of this algorithm it is possible to generate different kinematics covering most of practically interesting cases. One gets a possibility to introduce different sets of integration variables. As a result different sets of kinematical singularities can be regularized. To smooth sharp peaks some regularization formulas and procedures are used covering most typical cases. The algorithm is realized in the package CompHEP created for automatic calculations of collision and decay processes.