From ce05bae45170f65a62435871604c0e4650e162bb Mon Sep 17 00:00:00 2001 From: TjarkMiener Date: Tue, 9 Mar 2021 16:35:03 +0100 Subject: [PATCH] added factor for quarks --- exec/jointLklDM.cc | 3 +++ scripts/computeCLBands.C | 3 +++ scripts/jointLklDM.C | 3 +++ 3 files changed, 9 insertions(+) diff --git a/exec/jointLklDM.cc b/exec/jointLklDM.cc index 7adec42c..9bb603e2 100644 --- a/exec/jointLklDM.cc +++ b/exec/jointLklDM.cc @@ -1152,6 +1152,7 @@ void compute_branonBR(Float_t &mass, Int_t &nChannels, TString *channelval, Doub Float_t particle_mass[9] = {4.18,1.28,173.1,0.511e-3,0.106,1.7768,80.39,91.19,125.0}; // Classification of the elementary particles TString dirac_fermions = "cc_tt_bb_ee_mumu_tautau"; + TString quarks = "cc_tt_bb"; TString gauge_bosons = "WW_ZZ"; TString scalar_bosons = "hh"; // Loop over the channels and compute their branching ratios @@ -1168,6 +1169,8 @@ void compute_branonBR(Float_t &mass, Int_t &nChannels, TString *channelval, Doub ann_crosssection[iChannel] = (mass*mass)/(64. * TMath::Pi()*TMath::Pi()) * (4. * TMath::Power(mass,4) - 4. * mass*mass * particle_mass[iChannel]*particle_mass[iChannel] + 3. * TMath::Power(particle_mass[iChannel],4)) * TMath::Sqrt(1-((particle_mass[iChannel]*particle_mass[iChannel])/(mass*mass))); else if (scalar_bosons.Contains(particle_type[iChannel])) ann_crosssection[iChannel] = (mass*mass)/(32. * TMath::Pi()*TMath::Pi()) * TMath::Power((2.* mass*mass + particle_mass[iChannel]*particle_mass[iChannel]),2) * TMath::Sqrt(1-((particle_mass[iChannel]*particle_mass[iChannel])/(mass*mass))); + // Quarks with a factor of 3 (because of the color of the quarks) + if (quarks.Contains(particle_type[iChannel])) ann_crosssection[iChannel] *= 3.; // WW with a factor 2 (because the W is complex) if(!particle_type[iChannel].CompareTo("WW",TString::kIgnoreCase)) ann_crosssection[iChannel] *= 2.; // hh with a factor 1/2 (because the Higgs is real) diff --git a/scripts/computeCLBands.C b/scripts/computeCLBands.C index b27c9871..e8946e94 100644 --- a/scripts/computeCLBands.C +++ b/scripts/computeCLBands.C @@ -793,6 +793,7 @@ void compute_branonBR(Float_t &mass, Int_t &nChannels, TString *channelval, Doub Float_t particle_mass[9] = {4.18,1.28,173.1,0.511e-3,0.106,1.7768,80.39,91.19,125.0}; // Classification of the elementary particles TString dirac_fermions = "cc_tt_bb_ee_mumu_tautau"; + TString quarks = "cc_tt_bb"; TString gauge_bosons = "WW_ZZ"; TString scalar_bosons = "hh"; // Loop over the channels and compute their branching ratios @@ -809,6 +810,8 @@ void compute_branonBR(Float_t &mass, Int_t &nChannels, TString *channelval, Doub ann_crosssection[iChannel] = (mass*mass)/(64. * TMath::Pi()*TMath::Pi()) * (4. * TMath::Power(mass,4) - 4. * mass*mass * particle_mass[iChannel]*particle_mass[iChannel] + 3. * TMath::Power(particle_mass[iChannel],4)) * TMath::Sqrt(1-((particle_mass[iChannel]*particle_mass[iChannel])/(mass*mass))); else if (scalar_bosons.Contains(particle_type[iChannel])) ann_crosssection[iChannel] = (mass*mass)/(32. * TMath::Pi()*TMath::Pi()) * TMath::Power((2.* mass*mass + particle_mass[iChannel]*particle_mass[iChannel]),2) * TMath::Sqrt(1-((particle_mass[iChannel]*particle_mass[iChannel])/(mass*mass))); + // Quarks with a factor of 3 (because of the color of the quarks) + if (quarks.Contains(particle_type[iChannel])) ann_crosssection[iChannel] *= 3.; // WW with a factor 2 (because the W is complex) if(!particle_type[iChannel].CompareTo("WW",TString::kIgnoreCase)) ann_crosssection[iChannel] *= 2.; // hh with a factor 1/2 (because the Higgs is real) diff --git a/scripts/jointLklDM.C b/scripts/jointLklDM.C index bef3544c..aee871c2 100644 --- a/scripts/jointLklDM.C +++ b/scripts/jointLklDM.C @@ -1226,6 +1226,7 @@ void compute_branonBR(Float_t &mass, Int_t &nChannels, TString *channelval, Doub Float_t particle_mass[9] = {4.18,1.28,173.1,0.511e-3,0.106,1.7768,80.39,91.19,125.0}; // Classification of the elementary particles TString dirac_fermions = "cc_tt_bb_ee_mumu_tautau"; + TString quarks = "cc_tt_bb"; TString gauge_bosons = "WW_ZZ"; TString scalar_bosons = "hh"; // Loop over the channels and compute their branching ratios @@ -1242,6 +1243,8 @@ void compute_branonBR(Float_t &mass, Int_t &nChannels, TString *channelval, Doub ann_crosssection[iChannel] = (mass*mass)/(64. * TMath::Pi()*TMath::Pi()) * (4. * TMath::Power(mass,4) - 4. * mass*mass * particle_mass[iChannel]*particle_mass[iChannel] + 3. * TMath::Power(particle_mass[iChannel],4)) * TMath::Sqrt(1-((particle_mass[iChannel]*particle_mass[iChannel])/(mass*mass))); else if (scalar_bosons.Contains(particle_type[iChannel])) ann_crosssection[iChannel] = (mass*mass)/(32. * TMath::Pi()*TMath::Pi()) * TMath::Power((2.* mass*mass + particle_mass[iChannel]*particle_mass[iChannel]),2) * TMath::Sqrt(1-((particle_mass[iChannel]*particle_mass[iChannel])/(mass*mass))); + // Quarks with a factor of 3 (because of the color of the quarks) + if (quarks.Contains(particle_type[iChannel])) ann_crosssection[iChannel] *= 3.; // WW with a factor 2 (because the W is complex) if(!particle_type[iChannel].CompareTo("WW",TString::kIgnoreCase)) ann_crosssection[iChannel] *= 2.; // hh with a factor 1/2 (because the Higgs is real)