diff --git a/ECFAHiggsFactories/ttbar/ecms.dat b/ECFAHiggsFactories/ttbar/ecms.dat new file mode 100644 index 0000000..c4a9bcd --- /dev/null +++ b/ECFAHiggsFactories/ttbar/ecms.dat @@ -0,0 +1 @@ +ecms: [350, 400, 500, 550] diff --git a/ECFAHiggsFactories/ttbar/figures/Top.png b/ECFAHiggsFactories/ttbar/figures/Top.png new file mode 100644 index 0000000..5d1b49c Binary files /dev/null and b/ECFAHiggsFactories/ttbar/figures/Top.png differ diff --git a/ECFAHiggsFactories/ttbar/figures/TopProfile.png b/ECFAHiggsFactories/ttbar/figures/TopProfile.png new file mode 100644 index 0000000..0d8b609 Binary files /dev/null and b/ECFAHiggsFactories/ttbar/figures/TopProfile.png differ diff --git a/ECFAHiggsFactories/ttbar/figures/TopRMS.png b/ECFAHiggsFactories/ttbar/figures/TopRMS.png new file mode 100644 index 0000000..05eccc3 Binary files /dev/null and b/ECFAHiggsFactories/ttbar/figures/TopRMS.png differ diff --git a/ECFAHiggsFactories/ttbar/results/Madgraph/Top_350.root b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_350.root new file mode 100644 index 0000000..1a55be8 Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_350.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Madgraph/Top_400.root b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_400.root new file mode 100644 index 0000000..8f3407a Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_400.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Madgraph/Top_500.root b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_500.root new file mode 100644 index 0000000..37637fe Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_500.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Madgraph/Top_550.root b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_550.root new file mode 100644 index 0000000..9814a4b Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_550.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Sherpa/Top_350.root b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_350.root new file mode 100644 index 0000000..f24c0ce Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_350.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Sherpa/Top_400.root b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_400.root new file mode 100644 index 0000000..93c5cc9 Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_400.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Sherpa/Top_500.root b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_500.root new file mode 100644 index 0000000..bf941a0 Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_500.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Sherpa/Top_550.root b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_550.root new file mode 100644 index 0000000..c7e5837 Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_550.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Whizard/Top_350.root b/ECFAHiggsFactories/ttbar/results/Whizard/Top_350.root new file mode 100644 index 0000000..23b0299 Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Whizard/Top_350.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Whizard/Top_400.root b/ECFAHiggsFactories/ttbar/results/Whizard/Top_400.root new file mode 100644 index 0000000..c8dbb2a Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Whizard/Top_400.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Whizard/Top_500.root b/ECFAHiggsFactories/ttbar/results/Whizard/Top_500.root new file mode 100644 index 0000000..4abd748 Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Whizard/Top_500.root differ diff --git a/ECFAHiggsFactories/ttbar/results/Whizard/Top_550.root b/ECFAHiggsFactories/ttbar/results/Whizard/Top_550.root new file mode 100644 index 0000000..920e299 Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/Whizard/Top_550.root differ diff --git a/ECFAHiggsFactories/ttbar/results/XsectionSummary.root b/ECFAHiggsFactories/ttbar/results/XsectionSummary.root new file mode 100644 index 0000000..32735ae Binary files /dev/null and b/ECFAHiggsFactories/ttbar/results/XsectionSummary.root differ diff --git a/ECFAHiggsFactories/ttbar/results/xsectionSummary.dat b/ECFAHiggsFactories/ttbar/results/xsectionSummary.dat new file mode 100644 index 0000000..e10ed2f --- /dev/null +++ b/ECFAHiggsFactories/ttbar/results/xsectionSummary.dat @@ -0,0 +1,24 @@ +Top_350 sqrts= 350: +MadGraph5_aMC@NLO 0.25291 +- 3.81428e-05 pb +Pythia 0.261514 +- 0.00136477 pb +SHERPA-MC 0.253432 +- 4.40095e-06 pb +Whizard 0.252927 +- 3.38632e-06 pb + +Top_400 sqrts= 400: +MadGraph5_aMC@NLO 0.6231 +- 0.000291671 pb +Pythia 0.680722 +- 0.00334057 pb +SHERPA-MC 0.624237 +- 3.38092e-05 pb +Whizard 0.623093 +- 2.97182e-05 pb + +Top_500 sqrts= 500: +MadGraph5_aMC@NLO 0.54843 +- 0.00040801 pb +Pythia 0.604343 +- 0.0030691 pb +SHERPA-MC 0.549603 +- 4.35292e-05 pb +Whizard 0.548597 +- 4.08029e-05 pb + +Top_550 sqrts= 550: +MadGraph5_aMC@NLO 0.48082 +- 0.00038718 pb +Pythia 0.535849 +- 0.00276845 pb +SHERPA-MC 0.481725 +- 4.17002e-05 pb +Whizard 0.480957 +- 3.9504e-05 pb + diff --git a/ECFAHiggsFactories/ttbar/ttbar.yaml b/ECFAHiggsFactories/ttbar/ttbar.yaml new file mode 100644 index 0000000..fb57de8 --- /dev/null +++ b/ECFAHiggsFactories/ttbar/ttbar.yaml @@ -0,0 +1,35 @@ +Generators: + - Sherpa + - Whizard + - Madgraph + - Pythia + +OutputFormat: hepmc3 +OutDir: Run-Cards +Events: 10000 +EventMode: unweighted + +Model: SM +ISRMode: 0 + +Processes: + Top: + Final: [6, -6] + Order: [2,0] + RandomSeed: 42 + SqrtS: 350 + +Sherpa: + Run: + EW_SCHEME: 3 + FRAGMENTATION: None + + +ParticleData: + 23: + mass: 91.1876 + width: 2.4952 + 6: + mass: 173 + width: 1.4915 + diff --git a/k4GeneratorsConfig/src/analyze2f.cxx b/k4GeneratorsConfig/src/analyze2f.cxx index 8e2c269..bf90689 100644 --- a/k4GeneratorsConfig/src/analyze2f.cxx +++ b/k4GeneratorsConfig/src/analyze2f.cxx @@ -103,6 +103,14 @@ int main(int argc, char** argv) ss << "Invariant Mass(" << pdgIDa << "," << pdgIDb << ")"; TH1D* mpdgapdgb = new TH1D("mpdgapdgb",ss.str().c_str(),1000, 0., sqrts); + ss.clear(); ss.str(""); + ss << "PT(" << pdgIDa << "," << pdgIDb << ")"; + TH1D* ptpdgapdgb = new TH1D("ptpdgapdgb",ss.str().c_str(),1000, 0., sqrts); + + ss.clear(); ss.str(""); + ss << "PZ(" << pdgIDa << "," << pdgIDb << ")"; + TH1D* pzpdgapdgb = new TH1D("pzpdgapdgb",ss.str().c_str(),1000, 0., sqrts); + // loop over the events for (size_t i = 0; i < reader->getEntries(podio::Category::Event); ++i) { if (i) @@ -127,6 +135,8 @@ int main(int argc, char** argv) if ( particleA && particleB ){ edm4hep::LorentzVectorM part = *particleA + *particleB; mpdgapdgb->Fill(part.mass()); + ptpdgapdgb->Fill(part.pt()); + pzpdgapdgb->Fill(part.pz()); } } // release memory after an event