JHUGen MELA  JHUGen v7.5.6, MELA v2.4.2
Matrix element calculations as used in JHUGen.
Classes | Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
RooSpinTwo Class Referenceabstract

#include <RooSpinTwo.h>

Inheritance diagram for RooSpinTwo:
Inheritance graph
[legend]
Collaboration diagram for RooSpinTwo:
Collaboration graph
[legend]

Classes

struct  modelCouplings
 

Public Member Functions

 RooSpinTwo ()
 
 RooSpinTwo (const char *name, const char *title, modelMeasurables const &_measurables, modelParameters const &_parameters, modelCouplings const &_couplings, RooSpin::VdecayType _Vdecay1=RooSpin::kVdecayType_Zll, RooSpin::VdecayType _Vdecay2=RooSpin::kVdecayType_Zll, TVar::VerbosityLevel verbosity_=TVar::ERROR)
 
 RooSpinTwo (const RooSpinTwo &other, const char *name=0)
 
virtual TObject * clone (const char *newname) const =0
 
virtual ~RooSpinTwo ()
 
virtual Double_t evaluate () const =0
 
virtual Int_t getAnalyticalIntegral (RooArgSet &allVars, RooArgSet &analVars, const char *rangeName=0) const =0
 
virtual Double_t analyticalIntegral (Int_t code, const char *rangeName=0) const =0
 
virtual void printParameters () const
 
- Public Member Functions inherited from RooSpin
 RooSpin ()
 
 RooSpin (const char *name, const char *title, modelMeasurables const &_measurables, modelParameters const &_parameters, RooSpin::VdecayType _Vdecay1=RooSpin::kVdecayType_Zll, RooSpin::VdecayType _Vdecay2=RooSpin::kVdecayType_Zll, TVar::VerbosityLevel verbosity_=TVar::ERROR)
 
 RooSpin (const RooSpin &other, const char *name=0)
 
virtual ~RooSpin ()
 
void setVerbosity (TVar::VerbosityLevel verbosity_)
 
virtual void setDecayModes (RooSpin::VdecayType Vdecay1_, RooSpin::VdecayType Vdecay2_)
 
virtual void getMVGamV (Double_t *mV=0, Double_t *gamV=0) const
 
virtual void getMVprimeGamVprime (Double_t *mV=0, Double_t *gamV=0) const
 
virtual void defaultIntegration ()
 
virtual void alwaysIntegrate (Int_t code=1)
 

Protected Member Functions

virtual void evaluatePolarizationTerms (std::vector< Double_t > &Axxyyterm, const Int_t code, bool isGammaV1=false, bool isGammaV2=false) const =0
 
virtual void calculateCi (std::vector< Double_t > &ciRe, std::vector< Double_t > &ciIm, bool isGammaV1=false, bool isGammaV2=false) const
 
virtual void calculateAmplitudes (Double_t &A00Re, Double_t &A00Im, Double_t &AppRe, Double_t &AppIm, Double_t &A0pRe, Double_t &A0pIm, Double_t &Ap0Re, Double_t &Ap0Im, Double_t &AmmRe, Double_t &AmmIm, Double_t &A0mRe, Double_t &A0mIm, Double_t &Am0Re, Double_t &Am0Im, Double_t &ApmRe, Double_t &ApmIm, Double_t &AmpRe, Double_t &AmpIm, bool isGammaV1=false, bool isGammaV2=false) const
 
- Protected Member Functions inherited from RooSpin
virtual void calculatePropagator (Double_t &propRe, Double_t &propIm, Double_t mass, Int_t propType=1) const
 
virtual void calculateVffGVGA (Double_t &gV, Double_t &gA, RooSpin::VdecayType Vdecay, bool isGamma=false) const
 
virtual void calculateVffR1R2 (Double_t &R1Val, Double_t &R2Val, bool isGammaV1=false, bool isGammaV2=false) const
 
virtual Double_t calculateAmplitudeScale (int VGammaVpmode1=0, int VGammaVpmode2=0) const
 
virtual void calculateVprimeffGVGA (Double_t &gV, Double_t &gA, int whichVprime) const
 
virtual void calculateVprimeffR1R2 (Double_t &R1Val, Double_t &R2Val) const
 
virtual void setProxies (modelMeasurables _measurables)
 
virtual void setProxy (RooRealProxy &proxy, RooAbsReal *objectPtr)
 
virtual Bool_t checkFundamentalType (const RooRealProxy &proxy) const
 
virtual Bool_t computeNeededAmplitude (int, int) const
 

Protected Attributes

RooRealProxy b1Val
 
RooRealProxy b2Val
 
RooRealProxy b3Val
 
RooRealProxy b4Val
 
RooRealProxy b5Val
 
RooRealProxy b6Val
 
RooRealProxy b7Val
 
RooRealProxy b8Val
 
RooRealProxy b9Val
 
RooRealProxy b10Val
 
RooRealProxy b1ValIm
 
RooRealProxy b2ValIm
 
RooRealProxy b3ValIm
 
RooRealProxy b4ValIm
 
RooRealProxy b5ValIm
 
RooRealProxy b6ValIm
 
RooRealProxy b7ValIm
 
RooRealProxy b8ValIm
 
RooRealProxy b9ValIm
 
RooRealProxy b10ValIm
 
RooRealProxy Lambda
 
RooRealProxy f_spinz1
 
RooRealProxy f_spinz2
 
- Protected Attributes inherited from RooSpin
RooRealProxy h1
 
RooRealProxy h2
 
RooRealProxy Phi
 
RooRealProxy m1
 
RooRealProxy m2
 
RooRealProxy m12
 
RooRealProxy hs
 
RooRealProxy Phi1
 
RooRealProxy Y
 
RooRealProxy mX
 
RooRealProxy gamX
 
RooRealProxy mW
 
RooRealProxy gamW
 
RooRealProxy mZ
 
RooRealProxy gamZ
 
RooRealProxy mWprime
 
RooRealProxy gamWprime
 
RooRealProxy mZprime
 
RooRealProxy gamZprime
 
RooRealProxy Sin2ThetaW
 
RooRealProxy vev
 
RooRealProxy gVprimeff_decay1_left
 
RooRealProxy gVprimeff_decay1_right
 
RooRealProxy gVprimeff_decay2_left
 
RooRealProxy gVprimeff_decay2_right
 
RooSpin::VdecayType Vdecay1
 
RooSpin::VdecayType Vdecay2
 
Int_t intCodeStart
 
TVar::VerbosityLevel verbosity
 

Additional Inherited Members

- Public Types inherited from RooSpin
enum  VdecayType {
  kVdecayType_Wany =-1, kVdecayType_GammaOnshell =0, kVdecayType_Zll =1, kVdecayType_Znn =2,
  kVdecayType_Zuu =3, kVdecayType_Zdd =4, kVdecayType_Zud =5
}
 
