In this paper we present a new description about statistical probability density distributions (pdfs) of Polar Mesospheric Clouds (PMC) and noctilucent clouds (NLC). The analysis is based on observations of maximum backscatter, ice mass density, ice particle radius, and number density of ice particles measured by the ALOMAR RMR-lidar for all NLC seasons from 2002 to 2016. From this data set we derive a new class of pdfs that describe the statistics of PMC/NLC events which is different from previously statistical methods using the approach of an exponential distribution commonly named g-distribution. The new analysis describes successfully the probability statistic of ALOMAR lidar data. It turns out that the former g-function description is a special case of our new approach. In general the new statistical function can be applied to many kinds of different PMC parameters, e.g. maximum backscatter, integrated backscatter, ice mass density, ice water content, ice particle radius, ice particle number density or albedo measured by satellites. As a main advantage the new method allows to connect different observational PMC distributions of lidar, and satellite data, and also to compare with distributions from ice model studies. In particular, the statistical distributions of different ice parameters can be compared with each other on the basis of a common assessment that facilitate, for example, trend analysis of PMC/NLC.