def testmatch(
    name name,
    v1 v1,
    v2 v2,
    tolerance tolerance =1e-3
)

h1 =  Histo1D(5, 0.0, 100.0);

h2 =  Histo1D(5, 0.0, 100.0);

def meanmatch =  testmatch("means",    h1.mean(),   h2.mean());

def stddevmatch =  testmatch("std devs", h1.stdDev(), h2.stdDev());

tuple ok =  (stddevmatch and meanmatch);

p1 =  Profile1D(5, 0.0, 100.0);

p2 =  Profile1D(5, 0.0, 100.0);