[u/mrichter/AliRoot.git] / EVGEN / AliGenLcLib.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/


// Library class for particle pt and y distributions used for 
// LambdaC simulations.
// To be used with AliGenParam.
//
// Author: Annalisa Mastroserio <Annalisa.Mastroserio@cern.ch>
//

#include <TPDGCode.h>
#include <TMath.h>
#include <TRandom.h>
#include <TString.h>

#include "AliGenLcLib.h"
#include "AliLog.h"

ClassImp(AliGenLcLib)


//---------------------------------------------
//        LambdaC 
//---------------------------------------------
  Int_t AliGenLcLib::IpLcPlus(TRandom *)
{
  //PDG code
  return 4122;
}

Int_t AliGenLcLib::IpLcMinus(TRandom *)
{
  //PDG code
  return -4122;
}

Double_t AliGenLcLib::PtLcFlat( const Double_t *, const Double_t *)
{
  // FLAT pt-distribution
  return 1; 
}

Double_t AliGenLcLib::PtLcExp( const Double_t *x, const Double_t *)
{
  //     pt-distribution
  //return x[0]*TMath::Exp(-x[0]/0.16); // distribution used in LHC11f1 for the anchor runs : 139441, 139510, 139511, 130513, 130514, 130517. 
  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.
}

Double_t AliGenLcLib::YLcFlat(const Double_t *,const Double_t *)
{
  //LambdaC             y-distribution 
  return 1;
}


typedef Double_t (*GenFunc)   (const Double_t*,  const Double_t*);
typedef Int_t    (*GenFuncIp) (TRandom *);

GenFunc AliGenLcLib::GetPt(Int_t iPID, const char * sForm) const
{
  // Return pointer to Pt parameterisation
  printf("PID: %i, form: %s \n",iPID,sForm);	
  TString type(sForm);
  GenFunc func;

  switch(iPID) {

  case kLcPlus:  
    if     (type=="FLAT")                                         {func=PtLcFlat; break;}
    else if(type=="EXP")                                          {func=PtLcExp; break;}
    else {
      AliFatal(Form("Unknown Pt distribution form: %s",sForm));   func=0;
    }

  case kLcMinus:  
    if     (type=="FLAT")                                         {func=PtLcFlat; break;}
    else if(type=="EXP")                                          {func=PtLcExp; break;}
    else {
      AliFatal(Form("Unknown Pt distribution form: %s",sForm));   func=0;
    }

  default : AliFatal(Form("Unknown particle type: %i",iPID));      func=0;
  }//switch

  return func;
}

GenFunc AliGenLcLib::GetY(Int_t iPID, const char *sForm) const
{
  AliDebug(1,Form("PID: %i, form: %s",iPID,sForm));	
  GenFunc func;
  switch (iPID) {

  case kLcPlus:                                                func=YLcFlat; break;
  case kLcMinus:                                               func=YLcFlat; break;

  default  : AliFatal(Form("Unknown particle type: %i",iPID));     func=0; break;

  }//switch
  return func;
}

GenFuncIp AliGenLcLib::GetIp(Int_t iPID, const char *sForm) const
{
  // Return pointer to particle type parameterisation
  AliDebug(1,Form("PID: %i, form: %s",iPID,sForm));   //////////	

  switch (iPID){

  case kLcPlus:                                                  return IpLcPlus;
  case kLcMinus:                                                 return IpLcMinus;

  default  : AliFatal(Form("Unknown particle type: %i",iPID));  return 0;
  }
}
Commit	Line	Data
feb24caa	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	{
0f8eca74	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	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.
feb24caa	61	}
feb24caa	62
0f8eca74	63	Double_t AliGenLcLib::YLcFlat(const Double_t ,const Double_t )
feb24caa	64	{
feb24caa	65	//LambdaC y-distribution
0f8eca74	66	return 1;
feb24caa	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	printf("PID: %i, form: %s \n",iPID,sForm);
	78	TString type(sForm);
	79	GenFunc func;
	80
	81	switch(iPID) {
	82
	83	case kLcPlus:
	84	if (type=="FLAT") {func=PtLcFlat; break;}
	85	else if(type=="EXP") {func=PtLcExp; break;}
	86	else {
	87	AliFatal(Form("Unknown Pt distribution form: %s",sForm)); func=0;
	88	}
	89
	90	case kLcMinus:
	91	if (type=="FLAT") {func=PtLcFlat; break;}
	92	else if(type=="EXP") {func=PtLcExp; break;}
	93	else {
	94	AliFatal(Form("Unknown Pt distribution form: %s",sForm)); func=0;
	95	}
	96
	97	default : AliFatal(Form("Unknown particle type: %i",iPID)); func=0;
	98	}//switch
	99
	100	return func;
	101	}
	102
	103	GenFunc AliGenLcLib::GetY(Int_t iPID, const char *sForm) const
	104	{
	105	AliDebug(1,Form("PID: %i, form: %s",iPID,sForm));
	106	GenFunc func;
	107	switch (iPID) {
	108
	109	case kLcPlus: func=YLcFlat; break;
	110	case kLcMinus: func=YLcFlat; break;
	111
	112	default : AliFatal(Form("Unknown particle type: %i",iPID)); func=0; break;
	113
	114	}//switch
	115	return func;
	116	}
	117
	118	GenFuncIp AliGenLcLib::GetIp(Int_t iPID, const char *sForm) const
	119	{
	120	// Return pointer to particle type parameterisation
	121	AliDebug(1,Form("PID: %i, form: %s",iPID,sForm)); //////////
	122
	123	switch (iPID){
	124
	125	case kLcPlus: return IpLcPlus;
	126	case kLcMinus: return IpLcMinus;
	127
	128	default : AliFatal(Form("Unknown particle type: %i",iPID)); return 0;
	129	}
	130	}