+ Int_t AliGenPHOSlib::IpChargedPion(TRandom *ran)
+{
+// particle composition pi+, pi0, pi-
+//
+
+ Float_t random = ran->Rndm();
+
+ if ( (2.*random) < 1. )
+ {
+ return 211 ;
+ }
+ else
+ {
+ return -211 ;
+ }
+}
+
+//End Pions
+//======================================================================
+// Pi 0 Flat Distribution
+// Transverse momentum distribution PtPi0Flat
+// Rapidity distribution YPi0Flat
+// Particle distribution IdPi0Flat 111 (pi0)
+//
+
+Double_t AliGenPHOSlib::PtPi0Flat(Double_t */*px*/, Double_t *)
+{
+// Pion transverse momentum flat distribution
+
+return 1;
+
+}
+
+Double_t AliGenPHOSlib::YPi0Flat( Double_t */*py*/, Double_t *)
+{
+
+// pion y-distribution
+//
+ return 1.;
+}
+
+ Int_t AliGenPHOSlib::IpPi0Flat(TRandom *)
+{
+
+// particle composition pi0
+//
+ return 111 ;
+}
+// End Pi0Flat