search

key4hep / k4GeneratorsConfig

0 stars 1 forks source link

New process settings #9

Closed apricePhy closed 3 months ago

apricePhy commented 3 months ago

Hi Dirk,

I just pushed a new branch "new-process-settings” to the repo.

In addition to moving the isr,bes,sqrts to the process section I have made some other changes.

We now have global settings and process settings. This means the user can set sqrts either globally, as before, where it will be same for every process. Or the can set it in the process itself. If both are set then the process will win. In principle, any global setting can now be altered on a process by process basis, but I have not checked them all.

The structure of the directory will now be /Generator/ProcName/Run_ProcName{_Physics}.dat where _Physics can represent {isr,bes,…}. For example, the following 

  Muon91.2:
     Final: [13, -13]
     Order: [2,0]

  Muon91.2:
     Final: [13, -13]
     Order: [2,0]
     ISRMode: 1

will generate a folder /Generator/ProcName/ with two runcards and two executables

I also introduced process-specific selectors. So now, we can have different phase-space cuts for different processes.

Finally I added the ability to set defaults e.g Initial is always [11,-11] unless changed.

Play around and let me know if you run into anything. Alan

dirkzerwas commented 3 months ago

In FermionPair yaml, there were some inconsistencies in the selector settings (PT<0.5*sqrts) and PT fixed for neutrinos <45GeV In Madgraph the "inputfile name" did not have an extension and some other things. I moved the "filename calculation" to the initialisation of the Madgraph object, The full chain runs now. Will look at the other generators now for the two fermion input.

dirkzerwas commented 3 months ago

Some problems with inconsistent input/output names in the execution step for whizard and Madgraph, looking into that, but almost there....

dirkzerwas commented 3 months ago

CI went through, but looking at the output:

Difermion_Polarized: Madgraph 3.847 +- 0.128913 pb

Muon350: Madgraph 1.72679 +- 0.00186216 pb Sherpa 1.75815 +- 0.00299617 pb

Muon91.2: Madgraph 1720.03 +- 0.457748 pb Sherpa 1719.1 +- 1.07764 pb

MuonNeutrino350: Whizard 0.231006 +- 2.67451e-05 pb Madgraph 0.23104 +- 2.8156e-05 pb Sherpa 0.230949 +- 5.30189e-06 pb

MuonNeutrino91.2: Madgraph 2837.15 +- 2.57407 pb Sherpa 2766.95 +- 2.41293 pb

Tau350: Sherpa 0.873691 +- 0.000275726 pb Madgraph 0.87366 +- 0.000245258 pb

Tau91.2: Madgraph 1394.86 +- 1.00714 pb Sherpa 1410.67 +- 1.28874 pb

Top1000: Madgraph 0.185935 +- 0.000201299 pb Sherpa 0.184721 +- 0.000179541 pb

Top350: Madgraph 0.093295 +- 5.28213e-05 pb Sherpa 0.0896861 +- 0.00011159 pb

ZH250: Sherpa 0.197514 +- 0.00230028 pb Madgraph 0.19851 +- 0.00054668 pb

ZH350: Sherpa 0.133193 +- 0.00169692 pb Madgraph 0.12669 +- 0.00045771 pb

ZHDecay250: Madgraph 1.72148e-06 +- 1.49455e-09 pb

ZHDecay350: Madgraph 9.33635e-07 +- 1.52704e-09 pb

Not all executions went through (mostly whizard) But I think we should merge

dirkzerwas commented 2 months ago

correcting the cuts logic (missing "and" in the whizard .sin cuts section improves somewhat, but still missing some processes)

Difermion_Polarized: Madgraph 3.847 +- 0.128913 pb

Muon350: Madgraph 1.72679 +- 0.00186216 pb Sherpa 1.7551 +- 0.00297811 pb

Muon91.2: Madgraph 1720.03 +- 0.457748 pb Whizard 1721.38 +- 2.26812 pb Sherpa 1720.85 +- 1.06796 pb

MuonNeutrino350: Whizard 0.231006 +- 2.67451e-05 pb Madgraph 0.23104 +- 2.8156e-05 pb Sherpa 0.230948 +- 5.29197e-06 pb

MuonNeutrino91.2: Madgraph 2837.15 +- 2.57407 pb Sherpa 2766.7 +- 2.40499 pb

Tau350: Sherpa 0.873726 +- 0.000275327 pb Madgraph 0.87366 +- 0.000245258 pb Whizard 0.873822 +- 0.000791967 pb

Tau91.2: Madgraph 1394.86 +- 1.00714 pb Sherpa 1410.67 +- 1.28874 pb

Top1000: Madgraph 0.185935 +- 0.000201299 pb Sherpa 0.184721 +- 0.000179541 pb

Top350: Madgraph 0.093295 +- 5.28213e-05 pb Sherpa 0.0896861 +- 0.00011159 pb

ZH250: Sherpa 0.197514 +- 0.00230028 pb Madgraph 0.19851 +- 0.00054668 pb

