]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - PWG0/multPbPb/run.sh
AliAnalysisMultPbCentralitySelector:
[u/mrichter/AliRoot.git] / PWG0 / multPbPb / run.sh
index 004fef69ae4c09c969e4c76a72618b237307aebf..8f22b10bd0cb43e6f66049b665dcf7bd19d36481 100755 (executable)
@@ -17,6 +17,20 @@ centrBin=0
 centrEstimator="V0M"
 runTriggerStudy=no
 customSuffix=""
+ntrackletsTrigger=50
+rejectBGV0Trigger=kFALSE
+useTrackCentralityCut=0
+trackMin=0
+trackMax=100
+dataDir=""
+mcDir=""
+vzMin=-10
+vzMax=10
+etaMin=-0.5
+etaMax=0.5
+npart=381.188
+weakFactor=-1
+useSingleBin=kTRUE
 
 give_help() {
 
@@ -30,7 +44,8 @@ Available options:
                                   0 local
                                   1 caf    
                                   2 grid    
-  -c                           Run the correction
+  -c <data,mc>                 Run the correction data and MC are names of the folders. 
+                               ./output/ is added automatically in front of the folder names
   -s                           Run the trigger study task (by default it runs the multiplicity analysis)
  Proof settings
   -w nworkers                  Set the number of worker nodes
@@ -41,21 +56,13 @@ Available options:
                                   grid or a text file with a ESD per line
                                 - caf mode: a dataset
                                 - grid mode: a directory on alien
+  -h                           This help
  Options specific to the multiplicity analysis
-  -b <bin>                     Set centrality bin [default=$centrBin]
-  -e <estimator>               Set centrality estimator [default=$centrEstimator]
-                               Available choiches:
-                                - V0M = V0 multiplicity
-                                - FMD = FMD raw multiplicity
-                                - TRK = N. of tracks
-                                - TKL = N. of tracklets
-                                - CL0 = N. of clusters in layer 0
-                                - V0MvsFMD = correlation between V0 and FMD
-                                - TKLvsV0 = correlation between tracklets and V0
-                                - ZEMvsZDC = correlation between ZEM and ZDC     
+  -l                           Run over all centrality bins
   -o <option>                  Misc option [default=$option]
                                Available options: 
                                 - SAVE:     Move results to a different output folder*
+                                - DCA:      Use DCA cut with global tracks
                                 - ITSsa:    Use ITSsa tracks
                                 - TPC:      Use TPC only tracks
                                 - NOMCKINE: Skip MC kinematics (runs way faster)
@@ -64,25 +71,90 @@ Available options:
   -m                           Use this to run on Monte Carlo
   -x  <suffix>                 Set a custom suffix in the histo manager        
   -g                           Debug mode
-  -h                           This help
+  == The following options are only valid if running on a single bin ==
+  -b <bin>                     Set centrality bin [default=$centrBin]
+  -e <estimator>               Set centrality estimator [default=$centrEstimator]
+                               Available choiches:
+                                - V0M = V0 multiplicity
+                                - FMD = FMD raw multiplicity
+                                - TRK = N. of tracks
+                                - TKL = N. of tracklets
+                                - CL0 = N. of clusters in layer 0
+                                - V0MvsFMD = correlation between V0 and FMD
+                                - TKLvsV0 = correlation between tracklets and V0
+                                - ZEMvsZDC = correlation between ZEM and ZDC     
+  -y <min,max>                 Select centrality based on "good tracks" rather than on centrality
+                               estimator [off by default]
+  -0 <min,max>                 Select centrality based on v0 amplitude range rather than on centrality
+                               estimator [off by default]
+  -2 <min,max>                 Select centrality based on SPD outer layer clusters  rather than on centrality
+                               estimator [off by default]
+
+ Options specific to the trigger study task
+  -k <ntracklets>              Max number of tracklets to fill eta and pt 
+                               distributions [default=$ntrackletsTrigger]
+  -v                           Reject BG with the V0
+ Options specific for the corrections
+  -z <zmin,zmax>               Change vertex Z range [default = $vzMin,$vzMax]
+  -a <etamin,etamax>           Change eta range [default = $etaMin,$etaMax]
+  -p <npart>                   Number of participants, used only for dNdeta/npart [default=$npart]
+  -k <weakFrac>                Scale ration secondaries from strangeness/all rec by this factor [default=$weakFactor]
 ENDOFGUIDE
 
 }
 
-while getopts "x:sr:cgmd:o:w:n:e:b:t:" opt; do
+while getopts "x:sr:c:gmd:o:w:n:e:b:t:k:vy:0:2:hz:a:l" opt; do
   case $opt in
     r)
       run=yes
       runmode=$OPTARG
       ;;      
+    l)
+      useSingleBin=kFALSE
+      ;;      
+    y)
+      useTrackCentralityCut=1
+      trackMin=${OPTARG%%,*}
+      trackMax=${OPTARG##*,}
+      ;;      
+    0)
+      useTrackCentralityCut=2
+      trackMin=${OPTARG%%,*}
+      trackMax=${OPTARG##*,}
+      ;;      
+    2)
+      useTrackCentralityCut=3
+      trackMin=${OPTARG%%,*}
+      trackMax=${OPTARG##*,}
+      ;;      
     x)
       customSuffix=$OPTARG
       ;;      
+    p)
+      npart=$OPTARG
+      ;;      
+    k)
+      ntrackletsTrigger=$OPTARG
+      weakFactor=$OPTARG
+      ;;      
     s)
       runTriggerStudy=yes
       ;;      
+    v)
+      rejectBGV0Trigger=kTRUE
+      ;;      
     c)
       correct=yes
+      dataDir="./output/${OPTARG%%,*}"
+      mcDir="./output/${OPTARG##*,}"
+      ;;      
+    z)
+      vzMin=${OPTARG%%,*}
+      vzMax=${OPTARG##*,}
+      ;;      
+    a)
+      etaMin=${OPTARG%%,*}
+      etaMax=${OPTARG##*,}
       ;;      
     g)
       debug=kTRUE;
@@ -140,13 +212,13 @@ if [ "$run" = "yes" ]
     then
     if [ "$runTriggerStudy" = "yes" ]
        then
-       root $ropt runTriggerStudy.C\(\"$dataset\",$nev,$offset,$debug,$runmode,$isMC,\"$option\",$workers\)
+       root $ropt runTriggerStudy.C\(\"$dataset\",$nev,$offset,$debug,$runmode,$isMC,$ntrackletsTrigger,$rejectBGV0Trigger,\"$option\",$workers\)
     else
-       root $ropt run.C\(\"$dataset\",$nev,$offset,$debug,$runmode,$isMC,$centrBin,\"$centrEstimator\",\"$option\",\"$customSuffix\",$workers\)
+       root $ropt run.C\(\"$dataset\",$nev,$offset,$debug,$runmode,$isMC,$centrBin,\"$centrEstimator\",$useTrackCentralityCut,$trackMin,$trackMax,\"$option\",\"$customSuffix\",$workers,$useSingleBin\)
     fi
 fi
 
 if [ "$correct" = "yes" ]
     then
-    echo "To be implemented"
+    root $ropt correct.C+\(\"$dataDir\",\"$mcDir\",$vzMin,$vzMax,$etaMin,$etaMax,$npart,$weakFactor\);
 fi