X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=EVGEN%2FAliGenPHOSlib.cxx;h=1585d3151bf2b0c7307293ed6c98c17f9b4e5a4e;hb=7892d33ba50e3c3d7bf46f9f6473cdc091725143;hp=516bd7a0b33233009089f9f9a76884f84a1e5bf2;hpb=88cb7938ca21d4a80991d4e7aa564008c29340f7;p=u%2Fmrichter%2FAliRoot.git diff --git a/EVGEN/AliGenPHOSlib.cxx b/EVGEN/AliGenPHOSlib.cxx index 516bd7a0b33..1585d3151bf 100644 --- a/EVGEN/AliGenPHOSlib.cxx +++ b/EVGEN/AliGenPHOSlib.cxx @@ -1,4 +1,3 @@ - /************************************************************************** * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. * * * @@ -36,6 +35,10 @@ // GPS @ SUBATECH, Nantes , France (October 1999) // http://www-subatech.in2p3.fr/~photons/subatech // martinez@subatech.in2p3.fr +// Additional particle species simulation options has been added: +// Charged Pion, Charged Kaons, KLong Proton, Anti-Proton, Neutron, +// Anti-Neutron --> Changes made by Gustavo Conesa in November 2004 +// Add flat Omega(782) distribution in Nov. 2010 by Renzhuo WAN //====================================================================== #include "TMath.h" @@ -52,7 +55,7 @@ ClassImp(AliGenPHOSlib) // Rapidity distribution YPion // Particle distribution IdPion 111, 211 and -211 (pi0, pi+ and pi-) // - Double_t AliGenPHOSlib::PtPion(Double_t *px, Double_t *) + Double_t AliGenPHOSlib::PtPion(const Double_t *px, const Double_t *) { // Pion transverse momentum distribtuion taken // from AliGenMUONlib class, version 3.01 of aliroot @@ -79,7 +82,7 @@ ClassImp(AliGenPHOSlib) y=kb*TMath::Exp(-sqrt(x*x+kxmpi2)/kt); return y*x; } - Double_t AliGenPHOSlib::YPion( Double_t *py, Double_t *) + Double_t AliGenPHOSlib::YPion( const Double_t *py, const Double_t *) { // // pion y-distribution @@ -94,28 +97,28 @@ ClassImp(AliGenPHOSlib) return ka*TMath::Exp(-ex); } - Int_t AliGenPHOSlib::IpPion(TRandom *ran) +Int_t AliGenPHOSlib::IpPion(TRandom */*ran*/) +{ +// particle composition pi+, pi0, pi- +// + + return 111 ; +} + Int_t AliGenPHOSlib::IpChargedPion(TRandom *ran) { // particle composition pi+, pi0, pi- // Float_t random = ran->Rndm(); - if ( (3.*random) < 1. ) - { - return 211 ; - } - else - { - if ( (3.*random) >= 2.) + if ( (2.*random) < 1. ) { - return -211 ; + return 211 ; + } + else + { + return -211 ; } - else - { - return 111 ; - } - } } //End Pions @@ -126,7 +129,7 @@ ClassImp(AliGenPHOSlib) // Particle distribution IdPi0Flat 111 (pi0) // -Double_t AliGenPHOSlib::PtPi0Flat(Double_t *px, Double_t *) +Double_t AliGenPHOSlib::PtPi0Flat(const Double_t */*px*/, const Double_t *) { // Pion transverse momentum flat distribution @@ -134,7 +137,7 @@ return 1; } -Double_t AliGenPHOSlib::YPi0Flat( Double_t *py, Double_t *) +Double_t AliGenPHOSlib::YPi0Flat( const Double_t */*py*/, const Double_t *) { // pion y-distribution @@ -163,13 +166,12 @@ Double_t AliGenPHOSlib::YPi0Flat( Double_t *py, Double_t *) // 7=>BARYONS-BARYONBARS // SCALING EN MASSE PAR RAPPORT A PTPI - // MASS 1=>PI, 2=>K, 3=>ETA, 4=>OMEGA, 5=>ETA',6=>PHI + // MASS 0=>PI, 1=>K, 2=>ETA, 3=>OMEGA, 4=>ETA',5=>PHI const Double_t khm[10] = {0.1396, 0.494, 0.547, 0.782, 0.957, 1.02, - // MASS 7=>BARYON-BARYONBAR + // MASS 6=>BARYON-BARYONBAR 0.938, 0. , 0., 0.}; // VALUE MESON/PI AT 5 GEV const Double_t kfmax[10]={1., 1., 1., 1., 1., 1., 1., 1., 1., 1.}; - np--; Double_t f5=TMath::Power(((sqrt(100.018215)+2.)/(sqrt(100.+khm[np]*khm[np])+2.0)),12.3); Double_t kfmax2=f5/kfmax[np]; // PIONS @@ -182,16 +184,16 @@ Double_t AliGenPHOSlib::YPi0Flat( Double_t *py, Double_t *) // End Scaling //============================================================================ // K A O N S - Double_t AliGenPHOSlib::PtKaon( Double_t *px, Double_t *) + Double_t AliGenPHOSlib::PtKaon( const Double_t *px, const Double_t *) { // kaon // pt-distribution //____________________________________________________________ - return PtScal(*px,2); // 2==> Kaon in the PtScal function + return PtScal(*px,1); // 1==> Kaon in the PtScal function } - Double_t AliGenPHOSlib::YKaon( Double_t *py, Double_t *) + Double_t AliGenPHOSlib::YKaon( const Double_t *py, const Double_t *) { // y-distribution //____________________________________________________________ @@ -230,20 +232,43 @@ Double_t AliGenPHOSlib::YPi0Flat( Double_t *py, Double_t *) } } } + + Int_t AliGenPHOSlib::IpChargedKaon(TRandom *ran) +{ +// particle composition +// + + Float_t random = ran->Rndm(); + + if (random < 0.5) { + return 321; // K+ + } else { + return -321; // K- + } + + +} +Int_t AliGenPHOSlib::IpKaon0L(TRandom *) +{ + // particle composition + // + + return 130; // K^0 long +} // End Kaons //============================================================================ //============================================================================ // E T A S - Double_t AliGenPHOSlib::PtEta( Double_t *px, Double_t *) + Double_t AliGenPHOSlib::PtEta( const Double_t *px, const Double_t *) { // etas // pt-distribution //____________________________________________________________ - return PtScal(*px,3); // 3==> Eta in the PtScal function + return PtScal(*px,2); // 2==> Eta in the PtScal function } - Double_t AliGenPHOSlib::YEta( Double_t *py, Double_t *) + Double_t AliGenPHOSlib::YEta( const Double_t *py, const Double_t *) { // y-distribution //____________________________________________________________ @@ -274,7 +299,7 @@ Double_t AliGenPHOSlib::YPi0Flat( Double_t *py, Double_t *) // Particle distribution IdEtaFlat 111 (pi0) // -Double_t AliGenPHOSlib::PtEtaFlat(Double_t *px, Double_t *) +Double_t AliGenPHOSlib::PtEtaFlat(const Double_t */*px*/, const Double_t *) { // Eta transverse momentum flat distribution @@ -282,7 +307,7 @@ Double_t AliGenPHOSlib::PtEtaFlat(Double_t *px, Double_t *) } -Double_t AliGenPHOSlib::YEtaFlat( Double_t *py, Double_t *) +Double_t AliGenPHOSlib::YEtaFlat( const Double_t */*py*/, const Double_t *) { // // pion y-distribution @@ -297,20 +322,20 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t *py, Double_t *) // return 221 ; } -// End Pi0Flat +// End EtaFlat //============================================================================ //============================================================================ // O M E G A S - Double_t AliGenPHOSlib::PtOmega( Double_t *px, Double_t *) + Double_t AliGenPHOSlib::PtOmega( const Double_t *px, const Double_t *) { // omegas // pt-distribution //____________________________________________________________ - return PtScal(*px,4); // 4==> Omega in the PtScal function + return PtScal(*px,3); // 3==> Omega in the PtScal function } - Double_t AliGenPHOSlib::YOmega( Double_t *py, Double_t *) + Double_t AliGenPHOSlib::YOmega( const Double_t *py, const Double_t *) { // y-distribution //____________________________________________________________ @@ -334,18 +359,51 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t *py, Double_t *) } // End Omega //============================================================================ +//====================================================================== +// Omega(782) Flat Distribution +// Transverse momentum distribution PtOmegaFlat +// Rapidity distribution YOmegaFlat +// Particle distribution IdOmegaFlat 223(0mega) +// + +Double_t AliGenPHOSlib::PtOmegaFlat(const Double_t */*px*/, const Double_t *) +{ +// omega transverse momentum flat distribution + +return 1; + +} + +Double_t AliGenPHOSlib::YOmegaFlat( const Double_t */*py*/, const Double_t *) +{ + +// omega y-distribution +// + return 1.; +} + + Int_t AliGenPHOSlib::IpOmegaFlat(TRandom *) +{ + +// particle composition omega +// + return 223 ; +} +// End OmegaFlat + + //============================================================================ // E T A P R I M E - Double_t AliGenPHOSlib::PtEtaprime( Double_t *px, Double_t *) + Double_t AliGenPHOSlib::PtEtaprime( const Double_t *px, const Double_t *) { // etaprime // pt-distribution //____________________________________________________________ - return PtScal(*px,5); // 5==> Etaprime in the PtScal function + return PtScal(*px,4); // 4==> Etaprime in the PtScal function } - Double_t AliGenPHOSlib::YEtaprime( Double_t *py, Double_t *) + Double_t AliGenPHOSlib::YEtaprime( const Double_t *py, const Double_t *) { // y-distribution //____________________________________________________________ @@ -371,16 +429,16 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t *py, Double_t *) //=================================================================== //============================================================================ // P H I S - Double_t AliGenPHOSlib::PtPhi( Double_t *px, Double_t *) + Double_t AliGenPHOSlib::PtPhi( const Double_t *px, const Double_t *) { // phi // pt-distribution //____________________________________________________________ - return PtScal(*px,6); // 6==> Phi in the PtScal function + return PtScal(*px,5); // 5==> Phi in the PtScal function } - Double_t AliGenPHOSlib::YPhi( Double_t *py, Double_t *) + Double_t AliGenPHOSlib::YPhi( const Double_t *py, const Double_t *) { // y-distribution //____________________________________________________________ @@ -406,16 +464,16 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t *py, Double_t *) //=================================================================== //============================================================================ // B A R Y O N S == protons, protonsbar, neutrons, and neutronsbars - Double_t AliGenPHOSlib::PtBaryon( Double_t *px, Double_t *) + Double_t AliGenPHOSlib::PtBaryon( const Double_t *px, const Double_t *) { // baryons // pt-distribution //____________________________________________________________ - return PtScal(*px,7); // 7==> Baryon in the PtScal function + return PtScal(*px,6); // 6==> Baryon in the PtScal function } - Double_t AliGenPHOSlib::YBaryon( Double_t *py, Double_t *) + Double_t AliGenPHOSlib::YBaryon( const Double_t *py, const Double_t *) { // y-distribution //____________________________________________________________ @@ -453,128 +511,188 @@ Double_t AliGenPHOSlib::YEtaFlat( Double_t *py, Double_t *) return -2112; // n bar } } +} + + Int_t AliGenPHOSlib::IpProton(TRandom *) +{ +// particle composition +// + return 2212; // p + +} + Int_t AliGenPHOSlib::IpAProton(TRandom *) +{ +// particle composition +// + return -2212; // p bar + +} + + Int_t AliGenPHOSlib::IpNeutron(TRandom *) +{ +// particle composition +// + return 2112; // n + +} + Int_t AliGenPHOSlib::IpANeutron(TRandom *) +{ +// particle composition +// + return -2112; // n + } // End Baryons //=================================================================== - -typedef Double_t (*GenFunc) (Double_t*, Double_t*); - GenFunc AliGenPHOSlib::GetPt(Int_t param, const char* tname) const +typedef Double_t (*GenFunc) (const Double_t*, const Double_t*); +GenFunc AliGenPHOSlib::GetPt(Int_t param, const char* /*tname*/) const { // Return pinter to pT parameterisation GenFunc func; switch (param) - { - case kPion: + { + case kPion: func=PtPion; break; - case kPi0Flat: + case kPi0Flat: func=PtPi0Flat; break; - case kKaon: + case kKaon: func=PtKaon; break; - case kEta: + case kEta: func=PtEta; break; - case kEtaFlat: + case kEtaFlat: func=PtEtaFlat; break; - case kOmega: + case kOmega: func=PtOmega; break; - case kEtaPrime: + case kOmegaFlat: + func=PtOmegaFlat; + break; + case kEtaPrime: func=PtEtaprime; break; - case kBaryon: + case kBaryon: func=PtBaryon; break; - default: + default: func=0; printf(" unknown parametrisationn"); - } + } return func; } - GenFunc AliGenPHOSlib::GetY(Int_t param, const char* tname) const +GenFunc AliGenPHOSlib::GetY(Int_t param, const char* /*tname*/) const { -// Return pointer to Y parameterisation - GenFunc func; - switch (param) + // Return pointer to Y parameterisation + GenFunc func; + switch (param) { case kPion: - func=YPion; - break; + func=YPion; + break; case kPi0Flat: - func=YPi0Flat; - break; + func=YPi0Flat; + break; case kKaon: - func=YKaon; - break; + func=YKaon; + break; case kEta: - func=YEta; - break; - case kEtaFlat: - func=YEtaFlat; - break; + func=YEta; + break; + case kEtaFlat: + func=YEtaFlat; + break; case kOmega: - func=YOmega; - break; + func=YOmega; + break; + case kOmegaFlat: + func=YOmegaFlat; + break; case kEtaPrime: - func=YEtaprime; - break; + func=YEtaprime; + break; case kPhi: - func=YPhi; - break; + func=YPhi; + break; case kBaryon: - func=YBaryon; - break; + func=YBaryon; + break; default: - func=0; - printf(" unknown parametrisationn"); + func=0; + printf(" unknown parametrisationn"); } - return func; + return func; } typedef Int_t (*GenFuncIp) (TRandom *); - GenFuncIp AliGenPHOSlib::GetIp(Int_t param, const char* tname) const +GenFuncIp AliGenPHOSlib::GetIp(Int_t param, const char* /*tname*/) const { -// Return pointer to particle composition - GenFuncIp func; - switch (param) + // Return pointer to particle composition + GenFuncIp func; + switch (param) { case kPion: - func=IpPion; - break; + func=IpPion; + break; + case kChargedPion: + func=IpChargedPion; + break; case kPi0Flat: - func=IpPi0Flat; - break; + func=IpPi0Flat; + break; case kKaon: - func=IpKaon; - break; + func=IpKaon; + break; + case kChargedKaon: + func=IpChargedKaon; + break; + case kKaon0L: + func=IpKaon0L; + break; case kEta: - func=IpEta; - break; + func=IpEta; + break; case kEtaFlat: - func=IpEtaFlat; - break; - + func=IpEtaFlat; + break; + case kOmega: - func=IpOmega; - break; + func=IpOmega; + break; + case kOmegaFlat: + func=IpOmegaFlat; + break; case kEtaPrime: - func=IpEtaprime; - break; + func=IpEtaprime; + break; case kPhi: - func=IpPhi; - break; + func=IpPhi; + break; case kBaryon: - func=IpBaryon; - break; + func=IpBaryon; + break; + case kProton: + func=IpProton; + break; + case kAProton: + func=IpAProton; + break; + case kNeutron: + func=IpNeutron; + break; + case kANeutron: + func=IpANeutron; + break; + default: - func=0; - printf(" unknown parametrisationn"); + func=0; + printf(" unknown parametrisationn"); } - return func; + return func; }