#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 $") | |
| Double_t | shwfunc (Double_t *x, Double_t *par) |
|
||||||||||||
|
|
|
||||||||||||
|
Definition at line 21 of file ANtpInfoObjectFillerNue.cxx. 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 }
|
1.3.9.1