]>
Commit | Line | Data |
---|---|---|
da0e9ce3 | 1 | //-------------------------------------------------------------------------- |
2 | // | |
3 | // Environment: | |
4 | // This software is part of the EvtGen package developed jointly | |
5 | // for the BaBar and CLEO collaborations. If you use all or part | |
6 | // of it, please give an appropriate acknowledgement. | |
7 | // | |
8 | // Copyright Information: See EvtGen/COPYRIGHT | |
9 | // Copyright (C) 1998 Caltech, UCSB | |
10 | // | |
11 | // Module: EvtGen/EvtAbsLineShape.hh | |
12 | // | |
13 | // Description: Class to keep the particle properties for | |
14 | // one particle | |
15 | // | |
16 | // Modification history: | |
17 | // | |
18 | // Lange March 10, 2001 Module created | |
19 | // | |
20 | //------------------------------------------------------------------------ | |
21 | ||
22 | #ifndef EVTABSLINESHAPE_HH | |
23 | #define EVTABSLINESHAPE_HH | |
24 | ||
25 | #include "EvtGenBase/EvtSpinType.hh" | |
26 | #include "EvtGenBase/EvtId.hh" | |
27 | #include <vector> | |
28 | ||
29 | class EvtId; | |
30 | ||
31 | class EvtAbsLineShape { | |
32 | ||
33 | public: | |
34 | ||
35 | EvtAbsLineShape(); | |
36 | EvtAbsLineShape(double mass, double width, double maxRange, EvtSpinType::spintype sp); | |
37 | virtual ~EvtAbsLineShape(); | |
38 | EvtAbsLineShape& operator=(const EvtAbsLineShape& x); | |
39 | EvtAbsLineShape(const EvtAbsLineShape& x); | |
40 | ||
41 | double getMass() {return _mass;} | |
42 | double getMassMin() {return _massMin;} | |
43 | double getMassMax() {return _massMax;} | |
44 | double getMaxRange() {return _maxRange;} | |
45 | double getWidth() {return _width;} | |
46 | EvtSpinType::spintype getSpinType() {return _spin;} | |
47 | virtual double rollMass(); | |
48 | virtual EvtAbsLineShape* clone(); | |
49 | ||
50 | void reSetMass(double mass) { _mass=mass;} | |
51 | void reSetWidth(double width) { _width=width;} | |
52 | void reSetMassMin(double mass) { _massMin=mass;} | |
53 | void reSetMassMax(double mass) { _massMax=mass;} | |
54 | virtual void reSetBlatt(double blatt) {}; | |
55 | void includeBirthFactor(bool yesno) { _includeBirthFact = yesno; } | |
56 | void includeDecayFactor(bool yesno) { _includeDecayFact = yesno; } | |
57 | void setPWForDecay( int spin, EvtId d1, EvtId d2) { | |
58 | _userSetPW.push_back(spin); | |
59 | _userSetPWD1.push_back(d1); | |
60 | _userSetPWD2.push_back(d2); | |
61 | } | |
62 | void setPWForBirthL( int spin, EvtId par, EvtId othD) { | |
63 | _userSetBirthPW.push_back(spin); | |
64 | _userSetBirthOthD.push_back(othD); | |
65 | _userSetBirthPar.push_back(par); | |
66 | } | |
67 | ||
68 | virtual double getRandMass(EvtId *parId, int nDaug, EvtId *dauId, EvtId *othDaugId,double maxMass, double *dauMasses); | |
69 | virtual double getMassProb(double mass, double massPar, int nDaug, double *massDau); | |
70 | ||
71 | void fixForSP8() { _applyFixForSP8=true;} | |
72 | protected: | |
73 | ||
74 | bool _includeDecayFact; | |
75 | bool _includeBirthFact; | |
76 | double _mass; | |
77 | double _massMin; | |
78 | double _massMax; | |
79 | double _width; | |
80 | double _maxRange; | |
81 | ||
82 | // allow for special cases where the default method of picking the | |
83 | //lowest allowed partial wave for a decay is not the right answer. | |
84 | // string is "<spin> <daughter1> <daughter2>" | |
85 | //new 9/12/2003 Lange | |
86 | std::vector<EvtId> _userSetPWD1,_userSetPWD2; | |
87 | std::vector<int> _userSetPW; | |
88 | ||
89 | // also do it for birth factors | |
90 | std::vector<EvtId> _userSetBirthPar,_userSetBirthOthD; | |
91 | std::vector<int> _userSetBirthPW; | |
92 | ||
93 | EvtSpinType::spintype _spin; | |
94 | ||
95 | bool _applyFixForSP8; | |
96 | }; | |
97 | ||
98 | #endif | |
99 |