--- trieste/pamVMC/include/PamVMCPrimaryGenerator.h	2009/03/04 12:51:14	1.1
+++ trieste/pamVMC/include/PamVMCPrimaryGenerator.h	2009/03/24 14:04:01	1.2
@@ -29,7 +29,7 @@
 
   Int_t fPDG;
   Double_t fX0, fY0, fZ0;
-  Double_t fTHETA, fPHI;
+  Double_t fTHETA, fPHI;  // RADIANTS
   Double_t fP0;
   Bool_t fGOOD;
 
@@ -41,7 +41,7 @@
     { cout<<"PRIMARY particle infromation:"<<endl;
       cout<<"Pdg="<<fPDG<<endl; //add name later
       cout<<"Position: "<<"("<<fX0<<","<<fY0<<","<<fZ0<<")"<<endl;
-      cout<<"P0, Theta, Phi: "<<fP0<<","<<fTHETA<<","<<fPHI<<endl;
+      cout<<"P0, Theta, Phi: "<<fP0<<","<<fTHETA*180./3.1415<<","<<fPHI*180./3.1415<<endl;
       cout<<"GOOD Single Track:"<<fGOOD<<endl;
     }
 
@@ -53,7 +53,8 @@
 class PamVMCPrimaryGenerator : public TObject
 {
   public:
-    PamVMCPrimaryGenerator(TVirtualMCStack* stack); 
+    PamVMCPrimaryGenerator(TVirtualMCStack* stack, UInt_t seed); 
+    //PamVMCPrimaryGenerator(UInt_t seed);
     PamVMCPrimaryGenerator();
 
     virtual ~PamVMCPrimaryGenerator();
@@ -83,15 +84,19 @@
     void GenSphDist(Double_t, Double_t, Double_t, Double_t, Double_t);
     
     void GenDirection(Double_t thetamin, Double_t thetamax, Double_t phimin, Double_t phimax){
-      SetDirection(frnd->Uniform(thetamin,thetamax),frnd->Uniform(phimin,phimax)); }
+	SetDirection(frnd->Uniform(thetamin,thetamax),frnd->Uniform(phimin,phimax)); }
 
     //flat spectra generator
     void GenSpe(Double_t PEmin, Double_t PEmax, Bool_t isEnergy=kFALSE); 
     //power law spectra, gamma - differential spectral index
     void GenSpe(Double_t PEmin, Double_t PEmax, Double_t gamma, Bool_t isEnergy=kFALSE); 
 
-	//spherical dist. from radius r, allow only those particles going to cross 3 TOFs
-	void GenSphDist(Double_t);
+    //spherical dist. from radius r, allow only those particles going to cross 3 TOFs
+    void GenSphDist(Double_t);
+
+
+    void GenSphPhiThe(Double_t , Double_t , Double_t , Double_t , 
+		      Double_t , Double_t ); // cecilia
 
     // get methods
     Int_t GetParticle(){ return fprim.fPDG; };
@@ -139,7 +144,7 @@
 
     TRandom3 *        frnd;
     
-  ClassDef(PamVMCPrimaryGenerator,1)  //PamVMCPrimaryGenerator
+  ClassDef(PamVMCPrimaryGenerator,2)  //PamVMCPrimaryGenerator
 };
 
 // inline functions