enum  {
  prime_h1 =2, prime_h2 =3, prime_hs =5, prime_Phi =7,
  prime_Phi1 =11, prime_m1 =13, prime_m2 =17, prime_m12 =19,
  prime_Y =23
}
 
- Static Protected Attributes inherited from RooSpin
static constexpr Double_t GeVunit =1e-2
 

Detailed Description

Definition at line 7 of file RooSpinTwo.h.

Constructor & Destructor Documentation

◆ RooSpinTwo() [1/3]

RooSpinTwo::RooSpinTwo ( )

Definition at line 7 of file RooSpinTwo.cc.

7 : RooSpin(){}

◆ RooSpinTwo() [2/3]

RooSpinTwo::RooSpinTwo ( const char *  name,
const char *  title,
modelMeasurables const &  _measurables,
modelParameters const &  _parameters,
modelCouplings const &  _couplings,
RooSpin::VdecayType  _Vdecay1 = RooSpin::kVdecayType_Zll,
RooSpin::VdecayType  _Vdecay2 = RooSpin::kVdecayType_Zll,
TVar::VerbosityLevel  verbosity_ = TVar::ERROR 
)

Definition at line 8 of file RooSpinTwo.cc.

15  : RooSpin(
16  name, title,
17  _measurables, _parameters,
18  _Vdecay1, _Vdecay2,
19  verbosity_
20 ),
21 
22 b1Val("b1Val", "b1Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_1][0])),
23 b2Val("b2Val", "b2Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_2][0])),
24 b3Val("b3Val", "b3Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_3][0])),
25 b4Val("b4Val", "b4Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_4][0])),
26 b5Val("b5Val", "b5Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_5][0])),
27 b6Val("b6Val", "b6Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_6][0])),
28 b7Val("b7Val", "b7Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_7][0])),
29 b8Val("b8Val", "b8Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_8][0])),
30 b9Val("b9Val", "b9Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_9][0])),
31 b10Val("b10Val", "b10Val", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_10][0])),
32 
33 b1ValIm("b1ValIm", "b1ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_1][1])),
34 b2ValIm("b2ValIm", "b2ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_2][1])),
35 b3ValIm("b3ValIm", "b3ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_3][1])),
36 b4ValIm("b4ValIm", "b4ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_4][1])),
37 b5ValIm("b5ValIm", "b5ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_5][1])),
38 b6ValIm("b6ValIm", "b6ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_6][1])),
39 b7ValIm("b7ValIm", "b7ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_7][1])),
40 b8ValIm("b8ValIm", "b8ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_8][1])),
41 b9ValIm("b9ValIm", "b9ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_9][1])),
42 b10ValIm("b10ValIm", "b10ValIm", this, (RooAbsReal&)*(_couplings.bList[gGRAVITON_VV_10][1])),
43 
44 Lambda("Lambda", "Lambda", this, (RooAbsReal&)*(_couplings.Lambda)),
45 
46 f_spinz1("f_spinz1", "f_spinz1", this, (RooAbsReal&)*(_couplings.f_spinz1)),
47 f_spinz2("f_spinz2", "f_spinz2", this, (RooAbsReal&)*(_couplings.f_spinz2))
48 {}

◆ RooSpinTwo() [3/3]

RooSpinTwo::RooSpinTwo ( const RooSpinTwo &  other,
const char *  name = 0 
)

Definition at line 50 of file RooSpinTwo.cc.

50  :
51 RooSpin(other, name),
52 
53 b1Val("b1Val", this, other.b1Val),
54 b2Val("a2Val", this, other.b2Val),
55 b3Val("b3Val", this, other.b3Val),
56 b4Val("b4Val", this, other.b4Val),
57 b5Val("b5Val", this, other.b5Val),
58 b6Val("b6Val", this, other.b6Val),
59 b7Val("b7Val", this, other.b7Val),
60 b8Val("b8Val", this, other.b8Val),
61 b9Val("b9Val", this, other.b9Val),
62 b10Val("b10Val", this, other.b10Val),
63 
64 b1ValIm("b1ValIm", this, other.b1ValIm),
65 b2ValIm("a2ValIm", this, other.b2ValIm),
66 b3ValIm("b3ValIm", this, other.b3ValIm),
67 b4ValIm("b4ValIm", this, other.b4ValIm),
68 b5ValIm("b5ValIm", this, other.b5ValIm),
69 b6ValIm("b6ValIm", this, other.b6ValIm),
70 b7ValIm("b7ValIm", this, other.b7ValIm),
71 b8ValIm("b8ValIm", this, other.b8ValIm),
72 b9ValIm("b9ValIm", this, other.b9ValIm),
73 b10ValIm("b10ValIm", this, other.b10ValIm),
74 
75 Lambda("Lambda", this, other.Lambda),
76 
77 f_spinz1("f_spinz1", this, other.f_spinz1),
78 f_spinz2("f_spinz2", this, other.f_spinz2)
79 {}

◆ ~RooSpinTwo()

virtual RooSpinTwo::~RooSpinTwo ( )
inlinevirtual

Definition at line 30 of file RooSpinTwo.h.

30 {}

Member Function Documentation

◆ analyticalIntegral()

virtual Double_t RooSpinTwo::analyticalIntegral ( Int_t  code,
const char *  rangeName = 0 
) const
pure virtual

Implements RooSpin.

Implemented in RooSpinTwo_7DComplex_ppHVV.

◆ calculateAmplitudes()

void RooSpinTwo::calculateAmplitudes ( Double_t &  A00Re,
Double_t &  A00Im,
Double_t &  AppRe,
Double_t &  AppIm,
Double_t &  A0pRe,
Double_t &  A0pIm,
Double_t &  Ap0Re,
Double_t &  Ap0Im,
Double_t &  AmmRe,
Double_t &  AmmIm,
Double_t &  A0mRe,
Double_t &  A0mIm,
Double_t &  Am0Re,
Double_t &  Am0Im,
Double_t &  ApmRe,
Double_t &  ApmIm,
Double_t &  AmpRe,
Double_t &  AmpIm,
bool  isGammaV1 = false,
bool  isGammaV2 = false 
) const
protectedvirtual

Definition at line 139 of file RooSpinTwo.cc.

145  {
146  Double_t m1_=m1; if (Vdecay1==RooSpin::kVdecayType_GammaOnshell) m1_=0;
147  Double_t m2_=m2; if (Vdecay2==RooSpin::kVdecayType_GammaOnshell) m2_=0;
148 
149  std::vector<Double_t> ciRe;
150  std::vector<Double_t> ciIm;
151  calculateCi(ciRe, ciIm, isGammaV1, isGammaV2);
152 
153  Double_t propV1Re=1, propV2Re=1, propHRe=1;
154  Double_t propV1Im=0, propV2Im=0, propHIm=0;
155  if (Vdecay1!=RooSpin::kVdecayType_GammaOnshell) calculatePropagator(propV1Re, propV1Im, m1_, (isGammaV1 ? 0 : 1));
156  if (Vdecay2!=RooSpin::kVdecayType_GammaOnshell) calculatePropagator(propV2Re, propV2Im, m2_, (isGammaV2 ? 0 : 1));
157  calculatePropagator(propHRe, propHIm, m12, 2);
158 
159  Double_t ampScale = calculateAmplitudeScale((int) isGammaV1, (int) isGammaV2)*pow(GeVunit/Lambda, 2);
160 
161  Double_t c1Re = ciRe.at(0);
162  Double_t c2Re = ciRe.at(1);
163  Double_t c3Re = ciRe.at(2);
164  Double_t c41Re = ciRe.at(3);
165  Double_t c42Re = ciRe.at(4);
166  Double_t c5Re = ciRe.at(5);
167  Double_t c6Re = ciRe.at(6);
168  Double_t c7Re = ciRe.at(7);
169  Double_t c1Im = ciIm.at(0);
170  Double_t c2Im = ciIm.at(1);
171  Double_t c3Im = ciIm.at(2);
172  Double_t c41Im = ciIm.at(3);
173  Double_t c42Im = ciIm.at(4);
174  Double_t c5Im = ciIm.at(5);
175  Double_t c6Im = ciIm.at(6);
176  Double_t c7Im = ciIm.at(7);
177 
178  Double_t eta1 = m1_ / m12;
179  Double_t eta2 = m2_ / m12;
180  Double_t eta1p2 = eta1*eta2;
181 
182  Double_t eta1sq = eta1*eta1;
183  Double_t eta2sq = eta2*eta2;
184  Double_t eta1p2sq = pow(eta1p2, 2);
185 
186  Double_t etas = (1. - eta1sq - eta2sq)/2.;
187  if (pow(eta1+eta2, 2)>1.) etas = -etas;
188  Double_t x = etas;
189  Double_t xsq = x*x;
190  Double_t xxp = (pow(etas, 2)-eta1p2sq);
191  if (xxp<0) xxp=0;
192 
193  Double_t A00Re_tmp=0, A00Im_tmp=0,
194  AppRe_tmp=0, AppIm_tmp=0, A0pRe_tmp=0, A0pIm_tmp=0, Ap0Re_tmp=0, Ap0Im_tmp=0,
195  AmmRe_tmp=0, AmmIm_tmp=0, A0mRe_tmp=0, A0mIm_tmp=0, Am0Re_tmp=0, Am0Im_tmp=0,
196  ApmRe_tmp=0, ApmIm_tmp=0, AmpRe_tmp=0, AmpIm_tmp=0;
197 
198  if (Vdecay1!=RooSpin::kVdecayType_GammaOnshell && Vdecay2!=RooSpin::kVdecayType_GammaOnshell){
199  A00Re_tmp =
200  pow(m12, 4)*sqrt(2./3.)*
201  (
202  c1Re*(
203  eta1p2sq * (xsq - eta1p2sq/4.)
204  - (pow(eta1, 4)+pow(eta2, 4)) * xsq/2.
205  + (pow(eta1, 8)+pow(eta2, 8))/8.
206  + xsq/2.
207  - (pow(eta1, 4)+pow(eta2, 4))/4.
208  + 1.0/8.
209  )
210  + c2Re*2.*xxp*(
211  (pow(eta1, 4) + pow(eta2, 4))
212  - 2.*(eta1p2sq + 2.*xxp)
213  - 1.
214  )
215  -c3Re*8.*pow(xxp, 2)
216  + c41Re*2.*xxp*(1.+eta1sq-eta2sq)
217  + c42Re*2.*xxp*(1.-eta1sq+eta2sq)
218  );
219 
220  A00Im_tmp =
221  pow(m12, 4)*sqrt(2./3.)*
222  (
223  c1Im*(
224  eta1p2sq * (xsq - eta1p2sq/4.)
225  - (pow(eta1, 4)+pow(eta2, 4)) * xsq/2.
226  + (pow(eta1, 8)+pow(eta2, 8))/8.
227  + xsq/2.
228  - (pow(eta1, 4)+pow(eta2, 4))/4.
229  + 1.0/8.
230  )
231  + c2Im*2.*xxp*(
232  (pow(eta1, 4) + pow(eta2, 4))
233  - 2.*(eta1p2sq + 2.*xxp)
234  - 1.
235  )
236  -c3Im*8.*pow(xxp, 2)
237  + c41Im*2.*xxp*(1.+eta1sq-eta2sq)
238  + c42Im*2.*xxp*(1.-eta1sq+eta2sq)
239  );
240  }
241 
242  //-----------------------------------------------------------------------
243  // No m1_/m2_ singularities in A--
244 
245  AmmRe_tmp =
246  pow(m12, 4)/sqrt(6.)*
247  (
248  c1Re/4.*(
249  1. + 4.*xxp - pow(eta1sq-eta2sq, 2)
250  )
251  + c2Re*8.*xxp
252  + c5Im*8.*pow(xxp, 1.5)
253  - c6Im*4.*sqrt(xxp)
254  );
255 
256  AmmIm_tmp =
257  pow(m12, 4)/sqrt(6.)*
258  (
259  c1Im/4.*(
260  1. + 4.*xxp - pow(eta1sq-eta2sq, 2)
261  )
262  + c2Im*8.*xxp
263  - c5Re*8.*pow(xxp, 1.5)
264  + c6Re*4.*sqrt(xxp)
265  );
266 
267  //-----------------------------------------------------------------------
268  // No m1_/m2_ singularities in A++
269 
270  AppRe_tmp =
271  pow(m12, 4)/sqrt(6.)*
272  (
273  c1Re/4.*(
274  1. + 4.*xxp - pow(eta1sq-eta2sq, 2)
275  )
276  + c2Re*8.*xxp
277  - c5Im*8.*pow(xxp, 1.5)
278  + c6Im*4.*sqrt(xxp)
279  );
280 
281  AppIm_tmp =
282  pow(m12, 4)/sqrt(6.)*
283  (
284  c1Im/4.*(
285  1. + 4.*xxp - pow(eta1sq-eta2sq, 2)
286  )
287  + c2Im*8.*xxp
288  + c5Re*8.*pow(xxp, 1.5)
289  - c6Re*4.*sqrt(xxp)
290  );
291 
292  //-----------------------------------------------------------------------
293  // No m1_/m2_ singularities in A+- or A-+
294 
295  AmpRe_tmp = pow(m12, 4)*(c1Re/4.*(1.+4.*xxp-pow(eta1sq-eta2sq, 2)));
296  AmpIm_tmp = pow(m12, 4)*(c1Im/4.*(1.+4.*xxp-pow(eta1sq-eta2sq, 2)));
297  if (Vdecay1!=RooSpin::kVdecayType_GammaOnshell) { AmmIm_tmp *= eta1; AmmRe_tmp *= eta1; AppIm_tmp *= eta1; AppRe_tmp *= eta1; AmpIm_tmp *= eta1; AmpRe_tmp *= eta1; } // Do these multiplications here...
298  if (Vdecay2!=RooSpin::kVdecayType_GammaOnshell) { AmmIm_tmp *= eta2; AmmRe_tmp *= eta2; AppIm_tmp *= eta2; AppRe_tmp *= eta2; AmpIm_tmp *= eta2; AmpRe_tmp *= eta2; }
299  ApmRe_tmp = AmpRe_tmp; // ...for this reason!
300  ApmIm_tmp = AmpIm_tmp;
301 
302  //-----------------------------------------------------------------------
303 
304  Double_t A0m_0p_m0_p0_c4factor = 2.*xxp;
305  Double_t A0m_0p_m0_p0_c7factor = 4.*pow(xxp, 1.5); // x+-i
306 
307  if (Vdecay1!=RooSpin::kVdecayType_GammaOnshell){
308  Double_t A0m_0p_c1factor = (
309  -(pow(eta1, 6)-pow(eta2, 6))/8.
310  + (eta1sq-eta2sq)*(3.*eta1p2sq + 4.*xxp)/8.
311  - pow(eta1sq-eta2sq, 2)/8.
312  + xxp/2.
313  + (1. + (eta1sq-eta2sq))/8.
314  );
315  Double_t A0m_0p_c6factor = sqrt(xxp)*(1.+eta1sq-eta2sq); // x+-i
316 
317  A0mRe_tmp =
318  pow(m12, 4)*
319  (
320  c1Re*A0m_0p_c1factor
321  + c42Re*A0m_0p_m0_p0_c4factor
322  - c6Im*A0m_0p_c6factor
323  - c7Im*A0m_0p_m0_p0_c7factor
324  );
325 
326  A0mIm_tmp =
327  pow(m12, 4)*
328  (
329  c1Im*A0m_0p_c1factor
330  + c42Im*A0m_0p_m0_p0_c4factor
331  + c6Re*A0m_0p_c6factor
332  + c7Re*A0m_0p_m0_p0_c7factor
333  );
334 
335  //-----------------------------------------------------------------------
336 
337  A0pRe_tmp =
338  pow(m12, 4)*
339  (
340  c1Re*A0m_0p_c1factor
341  + c42Re*A0m_0p_m0_p0_c4factor
342  + c6Im*A0m_0p_c6factor
343  + c7Im*A0m_0p_m0_p0_c7factor
344  );
345 
346  A0pIm_tmp =
347  pow(m12, 4)*
348  (
349  c1Im*A0m_0p_c1factor
350  + c42Im*A0m_0p_m0_p0_c4factor
351  - c6Re*A0m_0p_c6factor
352  - c7Re*A0m_0p_m0_p0_c7factor
353  );
354 
355  if (Vdecay2!=RooSpin::kVdecayType_GammaOnshell){
356  A0mIm_tmp *= eta2; A0mRe_tmp *= eta2;
357  A0pIm_tmp *= eta2; A0pRe_tmp *= eta2;
358  }
359  }
360 
361  //-----------------------------------------------------------------------
362 
363  if (Vdecay2!=RooSpin::kVdecayType_GammaOnshell){
364  Double_t Am0_p0_c1factor = (
365  (pow(eta1, 6)-pow(eta2, 6))/8.
366  - (eta1sq-eta2sq)*(3.*eta1p2sq + 4.*xxp)/8.
367  - pow(eta1sq-eta2sq, 2)/8.
368  + xxp/2.
369  + (1. - (eta1sq-eta2sq))/8.
370  );
371  Double_t Am0_p0_c6factor = sqrt(xxp)*(1.-eta1sq+eta2sq); // x+-i
372 
373  Am0Re_tmp =
374  pow(m12, 4)*
375  (
376  c1Re*Am0_p0_c1factor
377  + c41Re*A0m_0p_m0_p0_c4factor
378  - c6Im*Am0_p0_c6factor
379  - c7Im*A0m_0p_m0_p0_c7factor
380  );
381 
382  Am0Im_tmp =
383  pow(m12, 4)*
384  (
385  c1Im*Am0_p0_c1factor
386  + c41Im*A0m_0p_m0_p0_c4factor
387  + c6Re*Am0_p0_c6factor
388  + c7Re*A0m_0p_m0_p0_c7factor
389  );
390 
391  //-----------------------------------------------------------------------
392 
393  Ap0Re_tmp =
394  pow(m12, 4)*
395  (
396  c1Re*Am0_p0_c1factor
397  + c41Re*A0m_0p_m0_p0_c4factor
398  + c6Im*Am0_p0_c6factor
399  + c7Im*A0m_0p_m0_p0_c7factor
400  );
401 
402  Ap0Im_tmp =
403  pow(m12, 4)*
404  (
405  c1Im*Am0_p0_c1factor
406  + c41Im*A0m_0p_m0_p0_c4factor
407  - c6Re*Am0_p0_c6factor
408  - c7Re*A0m_0p_m0_p0_c7factor
409  );
410 
411  if (Vdecay1!=RooSpin::kVdecayType_GammaOnshell){
412  Am0Im_tmp *= eta1; Am0Re_tmp *= eta1;
413  Ap0Im_tmp *= eta1; Ap0Re_tmp *= eta1;
414  }
415  }
416 
417  //-----------------------------------------------------------------------
418 
419  A00Re_tmp *= ampScale;
420  AmmRe_tmp *= ampScale;
421  AppRe_tmp *= ampScale;
422  A0mRe_tmp *= ampScale;
423  A0pRe_tmp *= ampScale;
424  Am0Re_tmp *= ampScale;
425  Ap0Re_tmp *= ampScale;
426  AmpRe_tmp *= ampScale;
427  ApmRe_tmp *= ampScale;
428  A00Im_tmp *= ampScale;
429  AmmIm_tmp *= ampScale;
430  AppIm_tmp *= ampScale;
431  A0mIm_tmp *= ampScale;
432  A0pIm_tmp *= ampScale;
433  Am0Im_tmp *= ampScale;
434  Ap0Im_tmp *= ampScale;
435  AmpIm_tmp *= ampScale;
436  ApmIm_tmp *= ampScale;
437 
438  std::vector<Double_t> A00_reals, A00_imags; A00_reals.push_back(A00Re_tmp); A00_imags.push_back(A00Im_tmp); A00_reals.push_back(propV1Re); A00_imags.push_back(propV1Im); A00_reals.push_back(propV2Re); A00_imags.push_back(propV2Im); A00_reals.push_back(propHRe); A00_imags.push_back(propHIm); AnaMelaHelpers::multiplyComplexNumbers(A00_reals, A00_imags, A00Re, A00Im);
439  std::vector<Double_t> Amm_reals, Amm_imags; Amm_reals.push_back(AmmRe_tmp); Amm_imags.push_back(AmmIm_tmp); Amm_reals.push_back(propV1Re); Amm_imags.push_back(propV1Im); Amm_reals.push_back(propV2Re); Amm_imags.push_back(propV2Im); Amm_reals.push_back(propHRe); Amm_imags.push_back(propHIm); AnaMelaHelpers::multiplyComplexNumbers(Amm_reals, Amm_imags, AmmRe, AmmIm);
440  std::vector<Double_t> App_reals, App_imags; App_reals.push_back(AppRe_tmp); App_imags.push_back(AppIm_tmp); App_reals.push_back(propV1Re); App_imags.push_back(propV1Im); App_reals.push_back(propV2Re); App_imags.push_back(propV2Im); App_reals.push_back(propHRe); App_imags.push_back(propHIm); AnaMelaHelpers::multiplyComplexNumbers(App_reals, App_imags, AppRe, AppIm);
441  std::vector<Double_t> A0m_reals, A0m_imags; A0m_reals.push_back(A0mRe_tmp); A0m_imags.push_back(A0mIm_tmp); A0m_reals.push_back(propV1Re); A0m_imags.push_back(propV1Im); A0m_reals.push_back(propV2Re); A0m_imags.push_back(propV2Im); A0m_reals.push_back(propHRe); A0m_imags.push_back(propHIm); AnaMelaHelpers::multiplyComplexNumbers(A0m_reals, A0m_imags, A0mRe, A0mIm);
442  std::vector<Double_t> A0p_reals, A0p_imags; A0p_reals.push_back(A0pRe_tmp); A0p_imags.push_back(A0pIm_tmp); A0p_reals.push_back(propV1Re); A0p_imags.push_back(propV1Im); A0p_reals.push_back(propV2Re); A0p_imags.push_back(propV2Im); A0p_reals.push_back(propHRe); A0p_imags.push_back(propHIm); AnaMelaHelpers::multiplyComplexNumbers(A0p_reals, A0p_imags, A0pRe, A0pIm);
443  std::vector<Double_t> Am0_reals, Am0_imags; Am0_reals.push_back(Am0Re_tmp); Am0_imags.push_back(Am0Im_tmp); Am0_reals.push_back(propV1Re); Am0_imags.push_back(propV1Im); Am0_reals.push_back(propV2Re); Am0_imags.push_back(propV2Im); Am0_reals.push_back(propHRe); Am0_imags.push_back(propHIm); AnaMelaHelpers::multiplyComplexNumbers(Am0_reals, Am0_imags, Am0Re, Am0Im);
444  std::vector<Double_t> Ap0_reals, Ap0_imags; Ap0_reals.push_back(Ap0Re_tmp); Ap0_imags.push_back(Ap0Im_tmp); Ap0_reals.push_back(propV1Re); Ap0_imags.push_back(propV1Im); Ap0_reals.push_back(propV2Re); Ap0_imags.push_back(propV2Im); Ap0_reals.push_back(propHRe); Ap0_imags.push_back(propHIm); AnaMelaHelpers::multiplyComplexNumbers(Ap0_reals, Ap0_imags, Ap0Re, Ap0Im);
445  std::vector<Double_t> Amp_reals, Amp_imags; Amp_reals.push_back(AmpRe_tmp); Amp_imags.push_back(AmpIm_tmp); Amp_reals.push_back(propV1Re); Amp_imags.push_back(propV1Im); Amp_reals.push_back(propV2Re); Amp_imags.push_back(propV2Im); Amp_reals.push_back(propHRe); Amp_imags.push_back(propHIm); AnaMelaHelpers::multiplyComplexNumbers(Amp_reals, Amp_imags, AmpRe, AmpIm);
446  std::vector<Double_t> Apm_reals, Apm_imags; Apm_reals.push_back(ApmRe_tmp); Apm_imags.push_back(ApmIm_tmp); Apm_reals.push_back(propV1Re); Apm_imags.push_back(propV1Im); Apm_reals.push_back(propV2Re); Apm_imags.push_back(propV2Im); Apm_reals.push_back(propHRe); Apm_imags.push_back(propHIm); AnaMelaHelpers::multiplyComplexNumbers(Apm_reals, Apm_imags, ApmRe, ApmIm);
447 
448 
449  if (
450  A00Re!=A00Re || A00Im!=A00Im ||
451  AppRe!=AppRe || AppIm!=AppIm ||
452  AmmRe!=AmmRe || AmmIm!=AmmIm ||
453  A0pRe!=A0pRe || A0pIm!=A0pIm ||
454  A0mRe!=A0mRe || A0mIm!=A0mIm ||
455  Ap0Re!=Ap0Re || Ap0Im!=Ap0Im ||
456  Am0Re!=Am0Re || Am0Im!=Am0Im ||
457  ApmRe!=ApmRe || ApmIm!=ApmIm ||
458  AmpRe!=AmpRe || AmpIm!=AmpIm ||
459  (
460  A00Re==0 && A00Im==0 &&
461  AppRe==0 && AppIm==0 &&
462  AmmRe==0 && AmmIm==0 &&
463  A0pRe==0 && A0pIm==0 &&
464  A0mRe==0 && A0mIm==0 &&
465  Ap0Re==0 && Ap0Im==0 &&
466  Am0Re==0 && Am0Im==0 &&
467  ApmRe==0 && ApmIm==0 &&
468  AmpRe==0 && AmpIm==0
469  )
470  ){
471  MELAerr << "Some of the amplitudes are NaN or all are 0:" << endl;
472  MELAerr << "A00Re=" << A00Re << ", A00Im=" << A00Im << endl;
473  MELAerr << "AppRe=" << AppRe << ", AppIm=" << AppIm << endl;
474  MELAerr << "AmmRe=" << AmmRe << ", AmmIm=" << AmmIm << endl;
475  MELAerr << "A0pRe=" << A0pRe << ", A0pIm=" << A0pIm << endl;
476  MELAerr << "A0mRe=" << A0mRe << ", A0mIm=" << A0mIm << endl;
477  MELAerr << "Ap0Re=" << Ap0Re << ", Ap0Im=" << Ap0Im << endl;
478  MELAerr << "Am0Re=" << Am0Re << ", Am0Im=" << Am0Im << endl;
479  MELAerr << "ApmRe=" << ApmRe << ", ApmIm=" << ApmIm << endl;
480  MELAerr << "AmpRe=" << AmpRe << ", AmpIm=" << AmpIm << endl;
481 
482  MELAerr << "Possible causes are" << endl;
483  MELAerr << "m12=" << m12 << endl;
484  MELAerr << "m1=" << m1_ << endl;
485  MELAerr << "m2=" << m2_ << endl;
486  MELAerr << "x=" << x << endl;
487  MELAerr << "xxp=" << xxp << endl;
488 
489  MELAerr << "c1 = (" << c1Re << ", " << c1Im << ")" << endl;
490  MELAerr << "c2 = (" << c2Re << ", " << c2Im << ")" << endl;
491  MELAerr << "c3 = (" << c3Re << ", " << c3Im << ")" << endl;
492  MELAerr << "c41 = (" << c41Re << ", " << c41Im << ")" << endl;
493  MELAerr << "c42 = (" << c42Re << ", " << c42Im << ")" << endl;
494  MELAerr << "c5 = (" << c5Re << ", " << c5Im << ")" << endl;
495  MELAerr << "c6 = (" << c6Re << ", " << c6Im << ")" << endl;
496  MELAerr << "c7 = (" << c7Re << ", " << c7Im << ")" << endl;
497  }
498 
499  return;
500 }

◆ calculateCi()

void RooSpinTwo::calculateCi ( std::vector< Double_t > &  ciRe,
std::vector< Double_t > &  ciIm,
bool  isGammaV1 = false,
bool  isGammaV2 = false 
) const
protectedvirtual

Definition at line 81 of file RooSpinTwo.cc.

81  {
82  Double_t mV;
83  getMVGamV(&mV);
84 
85  Double_t m1_=m1; if (Vdecay1==RooSpin::kVdecayType_GammaOnshell) m1_=0;
86  Double_t m2_=m2; if (Vdecay2==RooSpin::kVdecayType_GammaOnshell) m2_=0;
87  Double_t m1sq = pow(m1_, 2);
88  Double_t m2sq = pow(m2_, 2);
89  Double_t mVsq = pow(mV, 2);
90  Double_t m12sq = pow(m12, 2);
91 
92  Double_t s = (m12sq-m1sq-m2sq)/2.;
93  if (m1_>m2_+m12 || m2_>m1_+m12) s = -s;
94  Double_t kappa = s/pow(Lambda, 2);
95 
96  if (!isGammaV1 && !isGammaV2 && !(Vdecay1==RooSpin::kVdecayType_GammaOnshell || Vdecay2==RooSpin::kVdecayType_GammaOnshell)){ // ZZ/WW
97  Double_t c1Re = 2.*(b1Val + b2Val*kappa*(1.+m1sq/s)*(1.+m2sq/s) + b5Val*mVsq/s); ciRe.push_back(c1Re);
98  Double_t c2Re = -0.5*b1Val + b3Val*kappa*(1.-(m1sq+m2sq)/(2*s)) + 2.*b4Val*kappa + b7Val*kappa*mVsq/s; ciRe.push_back(c2Re);
99  Double_t c3Re = -(b2Val/2.+b3Val+2.*b4Val)*kappa*m12sq/s; ciRe.push_back(c3Re);
100  Double_t c41Re = -b1Val - b2Val*kappa - (b2Val*m1sq+b3Val*m2sq+2.*b6Val*mVsq)*kappa/s; ciRe.push_back(c41Re);
101  Double_t c42Re = -b1Val - b2Val*kappa - (b2Val*m2sq+b3Val*m1sq+2.*b6Val*mVsq)*kappa/s; ciRe.push_back(c42Re);
102  Double_t c5Re = 2.*b8Val*kappa*(m12sq)/s; ciRe.push_back(c5Re);
103  Double_t c6Re = b9Val*kappa*mVsq/s; ciRe.push_back(c6Re);
104  Double_t c7Re = b10Val*m12sq*mVsq*pow(kappa/s, 2); ciRe.push_back(c7Re);
105 
106  Double_t c1Im = 2.*(b1ValIm + b2ValIm*kappa*(1.+m1sq/s)*(1.+m2sq/s) + b5ValIm*mVsq/s); ciIm.push_back(c1Im);
107  Double_t c2Im = -0.5*b1ValIm + b3ValIm*kappa*(1.-(m1sq+m2sq)/(2*s)) + 2.*b4ValIm*kappa + b7ValIm*kappa*mVsq/s; ciIm.push_back(c2Im);
108  Double_t c3Im = -(b2ValIm/2.+b3ValIm+2.*b4ValIm)*kappa*m12sq/s; ciIm.push_back(c3Im);
109  Double_t c41Im = -b1ValIm - b2ValIm*kappa - (b2ValIm*m1sq+b3ValIm*m2sq+2.*b6ValIm*mVsq)*kappa/s; ciIm.push_back(c41Im);
110  Double_t c42Im = -b1ValIm - b2ValIm*kappa - (b2ValIm*m2sq+b3ValIm*m1sq+2.*b6ValIm*mVsq)*kappa/s; ciIm.push_back(c42Im);
111  Double_t c5Im = 2.*b8ValIm*kappa*(m12sq)/s; ciIm.push_back(c5Im);
112  Double_t c6Im = b9ValIm*kappa*mVsq/s; ciIm.push_back(c6Im);
113  Double_t c7Im = b10ValIm*m12sq*mVsq*pow(kappa/s, 2); ciIm.push_back(c7Im);
114  }
115  //else if ((!isGammaV1 || !isGammaV2) && !(Vdecay1==RooSpin::kVdecayType_GammaOnshell && Vdecay2==RooSpin::kVdecayType_GammaOnshell)){ // ZGs/ZG
116  //???
117  //}
118  //else{ // GG/GGs
119  else{ // Z/G/Gs - G/Gs
120  Double_t c1Re = 2.*(b1Val + b2Val*kappa*(1.+m1sq/s)*(1.+m2sq/s)); ciRe.push_back(c1Re);
121  Double_t c2Re = -0.5*b1Val + b3Val*kappa*(1.-(m1sq+m2sq)/(2*s)) + 2.*b4Val*kappa; ciRe.push_back(c2Re);
122  Double_t c3Re = -(b2Val/2.+b3Val+2.*b4Val)*kappa*m12sq/s; ciRe.push_back(c3Re);
123  Double_t c41Re = -b1Val - b2Val*kappa - (b2Val*m1sq+b3Val*m2sq)*kappa/s; ciRe.push_back(c41Re);
124  Double_t c42Re = -b1Val - b2Val*kappa - (b2Val*m2sq+b3Val*m1sq)*kappa/s; ciRe.push_back(c42Re);
125  Double_t c5Re = 2.*b8Val*kappa*(m12sq)/s; ciRe.push_back(c5Re);
126  Double_t c6Re = 0; ciRe.push_back(c6Re);
127  Double_t c7Re = 0; ciRe.push_back(c7Re);
128 
129  Double_t c1Im = 2.*(b1ValIm + b2ValIm*kappa*(1.+m1sq/s)*(1.+m2sq/s)); ciIm.push_back(c1Im);
130  Double_t c2Im = -0.5*b1ValIm + b3ValIm*kappa*(1.-(m1sq+m2sq)/(2*s)) + 2.*b4ValIm*kappa; ciIm.push_back(c2Im);
131  Double_t c3Im = -(b2ValIm/2.+b3ValIm+2.*b4ValIm)*kappa*m12sq/s; ciIm.push_back(c3Im);
132  Double_t c41Im = -b1ValIm - b2ValIm*kappa - (b2ValIm*m1sq+b3ValIm*m2sq)*kappa/s; ciIm.push_back(c41Im);
133  Double_t c42Im = -b1ValIm - b2ValIm*kappa - (b2ValIm*m2sq+b3ValIm*m1sq)*kappa/s; ciIm.push_back(c42Im);
134  Double_t c5Im = 2.*b8ValIm*kappa*(m12sq)/s; ciIm.push_back(c5Im);
135  Double_t c6Im = 0; ciIm.push_back(c6Im);
136  Double_t c7Im = 0; ciIm.push_back(c7Im);
137  }
138 }

◆ clone()

virtual TObject* RooSpinTwo::clone ( const char *  newname) const
pure virtual

Implements RooSpin.

Implemented in RooSpinTwo_7DComplex_ppHVV.

◆ evaluate()

virtual Double_t RooSpinTwo::evaluate ( ) const
pure virtual

Implements RooSpin.

Implemented in RooSpinTwo_7DComplex_ppHVV.

◆ evaluatePolarizationTerms()

virtual void RooSpinTwo::evaluatePolarizationTerms ( std::vector< Double_t > &  Axxyyterm,
const Int_t  code,
bool  isGammaV1 = false,
bool  isGammaV2 = false 
) const
protectedpure virtual

Implemented in RooSpinTwo_7DComplex_ppHVV.

◆ getAnalyticalIntegral()

virtual Int_t RooSpinTwo::getAnalyticalIntegral ( RooArgSet &  allVars,
RooArgSet &  analVars,
const char *  rangeName = 0 
) const
pure virtual

Implements RooSpin.

Implemented in RooSpinTwo_7DComplex_ppHVV.

◆ printParameters()

void RooSpinTwo::printParameters ( ) const
virtual

Reimplemented from RooSpin.

Definition at line 502 of file RooSpinTwo.cc.

502  {
503  MELAout << "b1: (" << b1Val << ", " << b1ValIm << ")" << endl;
504  MELAout << "b2: (" << b2Val << ", " << b2ValIm << ")" << endl;
505  MELAout << "b3: (" << b3Val << ", " << b3ValIm << ")" << endl;
506  MELAout << "b4: (" << b4Val << ", " << b4ValIm << ")" << endl;
507  MELAout << "b5: (" << b5Val << ", " << b5ValIm << ")" << endl;
508  MELAout << "b6: (" << b6Val << ", " << b6ValIm << ")" << endl;
509  MELAout << "b7: (" << b7Val << ", " << b7ValIm << ")" << endl;
510  MELAout << "b8: (" << b8Val << ", " << b8ValIm << ")" << endl;
511  MELAout << "b9: (" << b9Val << ", " << b9ValIm << ")" << endl;
512  MELAout << "b10: (" << b10Val << ", " << b10ValIm << ")" << endl;
513  MELAout << "Lambda: " << Lambda << endl;
514  MELAout << "f_spinz1: " << f_spinz1 << endl;
515  MELAout << "f_spinz2: " << f_spinz2 << endl;
516 
517  RooSpin::printParameters();
518 }

Member Data Documentation

◆ b10Val

RooRealProxy RooSpinTwo::b10Val
protected

Definition at line 48 of file RooSpinTwo.h.

◆ b10ValIm

RooRealProxy RooSpinTwo::b10ValIm
protected

Definition at line 59 of file RooSpinTwo.h.

◆ b1Val

RooRealProxy RooSpinTwo::b1Val
protected

Definition at line 39 of file RooSpinTwo.h.

◆ b1ValIm

RooRealProxy RooSpinTwo::b1ValIm
protected

Definition at line 50 of file RooSpinTwo.h.

◆ b2Val

RooRealProxy RooSpinTwo::b2Val
protected

Definition at line 40 of file RooSpinTwo.h.

◆ b2ValIm

RooRealProxy RooSpinTwo::b2ValIm
protected

Definition at line 51 of file RooSpinTwo.h.

◆ b3Val

RooRealProxy RooSpinTwo::b3Val
protected

Definition at line 41 of file RooSpinTwo.h.

◆ b3ValIm

RooRealProxy RooSpinTwo::b3ValIm
protected

Definition at line 52 of file RooSpinTwo.h.

◆ b4Val

RooRealProxy RooSpinTwo::b4Val
protected

Definition at line 42 of file RooSpinTwo.h.

◆ b4ValIm

RooRealProxy RooSpinTwo::b4ValIm
protected

Definition at line 53 of file RooSpinTwo.h.

◆ b5Val

RooRealProxy RooSpinTwo::b5Val
protected

Definition at line 43 of file RooSpinTwo.h.

◆ b5ValIm

RooRealProxy RooSpinTwo::b5ValIm
protected

Definition at line 54 of file RooSpinTwo.h.

◆ b6Val

RooRealProxy RooSpinTwo::b6Val
protected

Definition at line 44 of file RooSpinTwo.h.

◆ b6ValIm

RooRealProxy RooSpinTwo::b6ValIm
protected

Definition at line 55 of file RooSpinTwo.h.

◆ b7Val

RooRealProxy RooSpinTwo::b7Val
protected

Definition at line 45 of file RooSpinTwo.h.

◆ b7ValIm

RooRealProxy RooSpinTwo::b7ValIm
protected

Definition at line 56 of file RooSpinTwo.h.

◆ b8Val

RooRealProxy RooSpinTwo::b8Val
protected

Definition at line 46 of file RooSpinTwo.h.

◆ b8ValIm

RooRealProxy RooSpinTwo::b8ValIm
protected

Definition at line 57 of file RooSpinTwo.h.

◆ b9Val

RooRealProxy RooSpinTwo::b9Val
protected

Definition at line 47 of file RooSpinTwo.h.

◆ b9ValIm

RooRealProxy RooSpinTwo::b9ValIm
protected

Definition at line 58 of file RooSpinTwo.h.

◆ f_spinz1

RooRealProxy RooSpinTwo::f_spinz1
protected

Definition at line 63 of file RooSpinTwo.h.

◆ f_spinz2

RooRealProxy RooSpinTwo::f_spinz2
protected

Definition at line 64 of file RooSpinTwo.h.

◆ Lambda

RooRealProxy RooSpinTwo::Lambda
protected

Definition at line 61 of file RooSpinTwo.h.

The documentation for this class was generated from the following files:
RooSpin::calculateAmplitudeScale
virtual Double_t calculateAmplitudeScale(int VGammaVpmode1=0, int VGammaVpmode2=0) const
Definition: RooSpin.cc:366
RooSpin::printParameters
virtual void printParameters() const
Definition: RooSpin.cc:397
RooSpinTwo::b5ValIm
RooRealProxy b5ValIm
Definition: RooSpinTwo.h:54
gGRAVITON_VV_10
@ gGRAVITON_VV_10
Definition: raw_couplings.txt:144
RooSpinTwo::f_spinz1
RooRealProxy f_spinz1
Definition: RooSpinTwo.h:63
RooSpinTwo::f_spinz2
RooRealProxy f_spinz2
Definition: RooSpinTwo.h:64
RooSpin::m12
RooRealProxy m12
Definition: RooSpin.h:113
RooSpin::calculatePropagator
virtual void calculatePropagator(Double_t &propRe, Double_t &propIm, Double_t mass, Int_t propType=1) const
Definition: RooSpin.cc:133
RooSpinTwo::b7Val
RooRealProxy b7Val
Definition: RooSpinTwo.h:45
RooSpinTwo::b2ValIm
RooRealProxy b2ValIm
Definition: RooSpinTwo.h:51
RooSpinTwo::b4ValIm
RooRealProxy b4ValIm
Definition: RooSpinTwo.h:53
gGRAVITON_VV_1
@ gGRAVITON_VV_1
Definition: raw_couplings.txt:135
RooSpin::m1
RooRealProxy m1
Definition: RooSpin.h:111
RooSpinTwo::b2Val
RooRealProxy b2Val
Definition: RooSpinTwo.h:40
RooSpinTwo::calculateCi
virtual void calculateCi(std::vector< Double_t > &ciRe, std::vector< Double_t > &ciIm, bool isGammaV1=false, bool isGammaV2=false) const
Definition: RooSpinTwo.cc:81
RooSpinTwo::b10ValIm
RooRealProxy b10ValIm
Definition: RooSpinTwo.h:59
AnaMelaHelpers::multiplyComplexNumbers
void multiplyComplexNumbers(std::vector< Double_t > const &reals, std::vector< Double_t > const &imags, Double_t &resRe, Double_t &resIm)
Definition: RooSpin.cc:7
RooSpinTwo::b4Val
RooRealProxy b4Val
Definition: RooSpinTwo.h:42
gGRAVITON_VV_9
@ gGRAVITON_VV_9
Definition: raw_couplings.txt:143
MELAStreamHelpers::MELAout
MELAOutputStreamer MELAout
RooSpinTwo::b7ValIm
RooRealProxy b7ValIm
Definition: RooSpinTwo.h:56
modparameters::kappa
complex(8), public kappa
Definition: mod_Parameters.F90:882
gGRAVITON_VV_2
@ gGRAVITON_VV_2
Definition: raw_couplings.txt:136
RooSpinTwo::b8ValIm
RooRealProxy b8ValIm
Definition: RooSpinTwo.h:57
RooSpin::Vdecay2
RooSpin::VdecayType Vdecay2
Definition: RooSpin.h:136
RooSpinTwo::b5Val
RooRealProxy b5Val
Definition: RooSpinTwo.h:43
RooSpinTwo::b9ValIm
RooRealProxy b9ValIm
Definition: RooSpinTwo.h:58
RooSpin::kVdecayType_GammaOnshell
@ kVdecayType_GammaOnshell
Definition: RooSpin.h:30
gGRAVITON_VV_4
@ gGRAVITON_VV_4
Definition: raw_couplings.txt:138
gGRAVITON_VV_5
@ gGRAVITON_VV_5
Definition: raw_couplings.txt:139
RooSpinTwo::b6Val
RooRealProxy b6Val
Definition: RooSpinTwo.h:44
RooSpinTwo::b3Val
RooRealProxy b3Val
Definition: RooSpinTwo.h:41
RooSpinTwo::b10Val
RooRealProxy b10Val
Definition: RooSpinTwo.h:48
RooSpin::RooSpin
RooSpin()
Definition: RooSpin.cc:36
RooSpinTwo::b1ValIm
RooRealProxy b1ValIm
Definition: RooSpinTwo.h:50
RooSpinTwo::Lambda
RooRealProxy Lambda
Definition: RooSpinTwo.h:61
RooSpin::m2
RooRealProxy m2
Definition: RooSpin.h:112
gGRAVITON_VV_3
@ gGRAVITON_VV_3
Definition: raw_couplings.txt:137
RooSpinTwo::b8Val
RooRealProxy b8Val
Definition: RooSpinTwo.h:46
MELAStreamHelpers::MELAerr
MELAOutputStreamer MELAerr
gGRAVITON_VV_8
@ gGRAVITON_VV_8
Definition: raw_couplings.txt:142
gGRAVITON_VV_6
@ gGRAVITON_VV_6
Definition: raw_couplings.txt:140
RooSpinTwo::b3ValIm
RooRealProxy b3ValIm
Definition: RooSpinTwo.h:52
RooSpin::GeVunit
static constexpr Double_t GeVunit
Definition: RooSpin.h:141
RooSpin::getMVGamV
virtual void getMVGamV(Double_t *mV=0, Double_t *gamV=0) const
Definition: RooSpin.cc:305
RooSpinTwo::b9Val
RooRealProxy b9Val
Definition: RooSpinTwo.h:47
RooSpin::Vdecay1
RooSpin::VdecayType Vdecay1
Definition: RooSpin.h:135
RooSpinTwo::b6ValIm
RooRealProxy b6ValIm
Definition: RooSpinTwo.h:55
gGRAVITON_VV_7
@ gGRAVITON_VV_7
Definition: raw_couplings.txt:141
RooSpinTwo::b1Val
RooRealProxy b1Val
Definition: RooSpinTwo.h:39