|
17 | 17 | /// \author Xinye Peng <xinye.peng@cern.ch> |
18 | 18 | /// \author Biao Zhang <biao.zhang@cern.ch> |
19 | 19 | /// \author Oleksii Lubynets <oleksii.lubynets@cern.ch> |
| 20 | +/// \author Phil Stahlhut <phil.lennart.stahlhut@cern.ch> |
20 | 21 |
|
21 | 22 | #if !defined(__CINT__) || defined(__CLING__) |
22 | 23 |
|
@@ -88,13 +89,13 @@ int runMassFitter(const TString& configFileName) |
88 | 89 | config.ParseStream(is); |
89 | 90 | fclose(configFile); |
90 | 91 |
|
91 | | - Bool_t isMc = config["IsMC"].GetBool(); |
92 | | - Bool_t writeSignalPar = config["WriteSignalPar"].GetBool(); |
93 | | - TString inputFileName = config["InFileName"].GetString(); |
94 | | - TString reflFileName = config["ReflFileName"].GetString(); |
| 92 | + Bool_t const isMc = config["IsMC"].GetBool(); |
| 93 | + Bool_t const writeSignalPar = config["WriteSignalPar"].GetBool(); |
| 94 | + TString const inputFileName = config["InFileName"].GetString(); |
| 95 | + TString const reflFileName = config["ReflFileName"].GetString(); |
95 | 96 | TString outputFileName = config["OutFileName"].GetString(); |
96 | | - TString particleName = config["Particle"].GetString(); |
97 | | - TString collisionSystem = config["CollisionSystem"].GetString(); |
| 97 | + TString const particleName = config["Particle"].GetString(); |
| 98 | + TString const collisionSystem = config["CollisionSystem"].GetString(); |
98 | 99 |
|
99 | 100 | std::vector<std::string> inputHistoName; |
100 | 101 | std::vector<std::string> promptHistoName; |
@@ -274,40 +275,30 @@ int runMassFitter(const TString& configFileName) |
274 | 275 | inputFile->Close(); |
275 | 276 |
|
276 | 277 | // define output histos |
277 | | - auto hRawYieldsSignal = new TH1D("hRawYieldsSignal", ";" + sliceVarName + "(" + sliceVarUnit + ");raw yield", |
278 | | - nSliceVarBins, sliceVarLimits.data()); |
279 | | - auto hRawYieldsSignalCounted = new TH1D("hRawYieldsSignalCounted", ";" + sliceVarName + "(" + sliceVarUnit + ");raw yield via bin count", |
280 | | - nSliceVarBins, sliceVarLimits.data()); |
281 | | - auto hRawYieldsMean = new TH1D( |
282 | | - "hRawYieldsMean", ";" + sliceVarName + "(" + sliceVarUnit + ");mean (GeV/#it{c}^{2})", |
283 | | - nSliceVarBins, sliceVarLimits.data()); |
284 | | - auto hRawYieldsSigma = new TH1D( |
285 | | - "hRawYieldsSigma", ";" + sliceVarName + "(" + sliceVarUnit + ");width (GeV/#it{c}^{2})", |
286 | | - nSliceVarBins, sliceVarLimits.data()); |
287 | | - auto hRawYieldsSecSigma = new TH1D( |
288 | | - "hRawYieldsSecSigma", ";" + sliceVarName + "(" + sliceVarUnit + ");width (GeV/#it{c}^{2})", |
289 | | - nSliceVarBins, sliceVarLimits.data()); |
290 | | - auto hRawYieldsFracDoubleGaus = new TH1D( |
291 | | - "hRawYieldsFracDoubleGaus", ";" + sliceVarName + "(" + sliceVarUnit + ");fraction of double gaussian", |
292 | | - nSliceVarBins, sliceVarLimits.data()); |
293 | | - auto* hRawYieldsSignificance = new TH1D( |
294 | | - "hRawYieldsSignificance", |
295 | | - ";" + sliceVarName + "(" + sliceVarUnit + ");significance (3#sigma)", nSliceVarBins, sliceVarLimits.data()); |
296 | | - auto* hRawYieldsSgnOverBkg = |
297 | | - new TH1D("hRawYieldsSgnOverBkg", ";" + sliceVarName + "(" + sliceVarUnit + ");S/B (3#sigma)", |
298 | | - nSliceVarBins, sliceVarLimits.data()); |
299 | | - auto* hRawYieldsBkg = |
300 | | - new TH1D("hRawYieldsBkg", ";" + sliceVarName + "(" + sliceVarUnit + ");Background (3#sigma)", |
301 | | - nSliceVarBins, sliceVarLimits.data()); |
302 | | - auto* hRawYieldsChiSquareBkg = |
303 | | - new TH1D("hRawYieldsChiSquareBkg", |
304 | | - ";" + sliceVarName + "(" + sliceVarUnit + ");#chi^{2}/#it{ndf}", nSliceVarBins, sliceVarLimits.data()); |
305 | | - auto* hRawYieldsChiSquareTotal = |
306 | | - new TH1D("hRawYieldsChiSquareTotal", |
307 | | - ";" + sliceVarName + "(" + sliceVarUnit + ");#chi^{2}/#it{ndf}", nSliceVarBins, sliceVarLimits.data()); |
308 | | - auto* hReflectionOverSignal = |
309 | | - new TH1D("hReflectionOverSignal", ";" + sliceVarName + "(" + sliceVarUnit + ");Refl/Signal", |
310 | | - nSliceVarBins, sliceVarLimits.data()); |
| 278 | + auto* hRawYieldsSignal = new TH1D("hRawYieldsSignal", |
| 279 | + ";" + sliceVarName + "(" + sliceVarUnit + ");raw yield", nSliceVarBins, sliceVarLimits.data()); |
| 280 | + auto* hRawYieldsSignalCounted = new TH1D("hRawYieldsSignalCounted", |
| 281 | + ";" + sliceVarName + "(" + sliceVarUnit + ");raw yield via bin count", nSliceVarBins, sliceVarLimits.data()); |
| 282 | + auto* hRawYieldsMean = new TH1D("hRawYieldsMean", |
| 283 | + ";" + sliceVarName + "(" + sliceVarUnit + ");mean (GeV/#it{c}^{2})", nSliceVarBins, sliceVarLimits.data()); |
| 284 | + auto* hRawYieldsSigma = new TH1D("hRawYieldsSigma", |
| 285 | + ";" + sliceVarName + "(" + sliceVarUnit + ");width (GeV/#it{c}^{2})", nSliceVarBins, sliceVarLimits.data()); |
| 286 | + auto* hRawYieldsSecSigma = new TH1D("hRawYieldsSecSigma", |
| 287 | + ";" + sliceVarName + "(" + sliceVarUnit + ");width (GeV/#it{c}^{2})", nSliceVarBins, sliceVarLimits.data()); |
| 288 | + auto* hRawYieldsFracDoubleGaus = new TH1D("hRawYieldsFracDoubleGaus", |
| 289 | + ";" + sliceVarName + "(" + sliceVarUnit + ");fraction of double gaussian", nSliceVarBins, sliceVarLimits.data()); |
| 290 | + auto* hRawYieldsSignificance = new TH1D("hRawYieldsSignificance", |
| 291 | + ";" + sliceVarName + "(" + sliceVarUnit + ");significance (3#sigma)", nSliceVarBins, sliceVarLimits.data()); |
| 292 | + auto* hRawYieldsSgnOverBkg = new TH1D("hRawYieldsSgnOverBkg", |
| 293 | + ";" + sliceVarName + "(" + sliceVarUnit + ");S/B (3#sigma)", nSliceVarBins, sliceVarLimits.data()); |
| 294 | + auto* hRawYieldsBkg = new TH1D("hRawYieldsBkg", |
| 295 | + ";" + sliceVarName + "(" + sliceVarUnit + ");Background (3#sigma)", nSliceVarBins, sliceVarLimits.data()); |
| 296 | + auto* hRawYieldsChiSquareBkg = new TH1D("hRawYieldsChiSquareBkg", |
| 297 | + ";" + sliceVarName + "(" + sliceVarUnit + ");#chi^{2}/#it{ndf}", nSliceVarBins, sliceVarLimits.data()); |
| 298 | + auto* hRawYieldsChiSquareTotal = new TH1D("hRawYieldsChiSquareTotal", |
| 299 | + ";" + sliceVarName + "(" + sliceVarUnit + ");#chi^{2}/#it{ndf}", nSliceVarBins, sliceVarLimits.data()); |
| 300 | + auto* hReflectionOverSignal = new TH1D("hReflectionOverSignal", |
| 301 | + ";" + sliceVarName + "(" + sliceVarUnit + ");Refl/Signal", nSliceVarBins, sliceVarLimits.data()); |
311 | 302 |
|
312 | 303 | const Int_t nConfigsToSave = 6; |
313 | 304 | auto* hFitConfig = new TH2F("hfitConfig", "Fit Configurations", nConfigsToSave, 0, 6, nSliceVarBins, sliceVarLimits.data()); |
|
0 commit comments