ATO-17 - adding special testdEdx option to the Config.C + paremetrized ntracks and...
authormivanov <marian.ivanov@cern.ch>
Fri, 20 Jun 2014 10:28:57 +0000 (12:28 +0200)
committermivanov <marian.ivanov@cern.ch>
Fri, 20 Jun 2014 10:28:57 +0000 (12:28 +0200)
test/testdEdx/Config.C
test/testdEdx/rec.C
test/testdEdx/runtest.sh
test/testdEdx/sim.C

index f0e5906..94b9603 100644 (file)
@@ -2,7 +2,7 @@
 
 enum PprRun_t 
 {
-    test50,
+  test50, testdEdx,
     kParam_8000,   kParam_4000,  kParam_2000, 
     kHijing_cent1, kHijing_cent2, 
     kHijing_per1,  kHijing_per2, kHijing_per3, kHijing_per4,  kHijing_per5,
@@ -23,7 +23,7 @@ enum PprRun_t
 };
 
 const char* pprRunName[] = {
-    "test50",
+  "test50", "testdEdx",
     "kParam_8000",   "kParam_4000",  "kParam_2000", 
     "kHijing_cent1", "kHijing_cent2", 
     "kHijing_per1",  "kHijing_per2", "kHijing_per3", "kHijing_per4",  
@@ -61,7 +61,7 @@ const char * pprTrigConfName[] = {
 // This part for configuration    
 
 //static PprRun_t srun = kPythia6;
-static PprRun_t srun = test50;
+static PprRun_t srun = testdEdx;
 static PprRad_t srad = kGluonRadiation;
 static AliMagF::BMap_t smag = AliMagF::k5kG;
 static Int_t    sseed = 12345; //Set 0 to use the current time
@@ -421,6 +421,24 @@ AliGenerator* GeneratorFactory(PprRun_t srun) {
        gGener=gener;
       }
       break;
+    case testdEdx:
+      {
+        // generator for the dEdx simulation study
+       gSystem->Getenv("TestdEdxNTracks");
+       Int_t ntracks =50;
+        if (gSystem->Getenv("TestdEdxNTracks")) ntracks= atoi(gSystem->Getenv("TestdEdxNTracks"));
+       comment = comment.Append(":HIJINGparam test N particles");
+       AliGenHIJINGpara *gener = new AliGenHIJINGpara(ntracks);
+       gener->SetMomentumRange(0, 999999.);
+       gener->SetPhiRange(0., 360.);
+       // Set pseudorapidity range from -8 to 8.
+       Float_t thmin = EtaToTheta(2);   // theta min. <---> eta max
+       Float_t thmax = EtaToTheta(-2);  // theta max. <---> eta min 
+       gener->SetThetaRange(thmin,thmax);
+       gGener=gener;
+      }
+      break;
+
     case kParam_8000:
       {
        comment = comment.Append(":HIJINGparam N=8000");
index bf8446d..ae01a08 100644 (file)
@@ -9,15 +9,11 @@ void rec() {
                          Form("local://%s",gSystem->pwd()));
   reco.SetRunPlaneEff(kTRUE);
 
-  reco.SetRunQA("ALL:ALL") ;
-  
-  reco.SetQARefDefaultStorage("local://$ALICE_ROOT/QAref") ;
-  
-  for (Int_t det = 0 ; det < AliQA::kNDET ; det++) {
-    reco.SetQACycles((AliQAv1::DETECTORINDEX_t)det, 999) ;
-    reco.SetQAWriteExpert((AliQAv1::DETECTORINDEX_t)det) ; 
-  }
   
+  reco.SetRunQA(":"); 
+  reco.SetRunGlobalQA(kFALSE);
+
+
   TStopwatch timer;
   timer.Start();
   reco.Run();
index 41ff65b..8ef1ec7 100755 (executable)
@@ -1,20 +1,19 @@
 #!/bin/sh
-# Root
-# export ROOTSYS=/afs/cern.ch/alice/library/root/new
-# export PATH=$ROOTSYS/bin:$PATH
-# export LD_LIBRARY_PATH=$ROOTSYS/lib:$LD_LIBRARY_PATH
-# AliRoot
-# export ALICE=/afs/cern.ch/alice/library
-# export ALICE_LEVEL=new
-# export ALICE_ROOT=$ALICE/$ALICE_LEVEL
-# export ALICE_TARGET=`$ROOTSYS/bin/root-config --arch`
-# export PATH=$ALICE_ROOT/bin/tgt_${ALICE_TARGET}:$PATH
-# export LD_LIBRARY_PATH=$ALICE_ROOT/lib/tgt_${ALICE_TARGET}:$LD_LIBRARY_PATH
-# Geant3
-# export LD_LIBRARY_PATH=$ALICE/geant3/lib/tgt_${ALICE_TARGET}:$LD_LIBRARY_PATH
+
+#
+# Here we should write waht is the aim of test
+#
+
+export TestdEdxNTracks=$1
+nevents=$2
+echo Running dEdx digitzer test job 
+echo NEvents $nevents
+echo NTracks per event  $TestdEdxNTracks
+
 
 rm -rf *.root *.dat *.log fort* hlt hough raw* recraw/*.root recraw/*.log
-aliroot -b -q $1sim.C      2>&1 | tee sim.log
+echo aliroot -b -q sim.C\($nevents\)     
+aliroot -b -q sim.C\($nevents\)      2>&1 | tee sim.log
 mv syswatch.log simwatch.log
 aliroot -b -q $1rec.C      2>&1 | tee rec.log
 mv syswatch.log recwatch.log
@@ -27,5 +26,20 @@ aliroot -b -q rec.C      2>&1 | tee rec.log
 mv syswatch.log ../rawwatch.log
 aliroot -b -q aod.C 2>&1 | tee aod.log
 
+exit;
+
 
+submitMultiplicityScan(){
+#
+# Here we submit the jobs for the simulation//reconstruction
+# Parameters:
+#   nbins to be investigated   (default 5)
+#   max multiplicity per event (default central PbPb - 15000)
+#   minimal total statistic of the tracks (default 4 central events - 60000 tracks)
+# Jobs will be submitted per event     
+#
+# For each setting new directory will be created - indicating muiltiplicity
+# dir<ntracks>/dir<eventNr>  
 
+  
+}
index 48eebe1..58d7ba8 100644 (file)
@@ -19,13 +19,7 @@ void sim(Int_t nev=20) {
   simulator.SetSpecificStorage("GRP/GRP/Data",
                               Form("local://%s",gSystem->pwd()));
   
-  simulator.SetRunQA("ALL:ALL") ; 
-  
-  simulator.SetQARefDefaultStorage("local://$ALICE_ROOT/QAref") ;
-
-  for (Int_t det = 0 ; det < AliQA::kNDET ; det++) {
-    simulator.SetQACycles((AliQAv1::DETECTORINDEX_t)det, nev+1) ;
-  }
+  simulator.SetRunQA(":") ; 
   
   TStopwatch timer;
   timer.Start();