Line data    Source code

       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           6 : ClassImp(AliGenLcLib)
      33             : 
      34             : 
      35             : //---------------------------------------------
      36             : //        LambdaC 
      37             : //---------------------------------------------
      38             :   Int_t AliGenLcLib::IpLcPlus(TRandom *)
      39             : {
      40             :   //PDG code
      41           0 :   return 4122;
      42             : }
      43             : 
      44             : Int_t AliGenLcLib::IpLcMinus(TRandom *)
      45             : {
      46             :   //PDG code
      47           0 :   return -4122;
      48             : }
      49             : 
      50             : Double_t AliGenLcLib::PtLcFlat( const Double_t *, const Double_t *)
      51             : {
      52             :   // FLAT pt-distribution
      53           0 :   return 1; 
      54             : }
      55             : 
      56             : Double_t AliGenLcLib::PtLcExp( const Double_t *x, const Double_t *)
      57             : {
      58             :   //     pt-distribution
      59             :   //return x[0]*TMath::Exp(-x[0]/0.16); // distribution used in LHC11f1 for the anchor runs : 139441, 139510, 139511, 130513, 130514, 130517. 
      60           0 :   return TMath::GammaDist(x[0],2,0,1.7); //distribution as in LHC11a10a of the prompt Lc whose daughters are in |eta|<0.9. Used for Lb as well.
      61             : }
      62             : 
      63             : Double_t AliGenLcLib::YLcFlat(const Double_t *,const Double_t *)
      64             : {
      65             :   //LambdaC             y-distribution 
      66           0 :   return 1;
      67             : }
      68             : 
      69             : 
      70             : 
      71             : typedef Double_t (*GenFunc)   (const Double_t*,  const Double_t*);
      72             : typedef Int_t    (*GenFuncIp) (TRandom *);
      73             : 
      74             : GenFunc AliGenLcLib::GetPt(Int_t iPID, const char * sForm) const
      75             : {
      76             :   // Return pointer to Pt parameterisation
      77           0 :   printf("PID: %i, form: %s \n",iPID,sForm);  
      78           0 :   TString type(sForm);
      79             :   GenFunc func;
      80             : 
      81           0 :   switch(iPID) {
      82             : 
      83             :   case kLcPlus:  
      84           0 :     if     (type=="FLAT")                                         {func=PtLcFlat; break;}
      85           0 :     else if(type=="EXP")                                          {func=PtLcExp; break;}
      86             :     else {
      87           0 :       AliFatal(Form("Unknown Pt distribution form: %s",sForm));   func=0;
      88             :     }
      89             : 
      90             :   case kLcMinus:  
      91           0 :     if     (type=="FLAT")                                         {func=PtLcFlat; break;}
      92           0 :     else if(type=="EXP")                                          {func=PtLcExp; break;}
      93             :     else {
      94           0 :       AliFatal(Form("Unknown Pt distribution form: %s",sForm));   func=0;
      95             :     }
      96             : 
      97           0 :   default : AliFatal(Form("Unknown particle type: %i",iPID));      func=0;
      98           0 :   }//switch
      99             : 
     100             :   return func;
     101           0 : }
     102             : 
     103             : GenFunc AliGenLcLib::GetY(Int_t iPID, const char *sForm) const
     104             : {
     105           0 :   AliDebug(1,Form("PID: %i, form: %s",iPID,sForm));   
     106             :   GenFunc func;
     107           0 :   switch (iPID) {
     108             : 
     109           0 :   case kLcPlus:                                                func=YLcFlat; break;
     110           0 :   case kLcMinus:                                               func=YLcFlat; break;
     111             : 
     112           0 :   default  : AliFatal(Form("Unknown particle type: %i",iPID));     func=0; break;
     113             : 
     114             :   }//switch
     115           0 :   return func;
     116             : }
     117             : 
     118             : GenFuncIp AliGenLcLib::GetIp(Int_t iPID, const char *sForm) const
     119             : {
     120             :   // Return pointer to particle type parameterisation
     121           0 :   AliDebug(1,Form("PID: %i, form: %s",iPID,sForm));   //////////      
     122             : 
     123           0 :   switch (iPID){
     124             : 
     125           0 :   case kLcPlus:                                                  return IpLcPlus;
     126           0 :   case kLcMinus:                                                 return IpLcMinus;
     127             : 
     128           0 :   default  : AliFatal(Form("Unknown particle type: %i",iPID));  return 0;
     129             :   }
     130           0 : }