+ case herwig:
+//********************************************
+// Example for Charm Production with Pythia *
+//********************************************
+ {
+ AliGenHerwig *gener = new AliGenHerwig(-1);
+// final state kinematic cuts
+ gener->SetMomentumRange(0,7000);
+ gener->SetPhiRange(0. ,360.);
+ gener->SetThetaRange(0., 180.);
+ gener->SetYRange(-10,10);
+ gener->SetPtRange(0,7000);
+// vertex position and smearing
+ gener->SetOrigin(0,0,0); // vertex position
+ gener->SetVertexSmear(kPerEvent);
+ gener->SetSigma(0,0,5.6); // Sigma in (X,Y,Z) (cm) on IP position
+// Beam momenta
+ gener->SetBeamMomenta(7000,7000);
+// Beams
+ gener->SetProjectile("P");
+ gener->SetTarget("P");
+// Structure function
+ gener->SetStrucFunc(kGRVHO);
+// Hard scatering
+ gener->SetPtHardMin(200);
+ gener->SetPtRMS(20);
+// Min bias
+ gener->SetProcess(8000);
+// No Tracking
+ gener->SetTrackingFlag(0);
+ gGener = gener;
+ }
+ break;
+
+ case mcatnlo:
+ {
+ AliGenHerwig *gener = new AliGenHerwig(-1);
+ gener->SetMomentumRange(0,7000);
+ gener->SetPhiRange(0. ,360.);
+ gener->SetThetaRange(0., 180.);
+ gener->SetYRange(-10,10);
+ gener->SetPtRange(0,7000);
+ gener->SetOrigin(0,0,0); // vertex position
+ gener->SetVertexSmear(kPerEvent);
+ gener->SetSigma(0,0,5.6); // Sigma in (X,Y,Z) (cm) on IP position
+ gener->SetBeamMomenta(7000,7000);
+ gener->SetProjectile("P");
+ gener->SetTarget("P");
+ gener->SetStrucFunc(kCTEQ5M);
+ gener->SetProcess(-1705);
+ gener->SetHardProcessFile("sb.events");
+ gener->SetEventListRange(0,1);
+ gener->SetTrackingFlag(0);
+ gGener = gener;
+ }
+ break;
+