Merging full workflow into Tanay's HiggsDNA

Don't merge this yet! Let's do some testing first before we merge to master

This should work out of the box (fingers crossed). This will have the full h->aa->4gamma workflow included in the base class. This should also not affect any corrections that you are using currently. That still needs to be tested. There may be bugs that I haven't caught yet, so I need to do some testing... I also included some nice helper scripts for you in the scripts folder.

higgs_dna/workflows/base.py

@@ -533,10 +533,10 @@ class HggBaseProcessor(processor.ProcessorABC):  # type: ignore
 
                 # Append increased statistics to original array
                 tmp_photons = awkward.concatenate((tmp_photons, artificial_stats_photons))
+        """
 
         else:
             tmp_photons = events.Photon
-        """
 
         original_photons = tmp_photons
         #logger.info("original photons length %i"%len(original_photons)) # kt
@@ -620,11 +620,6 @@ class HggBaseProcessor(processor.ProcessorABC):  # type: ignore
         variations.append('nominal')
         logger.debug(f"[systematics variations] {variations}")
 
-        if "EvtMix" in dataset_name:
-            assert len(events.Photon[awkward.num(events.Photon, axis=1) >= 2]) * 10 == len(photons_dct["nominal"][awkward.num(photons_dct["nominal"], axis=1) >= 2])
-        elif "EvtMix" not in dataset_name:
-            assert len(events.Photon[awkward.num(events.Photon, axis=1) >= 2]) == len(photons_dct["nominal"][awkward.num(photons_dct["nominal"], axis=1) >= 2])
-
         for variation in variations:
             photons, jets = photons_dct["nominal"], events.Jet
             if variation == "nominal":