|
54 | 54 |
|
55 | 55 | using namespace RooFit; |
56 | 56 |
|
57 | | -HFInvMassFitter::HFInvMassFitter() : mHistoInvMass(nullptr), |
58 | | - mFitOption("L,E"), |
59 | | - mMinMass(0), |
60 | | - mMaxMass(5), |
61 | | - mTypeOfBkgPdf(Expo), |
62 | | - mTypeOfSgnPdf(SingleGaus), |
63 | | - mTypeOfReflPdf(1), |
64 | | - mMassParticle(TDatabasePDG::Instance()->GetParticle("D0")->Mass()), |
65 | | - mMass(1.865), |
66 | | - mMassLowLimit(0), |
67 | | - mMassUpLimit(0), |
68 | | - mMassReflLowLimit(0), |
69 | | - mMassReflUpLimit(0), |
70 | | - mSecMass(1.969), |
71 | | - mSigmaSgn(0.012), |
72 | | - mSecSigma(0.006), |
73 | | - mNSigmaForSidebands(4.), |
74 | | - mNSigmaForSgn(3.), |
75 | | - mSigmaSgnErr(0.), |
76 | | - mSigmaSgnDoubleGaus(0.025), |
77 | | - mFixedMean(false), |
78 | | - mBoundMean(false), |
79 | | - mBoundReflMean(false), |
80 | | - mFixedSigma(false), |
81 | | - mFixedSigmaDoubleGaus(false), |
82 | | - mBoundSigma(false), |
83 | | - mSigmaValue(0.012), |
84 | | - mParamSgn(0.1), |
85 | | - mFracDoubleGaus(0.2), |
86 | | - mFixedRawYield(-1.), |
87 | | - mFixedFracDoubleGaus(false), |
88 | | - mRatioDoubleGausSigma(0.), |
89 | | - mFixedRatioDoubleGausSigma(false), |
90 | | - mReflOverSgn(0), |
91 | | - mEnableReflections(false), |
92 | | - mRawYield(0), |
93 | | - mRawYieldErr(0), |
94 | | - mRawYieldCounted(0), |
95 | | - mRawYieldCountedErr(0), |
96 | | - mBkgYield(0), |
97 | | - mBkgYieldErr(0), |
98 | | - mSignificance(0), |
99 | | - mSignificanceErr(0), |
100 | | - mChiSquareOverNdfTotal(0), |
101 | | - mChiSquareOverNdfBkg(0), |
102 | | - mFixReflOverSgn(false), |
103 | | - mRooMeanSgn(nullptr), |
104 | | - mRooSigmaSgn(nullptr), |
105 | | - mRooSecSigmaSgn(nullptr), |
106 | | - mRooFracDoubleGaus(nullptr), |
107 | | - mSgnPdf(nullptr), |
108 | | - mBkgPdf(nullptr), |
109 | | - mReflPdf(nullptr), |
110 | | - mRooNSgn(nullptr), |
111 | | - mRooNBkg(nullptr), |
112 | | - mRooNRefl(nullptr), |
113 | | - mTotalPdf(nullptr), |
114 | | - mInvMassFrame(nullptr), |
115 | | - mReflFrame(nullptr), |
116 | | - mReflOnlyFrame(nullptr), |
117 | | - mResidualFrame(nullptr), |
118 | | - mResidualFrameForCalculation(nullptr), |
119 | | - mWorkspace(nullptr), |
120 | | - mIntegralHisto(0), |
121 | | - mIntegralBkg(0), |
122 | | - mIntegralSgn(0), |
123 | | - mHistoTemplateRefl(nullptr), |
124 | | - mDrawBgPrefit(false), |
125 | | - mHighlightPeakRegion(false) |
126 | | -{ |
127 | | - // default constructor |
128 | | -} |
129 | | - |
130 | | -HFInvMassFitter::HFInvMassFitter(const TH1* histoToFit, double minValue, double maxValue, int fitTypeBkg, int fitTypeSgn) : mHistoInvMass(nullptr), |
131 | | - mFitOption("L,E"), |
132 | | - mMinMass(minValue), |
133 | | - mMaxMass(maxValue), |
134 | | - mTypeOfBkgPdf(fitTypeBkg), |
135 | | - mTypeOfSgnPdf(fitTypeSgn), |
136 | | - mTypeOfReflPdf(1), |
137 | | - mMassParticle(TDatabasePDG::Instance()->GetParticle("D0")->Mass()), |
138 | | - mMass(1.865), |
139 | | - mMassLowLimit(0), |
140 | | - mMassUpLimit(0), |
141 | | - mMassReflLowLimit(0), |
142 | | - mMassReflUpLimit(0), |
143 | | - mSecMass(1.969), |
144 | | - mSigmaSgn(0.012), |
145 | | - mSecSigma(0.006), |
146 | | - mNSigmaForSidebands(3.), |
147 | | - mNSigmaForSgn(3.), |
148 | | - mSigmaSgnErr(0.), |
149 | | - mSigmaSgnDoubleGaus(0.025), |
150 | | - mFixedMean(false), |
151 | | - mBoundMean(false), |
152 | | - mBoundReflMean(false), |
153 | | - mFixedSigma(false), |
154 | | - mFixedSigmaDoubleGaus(false), |
155 | | - mBoundSigma(false), |
156 | | - mSigmaValue(0.012), |
157 | | - mParamSgn(0.1), |
158 | | - mFracDoubleGaus(0.2), |
159 | | - mFixedRawYield(-1.), |
160 | | - mFixedFracDoubleGaus(false), |
161 | | - mRatioDoubleGausSigma(0.), |
162 | | - mFixedRatioDoubleGausSigma(false), |
163 | | - mReflOverSgn(0), |
164 | | - mEnableReflections(false), |
165 | | - mRawYield(0), |
166 | | - mRawYieldErr(0), |
167 | | - mRawYieldCounted(0), |
168 | | - mRawYieldCountedErr(0), |
169 | | - mBkgYield(0), |
170 | | - mBkgYieldErr(0), |
171 | | - mSignificance(0), |
172 | | - mSignificanceErr(0), |
173 | | - mChiSquareOverNdfTotal(0), |
174 | | - mChiSquareOverNdfBkg(0), |
175 | | - mFixReflOverSgn(false), |
176 | | - mRooMeanSgn(nullptr), |
177 | | - mRooSigmaSgn(nullptr), |
178 | | - mRooSecSigmaSgn(nullptr), |
179 | | - mRooFracDoubleGaus(nullptr), |
180 | | - mSgnPdf(nullptr), |
181 | | - mBkgPdf(nullptr), |
182 | | - mReflPdf(nullptr), |
183 | | - mRooNSgn(nullptr), |
184 | | - mRooNBkg(nullptr), |
185 | | - mRooNRefl(nullptr), |
186 | | - mTotalPdf(nullptr), |
187 | | - mInvMassFrame(nullptr), |
188 | | - mReflFrame(nullptr), |
189 | | - mReflOnlyFrame(nullptr), |
190 | | - mResidualFrame(nullptr), |
191 | | - mResidualFrameForCalculation(nullptr), |
192 | | - mWorkspace(nullptr), |
193 | | - mIntegralHisto(0), |
194 | | - mIntegralBkg(0), |
195 | | - mIntegralSgn(0), |
196 | | - mHistoTemplateRefl(nullptr), |
197 | | - mDrawBgPrefit(false), |
198 | | - mHighlightPeakRegion(false) |
| 57 | +HFInvMassFitter::HFInvMassFitter(const TH1* histoToFit, |
| 58 | + double minValue, |
| 59 | + double maxValue, |
| 60 | + int fitTypeBkg, |
| 61 | + int fitTypeSgn) : mHistoInvMass(nullptr), |
| 62 | + mFitOption("L,E"), |
| 63 | + mMinMass(minValue), |
| 64 | + mMaxMass(maxValue), |
| 65 | + mTypeOfBkgPdf(fitTypeBkg), |
| 66 | + mTypeOfSgnPdf(fitTypeSgn), |
| 67 | + mTypeOfReflPdf(1), |
| 68 | + mMassParticle(TDatabasePDG::Instance()->GetParticle("D0")->Mass()), |
| 69 | + mMass(1.865), |
| 70 | + mMassLowLimit(0), |
| 71 | + mMassUpLimit(0), |
| 72 | + mMassReflLowLimit(0), |
| 73 | + mMassReflUpLimit(0), |
| 74 | + mSecMass(1.969), |
| 75 | + mSigmaSgn(0.012), |
| 76 | + mSecSigma(0.006), |
| 77 | + mNSigmaForSidebands(3.), |
| 78 | + mNSigmaForSgn(3.), |
| 79 | + mSigmaSgnErr(0.), |
| 80 | + mSigmaSgnDoubleGaus(0.025), |
| 81 | + mFixedMean(false), |
| 82 | + mBoundMean(false), |
| 83 | + mBoundReflMean(false), |
| 84 | + mFixedSigma(false), |
| 85 | + mFixedSigmaDoubleGaus(false), |
| 86 | + mBoundSigma(false), |
| 87 | + mSigmaValue(0.012), |
| 88 | + mParamSgn(0.1), |
| 89 | + mFracDoubleGaus(0.2), |
| 90 | + mFixedRawYield(-1.), |
| 91 | + mFixedFracDoubleGaus(false), |
| 92 | + mRatioDoubleGausSigma(0.), |
| 93 | + mFixedRatioDoubleGausSigma(false), |
| 94 | + mReflOverSgn(0), |
| 95 | + mEnableReflections(false), |
| 96 | + mRawYield(0), |
| 97 | + mRawYieldErr(0), |
| 98 | + mRawYieldCounted(0), |
| 99 | + mRawYieldCountedErr(0), |
| 100 | + mBkgYield(0), |
| 101 | + mBkgYieldErr(0), |
| 102 | + mSignificance(0), |
| 103 | + mSignificanceErr(0), |
| 104 | + mChiSquareOverNdfTotal(0), |
| 105 | + mChiSquareOverNdfBkg(0), |
| 106 | + mFixReflOverSgn(false), |
| 107 | + mRooMeanSgn(nullptr), |
| 108 | + mRooSigmaSgn(nullptr), |
| 109 | + mRooSecSigmaSgn(nullptr), |
| 110 | + mRooFracDoubleGaus(nullptr), |
| 111 | + mSgnPdf(nullptr), |
| 112 | + mBkgPdf(nullptr), |
| 113 | + mReflPdf(nullptr), |
| 114 | + mRooNSgn(nullptr), |
| 115 | + mRooNBkg(nullptr), |
| 116 | + mRooNRefl(nullptr), |
| 117 | + mTotalPdf(nullptr), |
| 118 | + mInvMassFrame(nullptr), |
| 119 | + mReflFrame(nullptr), |
| 120 | + mReflOnlyFrame(nullptr), |
| 121 | + mResidualFrame(nullptr), |
| 122 | + mRatioFrame(nullptr), |
| 123 | + mResidualFrameForCalculation(nullptr), |
| 124 | + mWorkspace(nullptr), |
| 125 | + mIntegralHisto(0), |
| 126 | + mIntegralBkg(0), |
| 127 | + mIntegralSgn(0), |
| 128 | + mHistoTemplateRefl(nullptr), |
| 129 | + mDrawBgPrefit(false), |
| 130 | + mHighlightPeakRegion(false) |
199 | 131 | { |
200 | 132 | // standard constructor |
201 | 133 | mHistoInvMass = dynamic_cast<TH1*>(histoToFit->Clone(histoToFit->GetTitle())); |
|
0 commit comments