changes from gsi. Using mult if no centrality. testfilterbit 128
[u/mrichter/AliRoot.git] / PWGGA / GammaConv / macros / printCuts.C
CommitLineData
baf06eb6 1
2
3using namespace std;
4void printCuts(TString str){
5
6
7
8 TString s = str(0,1);
9 Int_t goodId= s.Atoi();
10 s = str(1,1);
11 Int_t v0FinderType= s.Atoi();
12 s = str(2,1);
13 Int_t eProbCut= s.Atoi();
14 s = str(3,1);
15 Int_t ededxSigmaCut= s.Atoi();
16 s = str(4,1);
17 Int_t pidedxSigmaCut= s.Atoi();
18 s = str(5,1);
19 Int_t piMomdedxSigmaCut= s.Atoi();
20 s = str(6,1);
21 Int_t chi2GammaCut= s.Atoi();
22 s = str(7,1);
23 Int_t singlePtCut= s.Atoi();
24 s = str(8,1);
25 Int_t clsTPCCut= s.Atoi();
26 s = str(9,1);
27 Int_t etaCut= s.Atoi();
28
29 /* cout<<"etaCut: "<<etaCut<<endl;
30 cout<<"clsTPCCut: "<<clsTPCCut<<endl;
31 cout<<"singlePtCut: "<<singlePtCut<<endl;
32 cout<<"chi2GammaCut: "<<chi2GammaCut<<endl;
33 cout<<"piMomdedxSigmaCut: "<<piMomdedxSigmaCut<<endl;
34 cout<<"pidedxSigmaCut: "<<pidedxSigmaCut <<endl;
35 cout<<"ededxSigmaCut: "<<ededxSigmaCut <<endl;
36 cout<<"eProbCut: "<< eProbCut<<endl;
37 cout<<"v0FinderType: "<<v0FinderType <<endl;
38 cout<<"goodId: "<<goodId <<endl;
39 */
40
41 if(goodId !=9){
42 cout<<"Analysis Cut Selection too short or does not start with 9"<<endl;
43 return;
44 }
45
46 switch (v0FinderType){
47 case 0: // on fly V0 finder
48 cout<<"Using the on fly v0 finder"<<endl;
49 break;
50 case 1: // offline V0 finder
51 cout<<"Using the offline v0 finder"<<endl;
52 break;
53 default:
54 return;
55 }
56 switch(eProbCut){
57 case 0: // 0.
58 cout<<"Prob electron is 0.000"<<endl;
59 // kGCprobElectron = 0.000;
60 break;
61 case 1: // 0.001
62 cout<<"Prob electron is 0.001"<<endl;
63 // kGCprobElectron = 0.001;
64 break;
65 case 2: // 0.01
66 cout<<"Prob electron is 0.01"<<endl;
67 // kGCprobElectron = 0.01;
68 break;
69 default:
70 return;
71 }
72
73 switch(ededxSigmaCut){
74 case 0: // -10,10
75 cout<<"Using sigma dedx [-10,10]"<<endl;
76 // kGCPIDnSigmaBelowElectronLine=-10;
77 // kGCPIDnSigmaAboveElectronLine=10;
78 break;
79 case 1: // -5,5
80 cout<<"Using sigma dedx [-5,5]"<<endl;
81 // kGCPIDnSigmaBelowElectronLine=-5;
82 // kGCPIDnSigmaAboveElectronLine=5;
83 break;
84 case 2: // -3,5
85 cout<<"Using sigma dedx [-3,5]"<<endl;
86 // kGCPIDnSigmaBelowElectronLine=-3;
87 // kGCPIDnSigmaAboveElectronLine=5;
88 break;
89 default:
90 return;
91 }
92
93 switch(pidedxSigmaCut){
94 case 0: // -10
95 cout<<"using pidedxsigmacut: -10"<<endl;
96 // kGCPIDnSigmaAbovePionLine=-10;
97 break;
98 case 1: // 0
99 cout<<"using pidedxsigmacut: 0"<<endl;
100 // kGCPIDnSigmaAbovePionLine=0;
101 break;
102 case 2: // 1
103 cout<<"using pidedxsigmacut: 1"<<endl;
104 // kGCPIDnSigmaAbovePionLine=1;
105 break;
106 default:
107 return;
108 }
109
110 switch(piMomdedxSigmaCut){
111 case 0: // 0.5 GeV
112 cout<<"piMomdedxSigmaCut: 0.5"<<endl;
113 // kGCPIDMinPnSigmaAbovePionLine=0.5;
114 break;
115 case 1: // 1. GeV
116 cout<<"piMomdedxSigmaCut: 1"<<endl;
117 // kGCPIDMinPnSigmaAbovePionLine=1.;
118 break;
119 case 2: // 1.5 GeV
120 cout<<"piMomdedxSigmaCut: 1.5"<<endl;
121 // kGCPIDMinPnSigmaAbovePionLine=1.5;
122 break;
123 default:
124 return;
125 }
126
127 switch(chi2GammaCut){
128 case 0: // 100
129 cout<<"chi2CutConversion = 100."<<endl;
130 // kGCchi2CutConversion = 100.;
131 break;
132 case 1: // 50
133 cout<<"chi2CutConversion = 50."<<endl;
134 // kGCchi2CutConversion = 50.;
135 break;
136 case 2: // 30
137 cout<<"chi2CutConversion = 30."<<endl;
138 // kGCchi2CutConversion = 30.;
139 break;
140 default:
141 return;
142 }
143
144 switch(singlePtCut){
145 case 0: // 0.050 GeV
146 cout<<"kGCsingleptCut = 0,050"<<endl;
147 // kGCsingleptCut = 0.050;
148 break;
149 case 1: // 0.100 GeV
150 cout<<"kGCsingleptCut = 0,100"<<endl;
151 // kGCsingleptCut = 0.100;
152 break;
153 case 2: // 0.150 GeV
154 cout<<"kGCsingleptCut = 0,150"<<endl;
155 // kGCsingleptCut = 0.150;
156 break;
157 case 3: // 0.200 GeV
158 cout<<"kGCsingleptCut = 0,200"<<endl;
159 // kGCsingleptCut = 0.200;
160 break;
161 default:
162 return;
163 }
164
165 switch(clsTPCCut){
166 case 0: // 0
167 cout<<"kGCminClsTPCCut = 0"<<endl;
168 // kGCminClsTPCCut= 0.;
169 break;
170 case 1: // 70
171 // kGCminClsTPCCut= 70.;
172 cout<<"kGCminClsTPCCut = 70"<<endl;
173 break;
174 case 2: // 80
175 // kGCminClsTPCCut= 80.;
176 cout<<"kGCminClsTPCCut = 80"<<endl;
177 break;
178 case 3: // 100
179 // kGCminClsTPCCut= 100.;
180 cout<<"kGCminClsTPCCut = 100"<<endl;
181 break;
182 default:
183 return;
184 }
185
186 switch(etaCut){
187 case 0: // 0.9
188 cout<<"eta 0.9"<<endl;
189 // kGCetaCut = 0.9;
190 // kGCLineCutZRSlope = tan(2*atan(exp(-kGCetaCut)));
191 break;
192 case 1: // 1.2
193 cout<<"eta 1.2"<<endl;
194 // kGCetaCut = 1.2;
195 // kGCLineCutZRSlope = tan(2*atan(exp(-kGCetaCut)));
196 break;
197 case 2: // 1.4
198 cout<<"eta 1.4"<<endl;
199 // kGCetaCut = 1.4;
200 // kGCLineCutZRSlope = tan(2*atan(exp(-kGCetaCut)));
201 break;
202 default:
203 return;
204 }
205
206
207}