#include "AliGenAfterBurnerFlow.h"
#include "AliGenCocktailAfterBurner.h"
#include "AliMC.h"
+#include "AliRun.h"
#include "AliCollisionGeometry.h"
#include "AliGenCocktailEntry.h"
// private function
//
- if(TMath::Abs(pdg)>fgkPDG){
+ if(TMath::Abs(pdg)>=fgkPDG){
Error("AliAfterBurnerFlow","Overflow");
return;
}
////////////////////////////////////////////////////////////////////////////////////////////////////
-Float_t AliGenAfterBurnerFlow::GetNpNorm(Int_t npart)
+Float_t AliGenAfterBurnerFlow::GetNpNorm(Int_t npart) const
{
//
// Calculate npart norm.
////////////////////////////////////////////////////////////////////////////////////////////////////
-Bool_t AliGenAfterBurnerFlow::IsPrimary(Int_t pdg)
+Bool_t AliGenAfterBurnerFlow::IsPrimary(Int_t pdg) const
{
- if(pdg>fgkPDG) return kFALSE;
+ if(pdg>=fgkPDG) return kFALSE;
return fIsPrim[pdg];
}
Double_t CalcAngle(Double_t phi, Double_t phi0, Double_t phiRP, Double_t v2, Double_t v1=0.)
{
+ // Calculate relative angle
Double_t phi1 = phi-(phi+2*v1*TMath::Sin(phi-phiRP)+v2*TMath::Sin(2*(phi-phiRP))-phi0)/
(1.+2*v1*TMath::Cos(phi-phiRP)+ 2*v2*TMath::Cos(2*(phi-phiRP)));
if(TMath::Abs(phi/phi1-1.)<0.00001) return phi1;
void AliGenAfterBurnerFlow::InitPrimaries()
{
+ // Init the primary particle list
for(Int_t i=0; i<fgkPDG; i++) fIsPrim[i]=kFALSE;
//mesons
// M.G. Poghosyan
// PWG2 meeting on 06.05.2008 and 03.06.2008
- Info("Generate()","");
if (0)
for(Int_t ii=0; ii<fCounter;ii++)
//exclude incoming protons in PYTHIA
if(particle->GetPdgCode()==21) continue;
- if(TMath::Abs(pdg)>fgkPDG) continue;
+ if(TMath::Abs(pdg)>=fgkPDG) continue;
// is particle primary?
if(!fIsPrim[TMath::Abs(pdg)]) continue;
void AliGenAfterBurnerFlow::Rotate(Int_t i, Double_t phi, Bool_t IsPrim)
{
+ // Rotation
TParticle* particle = fStack->Particle(i);
TLorentzVector momentum;