Initialize fSepctrum = 0 in ctors

[u/mrichter/AliRoot.git] / macros / Config_gener.C

diff --git a/macros/Config_gener.C b/macros/Config_gener.C
index 31117bd65d4f8499160b9aa7347092a9c138e279..22ed1e253f7ca1555b6ae799f321e93c1f7fa0bb 100644 (file)
--- a/macros/Config_gener.C
+++ b/macros/Config_gener.C
@@ -1,10 +1,9 @@
-enum gentype_t {hijing, hijing_p, gun, box, pythia, 
+enum gentype_t {hijing, hijingParam, gun, box, pythia, 
                param1, param2, param3, 
                cocktail, fluka, halo, ntuple, scan, doublescan};
 
 gentype_t gentype=param1;
 
-
 ntracks=1;
 
 void Config()
@@ -101,7 +100,6 @@ geant3->SetCUTS(cut,cut, cut, cut, cut, cut,  cut,  cut, cut,  cut, tofmax);
      gener->SetSigma(0,0,0);         //Sigma in (X,Y,Z) (cm) on IP position
      break;
  case hijing:
-      case hijing:
      AliGenHijing *gener = new AliGenHijing(-1);
 // centre of mass energy 
      gener->SetEnergyCMS(5500);
@@ -140,25 +138,49 @@ geant3->SetCUTS(cut,cut, cut, cut, cut, cut,  cut,  cut, cut,  cut, tofmax);
      gener->SetVertexSmear(kPerEvent); 
 // no tracking
      gener->SetTrackingFlag(0);
+     break;
+     
  case pythia:
 //********************************************
 // Example for Charm  Production with Pythia *
 //********************************************
-
-     AliGenPythia *gener = new AliGenPythia(ntracks);
+     AliGenPythia *gener = new AliGenPythia(-1);
+//   final state kinematic cuts
      gener->SetMomentumRange(0,999);
-     gener->SetPhiRange(0,360);
+     gener->SetPhiRange(-180,180);
      gener->SetThetaRange(0., 180.);
      gener->SetYRange(-10,10);
      gener->SetPtRange(0,100);
-     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->SetStrucFunc(DO_Set_1);
-     gener->SetProcess(mb); 
+//   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
+//   Structure function
+// DO_Set_1
+// GRV_LO
+// GRV_HO
+// MRS_D_minus
+// MRS_D0
+// MRS_G
+// CTEQ_2pM
+// CTEQ_4M
+     gener->SetStrucFunc(GRV_HO);
+// Select corection for nuclear structure functions
+//     gener->SetNuclei(208,208);
+//
+//   Process type
+//   charm, beauty, charm_unforced, beauty_unforced, jpsi, jpsi_chi, mb
+     gener->SetProcess(beauty);
+//   
+//   Pt transfer of the hard scattering
+     gener->SetPtHard(0.,5.);
+//   Decay type (semielectronic, semimuonic, nodecay)
+     gener->SetForceDecay(semielectronic);
+//   Centre of mass energy 
      gener->SetEnergyCMS(5500.);
-     break;
-     
+//   No Tracking 
+     gener->SetTrackingFlag(0);
+     break;              
 
  case param1:
 //*******************************************************