← Back to branch summary

~maddevelopers/mg5amcnlo/2.5.0

« back to all changes in this revision

Viewing changes to aloha/template_files/oxxxxx.cc

  • Committer: olivier Mattelaer
  • Date: 2016-08-22 20:43:28 UTC
  • mfrom: (262.13.38 plugin_mode)
  • Revision ID: olivier.mattelaer@uclouvain.be-20160822204328-898681196houupn4
merge the plugin_mode branch

Show diffs side-by-side

added added

removed removed

Lines of Context:
aloha/template_files/oxxxxx.cc
10
10
  fo[1] = complex<double>(p[1]*nsf,p[2]*nsf);
11
11
  nh = nhel*nsf;
12
12
  if (fmass != 0.000){
13
 
    pp = min(p[0],pow(pow(p[1],2)+pow(p[2],2)+pow(p[3],2),0.5));
 
13
    pp = min(p[0],sqrt((p[1]*p[1])+(p[2]*p[2])+(p[3]*p[3])));
14
14
    if (pp == 0.000){
15
 
      sqm[0] = pow(abs(fmass),0.5);
 
15
      sqm[0] = sqrt(std::abs(fmass));
16
16
      sqm[1] = Sgn(sqm[0],fmass); 
17
17
      ip = -((1-nh)/2) * nhel ;
18
18
      im = (1+nh)/2 * nhel;
19
 
      fo[2] = im *sqm[abs(ip)];
20
 
      fo[3] = ip*nsf*sqm[abs(ip)];
21
 
      fo[4] = im*nsf*sqm[abs(im)];
22
 
      fo[5] = ip*sqm[abs(im)];
 
19
      fo[2] = im *sqm[std::abs(ip)];
 
20
      fo[3] = ip*nsf*sqm[std::abs(ip)];
 
21
      fo[4] = im*nsf*sqm[std::abs(im)];
 
22
      fo[5] = ip*sqm[std::abs(im)];
23
23
    }
24
24
    else{
25
 
      pp = min(p[0],pow(pow(p[1],2)+pow(p[2],2)+pow(p[3],2),0.5));
 
25
      pp = min(p[0],sqrt((p[1]*p[1])+(p[2]*p[2])+(p[3]*p[3])));
26
26
      sf[0] = double(1+nsf+(1-nsf)*nh)*0.5;
27
27
      sf[1] = double(1+nsf-(1-nsf)*nh)*0.5;
28
 
      omega[0] = pow(p[0]+pp,0.5);
 
28
      omega[0] = sqrt(p[0]+pp);
29
29
      omega[1] = fmass/omega[0];
30
30
      ip = (1+nh)/2 ;
31
31
      im = (1-nh)/2 ;
32
32
      sfomeg[0] = sf[0]*omega[ip];
33
33
      sfomeg[1] = sf[1]*omega[im];
34
34
      pp3 = max(pp+p[3],0.00);
35
 
      chi[0] = complex<double>(pow(pp3*0.5/pp,0.5),0.00);
 
35
      chi[0] = complex<double>(sqrt(pp3*0.5/pp),0.00);
36
36
      if (pp3 == 0.00){ 
37
37
        chi[1] = complex<double>(-nh,0.00);
38
38
      }
39
39
      else{
40
 
        chi[1] = complex<double>(nh*p[1],-p[2])/pow(2.0*pp*pp3,0.5);
 
40
        chi[1] = complex<double>(nh*p[1],-p[2])/sqrt(2.0*pp*pp3);
41
41
      }
42
42
      fo[2] = sfomeg[1]*chi[im];
43
43
      fo[3] = sfomeg[1]*chi[ip];
50
50
      sqp0p3 = 0.00;
51
51
    }
52
52
    else{
53
 
      sqp0p3 = pow(max(p[0]+p[3],0.00),0.5)*nsf;
 
53
      sqp0p3 = sqrt(max(p[0]+p[3],0.00))*nsf;
54
54
    }
55
55
    chi[0] = complex<double>(sqp0p3,0.00);
56
56
    if(sqp0p3 == 0.000){ 
57
 
      chi[1] = complex<double>(-nhel,0.00)*pow(2.0*p[0],0.5);
 
57
      chi[1] = complex<double>(-nhel,0.00)*sqrt(2.0*p[0]);
58
58
    }
59
59
    else{
60
60
      chi[1] = complex<double>(nh*p[1],-p[2])/sqp0p3;