--- PamVMC/include/PamVMCPrimaryGenerator.h	2009/06/12 18:39:19	1.5
+++ PamVMC/include/PamVMCPrimaryGenerator.h	2010/09/15 07:05:42	1.6
@@ -88,6 +88,13 @@
     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); 
+    void GenSpe_Flat(Double_t PEmin, Double_t PEmax, Double_t gamma, Bool_t isEnergy=kFALSE); 
+    void GenSpe_3par(Double_t PEmin, Double_t PEmax, Double_t a, Double_t b, Double_t c);
+    void GenSphericalPhiThe(); // cecilia // flusso istropo // sets phi the
+    void GenSphPhiThe(Double_t xmin, Double_t xmax, Double_t ymin, Double_t ymax, 
+		      Double_t zmin, Double_t zmax); // flusso isotropo entro S2 S3
+    //                                                  sets position and phi the
+    
 
     // get methods
     Int_t GetParticle(){ return fprim.fPDG; };
@@ -127,6 +134,7 @@
     Double_t KinEToMomentum(Double_t E0);
     Double_t MomentumToRig(Double_t P0) { return P0/fcharge; };
     Double_t RigToMomentum(Double_t R0){ return R0*fcharge; };
+    Double_t function3par(Double_t xx, Double_t a, Double_t b, Double_t c){return  5.*pow((xx + b * exp(-c * sqrt(xx))),-a);};
 
     // data members
     TVirtualMCStack*  fStack;