X-Git-Url: http://git.uio.no/git/?a=blobdiff_plain;f=PYTHIA6%2FAliDecayerPythia.cxx;h=06ff91ce7a15db6cb496b4ed8737233ff587fe2b;hb=cbd58b9ca6f605fd8df95cbeb1726a7e90dd23ae;hp=d4f3432f8483a92d40facd1d289461287b397b9b;hpb=513879494112344e9d1ac377f58c69f30cf006e4;p=u%2Fmrichter%2FAliRoot.git diff --git a/PYTHIA6/AliDecayerPythia.cxx b/PYTHIA6/AliDecayerPythia.cxx index d4f3432f848..06ff91ce7a1 100644 --- a/PYTHIA6/AliDecayerPythia.cxx +++ b/PYTHIA6/AliDecayerPythia.cxx @@ -24,9 +24,11 @@ #include "AliDecayerPythia.h" #include "AliPythia.h" +#include "AliLog.h" #include #include #include +#include ClassImp(AliDecayerPythia) @@ -57,7 +59,11 @@ Bool_t AliDecayerPythia::fgInit = kFALSE; AliDecayerPythia::AliDecayerPythia(): fPythia(AliPythia::Instance()), fDecay(kAll), - fHeavyFlavour(kTRUE) + fHeavyFlavour(kTRUE), + fLongLived(kFALSE), + fPatchOmegaDalitz(0), + fDecayerExodus(0), + fPi0(1) { // Default Constructor for (Int_t i=0; i< 501; i++) fBraPart[i]= 1.; @@ -68,14 +74,20 @@ AliDecayerPythia::AliDecayerPythia(const AliDecayerPythia &decayer): AliDecayer(decayer), fPythia(0), fDecay(kAll), - fHeavyFlavour(kTRUE) + fHeavyFlavour(kTRUE), + fLongLived(kFALSE), + fPatchOmegaDalitz(0), + fDecayerExodus(0), + fPi0(1) { // Copy Constructor decayer.Copy(*this); + for (Int_t i = 0; i < 501; i++) fBraPart[i] = 0.; } void AliDecayerPythia::Init() { + // Initialisation // if (!fgInit) { @@ -102,44 +114,110 @@ void AliDecayerPythia::Init() } } //...Switch off decay of pi0, K0S, Lambda, Sigma+-, Xi0-, Omega-. - + +/* if (fDecay != kNeutralPion) { fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 0); } else { fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1); } +*/ + if (fPi0) fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1); - - fPythia->SetMDCY(fPythia->Pycomp(310) ,1,0); - fPythia->SetMDCY(fPythia->Pycomp(3122),1,0); - fPythia->SetMDCY(fPythia->Pycomp(3112),1,0); - fPythia->SetMDCY(fPythia->Pycomp(3212),1,0); - fPythia->SetMDCY(fPythia->Pycomp(3222),1,0); - fPythia->SetMDCY(fPythia->Pycomp(3312),1,0); - fPythia->SetMDCY(fPythia->Pycomp(3322),1,0); - fPythia->SetMDCY(fPythia->Pycomp(3334),1,0); - + Int_t isw = 0; + if (fLongLived) isw = 1; + + fPythia->SetMDCY(fPythia->Pycomp(310) ,1, isw); + fPythia->SetMDCY(fPythia->Pycomp(3122),1, isw); + fPythia->SetMDCY(fPythia->Pycomp(3112),1, isw); +// fPythia->SetMDCY(fPythia->Pycomp(3212),1, isw); // Sigma0 decays elem. + fPythia->SetMDCY(fPythia->Pycomp(3222),1, isw); + fPythia->SetMDCY(fPythia->Pycomp(3312),1, isw); + fPythia->SetMDCY(fPythia->Pycomp(3322),1, isw); + fPythia->SetMDCY(fPythia->Pycomp(3334),1, isw); + // .. Force decay channels ForceDecay(); } +void AliDecayerPythia::SwitchOffParticle(Int_t kf) +{ +//switch off decay for particle "kf" +fPythia->SetMDCY(fPythia->Pycomp(kf),1,0); +} + void AliDecayerPythia::Decay(Int_t idpart, TLorentzVector* p) { // Decay a particle // - Float_t energy = p->Energy(); - Float_t theta = p->Theta(); - Float_t phi = p->Phi(); - + Float_t energy = p->Energy(); + Float_t theta = p->Theta(); + Float_t phi = p->Phi(); + + if(!fDecayerExodus) { Lu1Ent(0, idpart, energy, theta, phi); - fPythia->GetPrimaries(); + } else { + + // EXODUS decayer + if(idpart == 111){ + fPythia->SetMDCY(fPythia->Pycomp(22) ,1, 0); + Lu1Ent(0, idpart, energy, theta, phi); + fPythia->PizeroDalitz(); + fPythia->SetMDCY(fPythia->Pycomp(22) ,1, 1); + } + else if(idpart == 221){ + fPythia->SetMDCY(fPythia->Pycomp(22) ,1, 0); + Lu1Ent(0, idpart, energy, theta, phi); + fPythia->EtaDalitz(); + fPythia->SetMDCY(fPythia->Pycomp(22) ,1, 1); + } + else if(idpart == 113){ + Lu1Ent(0, idpart, energy, theta, phi); + fPythia->RhoDirect(); + } + else if(idpart == 223){ + fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 0); + Lu1Ent(0, idpart, energy, theta, phi); + fPythia->OmegaDirect(); + fPythia->OmegaDalitz(); + fPythia->SetMDCY(fPythia->Pycomp(111) ,1, 1); + } + else if(idpart == 331){ + fPythia->SetMDCY(fPythia->Pycomp(22) ,1, 0); + Lu1Ent(0, idpart, energy, theta, phi); + fPythia->EtaprimeDalitz(); + fPythia->SetMDCY(fPythia->Pycomp(22) ,1, 1); + } + else if(idpart == 333){ + fPythia->SetMDCY(fPythia->Pycomp(221) ,1, 0); + Lu1Ent(0, idpart, energy, theta, phi); + fPythia->PhiDirect(); + fPythia->PhiDalitz(); + fPythia->SetMDCY(fPythia->Pycomp(221) ,1, 1); + } + else if(idpart == 443){ + Lu1Ent(0, idpart, energy, theta, phi); + fPythia->JPsiDirect(); + } + } + + fPythia->GetPrimaries(); } Int_t AliDecayerPythia::ImportParticles(TClonesArray *particles) { // Import the decay products // - return fPythia->ImportParticles(particles, "All"); + const Float_t kconvT=0.001/2.999792458e8; // mm/c to seconds conversion + const Float_t kconvL=1./10; // mm to cm conversion + int np = fPythia->ImportParticles(particles, "All"); + // pythia assigns decay time in mm/c, convert to seconds + for (int ip=1;ipAt(ip); + if (!prod) continue; + prod->SetProductionVertex(prod->Vx()*kconvL,prod->Vy()*kconvL,prod->Vz()*kconvL,kconvT*prod->T()); + } + return np; } @@ -151,7 +229,6 @@ void AliDecayerPythia::ForceDecay() // Decay_t decay=fDecay; fPythia->SetMSTJ(21,2); - if (decay == kNoDecayHeavy) return; // // select mode @@ -177,8 +254,7 @@ void AliDecayerPythia::ForceDecay() ForceParticleDecay( 5232, products1, mult1, 3); ForceParticleDecay( 5332, products1, mult1, 3); ForceParticleDecay( 100443, 443, 1); // Psi' -> J/Psi X - ForceParticleDecay( 443, 13, 2); // J/Psi -> mu+ mu- - + ForceParticleDecay( 443, 13, 2); // J/Psi -> mu+ mu- ForceParticleDecay( 411,13,1); // D+/- ForceParticleDecay( 421,13,1); // D0 ForceParticleDecay( 431,13,1); // D_s @@ -187,6 +263,25 @@ void AliDecayerPythia::ForceDecay() ForceParticleDecay( 4232,13,1); // Sigma_c ForceParticleDecay( 4332,13,1); // Omega_c break; + case kChiToJpsiGammaToMuonMuon: + products[0] = 443; + products[1] = 22; + mult[0] = 1; + mult[1] = 1; + ForceParticleDecay( 20443, products, mult, 2); // Chi_1c -> J/Psi Gamma + ForceParticleDecay( 445, products, mult, 2); // Chi_2c -> J/Psi Gamma + ForceParticleDecay( 443, 13, 2); // J/Psi -> mu+ mu- + break; + case kChiToJpsiGammaToElectronElectron: + products[0] = 443; + products[1] = 22; + mult[0] = 1; + mult[1] = 1; + ForceParticleDecay( 20443, products, mult, 2); // Chi_1c -> J/Psi Gamma + ForceParticleDecay( 445, products, mult, 2); // Chi_2c -> J/Psi Gamma + ForceParticleDecay( 443, 11, 2); // J/Psi -> e+ e- + break; + case kBSemiMuonic: ForceParticleDecay( 511,13,1); // B0 ForceParticleDecay( 521,13,1); // B+/- @@ -196,6 +291,15 @@ void AliDecayerPythia::ForceDecay() ForceParticleDecay( 5232,13,1); // Sigma_b ForceParticleDecay( 5332,13,1); // Omega_b break; + case kDSemiMuonic: + ForceParticleDecay( 411,13,1); // D0 + ForceParticleDecay( 421,13,1); // D+/- + ForceParticleDecay( 431,13,1); // D_s + ForceParticleDecay( 4122,13,1); // Lambda_c + ForceParticleDecay( 4132,13,1); // Xsi_c + ForceParticleDecay( 4232,13,1); // Sigma_c + ForceParticleDecay( 4332,13,1); // Omega_c + break; case kSemiMuonic: ForceParticleDecay( 411,13,1); // D+/- ForceParticleDecay( 421,13,1); // D0 @@ -223,6 +327,9 @@ void AliDecayerPythia::ForceDecay() ForceParticleDecay(100553,13,2);// Upsilon' ForceParticleDecay(200553,13,2);// Upsilon'' break; + case kJpsiDiMuon: + ForceParticleDecay( 443,13,2); // J/Psi + break; case kBSemiElectronic: ForceParticleDecay( 511,11,1); // B0 ForceParticleDecay( 521,11,1); // B+/- @@ -259,6 +366,14 @@ void AliDecayerPythia::ForceDecay() ForceParticleDecay(100553,11,2);// Upsilon' ForceParticleDecay(200553,11,2);// Upsilon'' break; + case kPsiPrimeJpsiDiElectron: + products[0] = 443; + products[1] = 211; + mult[0] = 1; + mult[1] = 2; + ForceParticleDecay( 100443, products, mult, 2, 1); + ForceParticleDecay( 443,11,2); + break; case kBJpsiDiMuon: products[0] = 443; @@ -293,6 +408,13 @@ void AliDecayerPythia::ForceDecay() ForceParticleDecay( 531,443,1); // B_s ForceParticleDecay( 5122,443,1); // Lambda_b break; + case kBJpsiUndecayed: + ForceParticleDecay( 511,443,1); // B0 + ForceParticleDecay( 521,443,1); // B+/- + ForceParticleDecay( 531,443,1); // B_s + ForceParticleDecay( 5122,443,1); // Lambda_b + fPythia->SetMDCY(fPythia->Pycomp(443),1,0); // switch-off J/psi + break; case kBPsiPrimeDiElectron: ForceParticleDecay( 511,100443,1); // B0 ForceParticleDecay( 521,100443,1); // B+/- @@ -343,6 +465,10 @@ void AliDecayerPythia::ForceDecay() break; case kOmega: ForceOmega(); + break; + case kLambda: + ForceLambda(); + break; case kAll: break; case kNoDecay: @@ -351,6 +477,49 @@ void AliDecayerPythia::ForceDecay() case kNoDecayHeavy: case kNeutralPion: break; + case kNoDecayBeauty: + SwitchOffBDecay(); + break; + case kElectronEM: + ForceParticleDecay( 111,11,1); // pi^0 + ForceParticleDecay( 221,11,1); // eta + ForceParticleDecay( 113,11,1); // rho + ForceParticleDecay( 223,11,1); // omega + ForceParticleDecay( 331,11,1); // etaprime + ForceParticleDecay( 333,11,1); // phi + ForceParticleDecay( 443,11,1); // jpsi + break; + case kDiElectronEM: + ForceParticleDecay( 111,11,2); // pi^0 + ForceParticleDecay( 221,11,2); // eta + ForceParticleDecay( 113,11,2); // rho + ForceParticleDecay( 223,11,2); // omega + ForceParticleDecay( 331,11,2); // etaprime + ForceParticleDecay( 333,11,2); // phi + ForceParticleDecay( 443,11,2); // jpsi + break; + case kGammaEM: + ForceParticleDecay( 111, 22, 1); // pi^0 + ForceParticleDecay( 221, 22, 1); // eta + ForceParticleDecay( 113, 22, 1); // rho + ForceParticleDecay( 223, 22, 1); // omega + ForceParticleDecay( 331, 22, 1); // etaprime + ForceParticleDecay( 333, 22, 1); // phi + ForceParticleDecay( 443, 22, 1); // jpsi + ForceParticleDecay( 3212, 22, 1); // Sigma 0 + ForceParticleDecay( 310, 22, 1); // K0s + ForceParticleDecay( 2224, 22, 1); // Delta++ + ForceParticleDecay( 2214, 22, 1); // Delta+ + ForceParticleDecay( 1114, 22, 1); // Delta- + ForceParticleDecay( 2114, 22, 1); // Delta0 + ForceParticleDecay( 213, 22, 1); // rho+ + ForceParticleDecay( -213, 22, 1); // rho- + ForceParticleDecay( 313, 22, 1); // K0star + + break; + case kBeautyUpgrade: + ForceBeautyUpgrade(); + break; } } @@ -366,6 +535,34 @@ void AliDecayerPythia::SwitchOffHeavyFlavour() for (Int_t i = 156; i <= 160; i++) fPythia->SetMDME(i, 1, 0); } +void AliDecayerPythia::ForceBeautyUpgrade() +{ + // + // Force dedicated decay channels of signals ineresting + // for the ITS upgrade (Lb, Lc, Xi_c, B) + // + + // Lb: 50% of them in Lc pi+ and the rest with a Lc in final state + if(gRandom->Rndm()<0.50) { + const Int_t prod3[3]={4122,211,0}; + Int_t mult3[3]={1,1,1}; + ForceParticleDecay(5122,prod3,mult3,3,1); + } else { + ForceParticleDecay( 5122, 4122, 1); + } + // Xi_c + ForceParticleDecay( 4232, 3312, 1); + // B+ -> D0pi+ + const Int_t prod[2]={421,211}; + Int_t mult[2]={1,1}; + ForceParticleDecay(521,prod,mult,2,1); + // B0 -> D*-pi+ + const Int_t prod2[2]={413,211}; + ForceParticleDecay(511,prod2,mult,2,1); + // force charm hadronic decays (D mesons and Lc) + ForceHadronicD(0); +} + void AliDecayerPythia::Lu1Ent(Int_t flag, Int_t idpart, Double_t mom, Double_t theta, Double_t phi) { @@ -412,6 +609,14 @@ void AliDecayerPythia::ForceHadronicD(Int_t optUse4Bodies) Int_t iDeltaPP = 2224; Int_t productsD[2] = {kProton, kPiPlus}, multD[2] = {1, 1}; ForceParticleDecay(iDeltaPP, productsD, multD, 2); + // for Lambda_c -> Lambda(1520) pi+ -> p K- pi+ + Int_t iLambda1520 = 3124; + Int_t productsL[2] = {kProton, kKMinus}, multL[2] = {1, 1}; + ForceParticleDecay(iLambda1520, productsL, multL, 2); + // for Lambda_c -> Lambda pi+ + Int_t iLambda=3122; + //for Lambda_c -> antiK0 p + Int_t iK0bar=-311; Int_t decayP1[kNHadrons][4] = @@ -432,25 +637,63 @@ void AliDecayerPythia::ForceHadronicD(Int_t optUse4Bodies) }; Int_t decayP3[kNHadrons][4] = { - {-1 , -1 , -1 , -1}, + {kPiPlus , iPhi , 0 , 0}, {kKMinus , kPiPlus, iRho0 , 0 }, {-1 , -1 , -1 , -1}, {-1 , -1 , -1 , -1}, {kProton , kKMinus, kPiPlus , 0} }; + // for Lambda_c -> Lambda_1520 pi+ -> p K- pi+, D0-> K*0 pi+ pi- -> K3pi + Int_t decayP4[kNHadrons][4] = + { + {iKstarbar0 , kKPlus , 0 , 0}, + {iKstarbar0 , kPiPlus , kPiMinus, 0}, + {-1 , -1 , -1 , -1}, + {-1 , -1 , -1 , -1}, + {iLambda1520 , kPiPlus , 0 , 0} + }; + // for Lambda_c -> Lambda pi+ + Int_t decayP5[kNHadrons][4] = + { + {-1 , -1 , -1 , -1}, + {-1 , -1 , -1 , -1}, + {-1 , -1 , -1 , -1}, + {-1 , -1 , -1 , -1}, + {iLambda , kPiPlus, 0 , 0} + }; + + // for Lambda_c -> K0bar p + Int_t decayP6[kNHadrons][4] = + { + {-1 , -1 , -1 , -1}, + {-1 , -1 , -1 , -1}, + {-1 , -1 , -1 , -1}, + {-1 , -1 , -1 , -1}, + {kProton , iK0bar, 0 , 0} + }; + if(optUse4Bodies==0){ for(Int_t iDau=0;iDau<4;iDau++){ decayP2[1][iDau]=-1; decayP3[1][iDau]=-1; + decayP4[1][iDau]=-1; } } for (Int_t ihadron = 0; ihadron < kNHadrons; ihadron++) { Int_t kc = fPythia->Pycomp(hadron[ihadron]); - fPythia->SetMDCY(kc,1,1); Int_t ifirst = fPythia->GetMDCY(kc,2); Int_t ilast = ifirst + fPythia->GetMDCY(kc,3)-1; + Double_t norm = 0.; + for (channel=ifirst; channel<=ilast;channel++) norm+=fPythia->GetBRAT(channel); + if (norm < 1.-1.e-12 || norm > 1.+1.e-12) { + char pName[16]; + fPythia->Pyname(hadron[ihadron],pName); + AliWarning(Form("Total branching ratio of %s (PDG code = %d) not equal to 1 (= %f)",pName,hadron[ihadron],norm)); + } + fBraPart[kc] = norm; + fPythia->SetMDCY(kc,1,1); for (channel = ifirst; channel <= ilast; channel++) { if (( @@ -471,6 +714,25 @@ void AliDecayerPythia::ForceHadronicD(Int_t optUse4Bodies) fPythia->GetKFDP(channel,3) == decayP3[ihadron][2] && fPythia->GetKFDP(channel,4) == decayP3[ihadron][3] && fPythia->GetKFDP(channel,5) == 0 + ) || ( + fPythia->GetKFDP(channel,1) == decayP4[ihadron][0] && + fPythia->GetKFDP(channel,2) == decayP4[ihadron][1] && + fPythia->GetKFDP(channel,3) == decayP4[ihadron][2] && + fPythia->GetKFDP(channel,4) == decayP4[ihadron][3] && + fPythia->GetKFDP(channel,5) == 0 + ) || ( + fPythia->GetKFDP(channel,1) == decayP5[ihadron][0] && + fPythia->GetKFDP(channel,2) == decayP5[ihadron][1] && + fPythia->GetKFDP(channel,3) == decayP5[ihadron][2] && + fPythia->GetKFDP(channel,4) == decayP5[ihadron][3] && + fPythia->GetKFDP(channel,5) == 0 + ) || ( + fPythia->GetKFDP(channel,1) == decayP6[ihadron][0] && + fPythia->GetKFDP(channel,2) == decayP6[ihadron][1] && + fPythia->GetKFDP(channel,3) == decayP6[ihadron][2] && + fPythia->GetKFDP(channel,4) == decayP6[ihadron][3] && + fPythia->GetKFDP(channel,5) == 0 + )) { @@ -480,6 +742,7 @@ void AliDecayerPythia::ForceHadronicD(Int_t optUse4Bodies) fBraPart[kc] -= fPythia->GetBRAT(channel); } // selected channel ? } // decay channels + if (norm > 0) fBraPart[kc] /= norm; } // hadrons } @@ -494,7 +757,16 @@ void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t product, Int_t m fPythia->SetMDCY(kc,1,1); Int_t ifirst=fPythia->GetMDCY(kc,2); Int_t ilast=ifirst+fPythia->GetMDCY(kc,3)-1; - fBraPart[kc] = 1; + Double_t norm = 0.; + for (Int_t channel=ifirst; channel<=ilast;channel++){ + norm+=fPythia->GetBRAT(channel); + } + if (norm < 1.-1.e-12 || norm > 1.+1.e-12) { + char pName[16]; + fPythia->Pyname(particle,pName); + AliWarning(Form("Total branching ratio of %s (PDG code = %d) not equal to 1 (= %f)",pName,particle,norm)); + } + fBraPart[kc] = norm; // // Loop over decay channels for (Int_t channel=ifirst; channel<=ilast;channel++) { @@ -502,12 +774,13 @@ void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t product, Int_t m fPythia->SetMDME(channel,1,1); } else { fPythia->SetMDME(channel,1,0); - fBraPart[kc]-=fPythia->GetBRAT(channel); + fBraPart[kc] -= fPythia->GetBRAT(channel); } } + if (norm > 0.) fBraPart[kc] /= norm; } -void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t* products, Int_t* mult, Int_t npart) +void AliDecayerPythia::ForceParticleDecay(Int_t particle, const Int_t* products, Int_t* mult, Int_t npart, Bool_t flag) { // // Force decay of particle into products with multiplicity mult @@ -516,7 +789,14 @@ void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t* products, Int_t fPythia->SetMDCY(kc,1,1); Int_t ifirst=fPythia->GetMDCY(kc,2); Int_t ilast=ifirst+fPythia->GetMDCY(kc,3)-1; - fBraPart[kc] = 1; + Double_t norm = 0.; + for (Int_t channel=ifirst; channel<=ilast;channel++) norm+=fPythia->GetBRAT(channel); + if (norm < 1.-1.e-12 || norm > 1.+1.e-12) { + char pName[16]; + fPythia->Pyname(particle,pName); + AliWarning(Form("Total branching ratio of %s (PDG code = %d) not equal to 1 (= %f)",pName,particle,norm)); + } + fBraPart[kc] = norm; // // Loop over decay channels for (Int_t channel = ifirst; channel <= ilast; channel++) { @@ -524,13 +804,17 @@ void AliDecayerPythia::ForceParticleDecay(Int_t particle, Int_t* products, Int_t for (Int_t i = 0; i < npart; i++) { nprod += (CountProducts(channel, products[i]) >= mult[i]); } - if (nprod) { + if ((nprod && !flag) || ((nprod == npart) && flag)) { fPythia->SetMDME(channel,1,1); - } else { + } else { // fPythia->SetMDME(channel,1,0); - fBraPart[kc]-=fPythia->GetBRAT(channel); + fBraPart[kc] -= fPythia->GetBRAT(channel); } } + if (norm > 0.) fBraPart[kc] /= norm; + + + } void AliDecayerPythia::DefineParticles() @@ -614,7 +898,37 @@ void AliDecayerPythia::ForceOmega() } // decay channels } +void AliDecayerPythia::ForceLambda() +{ + // Force Lambda -> p pi- + Int_t kc=fPythia->Pycomp(3122); + fPythia->SetMDCY(kc,1,1); + Int_t ifirst = fPythia->GetMDCY(kc,2); + Int_t ilast = ifirst + fPythia->GetMDCY(kc,3)-1; + for (Int_t channel = ifirst; channel <= ilast; channel++) { + if ( + fPythia->GetKFDP(channel,1) == kProton && + fPythia->GetKFDP(channel,2) == kPiMinus && + fPythia->GetKFDP(channel,3) == 0 + ) + { + fPythia->SetMDME(channel,1,1); + } else { + fPythia->SetMDME(channel,1,0); + } // selected channel ? + } // decay channels +} + +void AliDecayerPythia::SwitchOffBDecay() +{ +// Switch off B-decays + Int_t heavyB[]={511,521,531,5122,5132,5232,5332}; + for(int i=0;i<4;i++) + { + fPythia->SetMDCY(fPythia->Pycomp(heavyB[i]),1,0); + } +} Float_t AliDecayerPythia::GetPartialBranchingRatio(Int_t kf) { @@ -648,27 +962,6 @@ void AliDecayerPythia::ReadDecayTable() } -#ifdef never -void AliDecayerPythia::Streamer(TBuffer &R__b) -{ - // Stream an object of class AliDecayerPythia. - - if (R__b.IsReading()) { - Version_t R__v = R__b.ReadVersion(); if (R__v) { } - AliDecayer::Streamer(R__b); - (AliPythia::Instance())->Streamer(R__b); - R__b >> (Int_t&)fDecay; - R__b.ReadStaticArray(fBraPart); - } else { - R__b.WriteVersion(AliDecayerPythia::IsA()); - AliDecayer::Streamer(R__b); - R__b << fPythia; - R__b << (Int_t)fDecay; - R__b.WriteArray(fBraPart, 501); - } -} -#endif - void AliDecayerPythia::Copy(TObject &) const { //