test that t0 sampling is different from exp

Author: comane

validphys2/src/validphys/tests/test_pythonmakereplica.py

@@ -14,6 +14,7 @@
 from validphys.api import API
 from validphys.tests.conftest import DATA
 from validphys.tests.test_covmats import CORR_DATA
+from validphys.tests.conftest import PDF
 
 SEED = 123456
 
@@ -81,6 +82,29 @@ def test_pseudodata_seeding(data_config, dataset_inputs, use_cuts):
     np.testing.assert_allclose(rep_1, rep_2)
 
 
+@pytest.mark.parametrize("dataset_inputs", [DATA, CORR_DATA])
+def test_pseudodata_t0_and_not(data_config, dataset_inputs):
+    """
+    Check that using t0-sampling produces different pseudodata to not using
+    t0-sampling.
+    """
+    config = dict(data_config)
+    config["dataset_inputs"] = dataset_inputs
+    config["use_cuts"] = "internal"
+    config["replica_mcseed"] = SEED
+    config["use_t0"] = True
+    config["t0pdfset"] = PDF
+
+    config2 = deepcopy(config)
+    config2["use_t0_sampling"] = True
+
+    rep_1 = API.make_replica(**config)
+    rep_2 = API.make_replica(**config2)
+
+    # Check that the two replicas are different
+    assert not np.allclose(rep_1, rep_2)
+
+
 @pytest.mark.parametrize("use_cuts", ["nocuts", "internal"])
 @pytest.mark.parametrize("dataset_inputs", [DATA, CORR_DATA, SINGLE_SYS_DATASETS])
 def test_pseudodata_has_correct_ndata(data_config, dataset_inputs, use_cuts):