In vmctest:
authorivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Mar 2011 17:42:20 +0000 (17:42 +0000)
committerivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Wed, 9 Mar 2011 17:42:20 +0000 (17:42 +0000)
- Updating macros according to AliRoot/test
- Now geometry is built from the code also in Geant4 simulation
  and it is not read from a file

test/vmctest/gun/commonConfig.C
test/vmctest/gun/genGunConfig.C
test/vmctest/gun/rec.C
test/vmctest/gun/runtest.sh
test/vmctest/gun/sim.C
test/vmctest/ppbench/commonConfig.C
test/vmctest/ppbench/g3Config.C
test/vmctest/ppbench/g4Config.C
test/vmctest/ppbench/runtest.sh
test/vmctest/ppbench/sim.C

index 0050bbd..bde0a39 100644 (file)
@@ -18,7 +18,7 @@ const char * pprTrigConfName[] = {
 
 enum ConfigVersion_t {
     kConfigV0,  // default configuration  
-    kConfigV1   // configuration for LHC production
+    kConfigV1   // configuration for LHC 2010 production
 };
 
 // Options 
@@ -29,8 +29,7 @@ static TString comment;
 // Functions
 void  LoadPythia();
 
-void commonConfig(Bool_t setRootGeometry = kFALSE,
-                  ConfigVersion_t configVersion = kConfigV1)
+void commonConfig(ConfigVersion_t configVersion = kConfigV0)
 {
   cout << "Running commonConfig.C ... " << endl;
 
@@ -49,11 +48,6 @@ void commonConfig(Bool_t setRootGeometry = kFALSE,
   // ALICE steering object (AliRunLoader)
   //=======================================================================
 
-  // Set Root geometry file
-  if ( setRootGeometry ) {
-    AliSimulation::Instance()->SetGeometryFile("geometry.root");
-  }
-
   AliRunLoader* rl 
     = AliRunLoader::Open("galice.root",
                          AliConfig::GetDefaultEventFolderName(),
@@ -108,31 +102,9 @@ void commonConfig(Bool_t setRootGeometry = kFALSE,
   Int_t   iTRD   =  1;
   Int_t   iZDC   =  1;
   Int_t   iEMCAL =  1;
-  Int_t   iACORDE = 0;
+  Int_t   iACORDE = 1;
   Int_t   iVZERO =  1;
-/*
-  Int_t   iABSO  =  0;
-  Int_t   iDIPO  =  0;
-  Int_t   iFMD   =  0;
-  Int_t   iFRAME =  0;
-  Int_t   iHALL  =  0;
-  Int_t   iITS   =  0;
-  Int_t   iMAG   =  0;
-  Int_t   iMUON  =  0;
-  Int_t   iPHOS  =  0;
-  Int_t   iPIPE  =  0;
-  Int_t   iPMD   =  0;
-  Int_t   iHMPID =  0;
-  Int_t   iSHIL  =  0;
-  Int_t   iT0    =  0;
-  Int_t   iTOF   =  0;
-  Int_t   iTPC   =  1;
-  Int_t   iTRD   =  0;
-  Int_t   iZDC   =  0;
-  Int_t   iEMCAL =  0;
-  Int_t   iACORDE = 0;
-  Int_t   iVZERO =  0;
-*/
+
   rl->CdGAFile();
   //=================== Alice BODY parameters =============================
   AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
@@ -287,7 +259,7 @@ void commonConfig(Bool_t setRootGeometry = kFALSE,
   {
       //=================== EMCAL parameters ============================
     if ( configVersion == kConfigV0 ) 
-      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
+      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETEV1");
     else if ( configVersion == kConfigV1 )  
       AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_FIRSTYEAR");
   }
index 1e18e32..aef47b5 100644 (file)
@@ -15,7 +15,7 @@ AliGenerator* genGunConfig()
   // Event generator
   //=======================================================================
 
-  // The cocktail iitself
+  // The cocktail itself
 
   AliGenCocktail *gener = new AliGenCocktail();
   gener->SetPhiRange(0, 360);
@@ -175,6 +175,77 @@ AliGenerator* genGunConfig()
   gfmd3i->SetPhi(50);
   gener->AddGenerator(gfmd3i,"Pi- for FMD3i",1);
   
+  //VZERO C
+  AliGenFixed *gv0c=new AliGenFixed(1);
+  gv0c->SetPart(kPiPlus);
+  gv0c->SetMomentum(1.5);
+  gv0c->SetTheta(170);
+  gv0c->SetPhi(50);
+  gener->AddGenerator(gv0c,"Pi+ for V0C",1);
+  
+  //VZERO A
+  AliGenFixed *gv0a=new AliGenFixed(1);
+  gv0a->SetPart(kPiMinus);
+  gv0a->SetMomentum(1.5);
+  gv0a->SetTheta(1.5);
+  gv0a->SetPhi(70);
+  gener->AddGenerator(gv0a,"Pi- for V0A",1);
+
+
+  //PMD
+  AliGenFixed *gpmd=new AliGenFixed(1);
+  gpmd->SetPart(kGamma);
+  gpmd->SetMomentum(2);
+  gpmd->SetTheta(12.6);
+  gpmd->SetPhi(60);
+  gener->AddGenerator(gpmd,"Gamma for PMD",1);
+
+  //ZDC
+  AliGenFixed *gzdc1=new AliGenFixed(1);
+  gzdc1->SetPart(kProton);
+  gzdc1->SetMomentum(700);
+  gzdc1->SetTheta(0.6);
+  gzdc1->SetPhi(60);
+  gener->AddGenerator(gzdc1,"Proton for ZDC",1);
+
+  AliGenFixed *gzdc2=new AliGenFixed(1);
+  gzdc2->SetPart(kNeutron);
+  gzdc2->SetMomentum(500);
+  gzdc2->SetTheta(0.6);
+  gzdc2->SetPhi(60);
+  gener->AddGenerator(gzdc2,"Neutron for ZDC",1);
+
+  //T0
+  AliGenFixed *gt0=new AliGenFixed(1);
+  gt0->SetPart(kPiPlus);
+  gt0->SetMomentum(2);
+  gt0->SetTheta(5.1);
+  gt0->SetPhi(60);
+  gener->AddGenerator(gt0,"Pi+ for T0",1);
+
+  AliGenFixed *gt01=new AliGenFixed(1);
+  gt01->SetPart(kPiMinus);
+  gt01->SetMomentum(2);
+  gt01->SetTheta(5.1);
+  gt01->SetPhi(60);
+  gener->AddGenerator(gt01,"Pi- for T0",1);
+
+
+  //ACORDE
+  AliGenFixed *gacorde=new AliGenFixed(1);
+  gacorde->SetPart(kMuonPlus);
+  gacorde->SetMomentum(20);
+  gacorde->SetTheta(90.);
+  gacorde->SetPhi(90);
+  gener->AddGenerator(gacorde,"Muon+ for ACORDE",1);
+
+  AliGenFixed *gacorde1=new AliGenFixed(1);
+  gacorde1->SetPart(kMuonMinus);
+  gacorde1->SetMomentum(20);
+  gacorde1->SetTheta(90.);
+  gacorde1->SetPhi(90);
+  gener->AddGenerator(gacorde1,"Muon- for ACORDE",1);
+
   gener->Init();
   
   return gener;
index 918fef3..d0882ab 100644 (file)
@@ -12,6 +12,9 @@ void rec() {
   reco.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
   reco.SetSpecificStorage("GRP/GRP/Data",
                          Form("local://%s",gSystem->pwd()));
+  reco.SetRunPlaneEff(kTRUE);
+
+  reco.SetFractionFriends(1.);
 
   TStopwatch timer;
   timer.Start();
index f0a9c33..eb263cc 100755 (executable)
@@ -24,7 +24,6 @@ fi
 
 if [ "$RUNG4" = "1" ]; then 
   rm -rf *.root *.dat *.log fort* hlt hough raw* recraw/*.root recraw/*.log
-  cp g3/geometry.root .
   aliroot -b -q  sim.C\($NEVENTS,\""$G4CONFIG"\"\)  2>&1 | tee sim.log
   aliroot -b -q rec.C      2>&1 | tee rec.log
   rm -fr $G4OUTDIR
index 5adf59c..9cf3c90 100644 (file)
@@ -3,12 +3,12 @@
 // Macro for running simulation in test/vmctest/gun.
 // From test/gun. 
 
-void sim(Int_t nev=4, const TString& config) {
+void sim(Int_t nev=1, const TString& config) {
 
   AliSimulation simulator(config);
-  simulator.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");
+  simulator.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO ACORDE");
   simulator.SetMakeDigitsFromHits("ITS TPC");
-  simulator.SetWriteRawData("","raw.root",kTRUE);
+  simulator.SetWriteRawData("ALL","raw.root",kTRUE);
 
   simulator.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
   simulator.SetSpecificStorage("GRP/GRP/Data",
index 0050bbd..bde0a39 100644 (file)
@@ -18,7 +18,7 @@ const char * pprTrigConfName[] = {
 
 enum ConfigVersion_t {
     kConfigV0,  // default configuration  
-    kConfigV1   // configuration for LHC production
+    kConfigV1   // configuration for LHC 2010 production
 };
 
 // Options 
@@ -29,8 +29,7 @@ static TString comment;
 // Functions
 void  LoadPythia();
 
-void commonConfig(Bool_t setRootGeometry = kFALSE,
-                  ConfigVersion_t configVersion = kConfigV1)
+void commonConfig(ConfigVersion_t configVersion = kConfigV0)
 {
   cout << "Running commonConfig.C ... " << endl;
 
@@ -49,11 +48,6 @@ void commonConfig(Bool_t setRootGeometry = kFALSE,
   // ALICE steering object (AliRunLoader)
   //=======================================================================
 
-  // Set Root geometry file
-  if ( setRootGeometry ) {
-    AliSimulation::Instance()->SetGeometryFile("geometry.root");
-  }
-
   AliRunLoader* rl 
     = AliRunLoader::Open("galice.root",
                          AliConfig::GetDefaultEventFolderName(),
@@ -108,31 +102,9 @@ void commonConfig(Bool_t setRootGeometry = kFALSE,
   Int_t   iTRD   =  1;
   Int_t   iZDC   =  1;
   Int_t   iEMCAL =  1;
-  Int_t   iACORDE = 0;
+  Int_t   iACORDE = 1;
   Int_t   iVZERO =  1;
-/*
-  Int_t   iABSO  =  0;
-  Int_t   iDIPO  =  0;
-  Int_t   iFMD   =  0;
-  Int_t   iFRAME =  0;
-  Int_t   iHALL  =  0;
-  Int_t   iITS   =  0;
-  Int_t   iMAG   =  0;
-  Int_t   iMUON  =  0;
-  Int_t   iPHOS  =  0;
-  Int_t   iPIPE  =  0;
-  Int_t   iPMD   =  0;
-  Int_t   iHMPID =  0;
-  Int_t   iSHIL  =  0;
-  Int_t   iT0    =  0;
-  Int_t   iTOF   =  0;
-  Int_t   iTPC   =  1;
-  Int_t   iTRD   =  0;
-  Int_t   iZDC   =  0;
-  Int_t   iEMCAL =  0;
-  Int_t   iACORDE = 0;
-  Int_t   iVZERO =  0;
-*/
+
   rl->CdGAFile();
   //=================== Alice BODY parameters =============================
   AliBODY *BODY = new AliBODY("BODY", "Alice envelop");
@@ -287,7 +259,7 @@ void commonConfig(Bool_t setRootGeometry = kFALSE,
   {
       //=================== EMCAL parameters ============================
     if ( configVersion == kConfigV0 ) 
-      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETE");
+      AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_COMPLETEV1");
     else if ( configVersion == kConfigV1 )  
       AliEMCAL *EMCAL = new AliEMCALv2("EMCAL", "EMCAL_FIRSTYEAR");
   }
index 8e2b7a8..d8344e8 100644 (file)
@@ -12,7 +12,7 @@ void Config()
   // AliRoot setup
   //
   gROOT->LoadMacro("$ALICE_ROOT/test/vmctest/ppbench/commonConfig.C");
-  commonConfig(kFALSE);
+  commonConfig();
 
   // Load Geant3 + Geant3 VMC libraries
   //
index 7b8c3b3..ed44b3c 100644 (file)
@@ -13,7 +13,7 @@ void Config()
   // AliRoot setup
   //
   gROOT->LoadMacro("$ALICE_ROOT/test/vmctest/ppbench/commonConfig.C");
-  commonConfig(kTRUE);
+  commonConfig();
 
   // Load Geant4 + Geant4 VMC libraries
   //
@@ -36,7 +36,7 @@ void Config()
                                  true);
 //      = new TG4RunConfiguration("geomRootToGeant4",
 //                                "QGSP_BERT_EMV+optical", 
-//                                "specialCuts+specialControls+stackPopper+stepLimiter",
+//                                "specialCuts+stackPopper+stepLimiter",
 //                                 true);
       
 
index 77ec40e..a925c98 100755 (executable)
@@ -3,7 +3,7 @@
 
 # Before running this script, you should run rungen.sh first.
 
-NEVENTS=10
+NEVENTS=3
 G3CONFIG="$ALICE_ROOT/test/vmctest/ppbench/g3Config.C" 
 G4CONFIG="$ALICE_ROOT/test/vmctest/ppbench/g4Config.C" 
 G3OUTDIR=g3
@@ -28,7 +28,6 @@ fi
 
 if [ "$RUNG4" = "1" ]; then 
   rm -rf *.root *.dat *.log fort* hlt hough raw* recraw/*.root recraw/*.log
-  cp g3/geometry.root .
   aliroot -b -q  sim.C\($NEVENTS,\""$G4CONFIG"\"\)  2>&1 | tee sim.log
   mv syswatch.log simwatch.log
   aliroot -b -q rec.C      2>&1 | tee rec.log
index 71e17fe..deb8805 100644 (file)
@@ -3,18 +3,14 @@
 // Macro for running simulation in test/vmctest/ppbench.
 // From test/ppbench. 
 
-void sim(Int_t nev=4, const TString& config) {
+void sim(Int_t nev=3, const TString& config) {
   if (gSystem->Getenv("EVENT"))
    nev = atoi(gSystem->Getenv("EVENT")) ;   
   
   AliSimulation simulator(config);
   simulator.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");
   simulator.SetMakeDigitsFromHits("ITS TPC");
-  //simulator.SetWriteRawData("ALL","raw.root",kTRUE);
-  //simulator.SetMakeDigits("");
-  //simulator.SetMakeSDigits("");
-  //simulator.SetMakeDigitsFromHits("");
-  simulator.SetWriteRawData("","raw.root",kTRUE);
+  simulator.SetWriteRawData("ALL","raw.root",kTRUE); 
 
   simulator.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
   simulator.SetSpecificStorage("GRP/GRP/Data",