Line data Source code
1 : #include "f_Decay.h"
2 :
3 : namespace Tauolapp
4 : {
5 :
6 : void TauolaDecay(int sign_type, double *polx, double* poly,
7 : double * polz, double * poln){
8 :
9 0 : taupos_.npa=1; //tau position in particle list
10 0 : taupos_.npb=1;
11 :
12 0 : double pol[4]={0};
13 0 : dekay_(&sign_type, pol);
14 0 : *polx=pol[0];
15 0 : *poly=pol[1];
16 0 : *polz=pol[2];
17 0 : *poln=pol[3];
18 :
19 : // std::cout << "Polarimetric: "<<pol[0]<<","<<pol[1]
20 : // <<","<<pol[2]<<","<<pol[3]<<std::endl;
21 :
22 :
23 0 : }
24 :
25 : void TauolaWriteDecayToEventRecord(int sign_type){
26 0 : taupos_.npa=1; //tau position in particle list
27 0 : taupos_.npb=1;
28 :
29 0 : double pol[4]={0};
30 :
31 0 : sign_type+=10;
32 0 : dekay_(&sign_type, pol); //write to event record
33 :
34 0 : }
35 :
36 : } // namespace Tauolapp
|