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
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
|
// -*- C++ -*-
//
// ScalarWaveFunction.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 2003-2011 Peter Richardson, 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.
//
#ifndef ThePEG_ScalarWaveFunction_H
#define ThePEG_ScalarWaveFunction_H
//
// This is the declaration of the ScalarWaveFunction class.
#include "WaveFunctionBase.h"
#include <ThePEG/Helicity/ScalarSpinInfo.h>
#include <ThePEG/EventRecord/Particle.h>
#include <ThePEG/EventRecord/RhoDMatrix.h>
namespace ThePEG {
namespace Helicity {
/** \ingroup Helicity
* \author Peter Richardson
*
* This class is the base class for scalar wavefunctions for use in
* helicity amplitude calculations. The general approach
* is to use a similar philosophy to the FORTRAN HELAS code but with
* additional structure.
*
* This class stores the scalar wavefunction as a complex number and inherits
* from the WaveFunctionBase class for the storage of the particles
* momentum and type.
*
* @see WaveFunctionBase
*/
class ScalarWaveFunction : public WaveFunctionBase {
public:
/** @name Standard constructors and destructors. */
//@{
/**
* Constructor, set the momentum, direction and Wavefunction.
* @param p The momentum.
* @param part The ParticleData pointer
* @param wave The wavefunction.
* @param dir The direction of the particle.
*/
ScalarWaveFunction(const Lorentz5Momentum & p,tcPDPtr part,
Complex wave,Direction dir=intermediate)
: WaveFunctionBase(p,part,dir), _wf(wave)
{
assert(iSpin()==1);
}
/**
* Constructor,set the 5-momentum and zero the wavefunction.
* @param p The 5-momentum.
* @param part The ParticleData pointer.
* @param dir The direction of the particle.
*/
ScalarWaveFunction(const Lorentz5Momentum & p,tcPDPtr part,Direction dir)
: WaveFunctionBase(p,part,dir), _wf(1.0)
{
assert(iSpin()==1);
}
static void calculateWaveFunctions(RhoDMatrix & rho,
tPPtr, Direction) {
rho=RhoDMatrix(PDT::Spin0);
}
static void constructSpinInfo(tPPtr part,Direction, bool time) {
tScalarSpinPtr inspin;
if(part->spinInfo()) inspin=dynamic_ptr_cast<tScalarSpinPtr>(part->spinInfo());
if(inspin) return;
assert(!part->spinInfo());
ScalarSpinPtr temp = new_ptr(ScalarSpinInfo(part->momentum(),time));
part->spinInfo(temp);
}
/**
* Default constructor.
*/
ScalarWaveFunction() : WaveFunctionBase(), _wf(1.0) {}
/**
* Special for spin correlations
*/
ScalarWaveFunction(tPPtr p,Direction dir,bool time)
: WaveFunctionBase(p->momentum(), p->dataPtr(), dir), _wf(1.0)
{
assert(iSpin()==1);
constructSpinInfo(p,dir,time);
}
/**
* Return the wavefunction.
*/
const Complex & wave() const {return _wf;}
public:
void transform(const LorentzRotation & r) {
transformMomentum(r);
}
private:
/**
* Complex number to store the wavefunction.
*/
Complex _wf;
};
}
}
#endif
|
