In vmctest/production:
authorivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 6 Jul 2010 13:52:49 +0000 (13:52 +0000)
committerivana <ivana@f7af4fe6-9843-0410-8265-dc069ae4e863>
Tue, 6 Jul 2010 13:52:49 +0000 (13:52 +0000)
Updated macros to be able to run with AliRoot v4-18-Rev-22 (and 23)

test/vmctest/production/Config.C
test/vmctest/production/rec.C
test/vmctest/production/sim.C
test/vmctest/production/simrun.C

index 9c11eda..6c1001a 100644 (file)
@@ -146,23 +146,8 @@ void Config()
   gAlice->SetTriggerDescriptor("p-p");
   //  AliSimulation::Instance()->SetTriggerConfig("p-p");
 
-  //
-  // FIELD
-  //
-  AliMagF* field = 0x0;
-  if (mag == kNoField) {
-    comment = comment.Append(" | L3 field 0.0 T");
-    field = new AliMagF("Maps","Maps", 0., 0., AliMagF::k5kGUniform,AliMagF::kBeamTypepp, energy/2.0);
-  } else if (mag == k5kG) {
-    comment = comment.Append(" | L3 field 0.5 T");
-    field = new AliMagF("Maps","Maps", -1., -1., AliMagF::k5kG,        AliMagF::kBeamTypepp, energy/2.0);
-  }
-
   printf("\n \n Comment: %s \n \n", comment.Data());
-
-  TGeoGlobalMagField::Instance()->SetField(field);
-    
   rl->CdGAFile();
   
   Int_t iABSO  = 1;
@@ -257,7 +242,7 @@ void Config()
       //============================ TPC parameters =====================
 
         AliTPC *TPC = new AliTPCv2("TPC", "Default");
-        TPC->SetPrimaryIonisation(); // not used with Geant3
+       // TPC->SetPrimaryIonisation(); // not used with Geant3
     }
 
 
@@ -324,27 +309,10 @@ void Config()
     {
         //=================== PHOS parameters ===========================
 
-        AliPHOS *PHOS = new AliPHOSv1("PHOS", "noCPV");
-        //Set simulation parameters different from the default ones.
-        AliPHOSSimParam* simEmc = AliPHOSSimParam::GetInstance() ;
-  
-        // APD noise of warm (+20C) PHOS:
-        // a2 = a1*(Y1/Y2)*(M1/M2), where a1 = 0.012 is APD noise at -25C,
-        // Y1 = 4.3 photo-electrons/MeV, Y2 = 1.7 p.e/MeV - light yields at -25C and +20C,
-        // M1 = 50, M2 = 50 - APD gain factors chosen for t1 = -25C and t2 = +20C,
-        // Y = MeanLightYield*APDEfficiency.
-
-        Float_t apdNoise = 0.012*2.5; 
-        simEmc->SetAPDNoise(apdNoise);
+       AliPHOS *PHOS = new AliPHOSv1("PHOS", "noCPV_Modules123");
 
-        //Raw Light Yield at +20C
-        simEmc->SetMeanLightYield(18800);
-
-        //ADC channel width at +18C.
-        simEmc->SetADCchannelW(0.0125);
     }
 
