Separated group seed

Author: ecole41

validphys2/src/validphys/n3fit_data.py

@@ -110,6 +110,20 @@ def replica_luxseed(replica, luxseed):
     return replica_nnseed(replica, luxseed)
 
 
+def group_replica_mcseed(replica_mcseed, groups_dataset_inputs_loaded_cd_with_cuts):
+    names_for_salt = []
+    for loaded_cd in groups_dataset_inputs_loaded_cd_with_cuts:
+        if loaded_cd.legacy_names is None:
+            names_for_salt.append(loaded_cd.setname)
+        else:
+            names_for_salt.append(loaded_cd.legacy_names[0])
+    name_salt = "-".join(names_for_salt)
+
+    name_seed = int(hashlib.sha256(name_salt.encode()).hexdigest(), 16) % 10**8
+    res = name_seed + replica_mcseed
+    return res
+
+
 class _Masks(TupleComp):
     """Class holding the training validation mask for a group of datasets
     If the same group of dataset receives the same trvlseed then the mask

validphys2/src/validphys/pseudodata.py

@@ -122,7 +122,7 @@ def read_replica_pseudodata(fit, context_index, replica):
 
 def make_replica(
     groups_dataset_inputs_loaded_cd_with_cuts,
-    replica_mcseed,
+    group_replica_mcseed,
     dataset_inputs_sampling_covmat,
     sep_mult=False,
     genrep=True,
@@ -192,7 +192,7 @@ def make_replica(
         )
 
     # Set random seed to replica_mcseed - Would also like to change this seed for each group via 'groupname' (this can't yet be accessed here)
-    rng = np.random.default_rng(seed=replica_mcseed)
+    rng = np.random.default_rng(seed=group_replica_mcseed)
     # construct covmat
     covmat = dataset_inputs_sampling_covmat
     covmat_sqrt = sqrt_covmat(covmat)