* provided "as is" without express or implied warranty. *
**************************************************************************/
-/*
-$Log$
-Revision 1.21 2002/11/21 16:13:03 alibrary
-Removing AliMCProcess and AliMC
-
-Revision 1.20 2002/10/29 14:26:49 hristov
-Code clean-up (F.Carminati)
-
-Revision 1.19 2002/03/12 11:06:03 morsch
-Add particle status code to argument list of SetTrack(..).
-
-Revision 1.18 2002/02/18 19:23:00 hristov
-Put protection in the destructor
-
-Revision 1.17 2001/11/26 15:46:35 morsch
-Set default value for vertex truncation in constructor.
-
-Revision 1.16 2001/10/16 08:45:37 morsch
-All vertex truncation and event vertex moved here.
-
-Revision 1.15 2001/07/28 10:43:41 hristov
-Default argument initialised once
-
-Revision 1.14 2001/07/27 17:13:12 morsch
-- SetTrack, KeepTrack and SetHighWaterMark methods added to delegate to local
-stack or stack owned by AliRun. (Piotr Skowronski, A.M.)
-- Changes needed for afterburners (Piotr Skowronski)
-
-Revision 1.13 2001/05/16 14:57:22 alibrary
-New files for folders and Stack
-
-Revision 1.12 2001/02/02 11:12:50 morsch
-Add Vertex() method that allows to get vertex from merging manager, if needed.
-
-Revision 1.11 2001/01/26 19:58:48 hristov
-Major upgrade of AliRoot code
-
-Revision 1.10 2000/12/21 15:30:18 fca
-Correcting coding convention violations
-
-Revision 1.9 2000/10/04 10:08:01 fca
-Correction of minor typing mistakes
-
-Revision 1.8 2000/10/02 21:28:14 fca
-Removal of useless dependecies via forward declarations
-
-Revision 1.7 2000/07/12 08:56:25 fca
-Coding convention correction and warning removal
-
-Revision 1.6 2000/07/11 18:24:59 fca
-Coding convention corrections + few minor bug fixes
-
-Revision 1.5 2000/06/08 13:34:50 fca
-Better control of momentum range in GenBox
-
-Revision 1.4 1999/09/29 09:24:29 fca
-Introduction of the Copyright and cvs Log
+/* $Id$ */
-*/
-
-///////////////////////////////////////////////////////////////////
-// //
-// Generate the final state of the interaction as the input //
-// to the MonteCarlo //
-//
+//-----------------------------------------------------------------
+// Generate the final state of the interaction as the input
+// to the MonteCarlo
+// Author: A.Morsch
+//-----------------------------------------------------------------
//Begin_Html
/*
<img src="picts/AliGeneratorClass.gif">
//End_Html
// //
///////////////////////////////////////////////////////////////////
-#include "TGenerator.h"
+#include <TGenerator.h>
+#include "AliCollisionGeometry.h"
#include "AliConfig.h"
#include "AliGenerator.h"
-#include "AliCollisionGeometry.h"
#include "AliRun.h"
#include "AliStack.h"
ClassImp(AliGenerator)
-TGenerator* AliGenerator::fgMCEvGen=0;
//_______________________________________________________________________
AliGenerator::AliGenerator():
+ fMCEvGen(0),
fThetaMin(0),
fThetaMax(0),
fPhiMin(0),
fChildWeight(0),
fAnalog(0),
fVertexSmear(kNoSmear),
- fVertexSource(kExternal),
+ fVertexSource(kInternal),
fCutVertexZ(0),
fTrackIt(0),
fOrigin(3),
fOsigma(3),
fVertex(3),
- fEventVertex(0),
fStack(0),
fCollisionGeometry(0)
{
//_______________________________________________________________________
AliGenerator::AliGenerator(Int_t npart):
+ fMCEvGen(0),
fThetaMin(0),
fThetaMax(0),
fPhiMin(0),
fChildWeight(0),
fAnalog(0),
fVertexSmear(kNoSmear),
- fVertexSource(kExternal),
+ fVertexSource(kInternal),
fCutVertexZ(0),
fTrackIt(0),
fOrigin(3),
fOsigma(3),
fVertex(3),
- fEventVertex(0),
fStack(0),
fCollisionGeometry(0)
{
AliGenerator::AliGenerator(const AliGenerator &gen):
TNamed(gen),
AliRndm(gen),
+ fMCEvGen(0),
fThetaMin(0),
fThetaMax(0),
fPhiMin(0),
fChildWeight(0),
fAnalog(0),
fVertexSmear(kNoSmear),
- fVertexSource(kExternal),
+ fVertexSource(kInternal),
fCutVertexZ(0),
fTrackIt(0),
fOrigin(3),
fOsigma(3),
fVertex(3),
- fEventVertex(0),
fStack(0)
{
//
//
fOrigin.Set(0);
fOsigma.Set(0);
- if (fgMCEvGen) {
- delete fgMCEvGen;
- fgMCEvGen=0;
+ if (fMCEvGen) {
+ delete fMCEvGen;
+ fMCEvGen=0;
}
}
//
if (fVertexSource == kInternal) {
VertexInternal();
- } else {
+ } else if (fVertexSource == kContainer) {
+ ;
+ } else if (fVertexSource == kExternal) {
VertexExternal();
}
}
// Obtain calculated vertex
// Default is gaussian smearing
Float_t random[6];
- Rndm(random,6);
- for (Int_t j = 0; j<3 ; j++) {
- fVertex[j]=
- fOrigin[j]+fOsigma[j]*TMath::Cos(2*random[2*j]*TMath::Pi())*
- TMath::Sqrt(-2*TMath::Log(random[2*j+1]));
+ Float_t dv[3];
+ Int_t j;
+
+ dv[2] = 1.e10;
+ while(TMath::Abs(dv[2]) > fCutVertexZ*fOsigma[2]) {
+ Rndm(random,6);
+ for (j=0; j < 3; j++) {
+ dv[j] = fOsigma[j]*TMath::Cos(2*random[2*j]*TMath::Pi())*
+ TMath::Sqrt(-2*TMath::Log(random[2*j+1]));
+ }
}
+ for (j=0; j < 3; j++) fVertex[j] = fOrigin[j] + dv[j];
+
}
//_______________________________________________________________________
-void AliGenerator::SetTrack(Int_t done, Int_t parent, Int_t pdg,
+void AliGenerator::PushTrack(Int_t done, Int_t parent, Int_t pdg,
Float_t *pmom, Float_t *vpos, Float_t *polar,
Float_t tof, TMCProcess mech, Int_t &ntr,
Float_t weight, Int_t is)
{
+ //
+ // Loads one track on the stack
+ //
if (fStack)
- fStack->SetTrack(done, parent, pdg, pmom, vpos, polar, tof,
+ fStack->PushTrack(done, parent, pdg, pmom, vpos, polar, tof,
mech, ntr, weight, is);
else
- gAlice->SetTrack(done, parent, pdg, pmom, vpos, polar, tof,
+ gAlice->PushTrack(done, parent, pdg, pmom, vpos, polar, tof,
mech, ntr, weight, is);
}
//_______________________________________________________________________
-void AliGenerator::SetTrack(Int_t done, Int_t parent, Int_t pdg,
+void AliGenerator::PushTrack(Int_t done, Int_t parent, Int_t pdg,
Double_t px, Double_t py, Double_t pz, Double_t e,
Double_t vx, Double_t vy, Double_t vz, Double_t tof,
Double_t polx, Double_t poly, Double_t polz,
TMCProcess mech, Int_t &ntr, Float_t weight, Int_t is)
{
+ //
+ // Loads one track on the stack
+ //
if (fStack)
- fStack->SetTrack(done, parent, pdg, px, py, pz, e, vx, vy, vz, tof,
+ fStack->PushTrack(done, parent, pdg, px, py, pz, e, vx, vy, vz, tof,
polx, poly, polz, mech, ntr, weight, is);
else
- gAlice->SetTrack(done, parent, pdg, px, py, pz, e, vx, vy, vz, tof,
+ gAlice->PushTrack(done, parent, pdg, px, py, pz, e, vx, vy, vz, tof,
polx, poly, polz, mech, ntr, weight, is);
}
//_______________________________________________________________________
void AliGenerator:: KeepTrack(Int_t itrack)
{
+ //
+ // Declare a track permanent on the stack
+ //
if (fStack)
fStack->KeepTrack(itrack);
else
//_______________________________________________________________________
void AliGenerator:: SetHighWaterMark(Int_t nt)
{
+ //
+ // Internal function to set the maximum index used in the stack
+ //
if (fStack)
fStack->SetHighWaterMark(nt);
else
gAlice->SetHighWaterMark(nt);
}
+void AliGenerator::FinishRun()
+{
+ ;
+}