]> git.uio.no Git - u/mrichter/AliRoot.git/blame - STARLIGHT/starlight/include/.svn/text-base/gammagammasingle.h.svn-base
STARLIGHT code and interface
[u/mrichter/AliRoot.git] / STARLIGHT / starlight / include / .svn / text-base / gammagammasingle.h.svn-base
CommitLineData
da32329d
AM
1///////////////////////////////////////////////////////////////////////////
2//
3// Copyright 2010
4//
5// This file is part of starlight.
6//
7// starlight is free software: you can redistribute it and/or modify
8// it under the terms of the GNU General Public License as published by
9// the Free Software Foundation, either version 3 of the License, or
10// (at your option) any later version.
11//
12// starlight is distributed in the hope that it will be useful,
13// but WITHOUT ANY WARRANTY; without even the implied warranty of
14// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
15// GNU General Public License for more details.
16//
17// You should have received a copy of the GNU General Public License
18// along with starlight. If not, see <http://www.gnu.org/licenses/>.
19//
20///////////////////////////////////////////////////////////////////////////
21//
22// File and Version Information:
23// $Rev:: $: revision of last commit
24// $Author:: $: author of last commit
25// $Date:: $: date of last commit
26//
27// Description:
28//
29//
30//
31///////////////////////////////////////////////////////////////////////////
32
33
34#ifndef GAMMAGAMMASINGLE_H
35#define GAMMAGAMMASINGLE_H
36
37
38#include <vector>
39
40#include "starlightconstants.h"
41#include "readinluminosity.h"
42#include "beambeamsystem.h"
43#include "randomgenerator.h"
44#include "eventchannel.h"
45#include "starlightconfig.h"
46#ifdef ENABLE_PYTHIA
47#include "pythiadecayer.h"
48#endif
49
50class Gammagammasingle : public eventChannel
51{
52 public:
53 Gammagammasingle(beamBeamSystem& bbsystem);
54 ~Gammagammasingle();
55
56 void singleCrossSection();
57 starlightConstants::event produceEvent(int &ievent);
58
59 upcEvent produceEvent();
60
61 private:
62 double _sigmax[starlightLimits::MAXWBINS][starlightLimits::MAXYBINS];//=new double[500][500]; //decreased from 1000*1000; too big! causes fault!
63 double _sigmaSum;
64 double _wdelta; //Added 7/26/07 for passing sigmadelta to pickw
65 double _remainwd;// "
66 int _ivalwd; // "
67
68 void pickw(double &w);
69 void picky(double &y);
70
71 void parentMomentum(double w,double y,double &E,double &px,double &py,double&pz);
72 double pp(double E);
73 void twoBodyDecay(starlightConstants::particleTypeEnum &ipid,double E,double W,double px0,double py0,double pz0,double &px1,double &py1,double&pz1,double &px2,double &py2,/*double &py2,*/double &pz2,int &iFbadevent);
74 // void transform(double betax,double betay,double betaz,double &E,double &px,double &py,double &pz,int &iFbadevent);
75 void thephi(double W,double px,double py,double pz,double E,double &theta,double &phi);
76
77 double getMass();
78 double getWidth();
79 double getSpin();
80
81 starlightConstants::particleTypeEnum _GGsingInputpidtest;
82 int _GGsingInputnumw;
83 int _GGsingInputnumy;
84 double _GGsingInputGamma_em;
85#ifdef ENABLE_PYTHIA
86 pythiaDecayer _pyDecayer;
87#endif
88
89};
90
91
92#endif // GAMMAGAMMASINGLE_H