ANtpInfoObjectFillerNue.cxx File Reference

#include "TH1F.h"
#include "TF1.h"
#include "TCanvas.h"
#include "TClonesArray.h"
#include "TMinuit.h"
#include "StandardNtuple/NtpStRecord.h"
#include "CandNtupleSR/NtpSRRecord.h"
#include "CandNtupleSR/NtpSREvent.h"
#include "CandNtupleSR/NtpSRTrack.h"
#include "CandNtupleSR/NtpSRShower.h"
#include "CandNtupleSR/NtpSRStrip.h"
#include "MessageService/MsgService.h"
#include "NueAna/NueAnaTools/SntpHelpers.h"
#include "AnalysisNtuples/Module/ANtpInfoObjectFillerNue.h"
#include "AnalysisNtuples/ANtpNueInfo.h"
#include "AnalysisNtuples/ANtpDefaultValue.h"
#include "VertexFinder/NtpVtxFinder/NtpVtxFinder.h"

Go to the source code of this file.

Functions

 CVSID ("$Id: ANtpInfoObjectFillerNue.cxx,v 1.3 2008/11/10 15:03:30 boehm Exp $")
static Double_t shwfunc (Double_t *x, Double_t *par)

Function Documentation

CVSID ( "$Id: ANtpInfoObjectFillerNue.  cxx,
v 1.3 2008/11/10 15:03:30 boehm Exp $"   
)
static Double_t shwfunc ( Double_t *  x,
Double_t *  par 
) [static]

Definition at line 21 of file ANtpInfoObjectFillerNue.cxx.

References MuELoss::R.

Referenced by ANtpInfoObjectFillerNue::Reset().

00022 {
00023 
00024   // function to be fitted for em showers par[0]=a par[1]=E0
00025   Float_t R = 1.46676;
00026   Float_t xx=R*x[0];
00027 //  cout<<"In shwfunc "<<xx<<" par[0] "<<par[0]<<" par[1] "<<par[1]<<" par[2] "<<par[2]<<endl;
00028 
00029   Double_t lnf = TMath::Log(R*par[2]*par[1])+(par[0]-1)*TMath::Log(xx*par[1])-
00030      par[1]*xx-TMath::LnGamma(par[0]);
00031 
00032   Double_t f = exp(lnf);
00033 
00034 //  Double_t f = 1.46676*par[2]*TMath::Power(par[1],par[0])*
00035 //               TMath::Power(xx,par[0]-1)*
00036 //               TMath::Exp(-par[1]*(xx))/TMath::Gamma(par[0]);
00037   return f;
00038 }


Generated on 16 Jan 2018 for loon by  doxygen 1.6.1