]>
Commit | Line | Data |
---|---|---|
33de81b8 | 1 | AliGenerator* CreateGenerator(); |
2 | ||
3 | void fastGenPA(Int_t nev = 1, char* filename = "galice.root") | |
4 | { | |
5 | // | |
6 | // Construction | |
7 | // | |
8 | // Output file | |
9 | TFile* file = new TFile(filename, "recreate"); | |
10 | // Create stack | |
11 | AliStack* stack = new AliStack(10000); | |
12 | stack->MakeTree(0, filename); | |
13 | ||
14 | // Create Header | |
15 | AliHeader* header = new AliHeader(); | |
16 | // Create Header Tree | |
17 | TTree* treeE = new TTree("TE","Headers"); | |
18 | treeE->Branch("Header", "AliHeader", &header, 4000, 0); | |
19 | treeE->Write(); | |
20 | // | |
21 | // Create and Initialize Generator | |
22 | AliGenerator *gener = CreateGenerator(); | |
23 | AliPythia* pyth = AliPythia::Instance(); | |
24 | ||
25 | gener->Init(); | |
26 | gener->SetStack(stack); | |
27 | ||
28 | // | |
29 | // Event Loop | |
30 | // | |
31 | Int_t iev; | |
32 | ||
33 | for (iev = 0; iev < nev; iev++) { | |
34 | ||
35 | printf("\n \n Event number %d \n \n", iev); | |
36 | ||
37 | // Initialize event | |
38 | header->Reset(0,iev); | |
39 | stack->Reset(); | |
40 | stack->BeginEvent(iev); | |
41 | ||
42 | // Generate event | |
43 | gener->Generate(); | |
44 | ||
45 | // Analysis | |
46 | Int_t npart = stack->GetNprimary(); | |
47 | // printf("Analyse %d Particles\n", npart); | |
48 | for (Int_t part=0; part<npart; part++) { | |
49 | TParticle *MPart = stack->Particle(part); | |
50 | Int_t mpart = MPart->GetPdgCode(); | |
51 | printf("Particle %d\n", mpart); | |
52 | } | |
53 | ||
54 | // Finish event | |
55 | header->SetNprimary(stack->GetNprimary()); | |
56 | header->SetNtrack(stack->GetNtrack()); | |
57 | // I/O | |
58 | // | |
59 | // stack->FinishEvent(); | |
60 | // header->SetStack(stack); | |
61 | // treeE->Fill(); | |
62 | // (stack->TreeK())->Write(0,TObject::kOverwrite); | |
63 | } // event loop | |
64 | // | |
65 | // Termination | |
66 | // Generator | |
67 | printf("Calling Finish Run \n"); | |
68 | ||
69 | gener->FinishRun(); | |
70 | // Header | |
71 | treeE->Write(0,TObject::kOverwrite); | |
72 | delete treeE; treeE = 0; | |
73 | // Stack | |
74 | stack->FinishRun(); | |
75 | // Write file | |
76 | gener->Write(); | |
77 | file->Write(); | |
78 | } | |
79 | ||
80 | ||
81 | AliGenerator* CreateGenerator() | |
82 | { | |
83 | AliGenCocktail *gener = new AliGenCocktail(); | |
84 | gener->SetTrackingFlag(0); | |
85 | AliGenHijing *hijing = new AliGenHijing(-1); | |
86 | // centre of mass energy | |
87 | hijing->SetEnergyCMS(TMath::Sqrt(82./208.) * 14000.); | |
88 | // impact parameter range | |
89 | hijing->SetImpactParameterRange(0., 6.); | |
90 | // reference frame | |
91 | hijing->SetReferenceFrame("CMS"); | |
92 | hijing->SetBoostLHC(1); | |
93 | // projectile | |
94 | hijing->SetProjectile("P", 1, 1); | |
95 | hijing->SetTarget ("A", 208, 82); | |
96 | // tell hijing to keep the full parent child chain | |
97 | hijing->KeepFullEvent(); | |
98 | // enable jet quenching | |
99 | hijing->SetJetQuenching(4); | |
100 | // enable shadowing | |
101 | hijing->SetShadowing(1); | |
102 | // Don't track spectators | |
103 | hijing->SetSpectators(0); | |
104 | // kinematic selection | |
105 | hijing->SetSelectAll(0); | |
106 | // | |
107 | AliGenSlowNucleons* gray = new AliGenSlowNucleons(1); | |
108 | AliSlowNucleonModel* model = new AliSlowNucleonModelExp(); | |
109 | gray->SetSlowNucleonModel(model); | |
110 | gray->SetDebug(1); | |
111 | ||
112 | ||
113 | gener->AddGenerator(hijing,"Hijing pPb", 1); | |
114 | gener->AddGenerator(gray, "Gray Particles",1); | |
115 | ||
116 | return gener; | |
117 | } | |
118 | ||
119 | ||
120 | ||
121 | ||
122 | ||
123 | ||
124 | ||
125 | ||
126 | ||
127 | ||
128 | ||
129 |