Introduction of the D+ and Ds (and the -Xi_c case)
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 21 May 2012 15:54:18 +0000 (15:54 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Mon, 21 May 2012 15:54:18 +0000 (15:54 +0000)
Annalisa Mastroserio

EVGEN/AliGenITSULib.cxx
EVGEN/AliGenITSULib.h

index 3e00576..f36dbdf 100644 (file)
@@ -67,10 +67,12 @@ GenFunc AliGenITSULib::GetPt(Int_t iPID, const char * sForm) const
 
   switch(TMath::Abs(iPID)) {
 
-   case kLb   :    func=PtLbDist; break;  
-   case kLc   :    func=PtLcDist; break;  
-   case kXi_c :    func=PtLcDist; break;
-   case kB    :    func=PtLbDist; break;
+   case kLb    :    func=PtLbDist; break;  
+   case kLc    :    func=PtLcDist; break;  
+   case kXi_c  :    func=PtLcDist; break;
+   case kB     :    func=PtLbDist; break;
+   case kDs    :    func=PtLcDist; break;
+   case kDplus :    func=PtLcDist; break;
    default : AliError(Form("Unknown particle type: %i, Pt dist is 0",iPID));      func=0;
   } 
  }else {
@@ -85,7 +87,7 @@ GenFunc AliGenITSULib::GetY(Int_t iPID, const char *sForm) const
 {
  GenFunc func;
 
- if(TMath::Abs(iPID) != kLc && TMath::Abs(iPID) != kLb && TMath::Abs(iPID) != kXi_c && TMath::Abs(iPID) != kB) {
+ if(TMath::Abs(iPID) != kLc && TMath::Abs(iPID) != kLb && TMath::Abs(iPID) != kXi_c && TMath::Abs(iPID) != kB && TMath::Abs(iPID)!=kDplus && TMath::Abs(iPID)!=kDs) {
   AliError(Form("Unknown PID: %i, form: %s, returning 0",iPID,sForm));   //////////    
   func=0;
  } else { 
@@ -100,18 +102,23 @@ GenFuncIp AliGenITSULib::GetIp(Int_t iPID, const char *sForm) const
  // Return pointer to particle type parameterisation
  GenFuncIp id;
 
- if(TMath::Abs(iPID) != kLc && TMath::Abs(iPID) != kLb && TMath::Abs(iPID) != kXi_c && TMath::Abs(iPID) != kB) {
+ if(TMath::Abs(iPID) != kLc && TMath::Abs(iPID) != kLb && TMath::Abs(iPID) != kXi_c && TMath::Abs(iPID) != kB && TMath::Abs(iPID)!=kDplus && TMath::Abs(iPID)!=kDs) {
   AliError(Form("Unknown PID: %i, form: %s, return 0",iPID,sForm));   //////////       
   id = 0;
  } else {
   switch (iPID){
-   case kLc:                                                 return id=IpLcPlus;
-   case -kLc:                                                return id=IpLcMinus;
-   case kLb:                                                 return id=IpLb;
-   case -kLb:                                                return id=IpLbBar;
-   case kXi_c:                                               return id=IpXic;
-   case kB:                                                  return id=IpBPlus;
-   case -kB:                                                 return id=IpBMinus;
+   case kLc    :                                  return id=IpLcPlus;
+   case -kLc   :                                  return id=IpLcMinus;
+   case kLb    :                                  return id=IpLb;
+   case -kLb   :                                  return id=IpLbBar;
+   case kXi_c  :                                  return id=IpXic;
+   case -kXi_c :                                  return id=IpXicBar;
+   case kB     :                                  return id=IpBPlus;
+   case -kB    :                                  return id=IpBMinus;
+   case kDs    :                                  return id=IpDsPlus;
+   case -kDs   :                                  return id=IpDsMinus;
+   case kDplus :                                  return id=IpDPlus;
+   case -kDplus:                                  return id=IpDMinus;
    default  : AliFatal(Form("Unknown particle type: %i",iPID));  id=0;
   }
 
index a87da05..2a1d335 100644 (file)
@@ -17,7 +17,7 @@ class AliGenITSULib :public AliGenLib {
 
  public:
 
-  enum EPartId {kLb=5122,kLc=4122,kXi_c = 4232, kB = 521};
+  enum EPartId {kLb=5122,kLc=4122,kXi_c = 4232,kB = 521,kDs=431,kDplus=411};
 
   //Getters
     
@@ -27,13 +27,18 @@ class AliGenITSULib :public AliGenLib {
 
  private:
 
-  static Int_t IpLcPlus(TRandom * /*ran*/)  {return   (int)kLc;}
-  static Int_t IpLcMinus(TRandom * /*ran*/) {return  -(int)kLc;}
-  static Int_t IpLb(TRandom * /*ran*/)      {return   (int)kLb;}
-  static Int_t IpLbBar(TRandom * /*ran*/)   {return  -(int)kLb;}
-  static Int_t IpXic(TRandom * /*ran*/)     {return (int)kXi_c;}
-  static Int_t IpBPlus(TRandom * /*ran*/)   {return    (int)kB;}
-  static Int_t IpBMinus(TRandom * /*ran*/)  {return   -(int)kB;}
+  static Int_t IpLcPlus(TRandom * /*ran*/)  {return     (int)kLc;}
+  static Int_t IpLcMinus(TRandom * /*ran*/) {return    -(int)kLc;}
+  static Int_t IpLb(TRandom * /*ran*/)      {return     (int)kLb;}
+  static Int_t IpLbBar(TRandom * /*ran*/)   {return    -(int)kLb;}
+  static Int_t IpXic(TRandom * /*ran*/)     {return   (int)kXi_c;}
+  static Int_t IpXicBar(TRandom * /*ran*/)  {return  -(int)kXi_c;}
+  static Int_t IpBPlus(TRandom * /*ran*/)   {return      (int)kB;}
+  static Int_t IpBMinus(TRandom * /*ran*/)  {return     -(int)kB;}
+  static Int_t IpDsPlus(TRandom * /*ran*/)  {return     (int)kDs;}
+  static Int_t IpDsMinus(TRandom * /*ran*/) {return    -(int)kDs;}
+  static Int_t IpDPlus(TRandom * /*ran*/)   {return  (int)kDplus;}
+  static Int_t IpDMinus(TRandom * /*ran*/)  {return -(int)kDplus;}
 
   static Double_t PtFlat(const Double_t * /*px*/, const Double_t * /*dummy*/) {return 1;}
   static Double_t YFlat (const Double_t * /*py*/, const Double_t * /*dummy*/) {return 1;}