/usr/include/CLHEP/GenericFunctions/AnalyticConvolution.hh is in libclhep-dev 2.1.2.3-1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
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 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 | // -*- C++ -*-
// $Id: AnalyticConvolution.hh,v 1.2 2003/09/06 14:04:13 boudreau Exp $
// ---------------------------------------------------------------------------//
// This function-object makes analytic convolutions of a gaussian plus either //
// an exponential, or else the function exp * (1+/-cos) //
// The choice depends on which constructor is used to build the analytic //
// convolution and which arguments are used. //
// //
// Joe Boudreau, Petar Maksimovic, Hongquan Niu, Craig Blocker //
// //
// ---------------------------------------------------------------------------//
#ifndef _AnalyticConvolution_h_
#define _AnalyticConvolution_h_
#include "CLHEP/GenericFunctions/AbsFunction.hh"
#include "CLHEP/GenericFunctions/Parameter.hh"
#include <complex>
namespace Genfun {
class Gaussian;
class Exponential;
class Cosine;
/**
* @author
* @ingroup genfun
*/
class AnalyticConvolution: public AbsFunction {
FUNCTION_OBJECT_DEF(AnalyticConvolution)
public:
// Flag for mixed or unmixed:
enum Type {MIXED =0, // PDF for mixed events
UNMIXED =1, // PDF for unmixed events
SMEARED_EXP =2, // Exponential (convolve) Gaussian
SMEARED_COS_EXP=3, // Exponential * Cosine (convolve) Gaussian
SMEARED_SIN_EXP=4, // Exponential * Sine (convolve) Gaussian
SMEARED_NEG_EXP=5}; // Negative exponential (convolve) Gaussian
// Constructor
AnalyticConvolution(Type=SMEARED_EXP);
// Copy constructor
AnalyticConvolution(const AnalyticConvolution &right);
// Destructor:
virtual ~AnalyticConvolution();
// Retreive function value
virtual double operator ()(double argument) const;
virtual double operator ()(const Argument & a) const {return operator() (a[0]);}
// Frequency of oscillation
Parameter & frequency();
const Parameter & frequency() const;
// Lifetime of exponential:
Parameter & lifetime();
const Parameter & lifetime() const;
// Width of the gaussian:
Parameter & sigma();
const Parameter & sigma() const;
// The mean of the gaussian:
Parameter & offset();
const Parameter & offset() const;
private:
// These are for calculating mixing terms.
double pow(double x, int n) const ;
double erfc(double x) const ;
std::complex<double> nwwerf(std::complex<double> z) const;
// It is illegal to assign an adjustable constant
const AnalyticConvolution & operator=(const AnalyticConvolution &right);
Parameter _lifetime;
Parameter _frequency;
Parameter _sigma;
Parameter _offset;
Type _type;
};
} // namespace Genfun
#endif
|