* Concatenation of multiple taxonomy & albedo catalogs
** Merging of classes
** Priority of classes
** Cleaning included (e.g. S vs K in Eos)
** Split X into EMP using PDF = f( a, i, H )


* Loop on level of detail: Population, a, i  [open code for N dimension]
** Taxonomic completeness
** MPC completeness
** Compute Number, Masses, etc
** Save in [a,i,H] boxes for each taxonomic class
*** Maximum level of detail reached [used to mask unreliable information]
*** Number: MPC, Missing MPC, Taxo     
*** Mass
*** Surface & Volume
*** Completeness level (MPC + Taxo)


* Details of box processing
** For all the known objects: count them
** For all unknown, apply fraction of known


* Critical inputs
** Bins: a, i, H
** Ranges: a, i, H
** Albedo limits for EMP
** Minimum number and fraction of MPC to trigger analysis
** Minimum number and fraction of Taxo to trigger analysis
** Definition of populations


** Threshold in SNR & Value VS uncertainty

* Asteroid input catalog
** Identifier: # Name Designation
** Taxonomy: Scheme + Class + Source
** Mass + Diameter + Albedo
** Orbital elements
** Family membership

** WITH UNCERTAINTIES!!!!

* Acronyms
** PDF  = Probability Density Function