1 |
/*****************************************************************************
|
2 |
* Project: RooFit *
|
3 |
* *
|
4 |
* This code was autogenerated by RooClassFactory *
|
5 |
*****************************************************************************/
|
6 |
|
7 |
// Your description goes here...
|
8 |
|
9 |
#include "Riostream.h"
|
10 |
|
11 |
#include "RooSUSY3Pdf.h"
|
12 |
#include "RooAbsReal.h"
|
13 |
#include "RooAbsCategory.h"
|
14 |
#include "TMath.h"
|
15 |
|
16 |
ClassImp(RooSUSY3Pdf)
|
17 |
|
18 |
RooSUSY3Pdf::RooSUSY3Pdf(const char *name, const char *title,
|
19 |
RooAbsReal& _inv,
|
20 |
RooAbsReal& _MZ,
|
21 |
RooAbsReal& _M,
|
22 |
RooAbsReal& _m0) :
|
23 |
RooAbsPdf(name,title),
|
24 |
inv("inv","inv",this,_inv),
|
25 |
MZ("MZ","MZ",this,_MZ),
|
26 |
M("M","M",this,_M),
|
27 |
m0("m0","m0",this,_m0)
|
28 |
{
|
29 |
}
|
30 |
|
31 |
|
32 |
RooSUSY3Pdf::RooSUSY3Pdf(const RooSUSY3Pdf& other, const char* name) :
|
33 |
RooAbsPdf(other,name),
|
34 |
inv("inv",this,other.inv),
|
35 |
MZ("MZ",this,other.MZ),
|
36 |
M("s",this,other.M),
|
37 |
m0("m0",this,other.m0)
|
38 |
{
|
39 |
}
|
40 |
|
41 |
|
42 |
|
43 |
Double_t RooSUSY3Pdf::evaluate() const
|
44 |
{
|
45 |
// ENTER EXPRESSION IN TERMS OF VARIABLE ARGUMENTS HERE
|
46 |
//quadratic
|
47 |
Double_t value = inv*TMath::Sqrt(TMath::Power(inv,4) - inv*inv*(m0*m0+M*M) + m0*m0*M*M)/TMath::Power((inv*inv-MZ*MZ),2)*
|
48 |
(-2*TMath::Power(inv,4) + inv*inv*(m0*m0+2*M*M) + m0*m0*M*M);
|
49 |
if (inv>=MZ){value = 0;}
|
50 |
if(value>0) return value;
|
51 |
else return 0.;
|
52 |
}
|
53 |
|
54 |
|
55 |
|