+// Generator of charm
+
+ AliGenCorrHF *gencharm = new AliGenCorrHF(1, 4);
+ gencharm->SetMomentumRange(0,9999);
+ gencharm->SetForceDecay(kAll);
+ ratioccbar = sigmaccbar/sigmaReaction;
+ if (!gMC) gencharm->SetDecayer(fDecayer);
+ gencharm->Init();
+ AddGenerator(gencharm,"CorrHFCharm",ratioccbar);
+
+//------------------------------------------------------------------
+// Generator of beauty
+
+ AliGenCorrHF *genbeauty = new AliGenCorrHF(1, 5);
+ genbeauty->SetMomentumRange(0,9999);
+ genbeauty->SetForceDecay(kAll);
+ ratiobbbar = sigmabbbar/sigmaReaction;
+ if (!gMC) genbeauty->SetDecayer(fDecayer);
+ genbeauty->Init();
+ AddGenerator(genbeauty,"CorrHFBeauty",ratiobbbar);
+
+//-------------------------------------------------------------------