ZH350: Sherpa 0.133193 +- 0.00169692 pb Madgraph 0.12669 +- 0.00045771 pb

ZHDecay250: Whizard 0.238792 +- 0.000112203 pb Madgraph 1.72148e-06 +- 1.49455e-09 pb

ZHDecay350: Whizard 0.128139 +- 6.38171e-05 pb Madgraph 9.33635e-07 +- 1.52704e-09 pb

dirkzerwas commented 2 months ago

Improvement when correcting the ISR implementation in Whizard.py, still some trouble with top

Difermion_Polarized: Madgraph 3.847 +- 0.128913 pb

Muon350: Madgraph 1.72679 +- 0.00186216 pb Whizard 1.78481 +- 0.00846396 pb Sherpa 1.7551 +- 0.00297811 pb

Muon91.2: Madgraph 1720.03 +- 0.457748 pb Whizard 1721.38 +- 2.26812 pb Sherpa 1720.85 +- 1.06796 pb

MuonNeutrino350: Whizard 0.231006 +- 2.67451e-05 pb Madgraph 0.23104 +- 2.8156e-05 pb Sherpa 0.230948 +- 5.29197e-06 pb

MuonNeutrino91.2: Madgraph 2837.15 +- 2.57407 pb Whizard 2839.97 +- 19.8289 pb Sherpa 2766.7 +- 2.40499 pb

Tau350: Sherpa 0.873726 +- 0.000275327 pb Madgraph 0.87366 +- 0.000245258 pb Whizard 0.873822 +- 0.000791967 pb

Tau91.2: Madgraph 1394.86 +- 1.00714 pb Sherpa 1410.67 +- 1.28874 pb Whizard 1443.28 +- 9.70288 pb

Top1000: Madgraph 0.185935 +- 0.000201299 pb Sherpa 0.184721 +- 0.000179541 pb

Top350: Madgraph 0.093295 +- 5.28213e-05 pb Sherpa 0.0896861 +- 0.00011159 pb

ZH250: Sherpa 0.197514 +- 0.00230028 pb Whizard 0.198518 +- 0.00069897 pb Madgraph 0.19851 +- 0.00054668 pb

ZH350: Whizard 0.126825 +- 0.000348292 pb Sherpa 0.133193 +- 0.00169692 pb Madgraph 0.12669 +- 0.00045771 pb

ZHDecay250: Whizard 0.238792 +- 0.000112203 pb Madgraph 1.72148e-06 +- 1.49455e-09 pb

ZHDecay350: Whizard 0.128139 +- 6.38171e-05 pb Madgraph 9.33635e-07 +- 1.52704e-09 pb

dirkzerwas commented 2 months ago

After doing special gymnastics for the top quark mass and width settings:

Difermion_Polarized: Madgraph 3.847 +- 0.128913 pb

Muon350: Madgraph 1.72679 +- 0.00186216 pb Whizard 1.78481 +- 0.00846396 pb Sherpa 1.7551 +- 0.00297811 pb

Muon91.2: Madgraph 1720.03 +- 0.457748 pb Whizard 1721.38 +- 2.26812 pb Sherpa 1720.85 +- 1.06796 pb

MuonNeutrino350: Whizard 0.231006 +- 2.67451e-05 pb Madgraph 0.23104 +- 2.8156e-05 pb Sherpa 0.230948 +- 5.29197e-06 pb

MuonNeutrino91.2: Madgraph 2837.15 +- 2.57407 pb Whizard 2839.97 +- 19.8289 pb Sherpa 2766.7 +- 2.40499 pb

Tau350: Sherpa 0.873726 +- 0.000275327 pb Madgraph 0.87366 +- 0.000245258 pb Whizard 0.873822 +- 0.000791967 pb

Tau91.2: Madgraph 1394.86 +- 1.00714 pb Sherpa 1410.67 +- 1.28874 pb Whizard 1443.28 +- 9.70288 pb

Top1000: Madgraph 0.185935 +- 0.000201299 pb Whizard 0.186629 +- 0.000826694 pb Sherpa 0.184721 +- 0.000179541 pb

Top350: Madgraph 0.093295 +- 5.28213e-05 pb Whizard 0.0930868 +- 0.000821829 pb Sherpa 0.0896861 +- 0.00011159 pb

ZH250: Sherpa 0.197514 +- 0.00230028 pb Whizard 0.198518 +- 0.00069897 pb Madgraph 0.19851 +- 0.00054668 pb

ZH350: Whizard 0.126825 +- 0.000348292 pb Sherpa 0.133193 +- 0.00169692 pb Madgraph 0.12669 +- 0.00045771 pb

ZHDecay250: Whizard 0.238792 +- 0.000112203 pb Madgraph 1.72148e-06 +- 1.49455e-09 pb

ZHDecay350: Whizard 0.128139 +- 6.38171e-05 pb Madgraph 9.33635e-07 +- 1.52704e-09 pb