Charm pt-distribution corrected. More realistic y-distribution for pi and K.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 27 Mar 2001 11:01:04 +0000 (11:01 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 27 Mar 2001 11:01:04 +0000 (11:01 +0000)
EVGEN/AliGenMUONlib.cxx

index 2ec6a0563ff635cb48d30ea2fa953057d3ede44d..1b55475a3780bfd5b6b659373850dc51d8ae9744 100644 (file)
 
 /*
 $Log$
 
 /*
 $Log$
+Revision 1.12  2001/03/09 13:01:41  morsch
+- enum constants for paramterisation type (particle family) moved to AliGen*lib.h
+- use AliGenGSIlib::kUpsilon, AliGenPHOSlib::kEtaPrime to access the constants
+
 Revision 1.11  2000/11/30 07:12:50  alibrary
 Introducing new Rndm and QA classes
 
 Revision 1.11  2000/11/30 07:12:50  alibrary
 Introducing new Rndm and QA classes
 
@@ -77,13 +81,15 @@ Double_t AliGenMUONlib::PtPion(Double_t *px, Double_t *dummy)
 Double_t AliGenMUONlib::YPion( Double_t *py, Double_t *dummy)
 {
 // Pion y
 Double_t AliGenMUONlib::YPion( Double_t *py, Double_t *dummy)
 {
 // Pion y
+  Double_t y=TMath::Abs(*py);
+/*
   const Double_t ka    = 7000.;
   const Double_t kdy   = 4.;
   const Double_t ka    = 7000.;
   const Double_t kdy   = 4.;
-
-  Double_t y=TMath::Abs(*py);
-  //
   Double_t ex = y*y/(2*kdy*kdy);
   return ka*TMath::Exp(-ex);
   Double_t ex = y*y/(2*kdy*kdy);
   return ka*TMath::Exp(-ex);
+*/
+  return 1.16526e+04+y*-3.79886e+03+y*y*4.31130e+02;
+  
 }
 //                 particle composition
 //
 }
 //                 particle composition
 //
@@ -133,14 +139,16 @@ Double_t AliGenMUONlib::PtKaon( Double_t *px, Double_t *dummy)
 Double_t AliGenMUONlib::YKaon( Double_t *py, Double_t *dummy)
 {
 // Kaon y
 Double_t AliGenMUONlib::YKaon( Double_t *py, Double_t *dummy)
 {
 // Kaon y
+  Double_t y=TMath::Abs(*py);
+/*
   const Double_t ka    = 1000.;
   const Double_t kdy   = 4.;
   const Double_t ka    = 1000.;
   const Double_t kdy   = 4.;
-  
-
-  Double_t y=TMath::Abs(*py);
   //
   Double_t ex = y*y/(2*kdy*kdy);
   return ka*TMath::Exp(-ex);
   //
   Double_t ex = y*y/(2*kdy*kdy);
   return ka*TMath::Exp(-ex);
+*/
+
+  return 1.16526e+04+y*-3.79886e+03+y*y*4.31130e+02;
 }
 
 //                 particle composition
 }
 
 //                 particle composition
@@ -273,9 +281,10 @@ Double_t AliGenMUONlib::PtCharm( Double_t *px, Double_t *dummy)
 // Charm pT
   const Double_t kpt0 = 4.08;
   const Double_t kxn  = 9.40;
 // Charm pT
   const Double_t kpt0 = 4.08;
   const Double_t kxn  = 9.40;
+
   Double_t x=*px;
   //
   Double_t x=*px;
   //
-  Double_t pass1 = 1.+(x/kpt0)*(x/kpt0);
+  Double_t pass1 = 1.+(x/kpt0);
   return x/TMath::Power(pass1,kxn);
 }
 //                  y-distribution
   return x/TMath::Power(pass1,kxn);
 }
 //                  y-distribution