Changes for #92479: Request to update ZDC initialization in macros/ConfigRaw2012.C
[u/mrichter/AliRoot.git] / test / genkine / gen / fastGen.C
index 9be11ce1b4c57e2d2f7f3d13c7b458a3a295f1e5..d053c91812e2ecd78d1acdc9d76491b91cf08a03 100644 (file)
@@ -11,6 +11,7 @@
 #include <TDatabasePDG.h>
 #include <TParticle.h>
 #include <TArrayI.h>
+#include <TTree.h>
 
 #include "AliGenerator.h"
 #include "AliPDG.h"
@@ -28,7 +29,7 @@ void GetFinalDecayProducts(Int_t ind, AliStack & stack , TArrayI & ar);
 void fastGen(Int_t nev = 1, char* filename = "galice.root")
 {
   AliPDG::AddParticlesToPdgDataBase();
-  TDatabasePDG*  pPdg = TDatabasePDG::Instance();
+  TDatabasePDG::Instance();
  
 
 
@@ -71,15 +72,15 @@ void fastGen(Int_t nev = 1, char* filename = "galice.root")
   // the initialization of the generator includes reading of the decay table.
 
   AliPythia * py= AliPythia::Instance();
-  py->SetMDME(731,1,0); //forbid D*+->D+ + pi0
-  py->SetMDME(732,1,0);//forbid D*+->D+ + gamma
+  py->SetMDME(737,1,0); //forbid D*+->D+ + pi0
+  py->SetMDME(738,1,0);//forbid D*+->D+ + gamma
 
   // Forbid all D0 decays except D0->K- pi+
-  for(Int_t d=741; d<=756; d++){ 
+  for(Int_t d=747; d<=762; d++){ 
     py->SetMDME(d,1,0);
   }
-  // decay 757 is D0->K- pi+
-  for(Int_t d=758; d<=801; d++){ 
+  // decay 763 is D0->K- pi+
+  for(Int_t d=764; d<=807; d++){ 
     py->SetMDME(d,1,0);
   }
 
@@ -102,7 +103,6 @@ void fastGen(Int_t nev = 1, char* filename = "galice.root")
     //  Generate event
     Int_t nprim = 0;
     Int_t ntrial = 0;
-    Int_t minmult = 1000;
     Int_t ndstar = 0;
    
    
@@ -110,9 +110,9 @@ void fastGen(Int_t nev = 1, char* filename = "galice.root")
     //-------------------------------------------------------------------------------------
 
     while(!ndstar) {
-      // Selection of events with multiplicity
-      // bigger than "minmult"
+      // Selection of events with D*
       stack->Reset();
+      stack->ConnectTree(rl->TreeK());
       gener->Generate();
       ntrial++;
       nprim = stack->GetNprimary();
@@ -176,8 +176,6 @@ void fastGen(Int_t nev = 1, char* filename = "galice.root")
   //                         Termination
   //  Generator
   gener->FinishRun();
-  //  Stack
-  stack->FinishRun();
   //  Write file
   rl->WriteHeader("OVERWRITE");
   gener->Write();