-
     if (iPMD)
     {
         //=================== PMD parameters ============================
@@ -440,8 +408,19 @@ void Config()
   geant4->ProcessGeantCommand("/mcVerbose/composedPhysicsList 2");  
   geant4->ProcessGeantCommand("/mcTracking/skipNeutrino true");
 
-
-
+  /*
+  // Set PAI model for TPC (TPC_Ne-CO2-N-2)
+  geant4->ProcessGeantCommand("/mcPhysics/emModel/selectMedium 219");
+  geant4->ProcessGeantCommand("/mcPhysics/emModel/setElossModel PAI");
+  geant4->ProcessGeantCommand("/mcPhysics/emModel/setFluctModel PAI");
+  geant4->ProcessGeantCommand("/mcPhysics/emModel/setParticles  all");
+  
+  // Set PAI model for TRD (TRD_XeCO2)
+  geant4->ProcessGeantCommand("/mcPhysics/emModel/selectMedium 291");
+  geant4->ProcessGeantCommand("/mcPhysics/emModel/setElossModel PAI");
+  geant4->ProcessGeantCommand("/mcPhysics/emModel/setFluctModel PAI");
+  geant4->ProcessGeantCommand("/mcPhysics/emModel/setParticles  all");
+  */
 
 
  //
index 8cc101a..a0d4c7e 100644 (file)
@@ -1,26 +1,90 @@
-// $Id$\r
-\r
-void rec() {\r
-\r
-  AliReconstruction reco;\r
-// switch off cleanESD\r
-  reco.SetCleanESD(kFALSE);\r
-\r
-  \r
-  reco.SetWriteESDfriend();\r
-  reco.SetWriteAlignmentData();\r
-  reco.SetRunQA("kFALSE:kFALSE");\r
-  \r
-  reco.SetRunPlaneEff(kTRUE);\r
-  reco.SetDefaultStorage("local://$ALICE_ROOT/OCDB");\r
-  \r
-  reco.SetSpecificStorage("GRP/GRP/Data",\r
-                          Form("local://%s",gSystem->pwd()));\r
-\r
-\r
-  TStopwatch timer;\r
-  timer.Start();\r
-  reco.Run();\r
-  timer.Stop();\r
-  timer.Print();\r
-}\r
+// $Id$
+
+void rec() {
+
+  AliReconstruction reco;
+
+//
+// switch off cleanESD, write ESDfriends and Alignment data
+
+  reco.SetCleanESD(kFALSE);
+  reco.SetWriteESDfriend();
+  reco.SetWriteAlignmentData();
+
+//
+// ITS Efficiency and tracking errors
+
+  reco.SetRunPlaneEff(kTRUE);
+  reco.SetUseTrackingErrorsForAlignment("ITS");
+
+//
+// Residual OCDB
+
+  reco.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Residual/");
+
+
+//
+// Vertex from RAW OCDB
+
+  reco.SetSpecificStorage("GRP/Calib/MeanVertexTPC","alien://folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("GRP/Calib/MeanVertex","alien://folder=/alice/data/2010/OCDB");
+
+
+//
+// EMCAL from RAW OCDB
+
+  reco.SetSpecificStorage("EMCAL/Calib/Data","alien://Folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("EMCAL/Calib/Pedestals","alien://Folder=/alice/data/2010/OCDB");
+
+//
+// PHOS from RAW OCDB
+
+  reco.SetSpecificStorage("PHOS/Calib/EmcBadChannels","alien://Folder=/alice/data/2010/OCDB");
+
+//
+// SPD and SDD from RAW OCDB
+
+  reco.SetSpecificStorage("ITS/Calib/SPDDead","alien://folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("TRIGGER/SPD/PITConditions","alien://folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("ITS/Calib/SPDNoise","alien://folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("ITS/Calib/CalibSDD","alien://Folder=/alice/data/2010/OCDB");
+
+//
+// TRD from RAW OCDB
+
+  reco.SetSpecificStorage("TRD/Calib/ChamberStatus","alien://folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("TRD/Calib/PadStatus","alien://folder=/alice/data/2010/OCDB");
+
+//
+// TPC from RAW OCDB
+
+  reco.SetSpecificStorage("TPC/Calib/PadGainFactor","alien://folder=/alice/data/2010/OCDB");
+
+//
+// V0 from RAW OCDB
+
+  reco.SetSpecificStorage("VZERO/Trigger/Data","alien://folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("VZERO/Calib/RecoParam","alien://folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("VZERO/Calib/Data","alien://folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("VZERO/Calib/TimeSlewing","alien://folder=/alice/data/2010/OCDB");
+  reco.SetSpecificStorage("VZERO/Calib/TimeDelays","alien://folder=/alice/data/2010/OCDB");
+
+//
+// TOF from RAW OCDB
+
+//  reco.SetSpecificStorage("TOF/Calib/Status","alien://folder=/alice/data/2010/OCDB");
+
+//
+// QA off
+
+  reco.SetRunQA(":");
+
+//
+// The rest
+  TStopwatch timer;
+  timer.Start();
+  reco.Run();
+  timer.Stop();
+  timer.Print();
+}
index 83e8f1c..575d6b2 100644 (file)
@@ -1,24 +1,68 @@
-// $Id$\r
-\r
-void sim(Int_t nev=20) {\r
-\r
-  AliSimulation simulator;\r
-  simulator.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");\r
-  simulator.SetMakeDigitsFromHits("ITS TPC");\r
-  simulator.SetRunHLT("");\r
-  \r
-  // The raw data are not written due to the huge increase of the \r
-  // virtual memory in HLT\r
-  // simulator.SetWriteRawData("ALL","raw.root",kTRUE);\r
-  \r
-  simulator.SetDefaultStorage("local://$ALICE_ROOT/OCDB");\r
\r
-  simulator.SetSpecificStorage("GRP/GRP/Data",\r
-                                Form("local://%s",gSystem->pwd()));\r
-  \r
-  TStopwatch timer;\r
-  timer.Start();\r
-  simulator.Run(nev);\r
-  timer.Stop();\r
-  timer.Print();\r
-}\r
+// $Id$
+
+void sim(Int_t nev=20) {
+
+  AliSimulation simulator;
+  simulator.SetMakeSDigits("TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0 VZERO");
+  simulator.SetMakeDigitsFromHits("ITS TPC");
+  
+//
+// Ideal OCDB
+
+  simulator.SetDefaultStorage("alien://Folder=/alice/simulation/2008/v4-15-Release/Ideal/");
+
+//
+// Mean verrtex from RAW OCDB
+
+ simulator.SetSpecificStorage("GRP/Calib/MeanVertexSPD", "alien://folder=/alice/data/2010/OCDB");
+
+//
+// SDD from RAW OCDB
+
+ simulator.SetSpecificStorage("ITS/Calib/CalibSDD","alien://Folder=/alice/data/2010/OCDB");
+
+//
+// EMCAL from RAW OCDB
+
+  simulator.SetSpecificStorage("EMCAL/Calib/Data","alien://Folder=/alice/data/2010/OCDB");
+
+//
+// TRD from RAW OCDB
+
+  simulator.SetSpecificStorage("TRD/Calib/ChamberStatus","alien://folder=/alice/data/2010/OCDB");
+  simulator.SetSpecificStorage("TRD/Calib/PadStatus","alien://folder=/alice/data/2010/OCDB");
+
+//
+// V0 from RAW OCDB
+
+  simulator.SetSpecificStorage("VZERO/Trigger/Data","alien://folder=/alice/data/2010/OCDB");
+  simulator.SetSpecificStorage("VZERO/Calib/RecoParam","alien://folder=/alice/data/2010/OCDB");
+  simulator.SetSpecificStorage("VZERO/Calib/Data","alien://folder=/alice/data/2010/OCDB");
+  simulator.SetSpecificStorage("VZERO/Calib/TimeSlewing","alien://folder=/alice/data/2010/OCDB");
+  simulator.SetSpecificStorage("VZERO/Calib/TimeDelays","alien://folder=/alice/data/2010/OCDB");
+
+//
+// TOF from RAW OCDB
+
+//  simulator.SetSpecificStorage("TOF/Calib/Status","alien://folder=/alice/data/2010/OCDB");
+
+//
+// Read GRP Data from RAW
+
+  simulator.SetSpecificStorage("GRP/GRP/Data", "alien://Folder=/alice/data/2010/OCDB");
+
+//
+// Vertex and Mag.field from OCDB
+
+  simulator.UseVertexFromCDB();
+  simulator.UseMagFieldFromGRP();
+
+//
+// The rest
+
+  TStopwatch timer;
+  timer.Start();
+  simulator.Run(nev);
+  timer.Stop();
+  timer.Print();
+}
index aed538b..8e83339 100644 (file)
@@ -6,6 +6,9 @@
 // Set job and simulation variables as :\r
 // root.exe -b -q simrun.C  --run <x> --event <y> --process <kPythia6/kPhojet/kPythia6ATLAS_Flat/kPythia6D6T> --field <kNoField/k5kG> --energy <900/2360/10000> --physicslist <QGSP_BERT_CHIPS/CHIPS/QGSP_BERT_EMV>\r
 //\r
+// For LHC10c9 production test:\r
+// root.exe -b -q simrun.C  --run 119846 --event 1 --process kPythia6D6T --field k5kG --energy 7000 --physicslist QGSP_BERT_CHIPS_OPTICAL\r
+//\r
 // By E. Sicking, CERN\r
 {\r
   int nrun = 0;\r
@@ -53,6 +56,7 @@
   }\r
 \r
   seed = nrun * 100000 + nevent;\r
+  //seed = nrun + nevent;   // LHC10c9 production test \r
   sprintf(sseed,"%d",seed);\r
 \r
   if (seed==0) {\r