Add HH load area profile generator

src/egon/data/datasets/hh_demand_profiles.py

@@ -112,6 +112,9 @@
 Base = declarative_base()
 
 
+# Set random_seed as long as no global_seed exists (see #351).
+RANDOM_SEED = 42
+
 # Define mapping of zensus household categories to eurostat categories
 # - Adults living in househould type
 # - number of kids not included even if in housholdtype name
@@ -499,12 +502,12 @@ def create_missing_zensus_data(
         if difference > 0:
             # add to any row
             split = split.round()
-            random_row = split.sample()
+            random_row = split.sample(random_state=RANDOM_SEED)
             split[random_row.index] = random_row + difference
         elif difference < 0:
             # subtract only from rows > 0
             split = split.round()
-            random_row = split[split > 0].sample()
+            random_row = split[split > 0].sample(random_state=RANDOM_SEED)
             split[random_row.index] = random_row + difference
         else:
             split = split.round()
@@ -923,6 +926,7 @@ def get_cell_demand_profile_ids(df_cell, pool_size):
     # np.random.default_rng().integers(low=0, high=pool_size[hh_type], size=sq) instead of random.sample
     # use random.choices() if with replacement
     # list of sample ids per hh_type in cell
+    random.seed(RANDOM_SEED)
     cell_profile_ids = [
         (hh_type, random.sample(range(pool_size[hh_type]), k=sq))
         if pool_size[hh_type] >= sq