1 #ifndef MMPULSESHAPEFUNCTION_H
2 #define MMPULSESHAPEFUNCTION_H
15 double pulse_shape(
double *time,
double *par) {
20 double func_electronic,dtsbeta,variable,puiss;
22 double fFunc_max = par[0];
23 double fAlpha = par[1];
24 double fTime_max = par[2];
25 double fBeta = par[3];
27 double albet = fAlpha*fBeta ;
29 if( albet <= 0 )
return( (
double)0. );
30 double dt = time[0]-fTime_max ;
33 variable=1.+dt/albet ;
34 puiss=TMath::Power(variable,fAlpha);
35 func_electronic=fFunc_max*puiss*TMath::Exp(-dtsbeta);
37 else func_electronic = 0. ;
39 return func_electronic ;
Definition: MMPulseShapeFunction.h:9