Rotation matrix in BuildGeometry has been changed to rotx999
[u/mrichter/AliRoot.git] / EVGEN / AliGenPythia.h
CommitLineData
a0134e39 1#ifndef ALIGENPYTHIA_H
2#define ALIGENPYTHIA_H
3da30618 3/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
5
6/* $Id$ */
7
f87cfe57 8
fe4da5cc 9#include "AliGenerator.h"
f1a48a38 10#include "AliPythia.h"
11#include "AliDecayer.h"
f87cfe57 12#include <TArrayI.h>
13
14class AliPythia;
15class TParticle;
fe4da5cc 16
17class AliGenPythia : public AliGenerator
18{
fe4da5cc 19 public:
20 AliGenPythia();
21 AliGenPythia(Int_t npart);
f87cfe57 22 AliGenPythia(const AliGenPythia &Pythia);
fe4da5cc 23 virtual ~AliGenPythia();
24 virtual void Generate();
25 virtual void Init();
26 // select process type
f1a48a38 27 virtual void SetProcess(Process_t proc = kPyCharm) {fProcess = proc;}
fe4da5cc 28 // select structure function
f1a48a38 29 virtual void SetStrucFunc(StrucFunc_t func = kGRV_HO) {fStrucFunc = func;}
fe4da5cc 30 // select pt of hard scattering
f1a48a38 31 virtual void SetPtHard(Float_t ptmin = 0, Float_t ptmax = 1.e10)
32 {fPtHardMin = ptmin; fPtHardMax = ptmax; }
fe4da5cc 33 // set centre of mass energy
f1a48a38 34 virtual void SetEnergyCMS(Float_t energy = 5500) {fEnergyCMS = energy;}
fe4da5cc 35 // force decay type
f1a48a38 36 virtual void SetForceDecay(Decay_t decay = kSemiMuonic) {fForceDecay = decay;}
811826d8 37 // treat protons as inside nuclei
38 virtual void SetNuclei(Int_t a1, Int_t a2);
fe4da5cc 39 // get cross section of process
40 virtual Float_t GetXsection() {return fXsection;}
a8228d85 41 // Check PDG code
f1a48a38 42 virtual Int_t CheckPDGCode(Int_t pdgcode);
9faa0850 43 virtual void FinishRun();
44
f87cfe57 45 // Assignment Operator
46 AliGenPythia & operator=(const AliGenPythia & rhs);
a0134e39 47 protected:
48 Process_t fProcess; // Process type
49 StrucFunc_t fStrucFunc; // Structure Function
50 Decay_t fForceDecay; // Decay channel are forced
51 Float_t fEnergyCMS; // Centre of mass energy
52 Float_t fKineBias; // Bias from kinematic selection
53 Int_t fTrials; // Number of trials
54 TArrayI fParentSelect; // Parent particles to be selected
55 TArrayI fChildSelect; // Decay products to be selected
56 Float_t fXsection; // Cross-section
a8a6107b 57 AliPythia *fPythia; //! Pythia
a0134e39 58 Float_t fPtHardMin; // lower pT-hard cut
59 Float_t fPtHardMax; // higher pT-hard cut
811826d8 60 Int_t fNucA1; // mass number nucleus side 1
61 Int_t fNucA2; // mass number nucleus side 2
f1a48a38 62 Bool_t fFullEvent; // Write Full event if true
c78d7f8c 63 AliDecayer *fDecayer; // ! pointer to the decayer instance
a0134e39 64 private:
65 // check if particle is selected as parent particle
66 Bool_t ParentSelected(Int_t ip);
67 // check if particle is selected as child particle
68 Bool_t ChildSelected(Int_t ip);
69 // all kinematic selection cuts go here
70 Bool_t KinematicSelection(TParticle *particle);
71 // adjust the weight from kinematic cuts
72 void AdjustWeights();
73
f87cfe57 74 ClassDef(AliGenPythia,1) // AliGenerator interface to Pythia
fe4da5cc 75};
76#endif
77
78
79
80
81