--- trieste/pamVMC/run_g4_testTRK.C	2009/03/04 12:51:04	1.1
+++ trieste/pamVMC/run_g4_testTRK.C	2009/03/24 14:03:55	1.2
@@ -1,3 +1,4 @@
+
 // $Id: run_g4.C,v 1.0 2007/06/01 
 //
 // Macro for running G4PAM_VMC with Geant4. 
@@ -27,6 +28,7 @@
     PamVMCApplication* appl 
       = new PamVMCApplication("PAMG4_VMC", "PAMELA GEANT4 VMC application");
     
+
     //appl->GetPrimaryGenerator()->SetNofPrimaries(1);  
         appl->SetVerboseLevel(0);
       
@@ -37,6 +39,11 @@
     // Initialize MC
     appl->InitMC(configMacro);
 
+    // Cecilia: set seed for random generator  
+    TRandom3 *rndgen = new TRandom3(123);
+    ((TGeant4*)gMC)->SetRandom(rndgen);
+
+
     // Customise Geant4 setting after initialization:
     // Physics list
 //     gMC->SetCut("CUTGAM", 1.0e-6);
@@ -53,27 +60,34 @@
     
     ((TGeant4*)gMC)->ProcessGeantMacro("g4config2.in");
     // Visualization settings
-     ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in");
+    ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in");
  
    ((TGeant4*)gMC)->SetMaxNStep(50000);
    //appl->GetPrimaryGenerator()->SetParticle(1000060120);
     appl->GetPrimaryGenerator()->SetParticle(kElectron);
     //appl->GetPrimaryGenerator()->SetDirection(0.,0.);
-     //appl->GetPrimaryGenerator()->SetKinEnergy(10.);
+    appl->GetPrimaryGenerator()->SetKinEnergy(10.);
      //appl->GetPrimaryGenerator()->SetRigidity(1.0); //GV 
     //appl->GetPrimaryGenerator()->SetPosition(1.,1.,130.);
+
     // Run MC
-    for (Int_t i=0; i<3; i++){
-      cout<<"Event "<<i<<endl;
-     appl->GetPrimaryGenerator()->SetKinEnergy(10.);
-     if (i==0) appl->GetPrimaryGenerator()->SetKinEnergy(0.01);
-       appl->GetPrimaryGenerator()->GenSphDist(120.);
-       //appl->GetPrimaryGenerator()->GenSphDist(120.,0.,90.,0.,360.);
-      // appl->GetPrimaryGenerator()->GenPosition(-22.4,22.4, -18.5, 18.5, 110.0,110.0);
-      //appl->GetPrimaryGenerator()->GenDirection(0.,TMath::Pi()/2.,0.,2*TMath::Pi());
-      //appl->GetPrimaryGenerator()->GenSpe(1.,5.,-2.75);
+    for (Int_t i=0; i<2; i++){
+      cout<<"in Loop: Event nr."<<i<<" --<<--<<--<<--<<--<<--<<--"<<endl;
+      // ritabrata 04.03.2009
+      //appl->GetPrimaryGenerator()->SetKinEnergy(0.1);
+      //if (i==0) appl->GetPrimaryGenerator()->SetKinEnergy(0.01);
+      //appl->GetPrimaryGenerator()->GenSphDist(120.);
+      // endritabrata
+
+      // cecilia
+      //appl->GetPrimaryGenerator()->GenPosition(-22.4,22.4, -18.5, 18.5, 110.0,110.0);
+      // appl->GetPrimaryGenerator()->GenPosition(-8.,8., -6., 6., 110.0,110.0);
+      appl->GetPrimaryGenerator()->GenSphPhiThe();
+      //appl->GetPrimaryGenerator()->GenDirection(0.,0.1, 0.,2*TMath::Pi()); // 5 gradi
+      // appl->GetPrimaryGenerator()->GenSpe(7.,20.,-2.75);
       //  appl->GetPrimaryGenerator()->Print();
-      // if (i==9) ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in");
+      // if (i==9) ((TGeant4*)gMC)->ProcessGeantMacro("g4vis.in"); 
+      appl->GetPrimaryGenerator()->Print();
       appl->RunMC(1);
 
     }