FASTSIM/fastGen.C

   1 AliGenerator*  CreateGenerator();
   2
   3 void fastGen(Int_t nev = 1, char* filename = "galice.root")
   4 {
   5 //  Runloader
   6
   7     AliRunLoader* rl = AliRunLoader::Open("galice.root","FASTRUN","recreate");
   8
   9     rl->SetCompressionLevel(2);
  10     rl->SetNumberOfEventsPerFile(nev);
  11     rl->LoadKinematics("RECREATE");
  12     rl->MakeTree("E");
  13     gAlice->SetRunLoader(rl);
  14
  15 //  Create stack
  16     rl->MakeStack();
  17     AliStack* stack      = rl->Stack();
  18
  19 //  Header
  20     AliHeader* header = rl->GetHeader();
  21 //
  22 //  Create and Initialize Generator
  23     AliGenerator *gener = CreateGenerator();
  24     gener->Init();
  25     gener->SetStack(stack);
  26
  27 //
  28 //                        Event Loop
  29 //
  30     Int_t iev;
  31
  32     for (iev = 0; iev < nev; iev++) {
  33
  34         printf("\n \n Event number %d \n \n", iev);
  35
  36 //  Initialize event
  37         header->Reset(0,iev);
  38         rl->SetEventNumber(iev);
  39         stack->Reset();
  40         rl->MakeTree("K");
  41 //      stack->ConnectTree();
  42
  43 //  Generate event
  44         gener->Generate();
  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         rl->TreeE()->Fill();
  62         rl->WriteKinematics("OVERWRITE");
  63
  64     } // event loop
  65 //
  66 //                         Termination
  67 //  Generator
  68     gener->FinishRun();
  69 //  Write file
  70     rl->WriteHeader("OVERWRITE");
  71     gener->Write();
  72     rl->Write();
  73
  74 }
  75
  76
  77 AliGenerator*  CreateGenerator()
  78 {
  79     gener = new AliGenPythia(1);
  80 //
  81 //
  82 //   vertex position and smearing
  83     gener->SetVertexSmear(kPerEvent);
  84 //   structure function
  85     gener->SetStrucFunc(kGRVHO);
  86 //   charm, beauty, charm_unforced, beauty_unforced, jpsi, jpsi_chi, mb
  87     gener->SetProcess(kPyJets);
  88 //   Centre of mass energy
  89     gener->SetEnergyCMS(5500.);
  90 //   Pt transfer of the hard scattering
  91     gener->SetPtHard(50.,50.2);
  92 //   Initialize generator
  93     return gener;
  94 }
  95
  96
  97
  98
  99
 100
 101
 102
 103
 104
 105
 106