StRoot: StEEmcPool/StEEmcPointMaker/eeSinglePeak.cxx Source File

00001 
00012 #include "eeSinglePeak.h"
00013 #include <TMath.h>
00014 #include <iostream>
00015 
00016 Double_t eeSinglePeak ( Double_t *X, Double_t *P ) 
00017 {
00018   // Shower-shape at shower-max: two gaussians with constrained means
00019   // and widths.
00020 
00021   Double_t x = X[0];
00022   Double_t mean = P[1];
00023 
00024   // Coefficient of first gaussian
00025   Double_t a = 1.0-P[3];
00026 
00027   // Coefficient of second gaussian
00028   Double_t b = P[3];
00029 
00030   // Width of first gaussian
00031   Double_t w1 = P[2];
00032   
00033   // Width of second gaussian
00034   Double_t w2 = P[2]*P[4];
00035 
00036   // Sum of two (normalized) gaussians
00037   Double_t g =
00038     a * TMath::Gaus( x, mean, w1, kTRUE ) + 
00039     b * TMath::Gaus( x, mean, w2, kTRUE );
00040           
00041   return P[0] * g;
00042   
00043 };
00044 
00045 Double_t eeDoublePeak ( Double_t *X, Double_t *P )
00046 {
00047   
00048   Double_t P1[] = { P[0], P[1], P[2], P[3], P[4] };
00049   Double_t P2[] = { P[5], P[6], P[7], P[8], P[9] };
00050 
00051   return eeSinglePeak( X, P1 ) + eeSinglePeak( X, P2 );
00052 
00053 }