// Decays are performed using Pythia.
// andreas.morsch@cern.ch
-#include "AliGenParam.h"
-#include "AliDecayer.h"
-#include "AliGenMUONlib.h"
-#include "AliRun.h"
-#include <TParticle.h>
-#include <TParticlePDG.h>
+#include <TCanvas.h>
#include <TDatabasePDG.h>
+#include <TF1.h>
+#include <TH1F.h>
#include <TLorentzVector.h>
+#include <TMath.h>
+#include <TParticle.h>
+#include <TParticlePDG.h>
+#include <TROOT.h>
#include <TVirtualMC.h>
-#include <TF1.h>
-#include <TCanvas.h>
-#include <TH1.h>
-#include <TMath.h>
+#include "AliDecayer.h"
+#include "AliGenMUONlib.h"
+#include "AliGenParam.h"
#include "AliMC.h"
+#include "AliRun.h"
ClassImp(AliGenParam)
//End_Html
//____________________________________________________________
-AliGenParam::AliGenParam()
+ AliGenParam::AliGenParam():
+ fPtParaFunc(0),
+ fYParaFunc(0),
+ fIpParaFunc(0),
+ fPtPara(0),
+ fYPara(0),
+ fParam(0),
+ fdNdy0(0.),
+ fYWgt(0.),
+ fPtWgt(0.),
+ fBias(0.),
+ fTrials(0),
+ fDeltaPt(0.01),
+ fDecayer(0)
{
-// Deafault constructor
- fPtPara = 0;
- fYPara = 0;
- fParam = 0;
- fAnalog = kAnalog;
- SetDeltaPt();
- fDecayer = 0;
-
-
+// Default constructor
}
//____________________________________________________________
-AliGenParam::AliGenParam(Int_t npart, AliGenLib * Library, Int_t param, char* tname):AliGenMC(npart)
+AliGenParam::AliGenParam(Int_t npart, AliGenLib * Library, Int_t param, char* tname)
+ :AliGenMC(npart),
+ fPtParaFunc(Library->GetPt(param, tname)),
+ fYParaFunc (Library->GetY (param, tname)),
+ fIpParaFunc(Library->GetIp(param, tname)),
+ fPtPara(0),
+ fYPara(0),
+ fParam(param),
+ fdNdy0(0.),
+ fYWgt(0.),
+ fPtWgt(0.),
+ fBias(0.),
+ fTrials(0),
+ fDeltaPt(0.01),
+ fDecayer(0)
{
// Constructor using number of particles parameterisation id and library
fName = "Param";
fTitle= "Particle Generator using pT and y parameterisation";
-
- fPtParaFunc = Library->GetPt(param, tname);
- fYParaFunc = Library->GetY (param, tname);
- fIpParaFunc = Library->GetIp(param, tname);
-
- fPtPara = 0;
- fYPara = 0;
- fParam = param;
fAnalog = kAnalog;
SetForceDecay();
- SetDeltaPt();
}
//____________________________________________________________
-AliGenParam::AliGenParam(Int_t npart, Int_t param, const char* tname, const char* name):AliGenMC(npart)
+AliGenParam::AliGenParam(Int_t npart, Int_t param, const char* tname, const char* name):
+ AliGenMC(npart),
+ fPtParaFunc(0),
+ fYParaFunc (0),
+ fIpParaFunc(0),
+ fPtPara(0),
+ fYPara(0),
+ fParam(param),
+ fdNdy0(0.),
+ fYWgt(0.),
+ fPtWgt(0.),
+ fBias(0.),
+ fTrials(0),
+ fDeltaPt(0.01),
+ fDecayer(0)
{
// Constructor using parameterisation id and number of particles
//
- fName = name;
- fTitle= "Particle Generator using pT and y parameterisation";
+ fName = name;
+ fTitle= "Particle Generator using pT and y parameterisation";
AliGenLib* pLibrary = new AliGenMUONlib();
-
fPtParaFunc = pLibrary->GetPt(param, tname);
fYParaFunc = pLibrary->GetY (param, tname);
fIpParaFunc = pLibrary->GetIp(param, tname);
- fPtPara = 0;
- fYPara = 0;
- fParam = param;
fAnalog = kAnalog;
fChildSelect.Set(5);
for (Int_t i=0; i<5; i++) fChildSelect[i]=0;
SetChildPtRange();
SetChildPhiRange();
SetChildThetaRange();
- SetDeltaPt();
}
//____________________________________________________________
Double_t (*PtPara) (Double_t*, Double_t*),
Double_t (*YPara ) (Double_t* ,Double_t*),
Int_t (*IpPara) (TRandom *))
- :AliGenMC(npart)
+ :AliGenMC(npart),
+
+ fPtParaFunc(PtPara),
+ fYParaFunc(YPara),
+ fIpParaFunc(IpPara),
+ fPtPara(0),
+ fYPara(0),
+ fParam(param),
+ fdNdy0(0.),
+ fYWgt(0.),
+ fPtWgt(0.),
+ fBias(0.),
+ fTrials(0),
+ fDeltaPt(0.01),
+ fDecayer(0)
{
// Constructor
// Gines Martinez 1/10/99
fName = "Param";
fTitle= "Particle Generator using pT and y parameterisation";
- fPtParaFunc = PtPara;
- fYParaFunc = YPara;
- fIpParaFunc = IpPara;
-//
- fPtPara = 0;
- fYPara = 0;
- fParam = param;
fAnalog = kAnalog;
fChildSelect.Set(5);
for (Int_t i=0; i<5; i++) fChildSelect[i]=0;
SetChildPtRange();
SetChildPhiRange();
SetChildThetaRange();
- SetDeltaPt();
}
AliGenParam::AliGenParam(const AliGenParam & Param)
- :AliGenMC(Param)
+ :AliGenMC(Param),
+ fPtParaFunc(0),
+ fYParaFunc(0),
+ fIpParaFunc(0),
+ fPtPara(0),
+ fYPara(0),
+ fParam(0),
+ fdNdy0(0.),
+ fYWgt(0.),
+ fPtWgt(0.),
+ fBias(0.),
+ fTrials(0),
+ fDeltaPt(0.01),
+ fDecayer(0)
{
// Copy constructor
Param.Copy(*this);
if (fPtPara) fPtPara->Delete();
fPtPara = new TF1(name, fPtParaFunc, fPtMin, fPtMax,0);
+ gROOT->GetListOfFunctions()->Remove(fPtPara);
// Set representation precision to 10 MeV
Int_t npx= Int_t((fPtMax - fPtMin) / fDeltaPt);
sprintf(name, "y-parameterisation for %s", GetName());
if (fYPara) fYPara->Delete();
fYPara = new TF1(name, fYParaFunc, fYMin, fYMax, 0);
+ gROOT->GetListOfFunctions()->Remove(fYPara);
+
sprintf(name, "pt-for-%s", GetName());
TF1 ptPara(name ,fPtParaFunc, 0, 15, 0);
fdNdy0=fYParaFunc(&y1,&y2);
//
// Integral over generation region
- Float_t intYS = yPara.Integral(fYMin, fYMax);
- Float_t intPt0 = ptPara.Integral(0,15);
- Float_t intPtS = ptPara.Integral(fPtMin,fPtMax);
+ Float_t intYS = yPara.Integral(fYMin, fYMax,(Double_t*) 0x0,1.e-6);
+ Float_t intPt0 = ptPara.Integral(0,15,(Double_t *) 0x0,1.e-6);
+ Float_t intPtS = ptPara.Integral(fPtMin,fPtMax,(Double_t*) 0x0,1.e-6);
Float_t phiWgt=(fPhiMax-fPhiMin)/2./TMath::Pi();
if (fAnalog == kAnalog) {
fYWgt = intYS/fdNdy0;
// Normalisation for selected kinematic region
//
Float_t ratio =
- fPtPara->Integral(ptMin,ptMax) / fPtPara->Integral( fPtPara->GetXmin(), fPtPara->GetXmax()) *
- fYPara->Integral(yMin,yMax)/fYPara->Integral(fYPara->GetXmin(),fYPara->GetXmax()) *
+ fPtPara->Integral(ptMin,ptMax,(Double_t *)0,1.e-6) / fPtPara->Integral( fPtPara->GetXmin(), fPtPara->GetXmax(),(Double_t *)0,1.e-6) *
+ fYPara->Integral(yMin,yMax,(Double_t *)0,1.e-6)/fYPara->Integral(fYPara->GetXmin(),fYPara->GetXmax(),(Double_t *)0,1.e-6) *
(phiMax-phiMin)/360.;
return TMath::Abs(ratio);
}