JHUGen MELA  v2.4.1
Matrix element calculations as used in JHUGen. MELA is an important tool that was used for the Higgs boson discovery and for precise measurements of its structure and interactions. Please see the website https://spin.pha.jhu.edu/ and papers cited there for more details, and kindly cite those papers when using this code.
ScalarPdfFactory_HVV.cc
Go to the documentation of this file.
1 #include "ScalarPdfFactory_HVV.h"
2 
3 
4 ScalarPdfFactory_HVV::ScalarPdfFactory_HVV(RooSpinZero::modelMeasurables const& measurables_, bool acceptance_, RooSpin::VdecayType V1decay_, RooSpin::VdecayType V2decay_, Bool_t OnshellH_) :
5 ScalarPdfFactory(measurables_, acceptance_, V1decay_, V2decay_, OnshellH_)
6 {
7  measurables.Y=0;
8  initAcceptanceParams();
9  makeParamsConst(true);
10  initPDF();
11 }
12 ScalarPdfFactory_HVV::ScalarPdfFactory_HVV(
13  RooSpinZero::modelMeasurables const& measurables_,
14  double gRatio_[4][8], double gZGsRatio_[4][1], double gGsGsRatio_[3][1], double gVVpRatio_[1][1], double gVpVpRatio_[1][1],
15  bool pmf_applied_, bool acceptance_, RooSpin::VdecayType V1decay_, RooSpin::VdecayType V2decay_, Bool_t OnshellH_
16 ) : ScalarPdfFactory(measurables_, gRatio_, gZGsRatio_, gGsGsRatio_, gVVpRatio_, gVpVpRatio_, pmf_applied_, acceptance_, V1decay_, V2decay_, OnshellH_)
17 {
18  measurables.Y=0;
19  initAcceptanceParams();
20  makeParamsConst(true);
21  initPDF();
22 }
23 
24 ScalarPdfFactory_HVV::~ScalarPdfFactory_HVV(){
25  destroyPDF();
26  destroyAcceptanceParams();
27 }
28 
29 void ScalarPdfFactory_HVV::initAcceptanceParams(){
30  if (acceptance){
31  accepParams.aPhi = new RooRealVar("aPhi", "aPhi", 1.);
32  accepParams.bPhi = new RooRealVar("bPhi", "bPhi", 4.88199e-03);
33  accepParams.cPhi = new RooRealVar("cPhi", "cPhi", 3.69579e-02);
34  accepParams.dPhi = new RooRealVar("dPhi", "dPhi", 0.);
35  accepParams.ePhi = new RooRealVar("ePhi", "ePhi", 0.);
36 
37  accepParams.aPhi1 = new RooRealVar("aPhi1", "aPhi1", 1.);
38  accepParams.bPhi1 = new RooRealVar("bPhi1", "bPhi1", -1.27958e-02);
39  accepParams.cPhi1 = new RooRealVar("cPhi1", "cPhi1", -1.64892e-01);
40  accepParams.dPhi1 = new RooRealVar("dPhi1", "dPhi1", 0.);
41  accepParams.ePhi1 = new RooRealVar("ePhi1", "ePhi1", 0.);
42 
43  accepParams.aH1 = new RooRealVar("aH1", "aH1", 1.);
44  accepParams.bH1 = new RooRealVar("bH1", "bH1", 2.64540e-02);
45  accepParams.cH1 = new RooRealVar("cH1", "cH1", 0.);
46  accepParams.dH1 = new RooRealVar("dH1", "dH1", 0.);
47  accepParams.eH1 = new RooRealVar("eH1", "eH1", 0.);
48 
49  accepParams.aH2 = new RooRealVar("aH2", "aH2", 1.);
50  accepParams.bH2 = new RooRealVar("bH2", "bH2", -3.73167e-01);
51  accepParams.cH2 = new RooRealVar("cH2", "cH2", 0.);
52  accepParams.dH2 = new RooRealVar("dH2", "dH2", 0.);
53  accepParams.eH2 = new RooRealVar("eH2", "eH2", 0.);
54 
55  accepParams.aHs = new RooRealVar("aHs", "aHs", 1.);
56  accepParams.bHs = new RooRealVar("bHs", "bHs", -1.55528e-01);
57  accepParams.cHs = new RooRealVar("cHs", "cHs", 0.);
58  accepParams.dHs = new RooRealVar("dHs", "dHs", 0.);
59  accepParams.eHs = new RooRealVar("eHs", "eHs", 0.);
60 
61  accepParams.aM1 = new RooRealVar("aM1", "aM1", 1.);
62  accepParams.bM1 = new RooRealVar("bM1", "bM1", -1.26554e-02);
63  accepParams.cM1 = new RooRealVar("cM1", "cM1", 3.13526e-05);
64  accepParams.dM1 = new RooRealVar("dM1", "dM1", 0.);
65 
66  accepParams.aM2 = new RooRealVar("aM2", "aM2", 1.);
67  accepParams.bM2 = new RooRealVar("bM2", "bM2", 5.75519e-04);
68  accepParams.cM2 = new RooRealVar("cM2", "cM2", -7.74696e-05);
69  accepParams.dM2 = new RooRealVar("dM2", "dM2", 0.);
70  }
71  else{
72  accepParams.aPhi = new RooRealVar("aPhi", "aPhi", 1.);
73  accepParams.bPhi = new RooRealVar("bPhi", "bPhi", 0.);
74  accepParams.cPhi = new RooRealVar("cPhi", "cPhi", 0.);
75  accepParams.dPhi = new RooRealVar("dPhi", "dPhi", 0.);
76  accepParams.ePhi = new RooRealVar("ePhi", "ePhi", 0.);
77 
78  accepParams.aPhi1 = new RooRealVar("aPhi1", "aPhi1", 1.);
79  accepParams.bPhi1 = new RooRealVar("bPhi1", "bPhi1", 0.);
80  accepParams.cPhi1 = new RooRealVar("cPhi1", "cPhi1", 0.);
81  accepParams.dPhi1 = new RooRealVar("dPhi1", "dPhi1", 0.);
82  accepParams.ePhi1 = new RooRealVar("ePhi1", "ePhi1", 0.);
83 
84  accepParams.aH1 = new RooRealVar("aH1", "aH1", 1.);
85  accepParams.bH1 = new RooRealVar("bH1", "bH1", 0.);
86  accepParams.cH1 = new RooRealVar("cH1", "cH1", 0.);
87  accepParams.dH1 = new RooRealVar("dH1", "dH1", 0.);
88  accepParams.eH1 = new RooRealVar("eH1", "eH1", 0.);
89 
90  accepParams.aH2 = new RooRealVar("aH2", "aH2", 1.);
91  accepParams.bH2 = new RooRealVar("bH2", "bH2", 0.);
92  accepParams.cH2 = new RooRealVar("cH2", "cH2", 0.);
93  accepParams.dH2 = new RooRealVar("dH2", "dH2", 0.);
94  accepParams.eH2 = new RooRealVar("eH2", "eH2", 0.);
95 
96  accepParams.aHs = new RooRealVar("aHs", "aHs", 1.);
97  accepParams.bHs = new RooRealVar("bHs", "bHs", 0.);
98  accepParams.cHs = new RooRealVar("cHs", "cHs", 0.);
99  accepParams.dHs = new RooRealVar("dHs", "dHs", 0.);
100  accepParams.eHs = new RooRealVar("eHs", "eHs", 0.);
101 
102  accepParams.aM1 = new RooRealVar("aM1", "aM1", 0.);
103  accepParams.bM1 = new RooRealVar("bM1", "bM1", 0.);
104  accepParams.cM1 = new RooRealVar("cM1", "cM1", 0.);
105  accepParams.dM1 = new RooRealVar("dM1", "dM1", 0.);
106 
107  accepParams.aM2 = new RooRealVar("aM2", "aM2", 0.);
108  accepParams.bM2 = new RooRealVar("bM2", "bM2", 0.);
109  accepParams.cM2 = new RooRealVar("cM2", "cM2", 0.);
110  accepParams.dM2 = new RooRealVar("dM2", "dM2", 0.);
111  }
112 }
113 void ScalarPdfFactory_HVV::destroyAcceptanceParams(){
114  delete accepParams.aM1;
115  delete accepParams.bM1;
116  delete accepParams.cM1;
117  delete accepParams.dM1;
118  delete accepParams.aM2;
119  delete accepParams.bM2;
120  delete accepParams.cM2;
121  delete accepParams.dM2;
122  delete accepParams.aPhi;
123  delete accepParams.bPhi;
124  delete accepParams.cPhi;
125  delete accepParams.dPhi;
126  delete accepParams.ePhi;
127  delete accepParams.aPhi1;
128  delete accepParams.bPhi1;
129  delete accepParams.cPhi1;
130  delete accepParams.dPhi1;
131  delete accepParams.ePhi1;
132  delete accepParams.aH1;
133  delete accepParams.bH1;
134  delete accepParams.cH1;
135  delete accepParams.dH1;
136  delete accepParams.eH1;
137  delete accepParams.aH2;
138  delete accepParams.bH2;
139  delete accepParams.cH2;
140  delete accepParams.dH2;
141  delete accepParams.eH2;
142  delete accepParams.aHs;
143  delete accepParams.bHs;
144  delete accepParams.cHs;
145  delete accepParams.dHs;
146  delete accepParams.eHs;
147 }
148 
149 void ScalarPdfFactory_HVV::makeParamsConst(bool yesNo){
150  SpinPdfFactory::makeParamsConst(yesNo);
151  if (acceptance && !yesNo){
152  accepParams.aPhi->setConstant(kFALSE);
153  accepParams.bPhi->setConstant(kFALSE);
154  accepParams.cPhi->setConstant(kFALSE);
155  accepParams.dPhi->setConstant(kFALSE);
156  accepParams.ePhi->setConstant(kFALSE);
157 
158  accepParams.aPhi1->setConstant(kFALSE);
159  accepParams.bPhi1->setConstant(kFALSE);
160  accepParams.cPhi1->setConstant(kFALSE);
161  accepParams.dPhi1->setConstant(kFALSE);
162  accepParams.ePhi1->setConstant(kFALSE);
163 
164  accepParams.aH1->setConstant(kFALSE);
165  accepParams.bH1->setConstant(kFALSE);
166  accepParams.cH1->setConstant(kFALSE);
167  accepParams.dH1->setConstant(kFALSE);
168  accepParams.eH1->setConstant(kFALSE);
169 
170  accepParams.aH2->setConstant(kFALSE);
171  accepParams.bH2->setConstant(kFALSE);
172  accepParams.cH2->setConstant(kFALSE);
173  accepParams.dH2->setConstant(kFALSE);
174  accepParams.eH2->setConstant(kFALSE);
175 
176  accepParams.aHs->setConstant(kFALSE);
177  accepParams.bHs->setConstant(kFALSE);
178  accepParams.cHs->setConstant(kFALSE);
179  accepParams.dHs->setConstant(kFALSE);
180  accepParams.eHs->setConstant(kFALSE);
181 
182  accepParams.aM1->setConstant(kFALSE);
183  accepParams.bM1->setConstant(kFALSE);
184  accepParams.cM1->setConstant(kFALSE);
185  accepParams.dM1->setConstant(kFALSE);
186 
187  accepParams.aM2->setConstant(kFALSE);
188  accepParams.bM2->setConstant(kFALSE);
189  accepParams.cM2->setConstant(kFALSE);
190  accepParams.dM2->setConstant(kFALSE);
191  }
192  else{
193  accepParams.aPhi->setConstant(kTRUE);
194  accepParams.bPhi->setConstant(kTRUE);
195  accepParams.cPhi->setConstant(kTRUE);
196  accepParams.dPhi->setConstant(kTRUE);
197  accepParams.ePhi->setConstant(kTRUE);
198 
199  accepParams.aPhi1->setConstant(kTRUE);
200  accepParams.bPhi1->setConstant(kTRUE);
201  accepParams.cPhi1->setConstant(kTRUE);
202  accepParams.dPhi1->setConstant(kTRUE);
203  accepParams.ePhi1->setConstant(kTRUE);
204 
205  accepParams.aH1->setConstant(kTRUE);
206  accepParams.bH1->setConstant(kTRUE);
207  accepParams.cH1->setConstant(kTRUE);
208  accepParams.dH1->setConstant(kTRUE);
209  accepParams.eH1->setConstant(kTRUE);
210 
211  accepParams.aH2->setConstant(kTRUE);
212  accepParams.bH2->setConstant(kTRUE);
213  accepParams.cH2->setConstant(kTRUE);
214  accepParams.dH2->setConstant(kTRUE);
215  accepParams.eH2->setConstant(kTRUE);
216 
217  accepParams.aHs->setConstant(kTRUE);
218  accepParams.bHs->setConstant(kTRUE);
219  accepParams.cHs->setConstant(kTRUE);
220  accepParams.dHs->setConstant(kTRUE);
221  accepParams.eHs->setConstant(kTRUE);
222 
223  accepParams.aM1->setConstant(kTRUE);
224  accepParams.bM1->setConstant(kTRUE);
225  accepParams.cM1->setConstant(kTRUE);
226  accepParams.dM1->setConstant(kTRUE);
227 
228  accepParams.aM2->setConstant(kTRUE);
229  accepParams.bM2->setConstant(kTRUE);
230  accepParams.cM2->setConstant(kTRUE);
231  accepParams.dM2->setConstant(kTRUE);
232  }
233 }
234 void ScalarPdfFactory_HVV::setZZ4fOrdering(bool flag){ PDF->setZZ4fOrdering(flag); }
235 
236 void ScalarPdfFactory_HVV::initPDF(){
237  PDF = new RooSpinZero_7DComplex_withAccep_HVV(
238  "PDF", "PDF",
239  measurables,
240  parameters,
241  couplings,
242  accepParams,
243  V1decay,V2decay
244  );
245  PDF_base = (RooSpin*)PDF;
246 }
247 
248 
249 
250 
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::dH1
RooRealVar * dH1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:25
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::cH2
RooRealVar * cH2
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:29
ScalarPdfFactory_HVV::destroyPDF
virtual void destroyPDF()
Definition: ScalarPdfFactory_HVV.h:31
RooSpin::modelMeasurables::Y
RooAbsReal * Y
Definition: RooSpin.h:59
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::aHs
RooRealVar * aHs
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:32
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::bPhi
RooRealVar * bPhi
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:13
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::dPhi1
RooRealVar * dPhi1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:20
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::cM2
RooRealVar * cM2
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:43
SpinPdfFactory::PDF_base
RooSpin * PDF_base
Definition: SpinPdfFactory.h:43
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::cM1
RooRealVar * cM1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:39
ScalarPdfFactory_HVV.h
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::aH2
RooRealVar * aH2
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:27
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::dHs
RooRealVar * dHs
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:35
SpinPdfFactory::V1decay
RooSpin::VdecayType V1decay
Definition: SpinPdfFactory.h:38
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::bM1
RooRealVar * bM1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:38
ScalarPdfFactory_HVV::~ScalarPdfFactory_HVV
~ScalarPdfFactory_HVV()
Definition: ScalarPdfFactory_HVV.cc:24
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::cPhi1
RooRealVar * cPhi1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:19
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::cH1
RooRealVar * cH1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:24
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::cHs
RooRealVar * cHs
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:34
ScalarPdfFactory_HVV::accepParams
RooSpinZero_7DComplex_withAccep_HVV::accepParameters accepParams
Definition: ScalarPdfFactory_HVV.h:10
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::dM1
RooRealVar * dM1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:40
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::eH2
RooRealVar * eH2
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:31
RooSpinZero_7DComplex_withAccep_HVV::setZZ4fOrdering
void setZZ4fOrdering(Bool_t flag=true)
Definition: RooSpinZero_7DComplex_withAccep_HVV.cc:422
ScalarPdfFactory_HVV::initPDF
virtual void initPDF()
Definition: ScalarPdfFactory_HVV.cc:236
ScalarPdfFactory_HVV::ScalarPdfFactory_HVV
ScalarPdfFactory_HVV(RooSpin::modelMeasurables const &measurables_, bool acceptance_=false, RooSpin::VdecayType V1decay_=RooSpin::kVdecayType_Zll, RooSpin::VdecayType V2decay_=RooSpin::kVdecayType_Zll, Bool_t OnshellH_=true)
Definition: ScalarPdfFactory_HVV.cc:4
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::dM2
RooRealVar * dM2
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:44
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::aPhi1
RooRealVar * aPhi1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:17
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::bPhi1
RooRealVar * bPhi1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:18
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::dH2
RooRealVar * dH2
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:30
SpinPdfFactory::parameters
RooSpin::modelParameters parameters
Definition: SpinPdfFactory.h:12
SpinPdfFactory::V2decay
RooSpin::VdecayType V2decay
Definition: SpinPdfFactory.h:39
ScalarPdfFactory_HVV::destroyAcceptanceParams
virtual void destroyAcceptanceParams()
Definition: ScalarPdfFactory_HVV.cc:113
ScalarPdfFactory_HVV::makeParamsConst
void makeParamsConst(bool yesNo=true)
Definition: ScalarPdfFactory_HVV.cc:149
RooSpin
Definition: RooSpin.h:25
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::dPhi
RooRealVar * dPhi
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:15
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::bH1
RooRealVar * bH1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:23
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::aH1
RooRealVar * aH1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:22
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::bH2
RooRealVar * bH2
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:28
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::aPhi
RooRealVar * aPhi
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:12
SpinPdfFactory::makeParamsConst
virtual void makeParamsConst(bool yesNo)
Definition: SpinPdfFactory.cc:154
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::ePhi1
RooRealVar * ePhi1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:21
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::bHs
RooRealVar * bHs
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:33
ScalarPdfFactory_HVV::setZZ4fOrdering
void setZZ4fOrdering(bool flag=true)
Definition: ScalarPdfFactory_HVV.cc:234
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::cPhi
RooRealVar * cPhi
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:14
RooSpinZero_7DComplex_withAccep_HVV
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:7
ScalarPdfFactory::acceptance
bool acceptance
Definition: ScalarPdfFactory.h:90
ScalarPdfFactory::couplings
RooSpinZero::modelCouplings couplings
Definition: ScalarPdfFactory.h:12
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::eHs
RooRealVar * eHs
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:36
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::eH1
RooRealVar * eH1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:26
RooSpin::modelMeasurables
Definition: RooSpin.h:50
SpinPdfFactory::measurables
RooSpin::modelMeasurables measurables
Definition: SpinPdfFactory.h:11
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::ePhi
RooRealVar * ePhi
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:16
ScalarPdfFactory_HVV::PDF
RooSpinZero_7DComplex_withAccep_HVV * PDF
Definition: ScalarPdfFactory_HVV.h:25
ScalarPdfFactory_HVV::initAcceptanceParams
virtual void initAcceptanceParams()
Definition: ScalarPdfFactory_HVV.cc:29
RooSpin::VdecayType
VdecayType
Definition: RooSpin.h:28
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::aM2
RooRealVar * aM2
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:41
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::bM2
RooRealVar * bM2
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:42
ScalarPdfFactory
Definition: ScalarPdfFactory.h:10
RooSpinZero_7DComplex_withAccep_HVV::accepParameters::aM1
RooRealVar * aM1
Definition: RooSpinZero_7DComplex_withAccep_HVV.h:37
Generated by   doxygen 1.8.17