Merge pull request #21 from usarica/newVH

v2.0.6
cms-analysis · Jun 23, 2017 · 6b87cf3 · 6b87cf3
2 parents 4b6b313 + 4ce5c76
commit 6b87cf3
Show file tree

Hide file tree

Showing 23 changed files with 1,882 additions and 456 deletions.
diff --git a/MELA/data/pAvgSmooth_JHUGen_Had_WH_HSMHiggs_13TeV.root b/MELA/data/pAvgSmooth_JHUGen_Had_WH_HSMHiggs_13TeV.root
diff --git a/MELA/data/pAvgSmooth_JHUGen_Had_ZH_HSMHiggs_13TeV.root b/MELA/data/pAvgSmooth_JHUGen_Had_ZH_HSMHiggs_13TeV.root
diff --git a/MELA/data/pAvgSmooth_JHUGen_JJQCD_HSMHiggs_13TeV.root b/MELA/data/pAvgSmooth_JHUGen_JJQCD_HSMHiggs_13TeV.root
diff --git a/MELA/data/pAvgSmooth_JHUGen_JJVBF_HSMHiggs_13TeV.root b/MELA/data/pAvgSmooth_JHUGen_JJVBF_HSMHiggs_13TeV.root
diff --git a/MELA/data/pAvgSmooth_JHUGen_JQCD_HSMHiggs_13TeV.root b/MELA/data/pAvgSmooth_JHUGen_JQCD_HSMHiggs_13TeV.root
diff --git a/MELA/data/pAvgSmooth_JHUGen_ZZGG_HSMHiggs.root b/MELA/data/pAvgSmooth_JHUGen_ZZGG_HSMHiggs.root
diff --git a/MELA/data/pAvgSmooth_MCFM_Had_WH_HSMHiggs_13TeV.root b/MELA/data/pAvgSmooth_MCFM_Had_WH_HSMHiggs_13TeV.root
diff --git a/MELA/data/pAvgSmooth_MCFM_Had_WH_bkgZZ_13TeV.root b/MELA/data/pAvgSmooth_MCFM_Had_WH_bkgZZ_13TeV.root
diff --git a/MELA/data/pAvgSmooth_MCFM_Had_ZH_HSMHiggs_13TeV.root b/MELA/data/pAvgSmooth_MCFM_Had_ZH_HSMHiggs_13TeV.root
diff --git a/MELA/data/pAvgSmooth_MCFM_Had_ZH_bkgZZ_13TeV.root b/MELA/data/pAvgSmooth_MCFM_Had_ZH_bkgZZ_13TeV.root
diff --git a/MELA/data/pAvgSmooth_MCFM_JJQCD_bkgZZ_13TeV.root b/MELA/data/pAvgSmooth_MCFM_JJQCD_bkgZZ_13TeV.root
diff --git a/MELA/data/pAvgSmooth_MCFM_JJVBF_HSMHiggs_13TeV.root b/MELA/data/pAvgSmooth_MCFM_JJVBF_HSMHiggs_13TeV.root
diff --git a/MELA/data/pAvgSmooth_MCFM_JJVBF_bkgZZ_13TeV.root b/MELA/data/pAvgSmooth_MCFM_JJVBF_bkgZZ_13TeV.root
diff --git a/MELA/data/pAvgSmooth_MCFM_ZZGG_HSMHiggs.root b/MELA/data/pAvgSmooth_MCFM_ZZGG_HSMHiggs.root
diff --git a/MELA/data/pAvgSmooth_MCFM_ZZGG_bkgZZ.root b/MELA/data/pAvgSmooth_MCFM_ZZGG_bkgZZ.root
diff --git a/MELA/data/pAvgSmooth_MCFM_ZZQQB_bkgZZ.root b/MELA/data/pAvgSmooth_MCFM_ZZQQB_bkgZZ.root
diff --git a/MELA/interface/Mela.h b/MELA/interface/Mela.h
@@ -289,12 +289,16 @@ class Mela{
 
   /***** ME CONSTANT HANDLES *****/
   // Constants that vary with sqrts due to application of PDFs
+  //
+  MelaPConstant* pAvgSmooth_JHUGen_JQCD_HSMHiggs[TVar::nFermionMassRemovalSchemes-1];
   //
   MelaPConstant* pAvgSmooth_JHUGen_JJQCD_HSMHiggs[TVar::nFermionMassRemovalSchemes-1];
   //
   MelaPConstant* pAvgSmooth_JHUGen_JJVBF_HSMHiggs[TVar::nFermionMassRemovalSchemes-1];
-  //
-  MelaPConstant* pAvgSmooth_JHUGen_JQCD_HSMHiggs[TVar::nFermionMassRemovalSchemes-1];
+  //
+  MelaPConstant* pAvgSmooth_JHUGen_Had_ZH_HSMHiggs[TVar::nFermionMassRemovalSchemes-1];
+  //
+  MelaPConstant* pAvgSmooth_JHUGen_Had_WH_HSMHiggs[TVar::nFermionMassRemovalSchemes-1];
   //
   MelaPConstant* pAvgSmooth_MCFM_JJQCD_bkgZJets_2l2q;
   // Decay ME constants that do not use PDFs
@@ -307,13 +311,41 @@ class Mela{
   MelaPConstant* pAvgSmooth_MCFM_ZZGG_HSMHiggs_4e;
   MelaPConstant* pAvgSmooth_MCFM_ZZGG_HSMHiggs_2mu2e;
   //
+  MelaPConstant* pAvgSmooth_MCFM_JJVBF_HSMHiggs_4mu;
+  MelaPConstant* pAvgSmooth_MCFM_JJVBF_HSMHiggs_4e;
+  MelaPConstant* pAvgSmooth_MCFM_JJVBF_HSMHiggs_2mu2e;
+  //
+  MelaPConstant* pAvgSmooth_MCFM_Had_ZH_HSMHiggs_4mu;
+  MelaPConstant* pAvgSmooth_MCFM_Had_ZH_HSMHiggs_4e;
+  MelaPConstant* pAvgSmooth_MCFM_Had_ZH_HSMHiggs_2mu2e;
+  //
+  MelaPConstant* pAvgSmooth_MCFM_Had_WH_HSMHiggs_4mu;
+  MelaPConstant* pAvgSmooth_MCFM_Had_WH_HSMHiggs_4e;
+  MelaPConstant* pAvgSmooth_MCFM_Had_WH_HSMHiggs_2mu2e;
+  //
   MelaPConstant* pAvgSmooth_MCFM_ZZGG_bkgZZ_4mu;
   MelaPConstant* pAvgSmooth_MCFM_ZZGG_bkgZZ_4e;
   MelaPConstant* pAvgSmooth_MCFM_ZZGG_bkgZZ_2mu2e;
   //
   MelaPConstant* pAvgSmooth_MCFM_ZZQQB_bkgZZ_4mu;
   MelaPConstant* pAvgSmooth_MCFM_ZZQQB_bkgZZ_4e;
   MelaPConstant* pAvgSmooth_MCFM_ZZQQB_bkgZZ_2mu2e;
+  //
+  MelaPConstant* pAvgSmooth_MCFM_JJVBF_bkgZZ_4mu;
+  MelaPConstant* pAvgSmooth_MCFM_JJVBF_bkgZZ_4e;
+  MelaPConstant* pAvgSmooth_MCFM_JJVBF_bkgZZ_2mu2e;
+  //
+  MelaPConstant* pAvgSmooth_MCFM_Had_ZH_bkgZZ_4mu;
+  MelaPConstant* pAvgSmooth_MCFM_Had_ZH_bkgZZ_4e;
+  MelaPConstant* pAvgSmooth_MCFM_Had_ZH_bkgZZ_2mu2e;
+  //
+  MelaPConstant* pAvgSmooth_MCFM_Had_WH_bkgZZ_4mu;
+  MelaPConstant* pAvgSmooth_MCFM_Had_WH_bkgZZ_4e;
+  MelaPConstant* pAvgSmooth_MCFM_Had_WH_bkgZZ_2mu2e;
+  //
+  MelaPConstant* pAvgSmooth_MCFM_JJQCD_bkgZZ_4mu;
+  MelaPConstant* pAvgSmooth_MCFM_JJQCD_bkgZZ_4e;
+  MelaPConstant* pAvgSmooth_MCFM_JJQCD_bkgZZ_2mu2e;
 
 
   /*****************/
@@ -339,14 +371,17 @@ class Mela{
     TVar::MatrixElement me_,
     TVar::Production prod_,
     TVar::Process proc_,
-    const char* relpath,
-    const char* spname
+    TString relpath,
+    TString spname,
+    const bool useSqrts=false
     );
+  void deletePConstantHandle(MelaPConstant*& handle);
   void computeConstant(float& prob);
   void setConstant();
   float getConstant_JHUGenUndecayed();
   float getConstant_4l();
   float getConstant_2l2q();
+  float getConstant_FourFermionDecay(int decid);
 
 };
 

diff --git a/MELA/interface/MelaPConstant.h b/MELA/interface/MelaPConstant.h
@@ -7,6 +7,25 @@
 
 
 class MelaPConstant{
+protected:
+
+  TVar::MatrixElement processME;
+  TVar::Production processProd;
+  TVar::Process processProc;
+
+  TF1* fcnLow;
+  TF1* fcnHigh;
+  TSpline3* fcnMid;
+  TString fname;
+
+  void GetFcnFromFile(const char* path, const char* spname);
+
+  double GetHPropagator(const MelaIO* RcdME, const TVar::VerbosityLevel& verbosity)const;
+  double GetZPropagator(const MelaIO* RcdME, const bool restricted, const TVar::VerbosityLevel& verbosity)const;
+  double GetAssociatedVjjPropagator(const MelaIO* RcdME, const TVar::VerbosityLevel& verbosity)const;
+  double GetVDaughterCouplings(const MelaIO* RcdME, const TVar::VerbosityLevel& verbosity)const;
+  double GetAlphaSatMZ(const MelaIO* RcdME, const int p, const TVar::VerbosityLevel& verbosity)const;
+
 public:
 
   MelaPConstant(
@@ -19,19 +38,12 @@ class MelaPConstant{
 
   virtual ~MelaPConstant();
 
-  double Eval(MelaIO* RcdME, TVar::VerbosityLevel verbosity)const;
-
-private:
-
-  TVar::MatrixElement processME;
-  TVar::Production processProd;
-  TVar::Process processProc;
+  double Eval(const MelaIO* RcdME, TVar::VerbosityLevel verbosity)const;
 
-  TF1* fcnLow;
-  TF1* fcnHigh;
-  TSpline3* fcnMid;
+  bool IsValid(){ return bool(fcnMid!=0); }
 
-  void GetFcnFromFile(const char* path, const char* spname);
+  TString GetFileName(){ return fname; }
+  TString GetSplineName(){ TString sname=""; if (IsValid()) sname=fcnMid->GetName(); return sname; }
 
 };
 

diff --git a/MELA/interface/SuperAssociatedMELA.h b/MELA/interface/SuperAssociatedMELA.h