Herwig

examples/FermionProduction.yaml → examples/ChargedLeptonProduction.yaml

@@ -4,6 +4,7 @@ Generators:
   - Madgraph
   - KKMC
   - Pythia
+  - Herwig
 
 OutputFormat: edm4hep
 OutDir: Run-Cards
@@ -18,10 +19,6 @@ Processes:
   Muon91.2:
      Final: [13, -13]
 
-  MuonNeutrino91.2:
-     Final: [14, -14]
-     ISRMode: 1
-
   Tau91.2:
      Final: [15, -15]
      ISRMode: 1
@@ -31,10 +28,6 @@ Processes:
      SqrtS: 350
      ISRMode: 1
 
-  MuonNeutrino350:
-     Final: [14, -14]
-     SqrtS: 350
-
   Tau350:
      Final: [15, -15]
      SqrtS: 350

examples/NeutrinoProduction.yaml

@@ -0,0 +1,58 @@
+Generators:
+  - Sherpa
+  - Whizard
+  - Madgraph
+  - KKMC
+  - Pythia
+
+OutputFormat: edm4hep
+OutDir: Run-Cards
+Events: 10000
+EventMode: unweighted
+
+SqrtS: 91.2
+Model: SM
+ISRMode: 0
+
+Processes:
+
+  MuonNeutrino91.2:
+     Final: [14, -14]
+     ISRMode: 1
+
+  MuonNeutrino350:
+     Final: [14, -14]
+     SqrtS: 350
+
+Selectors:
+  Process:
+    Tau350:
+      PT:
+        Max: 175
+        Min: 20
+        Flavour: [-15,15]
+    Muon350:
+      PT:
+        Max: 175
+        Min: 20
+        Flavour: [-13,13]
+    Muon91.2:
+      PT:
+        Max: 45.6
+        Min: 20
+        Flavour: [-13,13]
+
+
+Sherpa:
+  Run:
+    EW_SCHEME: 3
+
+
+ParticleData:
+  23:
+    mass: 91.1876
+    width: 2.4952
+
+Analysis:
+  Tools: [key4HEP]
+

examples/ZHProduction.yaml

@@ -1,4 +1,5 @@
 Generators:
+  - Herwig
   - Sherpa
   - Whizard
   - Madgraph

k4GeneratorsConfig/src/eventGenerationCollections2Root.cxx

@@ -28,7 +28,7 @@ void k4GeneratorsConfig::eventGenerationCollections2Root::Init() {
   TGaxis::SetExponentOffset(-0.08, -0.12, "y");
 
   // define the generator colors (as offset)
-  std::vector<std::string> genNames = {"Madgraph", "Sherpa", "Whizard", "KKMC", "Pythia", "Babayaga"};
+  std::vector<std::string> genNames = {"Madgraph", "Sherpa", "Whizard", "KKMC", "Pythia", "Babayaga", "Herwig"};
   for (unsigned int i = 0; i < genNames.size(); i++) {
     m_generatorColorOffset[genNames[i]] = i;
   }
@@ -464,7 +464,7 @@ void k4GeneratorsConfig::eventGenerationCollections2Root::writeCrossSectionFigur
     mgRMS->GetYaxis()->SetLabelSize(0.1);
 
     // generate a name and write a png
-    label << m_dirname << "/" << m_processesList[iProc] << "wRMS.png";
+    label << m_dirname << "/" << m_processesList[iProc] << "wRMS.pdf";
     c1->Print(label.str().c_str());
     label.clear();
     label.str("");
@@ -503,7 +503,7 @@ void k4GeneratorsConfig::eventGenerationCollections2Root::writeCrossSectionFigur
     mgDelta->GetYaxis()->SetLabelSize(0.1);
 
     // generate a name and write a png
-    label << m_dirname << "/" << m_processesList[iProc] << "wDelta.png";
+    label << m_dirname << "/" << m_processesList[iProc] << "wDelta.pdf";
     c1->Print(label.str().c_str());
     label.clear();
     label.str("");
@@ -571,7 +571,7 @@ void k4GeneratorsConfig::eventGenerationCollections2Root::writeAnalysisHistosFig
       // now build the legend:
       TLegend* topLegend = new TLegend(0.1, 0.0, 0.4, 0.4);
       std::ostringstream title;
-      title << m_procSqrtsList[proc].first << " #sqrt{s} =" << m_procSqrtsList[proc].second << "GeV";
+      title << m_procSqrtsList[proc].first << " #sqrt{s} =" << m_procSqrtsList[proc].second << " GeV";
       topLegend->SetHeader(title.str().c_str(), "C");
       // we need to get the histo from the top
       TList* topPadPrimitives = topPad->GetListOfPrimitives();
@@ -626,7 +626,7 @@ void k4GeneratorsConfig::eventGenerationCollections2Root::writeAnalysisHistosFig
       // done, save the canvas
       label << m_dirname << "/" << m_procSqrtsList[proc].first
             << (unsigned int)(m_procSqrtsList[proc].second * m_EnergyUnitCnv) << m_cnvAnalysisHistosNames[proc][ihisto]
-            << ".png";
+            << ".pdf";
       m_cnvAnalysisHistos[proc][ihisto]->Print(label.str().c_str());
       label.clear();
       label.str("");

python/Generators/Babayaga.py

@@ -139,6 +139,10 @@ def getParticleOperator(self, pdg, prop):
         # not used by Babayaga
         pass
 
+    def getParticlePropertyUnit(self):
+        # not used
+        pass
+
     def getParticleProperty(self, attr):
         # not used by Babayaga
         pass

python/Generators/GeneratorBase.py

@@ -323,7 +323,7 @@ def resetGeneratorDatacard(self):
     def getGeneratorCommand(self,key,value):
         pass
 
-    def addOption2GeneratorDatacard(self,key,value,replace=True):
+    def addOption2GeneratorDatacard(self,key,value,unit="",replace=True):
         # check if the key is already defined in the datacard, then we take the last one (TBC):
         if replace is True:
             if key in self.__datacardContent:
@@ -440,14 +440,19 @@ def prepareParticles(self,add2Datacard=True, writeParticleHeader=False):
             for attr in particle:
                 value = particle[attr]
                 prop  = self.getParticleProperty(attr)
+                unit  = self.getParticlePropertyUnit()
                 # writing out
                 if prop is not None:
                     command = self.getParticleOperator(pdg,prop)
                     if add2Datacard is True:
-                        self.addOption2GeneratorDatacard(command, value,replace=False)
+                        self.addOption2GeneratorDatacard(command, value, unit, replace=False)
                     else:
                         self.replaceOptionInGeneratorDatacard(command,value)
 
+    @abstractmethod
+    def getParticlePropertyUnit(self):
+        pass
+
     @abstractmethod
     def getParticleProperty(self, attr):
         pass