-
Notifications
You must be signed in to change notification settings - Fork 510
Expand file tree
/
Copy pathCachingTF1.h
More file actions
52 lines (45 loc) · 1.7 KB
/
Copy pathCachingTF1.h
File metadata and controls
52 lines (45 loc) · 1.7 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
/// @brief Extension to ROOT::TF1 allowing to cache integral
/// @author Sandro Wenzel, sandro.wenzel@cern.ch
#ifndef ALICEO2_CACHINGTF1_H
#define ALICEO2_CACHINGTF1_H
#include <TF1.h>
namespace o2
{
namespace math_utils
{
class CachingTF1 : public TF1
{
///
/// Class extending TF1 with capability to store expensive
/// internal caches (integral, etc) when streaming out
/// This can immensely speed up the construction of the TF1
/// (e.g., for the purpose of random number generation from arbitrary distributions)
///
public:
using TF1::TF1;
~CachingTF1() override = default;
// get reading access to fIntegral member
std::vector<double> const& getIntegralVector() const { return fIntegral; }
private:
// in the original TF1 implementation, these members
// are marked transient; by simply introducing something that
// points to them they will now be written correctly to disc
std::vector<double>* mIntegralCache = &fIntegral;
std::vector<double>* mAlphaCache = &fAlpha;
std::vector<double>* mBetaCache = &fBeta;
std::vector<double>* mGammaCache = &fGamma;
ClassDefOverride(CachingTF1, 1);
};
} // namespace math_utils
} // namespace o2
#endif