This should work for pp.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 1 Jun 2004 09:19:22 +0000 (09:19 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 1 Jun 2004 09:19:22 +0000 (09:19 +0000)
DPMJET/ConfigDPM.C

index 71acf3e..93e7c8a 100644 (file)
@@ -17,11 +17,29 @@ void Config()
     gRandom->SetSeed(123456);
 
 
+
     // libraries required by geant321
     gSystem->Load("libgeant321");
+    gSystem->Load("libdpmjet");
+    gSystem->Load("libTDPMjet");
 
     new TGeant3("C++ Interface to Geant3");
 
+    AliRunLoader* rl=0x0;
+
+    cout<<"Config.C: Creating Run Loader ..."<<endl;
+    rl = AliRunLoader::Open("galice.root",
+                           AliConfig::GetDefaultEventFolderName(),
+                           "recreate");
+    if (rl == 0x0)
+    {
+       gAlice->Fatal("Config.C","Can not instatiate the Run Loader");
+       return;
+    }
+    rl->SetCompressionLevel(2);
+    rl->SetNumberOfEventsPerFile(100);
+    gAlice->SetRunLoader(rl);
+
     if (!gSystem->Getenv("CONFIG_FILE"))
     {
         TFile  *rootfile = new TFile("galice.root", "recreate");
@@ -31,8 +49,9 @@ void Config()
 
     TGeant3 *geant3 = (TGeant3 *) gMC;
     AliIonPDGCodes *PDGcodes = new AliIonPDGCodes();
+
     PDGcodes->AddParticlesToPdgDataBase();
-    geant3->DefineIons();
+//    geant3->DefineIons();
 
     //
     // Set External decayer
@@ -82,11 +101,11 @@ void Config()
     // ####  AliGenDPMjet generation   ######################################
     AliGenDPMjet *gener = new AliGenDPMjet(1);
     gener->SetBeamEnergy(2700.);
-    gener->SetEnergyCMS(5400.);
-    gener->SetProjectile(208,82);
-    gener->SetTarget(208,82);
+    gener->SetEnergyCMS(5500.);
+    gener->SetProjectile(1, 1);
+    gener->SetTarget(1, 1);
     gener->SetImpactParameterRange(12.,16.); 
-    gener->SetTrackingFlag(1);
+    gener->SetTrackingFlag(0);
     gener->Init();
     // 
     // Activate this line if you want the vertex smearing to happen
@@ -95,12 +114,11 @@ void Config()
     //gener->SetVertexSmear(perTrack); 
     // Field (L3 0.4 T)
     AliMagFMaps* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., 1);
-    rootfile->cd();
     gAlice->SetField(field);    
+    rl->CdGAFile();   
 
-
-    Int_t   iABSO  =  1;
-    Int_t   iDIPO  =  1;
+    Int_t   iABSO  =  0;
+    Int_t   iDIPO  =  0;
     Int_t   iFMD   =  0;
     Int_t   iFRAME =  0;
     Int_t   iHALL  =  0;
@@ -108,15 +126,15 @@ void Config()
     Int_t   iMAG   =  0;
     Int_t   iMUON  =  0;
     Int_t   iPHOS  =  0;
-    Int_t   iPIPE  =  1;
+    Int_t   iPIPE  =  0;
     Int_t   iPMD   =  0;
     Int_t   iRICH  =  0;
-    Int_t   iSHIL  =  1;
+    Int_t   iSHIL  =  0;
     Int_t   iSTART =  0;
     Int_t   iTOF   =  0;
     Int_t   iTPC   =  0;
     Int_t   iTRD   =  0;
-    Int_t   iZDC   =  1;
+    Int_t   iZDC   =  0;
     Int_t   iEMCAL =  0;
     Int_t   iCRT   =  0;
     Int_t   iVZERO =  0;