File: Amplitude.cc

package info (click to toggle)
thepeg 1.8.0-1
  • links: PTS, VCS
  • area: main
  • in suites: jessie, jessie-kfreebsd, wheezy
  • size: 9,312 kB
  • ctags: 11,509
  • sloc: cpp: 57,129; sh: 11,315; java: 3,212; lisp: 1,402; makefile: 830; ansic: 58; perl: 3
file content (64 lines) | stat: -rw-r--r-- 1,969 bytes parent folder | download
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
// -*- C++ -*-
//
// Amplitude.cc is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2011 Leif Lonnblad
//
// ThePEG is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the Amplitude class.
//

#include "Amplitude.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"

using namespace ThePEG;

AbstractNoPIOClassDescription<Amplitude> Amplitude::initAmplitude;
// Definition of the static class description member.

void Amplitude::Init() {
  static ClassDocumentation<Amplitude> documentation
    ("This the abstract class from which any amplitude class, associated with ",
     "a vertex inherits from.");
}



Complex Amplitude::overestimateValue( const tcPDVector & particles,
				      const vector<Lorentz5Momentum> & momenta, 
				      const vector<int> & helicities ) {
  return value(particles,momenta,helicities);
}


Complex Amplitude::value( const PVector & particles,
			  const vector<int> & helicities ) {
  tcPDVector dataParticles;
  vector<Lorentz5Momentum> momenta;
  for ( PVector::const_iterator cit = particles.begin();
	cit != particles.end(); ++cit ) {
    dataParticles.push_back( (*cit)->dataPtr() );
    momenta.push_back( (*cit)->momentum() );
  }
  return value(dataParticles,momenta,helicities);
}


Complex Amplitude::overestimateValue( const PVector & particles,
				      const vector<int> & helicities ) {
  tcPDVector dataParticles;
  vector<Lorentz5Momentum> momenta;
  for ( PVector::const_iterator cit = particles.begin();
	cit != particles.end(); ++cit ) {
    dataParticles.push_back( (*cit)->dataPtr() );
    momenta.push_back( (*cit)->momentum() );
  }
  return overestimateValue(dataParticles,momenta,helicities);
}