modified: rec.C == changed name for the RecoParam OCDB Specific storage
authormivanov <marian.ivanov@cern.ch>
Thu, 3 Jul 2014 19:14:44 +0000 (21:14 +0200)
committermivanov <marian.ivanov@cern.ch>
Thu, 3 Jul 2014 19:14:44 +0000 (21:14 +0200)
    modified:   sim.C   == generation of RecoParam OCDB entry added
    modified:   submitSimJobs.sh == echo running commands which are executed

git reset HEAD <file>..." to unstage)

test/testdEdx/rec.C
test/testdEdx/sim.C
test/testdEdx/submitSimJobs.sh

index 0b540cd..7e0900e 100644 (file)
@@ -1,3 +1,7 @@
+const char * recoStorage="local:///cvmfs/alice.gsi.de/alice/data/2010/OCDB";
+Int_t run=0;
+
+
 void ModifyRecoParam(TObjArray* recoArray, Bool_t useIonTail, Double_t crossTalkCorrection){
   //
   // Modify reco param - and store it in the OCDB in local directory
@@ -8,7 +12,7 @@ void ModifyRecoParam(TObjArray* recoArray, Bool_t useIonTail, Double_t crossTalk
     p->SetUseIonTailCorrection(useIonTail);
     p->SetCrosstalkCorrection(crossTalkCorrection);
   }
-  TString localStorage = "local://"+gSystem->GetFromPipe("pwd")+"/OCDB";
+  TString localStorage = "local://"+gSystem->GetFromPipe("pwd")+"/OCDBrec";
   AliCDBStorage*pocdbStorage = AliCDBManager::Instance()->GetStorage(localStorage.Data());  
   AliCDBMetaData *metaData= new AliCDBMetaData();
   AliCDBId*   id1=new AliCDBId("TPC/Calib/RecoParam/", man->GetRun(), man->GetRun());
@@ -22,8 +26,6 @@ void rec(Bool_t useIonTail, Double_t crossTalkCorrection) {
   // run reconstruction
   // Parameters:
   //   useIonTail - switch for using ion tail correction - OCDB entry will be overritten in  
-  Int_t run=0;
-  const char * recoStorage="local:///cvmfs/alice.gsi.de/alice/data/2010/OCDB";
   //
   // stard reco setting
   //
@@ -66,7 +68,5 @@ void rec(Bool_t useIonTail, Double_t crossTalkCorrection) {
   for (Int_t i=0; i<4; i++){
     AliTPCRecoParam* recoParam = ( AliTPCRecoParam*)recoArray->At(i);
     printf("ipar=%d\t%d\t%f\n",i,recoParam->GetUseIonTailCorrection(), recoParam->GetCrosstalkCorrection());
-  }
-
-
+  } 
 }
index 58d7ba8..2d2af14 100644 (file)
@@ -1,4 +1,30 @@
-void sim(Int_t nev=20) {
+//
+// parameter to take from config file
+//
+const char * recoStorage="local:///cvmfs/alice.gsi.de/alice/data/2010/OCDB";
+Int_t run=0;
+
+
+void ModifyRecoParam(TObjArray* recoArray, Bool_t useIonTail, Double_t crossTalkCorrection){
+  //
+  // Modify reco param - and store it in the OCDB in local directory
+  //
+  AliCDBManager * man  =  AliCDBManager::Instance();
+  for (Int_t i=0; i<4; i++){
+    AliTPCRecoParam* p = ( AliTPCRecoParam*)recoArray->At(i);
+    p->SetUseIonTailCorrection(useIonTail);
+    p->SetCrosstalkCorrection(crossTalkCorrection);
+  }
+  TString localStorage = "local://"+gSystem->GetFromPipe("pwd")+"/OCDBsim";
+  AliCDBStorage*pocdbStorage = AliCDBManager::Instance()->GetStorage(localStorage.Data());  
+  AliCDBMetaData *metaData= new AliCDBMetaData();
+  AliCDBId*   id1=new AliCDBId("TPC/Calib/RecoParam/", man->GetRun(), man->GetRun());
+  pocdbStorage->Put(recoArray, (*id1), metaData);
+  AliCDBManager::Instance()->SetSpecificStorage("TPC/Calib/RecoParam/",localStorage.Data());
+}
+
+
+void sim(Int_t nev, Bool_t useIonTail, Double_t crossTalkCorrection) {
   gSystem->Load("liblhapdf");
   gSystem->Load("libEGPythia6");
   gSystem->Load("libpythia6");
@@ -7,8 +33,17 @@ void sim(Int_t nev=20) {
   gSystem->Load("libTHijing");
   gSystem->Load("libgeant321");
 
-  if (gSystem->Getenv("EVENT"))
-   nev = atoi(gSystem->Getenv("EVENT")) ;   
+  AliCDBManager * man = AliCDBManager::Instance();
+  man->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
+  man->SetSpecificStorage("TPC/Calib/RecoParam/",recoStorage);
+  man->SetRun(run);
+  AliCDBEntry* e = man->Get("TPC/Calib/RecoParam/",run); // get default
+  // modify content
+  TObjArray* recoArray = (TObjArray*)e->GetObject();
+  ModifyRecoParam(recoArray, useIonTail, crossTalkCorrection);
+
+
+  if (gSystem->Getenv("EVENT")) nev = atoi(gSystem->Getenv("EVENT")) ;   
   
   AliSimulation simulator;
   simulator.SetMakeSDigits("ITS TPC TRD TOF PHOS HMPID EMCAL MUON FMD ZDC PMD T0");
@@ -16,8 +51,7 @@ void sim(Int_t nev=20) {
   simulator.SetWriteRawData("ALL","raw.root",kTRUE);
 
   simulator.SetDefaultStorage("local://$ALICE_ROOT/OCDB");
-  simulator.SetSpecificStorage("GRP/GRP/Data",
-                              Form("local://%s",gSystem->pwd()));
+  simulator.SetSpecificStorage("GRP/GRP/Data", Form("local://%s",gSystem->pwd()));
   
   simulator.SetRunQA(":") ; 
   
@@ -26,4 +60,16 @@ void sim(Int_t nev=20) {
   simulator.Run(nev);
   timer.Stop();
   timer.Print();
+ //
+  // Print the OCDB setup which we used
+  //
+  AliCDBEntry* ocdbEntry = man->Get("TPC/Calib/RecoParam/",run);
+  TObjArray* recoArray = (TObjArray*)ocdbEntry->GetObject();
+  for (Int_t i=0; i<4; i++){
+    AliTPCRecoParam* recoParam = ( AliTPCRecoParam*)recoArray->At(i);
+    printf("ipar=%d\t%d\t%f\n",i,recoParam->GetUseIonTailCorrection(), recoParam->GetCrosstalkCorrection());
+  } 
+
+
+   
 }
index ff2baad..cc9173d 100644 (file)
@@ -13,6 +13,7 @@
 # global variables to be set
 DATE=`date +%Y_%m_%d_%H`
 nEvents=1               #fixed to 1 to make an event by event job submission
+
 ###############################################################################
 submitIonTailXTalkScan()
 {
@@ -94,19 +95,26 @@ runSim(){
 
 
   ## main body of the simulation part
-  rm -rf *.root *.dat *.log fort* hlt hough raw* recraw/*.root recraw/*.log
-  echo aliroot -b -q sim.C\($nevents,$ionTail,$xTalk\)     
+  rm -rf *.root *.dat *.log fort* hlt hough raw* recraw/*.root recraw/*.log GRP* 
+  printf   "\n ======================================================================\n\n"
+  echo Running: aliroot -b -q sim.C\($nevents,$ionTail,$xTalk\)     
   aliroot -b -q sim.C\($nevents,$ionTail,$xTalk\)            2>&1 | tee sim.log
   mv syswatch.log simwatch.log
-  echo aliroot -b -q rec.C\($ionTail\,$xTalk\) 
+  printf   "\n ======================================================================\n\n"
+  echo Running: aliroot -b -q rec.C\($ionTail\,$xTalk\) 
   aliroot -b -q rec.C\($ionTail\,$xTalk\)    2>&1 | tee rec.log    
   mv syswatch.log recwatch.log
-
   ## OCDB entries to be dumped in human readable format
   source $ALICE_ROOT/PWGPP/CalibMacros/AliOCDBtoolkit.sh
+  printf   "\n ======================================================================\n\n"
+  echo Running: ocdbMakeTable AliESDs.root "ESD" OCDBrec.list
   ocdbMakeTable AliESDs.root "ESD" OCDBrec.list
+  printf   "\n ======================================================================\n\n"
+  echo Running: ocdbMakeTable galice.root MC OCDBsim.list
   ocdbMakeTable galice.root MC OCDBsim.list
   ocdbFileName=$(cat OCDBrec.list | grep "TPC/Calib/RecoParam" | gawk '{print $2"/"$3}' )
+  printf   "\n ======================================================================\n\n"
+  echo Running: dumpObject $ocdbFileName  "object" "XML" RecoParam
   dumpObject $ocdbFileName  "object" "XML" RecoParam
 
   return 1;