StRoot  1
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Groups Pages
Cut.cxx
1 
2 #include "Cut.h"
3 
4 #include <iostream>
5 
6 ClassImp(Cut)
7 
8 int parseNumber(int a, int b);
9 
10 //
11 // initiatlize
12 //
13 NchCentrality Cut::mZdcCtbCent[2] = {kTotal,kFive};
14 Int_t Cut::mFlowCent[2] = {-1,-1};
15 Int_t Cut::mHMinusCent[2] = {-1,-1};
16 Int_t Cut::mNchCent[2] = {-1,1};
17 
18 // set all to false as default.
19 bool Cut::mDoFlowCent=false;
20 bool Cut::mDoZdcCtbCent=false;
21 bool Cut::mDoHMinusCent=false;
22 bool Cut::mDoNchCent=false;
23 bool Cut::mDoNchCentKludge=false;
24 
25 bool Cut::mDoSpectraCent=false; // default is not to use the spectra centrality
26 
27 Float_t Cut::mVertexZ[2] = {-200,200};
28 Float_t Cut::mVertexZSkip = 0;
29 
30 Float_t Cut::mEta[2] = {-.5,.5};
31 Int_t Cut::mFitPts[2] = {0,45};
32 Float_t Cut::mSDcaGl[2] = {-3,3};
33 Float_t Cut::mDcaPr[2] = {0,3};
34 Float_t Cut::mDcaXYGl[2] = {-3,3};
35 Int_t Cut::mMcPts[2] = {10,99};
36 Float_t Cut::mEtaTight[2] = { -.2, .2 };
37 Float_t Cut::mHitAvoid=0;
38 Int_t Cut::mCut = 0;
39 
40 char Cut::mHalf = 0;
41 char Cut::mHitHalf = 0;
42 char Cut::mGeomHalf = 0;
43 
44 
45 
46 
47 void
48 Cut::SetHalf(char half)
49 {
50  switch(half){
51  case 'e' : mHalf = 'e'; mVertexZ[1] = 0; break;
52  case 'w' : mHalf = 'w'; mVertexZ[0] = 0; break;
53  case 's' : mHalf = 's'; break; // enforce same side
54  default : cerr << "Wrong half " << half << endl; exit(1);
55  }
56 }
57 
58 void
59 Cut::SetHitHalf(char half)
60 {
61  switch(half){
62  case 'e' : mHitHalf = 'e'; break;
63  case 'w' : mHitHalf = 'w'; break;
64  case 's' : mHitHalf = 's'; break;
65  default : cerr << "Wrong hit half " << half << endl; exit(1);
66  }
67 }
68 
69 void
70 Cut::SetGeomHalf(char half)
71 {
72  switch(half){
73  case 'e' : mGeomHalf = 'e'; break;
74  case 'w' : mGeomHalf = 'w'; break;
75  case 's' : mGeomHalf = 's'; break;
76  default : cerr << "Wrong geom half " << half << endl; exit(1);
77  }
78 }
79 
80 /*
81 
82  Cuts defined as :
83  a. centrality
84  9=minbias
85  8=central
86  7=peripheral
87  6=central (version 2 (uses 5% hminus cut for other))
88  5=peripheral kludge (use 50-80 for corrections)
89  4=10-20%
90  3=20-30%
91  2=30-40%
92  1=40-60%
93  b. which centrality definition for spectra
94  9=flow
95  8=zdc-ctb
96  7=hminus
97  6=nch
98  c. which centrality definition for corrections/other
99  9=flow
100  8=zdc-ctb
101  7=hminus
102  6=nch
103  d. vertex
104  9=-200,200
105  8=-75,75
106  7=-95,95
107  6=-30,30
108  5=0,30
109  e. half?
110  9=none
111  8=hit avoid CM
112  7=hit avoid 10 cm around CM
113  f. fit pts
114  9=20-45
115  8=25
116  7=10
117  6=35 // really tight
118  5=24
119  4=23
120  3=30
121  g. dca 3d
122  9=1
123  8=1.2
124  7=3
125  h. eta
126  9=.5
127  8=.7
128  7= >0, <.5
129 */
130 /*
131 e.g. a=15796, b=5 : return = 1
132  a=5924, b=3 : return = 9, etc.
133 */
134 int parseNumber(int a, int b)
135 {
136  int div1=(int)pow(10,b);
137  int div2=(int)pow(10,b-1);
138 
139  int val=a%div1;
140  return (val==a) ? int(a/div2) : int(val/div2);
141 }
142 
143 void
144 Cut::SetCut(Int_t type)
145 {
146  if(type<10000000){
147  cout << "type " << type << " is too short" << endl;exit(1);
148  }
149 
150  mCut=type;
151 
152  int centType=parseNumber(type,8);
153  int centDefSpectra=parseNumber(type,7);
154  int centDefOther=parseNumber(type,6);
155  int vertex=parseNumber(type,5);
156  int half=parseNumber(type,4);
157  int fitPts=parseNumber(type,3);
158  int dca=parseNumber(type,2);
159  int eta=parseNumber(type,1);
160 
161  cout << "Set cut?" << endl
162  << "\tcentType=" << centType << endl
163  << "\tcentDefSpectra=" << centDefSpectra << endl
164  << "\tcentDefOther=" << centDefOther << endl
165  << "\tvertex=" << vertex << endl
166  << "\thalf=" << half << endl
167  << "\tfitPts=" << fitPts << endl
168  << "\tdca=" << dca << endl
169  << "\teta=" << eta << endl;
170 
171 
172  // list the peripheral definitions here for convenience
173  NchCentrality peripheralMin = kEighty,
174  peripheralMax = kSeventy;
175  NchCentrality kludgeMin = kEighty,
176  kludgeMax = kSixty;
177  //NchCentrality kludge2Min = kEighty,
178  // kludge2Max = kSixty;
179 
180 
181  switch(centType){
182  case 9: // minbias
183  if(mDoSpectraCent){
184  if(centDefSpectra==9){ // flow
185  mFlowCent[0]=0; mFlowCent[1]=9; mDoFlowCent=true;
186  }
187  else if(centDefSpectra==8){ // zdc-ctb
188  mZdcCtbCent[0]=kTotal; mZdcCtbCent[1]=kFive; mDoZdcCtbCent=true;
189  }
190  else if(centDefSpectra==7){ // hminus
191  mHMinusCent[0]=kTotal; mHMinusCent[1]=kFive; mDoHMinusCent=true;
192  }
193  else if(centDefSpectra==6){ // nch
194  mNchCent[0]=kTotal; mNchCent[1]=kFive; mDoNchCent=true;
195  }
196  else if(centDefSpectra==5){ // nch
197  mNchCent[0]=kTotal; mNchCent[1]=kFive; mDoNchCentKludge=true;
198  }
199  else{
200  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
201  }
202  }
203  else{ // not spectra centrality
204  if(centDefOther==9){ // flow
205  mFlowCent[0]=0; mFlowCent[1]=9; mDoFlowCent=true;
206  }
207  else if(centDefOther==8){ // zdc-ctb
208  mZdcCtbCent[0]=kTotal; mZdcCtbCent[1]=kFive; mDoZdcCtbCent=true;
209  }
210  else if(centDefOther==7){ // hminus
211  mHMinusCent[0]=kTotal; mHMinusCent[1]=kFive; mDoHMinusCent=true;
212  }
213  else if(centDefOther==6){ // nch
214  mNchCent[0]=kTotal; mNchCent[1]=kFive; mDoNchCent=true;
215  }
216  else if(centDefOther==5){ // nch
217  mNchCent[0]=kTotal; mNchCent[1]=kFive; mDoNchCentKludge=true;
218  }
219  else{
220  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
221  }
222  }
223  break;
224  case 8: // central
225  if(mDoSpectraCent){
226  if(centDefSpectra==9){ // flow
227  mFlowCent[0]=8; mFlowCent[1]=9; mDoFlowCent=true;
228  }
229  else if(centDefSpectra==8){ // zdc-ctb
230  mZdcCtbCent[0]=kFive; mZdcCtbCent[1]=kFive; mDoZdcCtbCent=true;
231  }
232  else if(centDefSpectra==7){ // hminus
233  mHMinusCent[0]=kFive; mHMinusCent[1]=kFive; mDoHMinusCent=true;
234  }
235  else if(centDefSpectra==6){ // nch
236  mNchCent[0]=kFive; mNchCent[1]=kFive; mDoNchCent=true;
237  }
238 
239  else{
240  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
241  }
242  }
243  else{
244  if(centDefOther==9){ // flow
245  mFlowCent[0]=8; mFlowCent[1]=9; mDoFlowCent=true;
246  }
247  else if(centDefOther==8){ // zdc-ctb
248  mZdcCtbCent[0]=kTen; mZdcCtbCent[1]=kFive; mDoZdcCtbCent=true;
249  }
250  else if(centDefOther==7){ // hminus
251  mHMinusCent[0]=kTen; mHMinusCent[1]=kFive; mDoHMinusCent=true;
252  }
253  else if(centDefOther==6){ // nch
254  mNchCent[0]=kTen; mNchCent[1]=kFive; mDoNchCent=true;
255  }
256  else{
257  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
258  }
259  }
260  break;
261  case 7: // peripheral
262  if(mDoSpectraCent){
263  if(centDefSpectra==9){ // flow
264  mFlowCent[0]=1; mFlowCent[1]=3; mDoFlowCent=true;
265  }
266  else if(centDefSpectra==8){ // zdc-ctb
267  mZdcCtbCent[0]=peripheralMin;
268  mZdcCtbCent[1]=peripheralMax; mDoZdcCtbCent=true;
269  }
270  else if(centDefSpectra==7){ // hminus
271  mHMinusCent[0]=peripheralMin;
272  mHMinusCent[1]=peripheralMax; mDoHMinusCent=true;
273  }
274  else if(centDefSpectra==6){ // nch
275  mNchCent[0]=peripheralMin;
276  mNchCent[1]=peripheralMax; mDoNchCent=true;
277  }
278  else if(centDefSpectra==5){ // nch
279  mNchCent[0]=peripheralMin;
280  mNchCent[1]=peripheralMax; mDoNchCentKludge=true;
281  }
282  else{
283  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
284  }
285  }
286  else{
287  if(centDefOther==9){ // flow
288  mFlowCent[0]=1; mFlowCent[1]=3; mDoFlowCent=true;
289  }
290  else if(centDefOther==8){ // zdc-ctb
291  mZdcCtbCent[0]=peripheralMin;
292  mZdcCtbCent[1]=peripheralMax; mDoZdcCtbCent=true;
293  }
294  else if(centDefOther==7){ // hminus
295  mHMinusCent[0]=peripheralMin;
296  mHMinusCent[1]=peripheralMax; mDoHMinusCent=true;
297  }
298  else if(centDefOther==6){ // nch
299  mNchCent[0]=peripheralMin;
300  mNchCent[1]=peripheralMax; mDoNchCent=true;
301  }
302  else if(centDefOther==5){ // nch
303  mNchCent[0]=peripheralMin;
304  mNchCent[1]=peripheralMax; mDoNchCentKludge=true;
305  }
306  else{
307  cerr << "Unknown centrality definition for spectra: " << centDefSpectra<<endl; exit(1);
308  }
309  }
310  break;
311  case 6: // central
312  if(mDoSpectraCent){
313  if(centDefSpectra==9){ // flow
314  mFlowCent[0]=8; mFlowCent[1]=9; mDoFlowCent=true;
315  }
316  else if(centDefSpectra==8){ // zdc-ctb
317  mZdcCtbCent[0]=kFive; mZdcCtbCent[1]=kFive; mDoZdcCtbCent=true;
318  }
319  else if(centDefSpectra==7){ // hminus
320  mHMinusCent[0]=kFive; mHMinusCent[1]=kFive; mDoHMinusCent=true;
321  }
322  else if(centDefSpectra==6){ // nch
323  mNchCent[0]=kFive; mNchCent[1]=kFive; mDoNchCent=true;
324  }
325  else if(centDefSpectra==5){ // nch
326  mNchCent[0]=kFive; mNchCent[1]=kFive; mDoNchCentKludge=true;
327  }
328  else{
329  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
330  }
331  }
332  else{
333  if(centDefOther==9){ // flow
334  mFlowCent[0]=8; mFlowCent[1]=9; mDoFlowCent=true;
335  }
336  else if(centDefOther==8){ // zdc-ctb
337  mZdcCtbCent[0]=kFive; mZdcCtbCent[1]=kFive; mDoZdcCtbCent=true;
338  }
339  else if(centDefOther==7){ // hminus
340  mHMinusCent[0]=kFive; mHMinusCent[1]=kFive; mDoHMinusCent=true;
341  }
342  else if(centDefOther==6){ // nch
343  mNchCent[0]=kFive; mNchCent[1]=kFive; mDoNchCent=true;
344  }
345  else if(centDefOther==5){ // nch
346  mNchCent[0]=kFive; mNchCent[1]=kFive; mDoNchCentKludge=true;
347  }
348  else{
349  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
350  }
351  }
352  break;
353 
354  case 5: // peripheral kludge
355  if(mDoSpectraCent){
356  if(centDefSpectra==9){ // flow
357  mFlowCent[0]=1; mFlowCent[1]=3; mDoFlowCent=true;
358  }
359  else if(centDefSpectra==8){ // zdc-ctb
360  mZdcCtbCent[0]=peripheralMin;
361  mZdcCtbCent[1]=peripheralMax; mDoZdcCtbCent=true;
362  }
363  else if(centDefSpectra==7){ // hminus
364  mHMinusCent[0]=peripheralMin;
365  mHMinusCent[1]=peripheralMax; mDoHMinusCent=true;
366  }
367  else if(centDefSpectra==6){ // nch
368  mNchCent[0]=peripheralMin;
369  mNchCent[1]=peripheralMax; mDoNchCent=true;
370  }
371  else if(centDefSpectra==5){ // nch
372  mNchCent[0]=peripheralMin;
373  mNchCent[1]=peripheralMax; mDoNchCentKludge=true;
374  }
375  else{
376  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
377  }
378  }
379  else{
380  if(centDefOther==9){ // flow
381  mFlowCent[0]=1; mFlowCent[1]=3; mDoFlowCent=true;
382  }
383  else if(centDefOther==8){ // zdc-ctb
384  mZdcCtbCent[0]=kludgeMin;
385  mZdcCtbCent[1]=kludgeMax; mDoZdcCtbCent=true;
386  }
387  else if(centDefOther==7){ // hminus
388  mHMinusCent[0]=kludgeMin;
389  mHMinusCent[1]=kludgeMax; mDoHMinusCent=true;
390  }
391  else if(centDefOther==6){ // nch
392  mNchCent[0]=kludgeMin;
393  mNchCent[1]=kludgeMax; mDoNchCent=true;
394  }
395  else if(centDefOther==5){ // nch
396  mNchCent[0]=kludgeMin;
397  mNchCent[1]=kludgeMax; mDoNchCentKludge=true;
398  }
399  else{
400  cerr << "Unknown centrality definition for spectra: " << centDefSpectra<<endl; exit(1);
401  }
402  }
403  break;
404 
405  case 4: // 10-20%
406  if(mDoSpectraCent){
407  if(centDefSpectra==9){ // flow
408  mFlowCent[0]=6; mFlowCent[1]=7; mDoFlowCent=true;
409  }
410  else if(centDefSpectra==8){ // zdc-ctb
411  mZdcCtbCent[0]=kTwenty;
412  mZdcCtbCent[1]=kTwenty; mDoZdcCtbCent=true;
413  }
414  else if(centDefSpectra==7){ // hminus
415  mHMinusCent[0]=kTwenty;
416  mHMinusCent[1]=kTwenty; mDoHMinusCent=true;
417  }
418  else if(centDefSpectra==6){ // nch
419  mNchCent[0]=kTwenty;
420  mNchCent[1]=kTwenty; mDoNchCent=true;
421  }
422  else if(centDefSpectra==5){ // nch
423  mNchCent[0]=kTwenty;
424  mNchCent[1]=kTwenty; mDoNchCentKludge=true;
425  }
426  else{
427  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
428  }
429  }
430  else{
431  if(centDefOther==9){ // flow
432  mFlowCent[0]=6; mFlowCent[1]=7; mDoFlowCent=true;
433  }
434  else if(centDefOther==8){ // zdc-ctb
435  mZdcCtbCent[0]=kTwenty;
436  mZdcCtbCent[1]=kTwenty; mDoZdcCtbCent=true;
437  }
438  else if(centDefOther==7){ // hminus
439  mHMinusCent[0]=kTwenty;
440  mHMinusCent[1]=kTwenty; mDoHMinusCent=true;
441  }
442  else if(centDefOther==6){ // nch
443  mNchCent[0]=kTwenty;
444  mNchCent[1]=kTwenty; mDoNchCent=true;
445  }
446  else if(centDefOther==5){ // nch
447  mNchCent[0]=kTwenty;
448  mNchCent[1]=kTwenty; mDoNchCentKludge=true;
449  }
450  else{
451  cerr << "Unknown centrality definition for spectra: " << centDefSpectra<<endl; exit(1);
452  }
453  }
454  break;
455  case 3: // 20-30
456  if(mDoSpectraCent){
457  if(centDefSpectra==9){ // flow
458  mFlowCent[0]=4; mFlowCent[1]=5; mDoFlowCent=true;
459  }
460  else if(centDefSpectra==8){ // zdc-ctb
461  mZdcCtbCent[0]=kThirty;
462  mZdcCtbCent[1]=kThirty; mDoZdcCtbCent=true;
463  }
464  else if(centDefSpectra==7){ // hminus
465  mHMinusCent[0]=kThirty;
466  mHMinusCent[1]=kThirty; mDoHMinusCent=true;
467  }
468  else if(centDefSpectra==6){ // nch
469  mNchCent[0]=kThirty;
470  mNchCent[1]=kThirty; mDoNchCent=true;
471  }
472  else if(centDefSpectra==5){ // nch
473  mNchCent[0]=kThirty;
474  mNchCent[1]=kThirty; mDoNchCentKludge=true;
475  }
476  else{
477  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
478  }
479  }
480  else{
481  if(centDefOther==9){ // flow
482  mFlowCent[0]=1; mFlowCent[1]=3; mDoFlowCent=true;
483  }
484  else if(centDefOther==8){ // zdc-ctb
485  mZdcCtbCent[0]=kThirty;
486  mZdcCtbCent[1]=kThirty; mDoZdcCtbCent=true;
487  }
488  else if(centDefOther==7){ // hminus
489  mHMinusCent[0]=kThirty;
490  mHMinusCent[1]=kThirty; mDoHMinusCent=true;
491  }
492  else if(centDefOther==6){ // nch
493  mNchCent[0]=kThirty;
494  mNchCent[1]=kThirty; mDoNchCent=true;
495  }
496  else if(centDefOther==5){ // nch
497  mNchCent[0]=kThirty;
498  mNchCent[1]=kThirty; mDoNchCentKludge=true;
499  }
500  else{
501  cerr << "Unknown centrality definition for spectra: " << centDefSpectra<<endl; exit(1);
502  }
503  }
504  break;
505 
506  case 2: // 30-40
507  if(mDoSpectraCent){
508  if(centDefSpectra==9){ // flow
509  mFlowCent[0]=3; mFlowCent[1]=3; mDoFlowCent=true;
510  }
511  else if(centDefSpectra==8){ // zdc-ctb
512  mZdcCtbCent[0]=kForty;
513  mZdcCtbCent[1]=kForty; mDoZdcCtbCent=true;
514  }
515  else if(centDefSpectra==7){ // hminus
516  mHMinusCent[0]=kForty;
517  mHMinusCent[1]=kForty; mDoHMinusCent=true;
518  }
519  else if(centDefSpectra==6){ // nch
520  mNchCent[0]=kForty;
521  mNchCent[1]=kForty; mDoNchCent=true;
522  }
523  else if(centDefSpectra==5){ // nch
524  mNchCent[0]=kForty;
525  mNchCent[1]=kForty; mDoNchCentKludge=true;
526  }
527  else{
528  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
529  }
530  }
531  else{
532  if(centDefOther==9){ // flow
533  mFlowCent[0]=1; mFlowCent[1]=3; mDoFlowCent=true;
534  }
535  else if(centDefOther==8){ // zdc-ctb
536  mZdcCtbCent[0]=kForty;
537  mZdcCtbCent[1]=kForty; mDoZdcCtbCent=true;
538  }
539  else if(centDefOther==7){ // hminus
540  mHMinusCent[0]=kForty;
541  mHMinusCent[1]=kForty; mDoHMinusCent=true;
542  }
543  else if(centDefOther==6){ // nch
544  mNchCent[0]=kForty;
545  mNchCent[1]=kForty; mDoNchCent=true;
546  }
547  else if(centDefOther==5){ // nch
548  mNchCent[0]=kForty;
549  mNchCent[1]=kForty; mDoNchCentKludge=true;
550  }
551  else{
552  cerr << "Unknown centrality definition for spectra: " << centDefSpectra<<endl; exit(1);
553  }
554  }
555  break;
556 
557  case 1: // 40-60
558  if(mDoSpectraCent){
559  if(centDefSpectra==9){ // flow
560  mFlowCent[0]=2; mFlowCent[1]=2; mDoFlowCent=true;
561  }
562  else if(centDefSpectra==8){ // zdc-ctb
563  mZdcCtbCent[0]=kSixty;
564  mZdcCtbCent[1]=kFifty; mDoZdcCtbCent=true;
565  }
566  else if(centDefSpectra==7){ // hminus
567  mHMinusCent[0]=kSixty;
568  mHMinusCent[1]=kFifty; mDoHMinusCent=true;
569  }
570  else if(centDefSpectra==6){ // nch
571  mNchCent[0]=kSixty;
572  mNchCent[1]=kFifty; mDoNchCent=true;
573  }
574  else if(centDefSpectra==5){ // nch
575  mNchCent[0]=kSixty;
576  mNchCent[1]=kFifty; mDoNchCentKludge=true;
577  }
578  else{
579  cerr << "unknown centrality definition for spectra " << centDefSpectra<<endl; exit(1);
580  }
581  }
582  else{
583  if(centDefOther==9){ // flow
584  mFlowCent[0]=1; mFlowCent[1]=3; mDoFlowCent=true;
585  }
586  else if(centDefOther==8){ // zdc-ctb
587  mZdcCtbCent[0]=kSixty;
588  mZdcCtbCent[1]=kFifty; mDoZdcCtbCent=true;
589  }
590  else if(centDefOther==7){ // hminus
591  mHMinusCent[0]=kSixty;
592  mHMinusCent[1]=kFifty; mDoHMinusCent=true;
593  }
594  else if(centDefOther==6){ // nch
595  mNchCent[0]=kSixty;
596  mNchCent[1]=kFifty; mDoNchCent=true;
597  }
598  else if(centDefOther==5){ // nch
599  mNchCent[0]=kSixty;
600  mNchCent[1]=kFifty; mDoNchCentKludge=true;
601  }
602  else{
603  cerr << "Unknown centrality definition for spectra: " << centDefSpectra<<endl; exit(1);
604  }
605  }
606  break;
607 
608  default:
609  cerr << "Unknown centrality type: " << centType << endl; exit(1);
610  }
611 
612  switch(vertex){
613  case 9:
614  mVertexZ[0]=-200; mVertexZ[1]=200; break;
615  case 8:
616  mVertexZ[0]=-75; mVertexZ[1]=75; break;
617  case 7:
618  mVertexZ[0]=-95; mVertexZ[1]=95; break;
619  case 6:
620  mVertexZ[0]=-30; mVertexZ[1]=30; break;
621  case 5:
622  mVertexZ[0]=0; mVertexZ[1]=30; break;
623  default:
624  cerr << "Unknown vertex type: " << vertex << endl; exit(1);
625  }
626 
627  switch(half){
628  case 9:
629  break;
630  case 8:
631  SetHitHalf('s'); break;
632  case 7:
633  SetHitAvoid(10); break;
634  default:
635  cerr << "Unknown half type : " << half << endl; exit(1);
636  }
637 
638  switch(fitPts){
639  case 9:
640  mFitPts[0]=20; mFitPts[1]=99; break;
641  case 8:
642  mFitPts[0]=25; mFitPts[1]=99; break;
643  case 7:
644  mFitPts[0]=10; mFitPts[1]=99; break;
645  case 6:
646  mFitPts[0]=35; mFitPts[1]=99; break;
647  case 5:
648  mFitPts[0]=24; mFitPts[1]=99; break;
649  case 4:
650  mFitPts[0]=23; mFitPts[1]=99; break;
651  case 3:
652  mFitPts[0]=30; mFitPts[1]=99; break;
653  default:
654  cerr << "Unknown fit pts type : " << fitPts << endl; exit(1);
655  }
656 
657  switch(dca){
658  case 9:
659  mSDcaGl[0]=-1; mSDcaGl[1]=1; break;
660  case 8:
661  mSDcaGl[0]=-1.2; mSDcaGl[1]=1.2; break;
662  case 7:
663  mSDcaGl[0]=-3; mSDcaGl[1]=3; break;
664 
665  default:
666  cerr << "Unknown dca type : " << half << endl; exit(1);
667  }
668 
669  switch(eta){
670  case 9:
671  mEta[0]=-.5; mEta[1]=.5; break;
672  case 8:
673  mEta[0]=-.7; mEta[1]=.7; break;
674  case 7:
675  mEta[0]=0; mEta[1]=0.5; break;
676  default:
677  cerr << "Unknown dca type : " << half << endl; exit(1);
678  }
679 
680 
681 }
682 
683 
684 void
685 Cut::ShowCuts()
686 {
687 
688  cout << "******************************************************" << endl;
689  cout << "Cut::ShowCuts()" << endl;
690  cout << "Cut=" << mCut << endl;
691  cout << "event cuts:" << endl;
692 
693  if(mDoSpectraCent)
694  cout << "\tWill use spectra centrality definition" << endl;
695  else
696  cout << "\tWill use 'other' centrality definition" << endl;
697 
698  if(mDoFlowCent)
699  cout << "\tflow cent : " << mFlowCent[0] << " -- " << mFlowCent[1] << endl;
700  else if(mDoZdcCtbCent)
701  cout << "\tzdc-ctb cent : " << (int) mZdcCtbCent[0] << " -- "
702  << (int) mZdcCtbCent[1] << endl;
703  else if(mDoHMinusCent)
704  cout << "\th minus cent : " << (int) mHMinusCent[0] << " -- "
705  << (int) mHMinusCent[1] << endl;
706  else if(mDoNchCent)
707  cout << "\th nCh cent : " << (int) mNchCent[0] << " -- "
708  << (int) mNchCent[1] << endl;
709  else if(mDoNchCentKludge){
710  cout << "\tDoing nCh KLUDGE" << endl;
711  cout << "\th nCh cent : " << (int) mNchCent[0] << " -- "
712  << (int) mNchCent[1] << endl;
713  }
714  else{
715  cout << "\tUnknown centrality definition. Goodbye" << endl; exit(1);
716  }
717  if(mDoNchCentKludge)
718  cout << "\tDoing nCh KLUDGE" << endl;
719 
720  cout << "\tvertex z : " << mVertexZ[0] << " -- " << mVertexZ[1] << endl;
721  if(mVertexZSkip){
722  cout << "\tvtx z skip: "<< mVertexZSkip << endl;
723  }
724  if(mHalf){
725  cout << "\twill process half " << "'" <<mHalf <<"'" << endl;
726  }
727  if(mGeomHalf){
728  cout << "\twill process geom half " << "'" << mGeomHalf << "'" << endl;
729  }
730 
731  cout << "track cuts:" << endl
732  << "\teta : " << mEta[0] << " -- " << mEta[1] << endl
733  << "\tfit pts : " << mFitPts[0] << " -- " << mFitPts[1] << endl
734  << "\tsdca : " << mSDcaGl[0] << " -- " << mSDcaGl[1] << endl
735  << "\tmc pts : " << mMcPts[0] << " -- " << mMcPts[1] << endl
736  << "\teta tight : " << mEtaTight[0] << "-- " << mEtaTight[1] << endl
737  << endl;
738  if(mHitAvoid){
739  cout << "\ttracks with hits outside of " << mHitAvoid << " cm" <<endl;
740  }
741  if(mHitHalf){
742  cout << "\twill process tracks with hits on half " << mHitHalf << endl;
743  }
744  cout << "******************************************************" << endl;
745 
746 
747 
748 }
Cut
Definition: Cut.h:18
Generated by   doxygen 1.8.5