ref ttbar

key4hep · Sep 26, 2024 · ac528c3 · ac528c3
1 parent f9d7ecb
commit ac528c3
Show file tree

Hide file tree

Showing 20 changed files with 70 additions and 0 deletions.
diff --git a/ECFAHiggsFactories/ttbar/ecms.dat 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
diff --git a/ECFAHiggsFactories/ttbar/figures/TopProfile.png b/ECFAHiggsFactories/ttbar/figures/TopProfile.png
diff --git a/ECFAHiggsFactories/ttbar/figures/TopRMS.png b/ECFAHiggsFactories/ttbar/figures/TopRMS.png
diff --git a/ECFAHiggsFactories/ttbar/results/Madgraph/Top_350.root b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_350.root
diff --git a/ECFAHiggsFactories/ttbar/results/Madgraph/Top_400.root b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_400.root
diff --git a/ECFAHiggsFactories/ttbar/results/Madgraph/Top_500.root b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_500.root
diff --git a/ECFAHiggsFactories/ttbar/results/Madgraph/Top_550.root b/ECFAHiggsFactories/ttbar/results/Madgraph/Top_550.root
diff --git a/ECFAHiggsFactories/ttbar/results/Sherpa/Top_350.root b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_350.root
diff --git a/ECFAHiggsFactories/ttbar/results/Sherpa/Top_400.root b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_400.root
diff --git a/ECFAHiggsFactories/ttbar/results/Sherpa/Top_500.root b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_500.root
diff --git a/ECFAHiggsFactories/ttbar/results/Sherpa/Top_550.root b/ECFAHiggsFactories/ttbar/results/Sherpa/Top_550.root
diff --git a/ECFAHiggsFactories/ttbar/results/Whizard/Top_350.root b/ECFAHiggsFactories/ttbar/results/Whizard/Top_350.root
diff --git a/ECFAHiggsFactories/ttbar/results/Whizard/Top_400.root b/ECFAHiggsFactories/ttbar/results/Whizard/Top_400.root
diff --git a/ECFAHiggsFactories/ttbar/results/Whizard/Top_500.root b/ECFAHiggsFactories/ttbar/results/Whizard/Top_500.root
diff --git a/ECFAHiggsFactories/ttbar/results/Whizard/Top_550.root b/ECFAHiggsFactories/ttbar/results/Whizard/Top_550.root
diff --git a/ECFAHiggsFactories/ttbar/results/XsectionSummary.root b/ECFAHiggsFactories/ttbar/results/XsectionSummary.root
diff --git a/ECFAHiggsFactories/ttbar/results/xsectionSummary.dat 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
@@ -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
@@ -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