I/O enabled.

[u/mrichter/AliRoot.git] / FASTSIM / fastGen.C

diff --git a/FASTSIM/fastGen.C b/FASTSIM/fastGen.C
index 35b1bb9bc89ddd38c9908402ea0b6d4d2648253e..e68598298892fa4a98a17507d3a46bc79a329437 100644 (file)
--- a/FASTSIM/fastGen.C
+++ b/FASTSIM/fastGen.C
@@ -34,6 +34,7 @@ void fastGen(Int_t nev = 1, char* filename = "galice.root")
        
 //  Initialize event
        header->Reset(0,iev);
        
 //  Initialize event
        header->Reset(0,iev);

+       stack->Reset();

        stack->BeginEvent(iev);
 
 //  Generate event
        stack->BeginEvent(iev);
 
 //  Generate event


@@ -52,12 +53,10 @@ void fastGen(Int_t nev = 1, char* filename = "galice.root")
        header->SetNtrack(stack->GetNtrack());  
 //      I/O
 //     
        header->SetNtrack(stack->GetNtrack());  
 //      I/O
 //     

-//     stack->FinishEvent();
-//     header->SetStack(stack);
-//     treeE->Fill();
-//     Reset stack
-
-       stack->Reset();
+       stack->FinishEvent();
+       header->SetStack(stack);
+       treeE->Fill();
+       (stack->TreeK())->Write(0,TObject::kOverwrite);

     } // event loop
 //
 //                         Termination
     } // event loop
 //
 //                         Termination


@@ -86,9 +85,9 @@ AliGenerator*  CreateGenerator()
 //   charm, beauty, charm_unforced, beauty_unforced, jpsi, jpsi_chi, mb
     gener->SetProcess(kPyJets);
 //   Centre of mass energy 
 //   charm, beauty, charm_unforced, beauty_unforced, jpsi, jpsi_chi, mb
     gener->SetProcess(kPyJets);
 //   Centre of mass energy 

-    gener->SetEnergyCMS(14000.);
+    gener->SetEnergyCMS(5500.);

 //   Pt transfer of the hard scattering
 //   Pt transfer of the hard scattering

-    gener->SetPtHard(5.,5.1);
+    gener->SetPtHard(50.,50.2);

 //   Initialize generator    
     return gener;
 }
 //   Initialize generator    
     return gener;
 }