message: Allow to have radiative decays for Upsilon as well
diff: @@ -19,34 +19,32 @@ AliGenerator* GenParamCustomRadiativeDecay()
  AliGenCocktail *gener = new AliGenCocktail();
  gener->UsePerEventRates();
  - AliGenParam *jpsiGener = new AliGenParam(1,-1,PtCustom,YCustom,V2Custom,IpCustom);
  - jpsiGener->SetMomentumRange(0,1e6);
  - jpsiGener->SetPtRange(0,999.);
  - jpsiGener->SetYRange(-4.2, -2.3);
  - jpsiGener->SetPhiRange(0., 360.);
  + AliGenParam *quarkoniaGener = new AliGenParam(1,-1,PtCustom,YCustom,V2Custom,IpCustom);
  + quarkoniaGener->SetMomentumRange(0,1e6);
  + quarkoniaGener->SetPtRange(0,999.);
  + quarkoniaGener->SetYRange(-4.2, -2.3);
  + quarkoniaGener->SetPhiRange(0., 360.);
  - jpsiGener->SetTrackingFlag(1);
  - jpsiGener->SetForceDecay(kNoDecay); // Jpsi particles decay are switched-off for Pythia
  + quarkoniaGener->SetTrackingFlag(1);
  + quarkoniaGener->SetForceDecay(kNoDecay); // quarkonia particles decay are switched-off for Pythia
  // evtGen (for radiative decays)
  AliGenEvtGen *evtGener = new AliGenEvtGen();
  evtGener->SetForceDecay(kDiMuon);
  - evtGener->SetParticleSwitchedOff(AliGenEvtGen::kCharmPart);
  + evtGener->SetParticleSwitchedOff(AliGenEvtGen::kAllPart); // will be used by both J/psi and Upsilon
  - // add Jpsi generator to cocktail
  - gener->AddGenerator(jpsiGener,"Jpsi",1);
  + // add quarkonia generator to cocktail
  + gener->AddGenerator(quarkoniaGener,"quarkonia",1);
  // add EvtGen generator to cocktail
  gener->AddGenerator(evtGener,"EvtGen",1);
  return gener;
  -
  - return gener;
  }
  //-------------------------------------------------------------------------
  Int_t IpCustom( TRandom *)
  {
  - // particle to simulate (e.g. 443 for J/psi)
  + // particle to simulate (e.g. 443 for J/psi, 553 for Upsilon)
  return VAR_GENPARAMCUSTOM_PDGPARTICLECODE;
  }
  @@ -60,7 +58,8 @@ Double_t PtCustom( const Double_t *px, const Double_t */*dummy*/ )
  p1 = VAR_GENPARAMCUSTOM_PT_P1; // 18.05;
  p2 = VAR_GENPARAMCUSTOM_PT_P2; // 2.05;
  p3 = VAR_GENPARAMCUSTOM_PT_P3; // 3.34;
  - return p0 *x / TMath::Power( p1 + TMath::Power(x,p2), p3 );
  +// return p0 *x / TMath::Power( p1 + TMath::Power(x,p2), p3 );
  + return p0 *x / TMath::Power( 1. + TMath::Power(x/p1,p2), p3 );
  }
  //-------------------------------------------------------------------------
  @@ -71,7 +70,8 @@ Double_t YCustom( const Double_t *py, const Double_t */*dummy*/ )
  Float_t p0,p1;
  p0 = VAR_GENPARAMCUSTOM_Y_P0; // 4.08e5;
  p1 = VAR_GENPARAMCUSTOM_Y_P1; // 7.1e4;
  - return p0 + p1*y;
  +// return p0 + p1*y;
  + return p0*TMath::Exp(-(1./2.)*TMath::Power(((y)/p1),2));
  }
  //-------------------------------------------------------------------------
  