EVGEN/AliGenLcLib.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16
  17 // Library class for particle pt and y distributions used for
  18 // LambdaC simulations.
  19 // To be used with AliGenParam.
  20 //
  21 // Author: Annalisa Mastroserio <Annalisa.Mastroserio@cern.ch>
  22 //
  23
  24 #include <TPDGCode.h>
  25 #include <TMath.h>
  26 #include <TRandom.h>
  27 #include <TString.h>
  28
  29 #include "AliGenLcLib.h"
  30 #include "AliLog.h"
  31
  32 ClassImp(AliGenLcLib)
  33
  34
  35 //---------------------------------------------
  36 //        LambdaC
  37 //---------------------------------------------
  38   Int_t AliGenLcLib::IpLcPlus(TRandom *)
  39 {
  40   //PDG code
  41   return 4122;
  42 }
  43
  44 Int_t AliGenLcLib::IpLcMinus(TRandom *)
  45 {
  46   //PDG code
  47   return -4122;
  48 }
  49
  50 Double_t AliGenLcLib::PtLcFlat( const Double_t *, const Double_t *)
  51 {
  52   // FLAT pt-distribution
  53   return 1;
  54 }
  55
  56 Double_t AliGenLcLib::PtLcExp( const Double_t *x, const Double_t *)
  57 {
  58   //     EXP  pt-distribution
  59   return x[0]*TMath::Exp(-x[0]/0.17);
  60 }
  61
  62 Double_t AliGenLcLib::YLcFlat(const Double_t */*x*/,const Double_t *)
  63 {
  64   //LambdaC             y-distribution
  65   return 5;
  66 }
  67
  68
  69
  70 typedef Double_t (*GenFunc)   (const Double_t*,  const Double_t*);
  71 typedef Int_t    (*GenFuncIp) (TRandom *);
  72
  73 GenFunc AliGenLcLib::GetPt(Int_t iPID, const char * sForm) const
  74 {
  75   // Return pointer to Pt parameterisation
  76   printf("PID: %i, form: %s \n",iPID,sForm);
  77   TString type(sForm);
  78   GenFunc func;
  79
  80   switch(iPID) {
  81
  82   case kLcPlus:
  83     if     (type=="FLAT")                                         {func=PtLcFlat; break;}
  84     else if(type=="EXP")                                          {func=PtLcExp; break;}
  85     else {
  86       AliFatal(Form("Unknown Pt distribution form: %s",sForm));   func=0;
  87     }
  88
  89   case kLcMinus:
  90     if     (type=="FLAT")                                         {func=PtLcFlat; break;}
  91     else if(type=="EXP")                                          {func=PtLcExp; break;}
  92     else {
  93       AliFatal(Form("Unknown Pt distribution form: %s",sForm));   func=0;
  94     }
  95
  96   default : AliFatal(Form("Unknown particle type: %i",iPID));      func=0;
  97   }//switch
  98
  99   return func;
 100 }
 101
 102 GenFunc AliGenLcLib::GetY(Int_t iPID, const char *sForm) const
 103 {
 104   AliDebug(1,Form("PID: %i, form: %s",iPID,sForm));
 105   GenFunc func;
 106   switch (iPID) {
 107
 108   case kLcPlus:                                                func=YLcFlat; break;
 109   case kLcMinus:                                               func=YLcFlat; break;
 110
 111   default  : AliFatal(Form("Unknown particle type: %i",iPID));     func=0; break;
 112
 113   }//switch
 114   return func;
 115 }
 116
 117 GenFuncIp AliGenLcLib::GetIp(Int_t iPID, const char *sForm) const
 118 {
 119   // Return pointer to particle type parameterisation
 120   AliDebug(1,Form("PID: %i, form: %s",iPID,sForm));   //////////
 121
 122   switch (iPID){
 123
 124   case kLcPlus:                                                  return IpLcPlus;
 125   case kLcMinus:                                                 return IpLcMinus;
 126
 127   default  : AliFatal(Form("Unknown particle type: %i",iPID));  return 0;
 128   }
 129 }