Fixed decay modes
[u/mrichter/AliRoot.git] / test / genkine / gen / fastGen.C
index 9be11ce1b4c57e2d2f7f3d13c7b458a3a295f1e5..f02a8ec78061e6100cd27590e0a7aae15a148e98 100644 (file)
@@ -28,7 +28,7 @@ void GetFinalDecayProducts(Int_t ind, AliStack & stack , TArrayI & ar);
 void fastGen(Int_t nev = 1, char* filename = "galice.root")
 {
   AliPDG::AddParticlesToPdgDataBase();
 void fastGen(Int_t nev = 1, char* filename = "galice.root")
 {
   AliPDG::AddParticlesToPdgDataBase();
-  TDatabasePDG*  pPdg = TDatabasePDG::Instance();
+  TDatabasePDG::Instance();
  
 
 
  
 
 
@@ -71,15 +71,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();
   // 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+
 
   // 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);
   }
     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);
   }
 
     py->SetMDME(d,1,0);
   }
 
@@ -102,7 +102,6 @@ void fastGen(Int_t nev = 1, char* filename = "galice.root")
     //  Generate event
     Int_t nprim = 0;
     Int_t ntrial = 0;
     //  Generate event
     Int_t nprim = 0;
     Int_t ntrial = 0;
-    Int_t minmult = 1000;
     Int_t ndstar = 0;
    
    
     Int_t ndstar = 0;
    
    
@@ -110,9 +109,9 @@ void fastGen(Int_t nev = 1, char* filename = "galice.root")
     //-------------------------------------------------------------------------------------
 
     while(!ndstar) {
     //-------------------------------------------------------------------------------------
 
     while(!ndstar) {
-      // Selection of events with multiplicity
-      // bigger than "minmult"
+      // Selection of events with D*
       stack->Reset();
       stack->Reset();
+      stack->ConnectTree(rl->TreeK());
       gener->Generate();
       ntrial++;
       nprim = stack->GetNprimary();
       gener->Generate();
       ntrial++;
       nprim = stack->GetNprimary();
@@ -176,8 +175,6 @@ void fastGen(Int_t nev = 1, char* filename = "galice.root")
   //                         Termination
   //  Generator
   gener->FinishRun();
   //                         Termination
   //  Generator
   gener->FinishRun();
-  //  Stack
-  stack->FinishRun();
   //  Write file
   rl->WriteHeader("OVERWRITE");
   gener->Write();
   //  Write file
   rl->WriteHeader("OVERWRITE");
   gener->Write();