updated vertex selection

[u/mrichter/AliRoot.git] / MUON / Config.C

diff --git a/MUON/Config.C b/MUON/Config.C
index 435bfd4f00332b834a94ff6038adb4a26b8cba10..0b085c981f291bcc16604e79fc6ed2cdb301c072 100644 (file)
--- a/MUON/Config.C
+++ b/MUON/Config.C
@@ -29,8 +29,15 @@ void Config(char directory[100]="", char option[6]="param", const char* digitsto
   // Config file for MUON test
   //=====================================================================
   //  Libraries required by geant321
+  gSystem->Load("liblhapdf.so");      // Parton density functions
+  gSystem->Load("libpythia6.so");     // Pythia
   gSystem->Load("libgeant321.so");
+  gSystem->Load("libEG");
+  gSystem->Load("libEGPythia6");
+  gSystem->Load("libAliPythia6.so");  // ALICE specific implementations
+
   new TGeant3TGeo("C++ Interface to Geant3");
+
   //=======================================================================
   //  Create the output file    
   Text_t filename[100];
@@ -126,6 +133,19 @@ void Config(char directory[100]="", char option[6]="param", const char* digitsto
     gener->SetTrackingFlag(1);
     gener->Init();
   }
+  if (!strcmp(option,"paramJpsi")) {
+    AliGenParam *gener = new AliGenParam(1, AliGenMUONlib::kJpsi);
+    gener->SetMomentumRange(0,999);
+    gener->SetPtRange(0,100.);
+    gener->SetPhiRange(0., 360.);
+    gener->SetCutOnChild(1);
+    gener->SetChildPhiRange(0.,360.);
+    gener->SetChildThetaRange(171.0,178.0);
+    gener->SetOrigin(0,0,0);
+    gener->SetForceDecay(kDiMuon);
+    gener->SetTrackingFlag(1);
+    gener->Init();
+  }
   if (!strcmp(option,"hijing")) { //Hijing generator from ConfigPPR in macros
     AliGenHijing *gener = new AliGenHijing(-1);
     // centre of mass energy 
@@ -167,8 +187,7 @@ void Config(char directory[100]="", char option[6]="param", const char* digitsto
   }  
   //============================================================= 
   // Field (L3 0.5 T)
-  AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 1, 1., 10., AliMagFMaps::k5kG);
-  gAlice->SetField(field);
+  TGeoGlobalMagField::Instance()->SetField(new AliMagF("Maps","Maps", -1., -1, AliMagF::k5kG));
   //============================================================= 
   //=================== Alice BODY parameters =============================
   AliBODY *BODY = new AliBODY("BODY","Alice envelop");
@@ -203,8 +222,9 @@ void Config(char directory[100]="", char option[6]="param", const char* digitsto
   // Noise-only digits in tracker/trigger (0=no noise, 1=default (noise in tracker), 2=noise in tracker and trigger):
   //MUON->SetDigitizerWithNoise(kFALSE);
 
-  // Use fast raw data decoder
-  // MUON->SetFastDecoder(kTRUE);  
+  // Use non-high performance raw data decoder 
+  //MUON->SetFastTrackerDecoder(kFALSE);  
+  //MUON->SetFastTriggerDecoder(kFALSE);  
   
   //
   // If SetAlign, the detection elements transformations
@@ -214,8 +234,8 @@ void Config(char directory[100]="", char option[6]="param", const char* digitsto
   // To generate and read scaler trigger events in rawdata
   // MUON->SetTriggerScalerEvent();
   
-  // To switch on the tail effect
-  // MUON->SetTailEffect(kTRUE);
+  // To switch off the tail effect
+  // MUON->SetTailEffect(kFALSE);
 
   // If you want to play with builders, first reset the geometry builder,
   // and then add yours.