Working on the electron cut
[u/mrichter/AliRoot.git] / PWGLF / SPECTRA / ChargedHadrons / multPbPb / run.sh
1 #!/bin/bash
2
3 # defaults
4 isMC=kFALSE
5 run=no
6 correct=no
7 nev=-1
8 offset=0
9 debug=kFALSE
10 runmode=1
11 dataset=/alice/sim/LHC10f8a_130844
12 ropt="-l"
13 option="DCA,SAVE"
14 workers=26
15 analysismode=9; #SPD + field on
16 centrBin=-1
17 centrEstimator="V0M"
18 runTriggerStudy=no
19 customSuffix=""
20 ntrackletsTrigger=50
21 rejectBGV0Trigger=kFALSE
22 useTrackCentralityCut=0
23 trackMin=0
24 trackMax=100
25 dataDir=""
26 mcDir=""
27 vzMin=-10
28 vzMax=10
29 etaMin=-0.5
30 etaMax=0.5
31 npart=381.188
32 weakFactor=-1
33 useSingleBin=kTRUE
34 #OUTPATH=output.BAK2010
35 OUTPATH=output
36
37 give_help() {
38
39 cat <<ENDOFGUIDE
40 This scripts runs the mupliplicity analysis and the trigger study task
41
42 Available options:
43  Mode control, at least one of the following options should be used
44   -r <mode>                    Run the task
45                                Modes [default=$runmode]:
46                                   0 local
47                                   1 caf    
48                                   2 grid    (remeber to set run list)
49                                   3 prooflite
50   -c <data,mc>                 Run the correction data and MC are names of the folders. 
51                                ./output/ is added automatically in front of the folder names
52   -s                           Run the trigger study task (by default it runs the multiplicity analysis)
53  Proof settings
54   -w nworkers                  Set the number of worker nodes (0 == 1 worker per node)
55   -n <nev>                     Number of events to be analized
56  Grid Settings
57   -u <runs>                    Runs to be processed in the data folder
58  Misc
59   -d <dataset>                 Dataset or data collection (according to run mode) [default=$dataset]
60                                 - local mode: a single ESD file, an xml collection of files on 
61                                   grid or a text file with a ESD per line
62                                 - caf mode: a dataset
63                                 - grid mode: a directory on alien (don't forget the run list)
64   -h                           This help
65  Options specific to the multiplicity analysis
66   -l                           Run over all centrality bins
67   -o <option>                  Misc option [default=$option]
68                                Available options: 
69                                 - SAVE:     Move results to a different output folder*
70                                 - DCA:      Use DCA cut with global tracks
71                                 - ITSsa:    Use ITSsa tracks
72                                 - TPC:      Use TPC only tracks
73                                 - NOMCKINE: Skip MC kinematics (runs way faster)
74                                 - NoElectrons: apply cuts to reject electrons in the rec sample, 
75                                   exclude electrons by PDG code in the MC sample
76                                 * == can be used in trigger studies task
77   -t <option>                  Command line option for root [defaul=$ropt]
78   -m                           Use this to run on Monte Carlo
79   -x  <suffix>                 Set a custom suffix in the histo manager        
80   -g                           Debug mode
81   == The following options are only valid if running on a single bin ==
82   -b <bin>                     Set centrality bin [default=$centrBin]
83   -e <estimator>               Set centrality estimator [default=$centrEstimator]
84                                Available choiches:
85                                 - V0M = V0 multiplicity
86                                 - FMD = FMD raw multiplicity
87                                 - TRK = N. of tracks
88                                 - TKL = N. of tracklets
89                                 - CL0 = N. of clusters in layer 0
90                                 - V0MvsFMD = correlation between V0 and FMD
91                                 - TKLvsV0 = correlation between tracklets and V0
92                                 - ZEMvsZDC = correlation between ZEM and ZDC     
93   -y <min,max>                 Select centrality based on "good tracks" rather than on centrality
94                                estimator [off by default]
95   -0 <min,max>                 Select centrality based on v0 amplitude range rather than on centrality
96                                estimator [off by default]
97   -2 <min,max>                 Select centrality based on SPD outer layer clusters  rather than on centrality
98                                estimator [off by default]
99
100  Options specific to the trigger study task
101   -k <ntracklets>              Max number of tracklets to fill eta and pt 
102                                distributions [default=$ntrackletsTrigger]
103   -v                           Reject BG with the V0
104  Options specific for the corrections
105   -z <zmin,zmax>               Change vertex Z range [default = $vzMin,$vzMax]
106   -a <etamin,etamax>           Change eta range [default = $etaMin,$etaMax]
107   -p <npart>                   Number of participants, used only for dNdeta/npart [default=$npart]
108   -k <weakFrac>                Scale ration secondaries from strangeness/all rec by this factor [default=$weakFactor]
109   -b <bin>                     Set centrality bin to be corrected. Only valid if you processed multiple 
110                                bins at one (it changes the suffix of the multPbPbtracks.root file). It it's -1,
111                                a file without suffix is searched for. This options applyies both to the data and to the 
112                                MC file. [default=$centrBin]
113 ENDOFGUIDE
114
115 }
116
117 while getopts "x:sr:c:gmd:o:w:n:e:b:t:k:vy:0:2:hz:a:lp:u:" opt; do
118   case $opt in
119     r)
120       run=yes
121       runmode=$OPTARG
122       ;;      
123     l)
124       useSingleBin=kFALSE
125       ;;
126     u) 
127       runList=$OPTARG
128       ;;
129     y)
130       useTrackCentralityCut=1
131       trackMin=${OPTARG%%,*}
132       trackMax=${OPTARG##*,}
133       ;;      
134     0)
135       useTrackCentralityCut=2
136       trackMin=${OPTARG%%,*}
137       trackMax=${OPTARG##*,}
138       ;;      
139     2)
140       useTrackCentralityCut=3
141       trackMin=${OPTARG%%,*}
142       trackMax=${OPTARG##*,}
143       ;;      
144     x)
145       customSuffix=$OPTARG
146       ;;      
147     p)
148       npart=$OPTARG
149       ;;      
150     k)
151       ntrackletsTrigger=$OPTARG
152       weakFactor=$OPTARG
153       ;;      
154     s)
155       runTriggerStudy=yes
156       ;;      
157     v)
158       rejectBGV0Trigger=kTRUE
159       ;;      
160     c)
161       correct=yes
162       dataDir="./$OUTPATH/${OPTARG%%,*}"
163       mcDir="./$OUTPATH/${OPTARG##*,}"
164       ;;      
165     z)
166       vzMin=${OPTARG%%,*}
167       vzMax=${OPTARG##*,}
168       ;;      
169     a)
170       etaMin=${OPTARG%%,*}
171       etaMax=${OPTARG##*,}
172       ;;      
173     g)
174       debug=kTRUE;
175       ;;
176     m)
177       isMC=kTRUE
178       ;;
179     d)
180       dataset=$OPTARG
181      ;;
182     e)
183       centrEstimator=$OPTARG
184      ;;
185     b)
186       centrBin=$OPTARG
187      ;;
188     o)
189       option=$OPTARG
190       ;;
191     t)
192       ropt=$OPTARG
193       ;;
194     w) 
195       workers=$OPTARG
196       ;;
197     n) 
198       nev=$OPTARG
199       ;;
200     h)
201       give_help
202       exit 1
203       ;;
204     \?)
205       echo "Invalid option: -$OPTARG" >&2
206       give_help
207       exit 1
208       ;;
209     :)
210       echo "Option -$OPTARG requires an argument." >&2
211       give_help
212       exit 1
213       ;;
214   esac
215 done
216
217 if [ "$run" = "$correct" ]
218     then 
219     echo "One and only one option between -r and -c must be selected"
220     give_help
221     exit 1
222 fi
223
224
225 if [ "$run" = "yes" ]
226     then
227     if [ "$runTriggerStudy" = "yes" ]
228         then
229         root $ropt runTriggerStudy.C\(\"$dataset\",$nev,$offset,$debug,$runmode,$isMC,$ntrackletsTrigger,$rejectBGV0Trigger,\"$option\",$workers\)
230     else
231         root $ropt run.C\(\"$dataset\",$nev,$offset,$debug,$runmode,$isMC,$centrBin,\"$centrEstimator\",$useTrackCentralityCut,$trackMin,$trackMax,\"$option\",\"$customSuffix\",$workers,$useSingleBin,\"$runList\"\)
232     fi
233 fi
234
235 if [ "$correct" = "yes" ]
236     then
237     root $ropt correct.C+\(\"$dataDir\",\"$mcDir\",$vzMin,$vzMax,$etaMin,$etaMax,$npart,$weakFactor,$centrBin\);
238 fi