correcting printf format specifiers
[u/mrichter/AliRoot.git] / PYTHIA8 / pythia8.C
1 AliGenerator*  CreateGenerator();
2
3 void pythia8(Int_t nev = 1, char* filename = "galice.root")
4 {
5 //  Runloader
6     gSystem->Load("liblhapdf.so");    
7     gSystem->Load("libEGPythia6.so"); 
8     gSystem->Load("libpythia6.so");   
9     gSystem->Load("libAliPythia6.so");
10     gSystem->Load("libpythia8.so");   
11     gSystem->Load("libEGPythia8.so"); 
12     gSystem->Load("libAliPythia8.so");
13     
14     AliRunLoader* rl = AliRunLoader::Open("galice.root","FASTRUN","recreate");
15     
16     rl->SetCompressionLevel(2);
17     rl->SetNumberOfEventsPerFile(nev);
18     rl->LoadKinematics("RECREATE");
19     rl->MakeTree("E");
20     gAlice->SetRunLoader(rl);
21
22 //  Create stack
23     rl->MakeStack();
24     AliStack* stack      = rl->Stack();
25  
26 //  Header
27     AliHeader* header = rl->GetHeader();
28 //
29 //  Create and Initialize Generator
30     AliGenerator *gener = CreateGenerator();
31     gener->Init();
32     gener->SetStack(stack);
33     
34 //
35 //                        Event Loop
36 //
37     Int_t iev;
38      
39     for (iev = 0; iev < nev; iev++) {
40
41         printf("\n \n Event number %d \n \n", iev);
42         
43 //  Initialize event
44         header->Reset(0,iev);
45         rl->SetEventNumber(iev);
46         stack->Reset();
47         rl->MakeTree("K");
48 //      stack->ConnectTree();
49     
50 //  Generate event
51         gener->Generate();
52 //  Analysis
53         Int_t npart = stack->GetNprimary();
54         printf("Analyse %d Particles\n", npart);
55         for (Int_t part=0; part<npart; part++) {
56             TParticle *MPart = stack->Particle(part);
57             Int_t mpart  = MPart->GetPdgCode();
58 //          printf("Particle %d\n", mpart);
59         }
60         
61 //  Finish event
62         header->SetNprimary(stack->GetNprimary());
63         header->SetNtrack(stack->GetNtrack());  
64 //      I/O
65 //      
66         stack->FinishEvent();
67         header->SetStack(stack);
68         rl->TreeE()->Fill();
69         rl->WriteKinematics("OVERWRITE");
70
71     } // event loop
72 //
73 //                         Termination
74 //  Generator
75     gener->FinishRun();
76 //  Write file
77     rl->WriteHeader("OVERWRITE");
78     gener->Write();
79     rl->Write();
80     
81 }
82
83
84 AliGenerator*  CreateGenerator()
85 {
86     AliGenPythiaPlus* gener = new AliGenPythiaPlus(AliPythia8::Instance());
87 //
88 //
89     gener->SetProcess(kPyMb);
90 //   Centre of mass energy 
91     gener->SetEnergyCMS(14000.);
92 //   Initialize generator    
93     gener->SetEventListRange(-1, 10)
94     return gener;
95 }
96
97
98
99
100
101
102
103