//------------------------------------------------------------------------------// EvtGen is a particle decay simulator specifically designed for the needs of B-physics studies. Complete manual of EvtGen is available at address: http://robbep.home.cern.ch/robbep/EvtGen/GuideEvtGen.pdf //------------------------------------------------------------------------------// [14/12/2009] 1)This module contains: - EvtGen original code organized in three directories: EvtGenBase-EvtGenModels-EvtGen. EvtGenModels directory contains all models avilable in EvtGen to decay particles. - The interface classes AliGenEvtGen and AliDecayerEvtGen for AliRoot users: the methods to decay particles are in AliDecayerEvtGen which represents the implementation of AliDecayer using EvtGen package 2)EvtGen required libphotos.so to generate Final State Radiation. Photos code is located in PHOTOS directory. 3)The implemented cases to force beauty hadrons interest those channels: B->J/psi+X - B->J/psi+X,J/psi->e-e+ - B->J/psi+X,J/psi->mu-mu+ - B->e+X. Those are decayed by exclusive models of EvtGen, and all parameters and models used to decay them are not modified with respect to the official EvtGen release. Some of particles produced in the decay chain (i.e. Lamda_c from Lambda_b, Xi_c0 from Xi_b-, etc...) are decayed by Pythia with the configuration setting in the AliDecayerPythia. 4)Polarization is still to be implemented in the interface classes. It would be done soon. 5)To use EvtGen to decay beauty particles adjust those lines in the Config.C: //----- load libraries (after loading libpythia6.so) gSystem->Load("libphotos.so"); gSystem->Load("libEvtGenBase"); gSystem->Load("libEvtGenModels"); gSystem->Load("libEvtGen"); gSystem->Load("libTEvtGen"); //----- declare an AliGenCocktail AliGenCocktail *generCock=new AliGenCocktail(); generCock->UsePerEventRates(); //----- declare Pythia configuration: switch-off beauty decays in Pythia AliGenPythia *pythia= ..... .... pythia->SetForceDecay(kNoDecayBeauty); //----- declare EvtGen configuration and put the two generators in the cocktail AliGenEvtGen *gene = new AliGenEvtGen(); gene->SetForceDecay(kBJpsiDiElectron); gene->SetParticleSwitchedOff(AliGenEvtGen::kBeautyPart); generCock->AddGenerator(pythia, "Pythia", 1.); generCock->AddGenerator(gene,"gene",1.);