Removing warnings (Panos)
[u/mrichter/AliRoot.git] / STEER / AliEventTagCuts.cxx
1 /**************************************************************************
2  * Author: Panos Christakoglou.                                           *
3  * Contributors are mentioned in the code where appropriate.              *
4  *                                                                        *
5  * Permission to use, copy, modify and distribute this software and its   *
6  * documentation strictly for non-commercial purposes is hereby granted   *
7  * without fee, provided that the above copyright notice appears in all   *
8  * copies and that both the copyright notice and this permission notice   *
9  * appear in the supporting documentation. The authors make no claims     *
10  * about the suitability of this software for any purpose. It is          *
11  * provided "as is" without express or implied warranty.                  *
12  **************************************************************************/
13
14 /* $Id$ */
15
16 //-----------------------------------------------------------------
17 //           AliEventTagCuts class
18 //   This is the class to deal with the event tag level cuts
19 //   Origin: Panos Christakoglou, UOA-CERN, Panos.Christakoglou@cern.ch
20 //-----------------------------------------------------------------
21
22 class AliLog;
23 class AliESD;
24
25 #include "AliEventTag.h"
26 #include "AliEventTagCuts.h"
27
28 ClassImp(AliEventTagCuts)
29
30
31 //___________________________________________________________________________
32 AliEventTagCuts::AliEventTagCuts() :
33   TObject(),
34   
35   fNParticipantsMin(-1), fNParticipantsMax(10000),
36   fNParticipantsFlag(kFALSE),
37   fImpactParamMin(-1.0), fImpactParamMax(1000.0),
38   fImpactParamFlag(kFALSE),
39
40   fVxMin(-1000.0), fVxMax(1000.0), 
41   fVxFlag(kFALSE),
42   fVyMin(-1000.0), fVyMax(1000.0),  
43   fVyFlag(kFALSE),
44   fVzMin(-1000.0), fVzMax(1000.0),
45   fVzFlag(kFALSE),
46   fPrimaryVertexFlag(1),
47   fPVFlag(kFALSE),
48
49   fPrimaryVertexZErrorMin(-10000.), fPrimaryVertexZErrorMax(10000.),
50   fPVzErrorFlag(kFALSE),
51   fTriggerMask(0),
52   fTriggerMaskFlag(kFALSE),
53   fTriggerCluster(0),
54   fTriggerClusterFlag(kFALSE),
55  
56   fZDCNeutron1EnergyMin(-1.0), fZDCNeutron1EnergyMax(100000.0),
57   fZDCNeutron1EnergyFlag(kFALSE),
58   fZDCProton1EnergyMin(-1.0), fZDCProton1EnergyMax(100000.0),
59   fZDCProton1EnergyFlag(kFALSE),
60   fZDCNeutron2EnergyMin(-1.0), fZDCNeutron2EnergyMax(100000.0),
61   fZDCNeutron2EnergyFlag(kFALSE),
62   fZDCProton2EnergyMin(-1.0), fZDCProton2EnergyMax(100000.0),
63   fZDCProton2EnergyFlag(kFALSE),
64   fZDCEMEnergyMin(-1.0), fZDCEMEnergyMax(100000.0),
65   fZDCEMEnergyFlag(kFALSE),
66   fT0VertexZMin(-10000.0), fT0VertexZMax(10000.0),  
67   fT0VertexZFlag(kFALSE),
68   fMultMin(0), fMultMax(100000),  
69   fMultFlag(kFALSE),
70   fPosMultMin(-1), fPosMultMax(100000),
71   fPosMultFlag(kFALSE),
72   fNegMultMin(-1), fNegMultMax(100000),
73   fNegMultFlag(kFALSE),
74   fNeutrMultMin(-1), fNeutrMultMax(100000),
75   fNeutrMultFlag(kFALSE),
76   fNV0sMin(-1), fNV0sMax(1000000),
77   fNV0sFlag(kFALSE),
78   fNCascadesMin(-1), fNCascadesMax(100000),
79   fNCascadesFlag(kFALSE),
80   fNKinksMin(-1), fNKinksMax(1000000),
81   fNKinksFlag(kFALSE),
82
83   fNPMDTracksMin(-1), fNPMDTracksMax(100000),
84   fNPMDTracksFlag(kFALSE),
85   fNFMDTracksMin(-1), fNFMDTracksMax(100000),
86   fNFMDTracksFlag(kFALSE),
87   fNPHOSClustersMin(-1), fNPHOSClustersMax(100000),
88   fNPHOSClustersFlag(kFALSE),
89   fNEMCALClustersMin(-1), fNEMCALClustersMax(100000),
90   fNEMCALClustersFlag(kFALSE),
91   fNJetCandidatesMin(-1), fNJetCandidatesMax(100000),
92   fNJetCandidatesFlag(kFALSE),
93
94   fTopJetEnergyMin(-1.0), 
95   fTopJetEnergyMinFlag(kFALSE),
96   fTopNeutralEnergyMin(-1.0), 
97   fTopNeutralEnergyMinFlag(kFALSE),
98   fNHardPhotonCandidatesMin(-1), fNHardPhotonCandidatesMax(100000),
99   fNHardPhotonCandidatesFlag(kFALSE),
100   fNChargedAbove1GeVMin(-1), fNChargedAbove1GeVMax(100000),
101   fNChargedAbove1GeVFlag(kFALSE),
102   fNChargedAbove3GeVMin(-1), fNChargedAbove3GeVMax(100000),
103   fNChargedAbove3GeVFlag(kFALSE),
104   fNChargedAbove10GeVMin(-1), fNChargedAbove10GeVMax(100000),
105   fNChargedAbove10GeVFlag(kFALSE),
106   fNMuonsAbove1GeVMin(-1), fNMuonsAbove1GeVMax(100000),
107   fNMuonsAbove1GeVFlag(kFALSE),
108   fNMuonsAbove3GeVMin(-1), fNMuonsAbove3GeVMax(100000),
109   fNMuonsAbove3GeVFlag(kFALSE),
110   fNMuonsAbove10GeVMin(-1), fNMuonsAbove10GeVMax(100000), 
111   fNMuonsAbove10GeVFlag(kFALSE),
112   fNElectronsAbove1GeVMin(-1), fNElectronsAbove1GeVMax(100000),
113   fNElectronsAbove1GeVFlag(kFALSE),
114   fNElectronsAbove3GeVMin(-1), fNElectronsAbove3GeVMax(100000),
115   fNElectronsAbove3GeVFlag(kFALSE),
116   fNElectronsAbove10GeVMin(-1), fNElectronsAbove10GeVMax(100000),
117   fNElectronsAbove10GeVFlag(kFALSE),
118   fNElectronsMin(-1), fNElectronsMax(100000),
119   fNElectronsFlag(kFALSE),
120   fNMuonsMin(-1), fNMuonsMax(100000),
121   fNMuonsFlag(kFALSE),
122   fNPionsMin(-1), fNPionsMax(100000),
123   fNPionsFlag(kFALSE),
124   fNKaonsMin(-1), fNKaonsMax(100000),
125   fNKaonsFlag(kFALSE),
126   fNProtonsMin(-1), fNProtonsMax(100000),
127   fNProtonsFlag(kFALSE),
128   fNLambdasMin(-1), fNLambdasMax(100000),
129   fNLambdasFlag(kFALSE),
130   fNPhotonsMin(-1), fNPhotonsMax(100000),
131   fNPhotonFlag(kFALSE),
132   fNPi0sMin(-1), fNPi0sMax(100000), 
133   fNPi0sFlag(kFALSE),
134   fNNeutronsMin(-1), fNNeutronsMax(100000), 
135   fNNeutronsFlag(kFALSE),
136   fNKaon0sMin(-1), fNKaon0sMax(100000), 
137   fNKaon0sFlag(kFALSE),
138   fTotalPMin(-1.0), fTotalPMax(1000000.0),
139   fTotalPFlag(kFALSE),
140   fMeanPtMin(-1.0), fMeanPtMax(100000.0),
141   fMeanPtFlag(kFALSE),
142   fTopPtMin(-1.0),
143   fTopPtMinFlag(kFALSE),
144   fTotalNeutralPMin(-1.0), fTotalNeutralPMax(1000000.0),   
145   fTotalNeutralPFlag(kFALSE),
146   fMeanNeutralPtMin(-1.0), fMeanNeutralPtMax(1000000.0), 
147   fMeanNeutralPtFlag(kFALSE),
148   fTopNeutralPtMin(-1.0), 
149   fTopNeutralPtMinFlag(kFALSE),
150   fEventPlaneAngleMin(-10000000.0), fEventPlaneAngleMax(10000000.0), 
151   fEventPlaneAngleFlag(kFALSE),
152   fHBTRadiiMin(-1.0), fHBTRadiiMax(100000.0), 
153   fHBTRadiiFlag(kFALSE)
154 {
155   //Default constructor which calls the Reset method.
156   Reset();
157 }
158
159 //___________________________________________________________________________
160 AliEventTagCuts::~AliEventTagCuts() {  
161   //Defaut destructor.
162 }
163
164 //___________________________________________________________________________
165 void AliEventTagCuts::Reset() {
166   //Sets dummy values to every private member.
167   fNParticipantsFlag = kFALSE;
168   fImpactParamFlag = kFALSE;
169
170   fVxFlag = kFALSE;
171   fVyFlag = kFALSE;
172   fVzFlag = kFALSE;
173   fPVFlag = kFALSE;
174   fPVzErrorFlag = kFALSE;
175
176   fTriggerMaskFlag = kFALSE;
177   fTriggerClusterFlag = kFALSE;
178
179   fZDCNeutron1EnergyFlag = kFALSE;
180   fZDCProton1EnergyFlag = kFALSE;
181   fZDCNeutron2EnergyFlag = kFALSE;
182   fZDCProton2EnergyFlag = kFALSE;
183   fZDCEMEnergyFlag = kFALSE;
184   fT0VertexZFlag = kFALSE;
185   fMultFlag = kFALSE;
186   fPosMultFlag = kFALSE;
187   fNegMultFlag = kFALSE;
188   fNeutrMultFlag = kFALSE;
189   fNV0sFlag = kFALSE;
190   fNCascadesFlag = kFALSE;
191   fNKinksFlag = kFALSE;
192
193   fNPMDTracksFlag = kFALSE;
194   fNFMDTracksFlag = kFALSE;
195   fNPHOSClustersFlag = kFALSE;
196   fNEMCALClustersFlag = kFALSE;
197   fNJetCandidatesFlag = kFALSE;
198
199   fTopJetEnergyMinFlag = kFALSE;
200   fNHardPhotonCandidatesFlag = kFALSE;
201   fTopNeutralEnergyMinFlag = kFALSE;
202   fNChargedAbove1GeVFlag = kFALSE;
203   fNChargedAbove3GeVFlag = kFALSE;
204   fNChargedAbove10GeVFlag = kFALSE;
205   fNMuonsAbove1GeVFlag = kFALSE;
206   fNMuonsAbove3GeVFlag = kFALSE;
207   fNMuonsAbove10GeVFlag = kFALSE;
208   fNElectronsAbove1GeVFlag = kFALSE;
209   fNElectronsAbove3GeVFlag = kFALSE;
210   fNElectronsAbove10GeVFlag = kFALSE;
211   fNElectronsFlag = kFALSE;
212   fNMuonsFlag = kFALSE;
213   fNPionsFlag = kFALSE;
214   fNKaonsFlag = kFALSE;
215   fNProtonsFlag = kFALSE;
216   fNLambdasFlag = kFALSE;
217   fNPhotonFlag = kFALSE;
218   fNPi0sFlag = kFALSE;
219   fNNeutronsFlag = kFALSE;
220   fNKaon0sFlag = kFALSE;
221   fTotalPFlag = kFALSE;
222   fMeanPtFlag = kFALSE;
223   fTopPtMinFlag = kFALSE;
224   fTotalNeutralPFlag = kFALSE;
225   fMeanNeutralPtFlag = kFALSE;
226   fTopNeutralPtMinFlag = kFALSE;
227   fEventPlaneAngleFlag = kFALSE;
228   fHBTRadiiFlag = kFALSE;
229   
230   fVxMin = -1000.0; fVxMax = 1000.0; 
231   fVyMin = -1000.0; fVyMax = 1000.0;  
232   fVzMin = -1000.0; fVzMax = 1000.0;
233   fNParticipantsMin = -1; fNParticipantsMax = 10000;
234   fImpactParamMin = -1.0; fImpactParamMax = 1000.0;
235   fPrimaryVertexFlag = 1;
236
237   fPrimaryVertexZErrorMin = -10000.; fPrimaryVertexZErrorMax = 10000.;
238   fTriggerMask = 0;
239   fTriggerCluster = 0;
240  
241   fZDCNeutron1EnergyMin = -1.0; fZDCNeutron1EnergyMax = 100000.0;
242   fZDCProton1EnergyMin = -1.0; fZDCProton1EnergyMax = 100000.0;
243   fZDCNeutron2EnergyMin = -1.0; fZDCNeutron2EnergyMax = 100000.0;
244   fZDCProton2EnergyMin = -1.0; fZDCProton2EnergyMax = 100000.0;
245   fZDCEMEnergyMin = -1.0; fZDCEMEnergyMax = 100000.0;
246   fT0VertexZMin = -10000.0; fT0VertexZMax = 10000.0;  
247   fMultMin = 0; fMultMax = 100000;  
248   fPosMultMin = -1; fPosMultMax = 100000;
249   fNegMultMin = -1; fNegMultMax = 100000;
250   fNeutrMultMin = -1; fNeutrMultMax = 100000;
251   fNV0sMin = -1; fNV0sMax = 1000000;
252   fNCascadesMin = -1; fNCascadesMax = 100000;
253   fNKinksMin = -1; fNKinksMax = 1000000;
254
255   fNPMDTracksMin = -1, fNPMDTracksMax = 100000;
256   fNFMDTracksMin = -1, fNFMDTracksMax = 100000;
257   fNPHOSClustersMin = -1, fNPHOSClustersMax = 100000;
258   fNEMCALClustersMin = -1, fNEMCALClustersMax = 100000;
259   fNJetCandidatesMin = -1, fNJetCandidatesMax = 100000;
260
261   fTopJetEnergyMin = -1.0; 
262   fNHardPhotonCandidatesMin = -1; fNHardPhotonCandidatesMax = 100000;
263   fTopNeutralEnergyMin = -1.0; 
264   fNChargedAbove1GeVMin = -1; fNChargedAbove1GeVMax = 100000;
265   fNChargedAbove3GeVMin = -1; fNChargedAbove3GeVMax = 100000;
266   fNChargedAbove10GeVMin = -1; fNChargedAbove10GeVMax = 100000;
267   fNMuonsAbove1GeVMin = -1; fNMuonsAbove1GeVMax = 100000;
268   fNMuonsAbove3GeVMin = -1; fNMuonsAbove3GeVMax = 100000;
269   fNMuonsAbove10GeVMin = -1; fNMuonsAbove10GeVMax = 100000; 
270   fNElectronsAbove1GeVMin = -1; fNElectronsAbove1GeVMax = 100000;
271   fNElectronsAbove3GeVMin = -1; fNElectronsAbove3GeVMax = 100000;
272   fNElectronsAbove10GeVMin = -1; fNElectronsAbove10GeVMax = 100000;
273   fNElectronsMin = -1; fNElectronsMax = 100000;
274   fNMuonsMin = -1; fNMuonsMax = 100000;
275   fNPionsMin = -1; fNPionsMax = 100000;
276   fNKaonsMin = -1; fNKaonsMax = 100000;
277   fNProtonsMin = -1; fNProtonsMax = 100000;
278   fNLambdasMin = -1; fNLambdasMax = 100000;
279   fNPhotonsMin = -1; fNPhotonsMax = 100000;
280   fNPi0sMin = -1; fNPi0sMax = 100000; 
281   fNNeutronsMin = -1; fNNeutronsMax = 100000; 
282   fNKaon0sMin = -1; fNKaon0sMax = 100000; 
283   fTotalPMin = -1.0; fTotalPMax = 1000000.0;
284   fMeanPtMin = -1.0; fMeanPtMax = 100000.0;
285   fTopPtMin = -1.0; fTotalNeutralPMin = -1.0;
286   fTotalNeutralPMax = 1000000.0;   
287   fMeanNeutralPtMin = -1.0; fMeanNeutralPtMax = 1000000.0; 
288   fTopNeutralPtMin = -1.0; 
289   fEventPlaneAngleMin = -10000000.0; fEventPlaneAngleMax = 10000000.0; 
290   fHBTRadiiMin = -1.0; fHBTRadiiMax = 100000.0; 
291 }
292
293 //___________________________________________________________________________
294 void AliEventTagCuts::SetPrimaryVertexXRange(Float_t low, Float_t high) {
295   //Sets the primary vertex x range 
296   //and the corresponding flag to kTRUE if the cut is used.
297   fVxMin = low;
298   fVxMax = high; 
299   fVxFlag = kTRUE;
300 }
301
302 //___________________________________________________________________________
303 void AliEventTagCuts::SetPrimaryVertexYRange(Float_t low, Float_t high) {
304   //Sets the primary vertex y range 
305   //and the corresponding flag to kTRUE if the cut is used.
306   fVyMin = low;
307   fVyMax = high; 
308   fVyFlag = kTRUE;
309 }
310
311 //___________________________________________________________________________
312 void AliEventTagCuts::SetPrimaryVertexZRange(Float_t low, Float_t high) {
313   //Sets the primary vertex z range 
314   //and the corresponding flag to kTRUE if the cut is used.
315   fVzMin = low;
316   fVzMax = high; 
317   fVzFlag = kTRUE;
318 }
319
320 //___________________________________________________________________________
321 void AliEventTagCuts::SetPrimaryVertexZErrorRange(Float_t low, Float_t high) {
322   //Sets the primary vertex z error range 
323   //and the corresponding flag to kTRUE if the cut is used.
324   fPrimaryVertexZErrorMin = low;
325   fPrimaryVertexZErrorMax = high; 
326   fPVzErrorFlag = kTRUE;
327 }
328
329 //___________________________________________________________________________
330 void AliEventTagCuts::SetTriggerMask(ULong64_t trmask) {
331   //Sets the trigger mask 
332   //and the corresponding flag to kTRUE if the cut is used.
333   fTriggerMask = trmask;
334   fTriggerMaskFlag = kTRUE;
335 }
336
337 //___________________________________________________________________________
338 void AliEventTagCuts::SetTriggerCluster(UChar_t trcluster) {
339   //Sets the trigger cluster 
340   //and the corresponding flag to kTRUE if the cut is used.
341   fTriggerCluster = trcluster;
342   fTriggerClusterFlag = kTRUE;
343 }
344
345 //___________________________________________________________________________
346 void AliEventTagCuts::SetMultiplicityRange(Int_t low, Int_t high) {
347   //Sets the primary multiplicity range 
348   //and the corresponding flag to kTRUE if the cut is used.
349   fMultMin = low;
350   fMultMax = high;
351   fMultFlag = kTRUE;
352 }
353
354 //___________________________________________________________________________
355 void AliEventTagCuts::SetNParticipantsRange(Int_t low, Int_t high) {
356   //Sets the number of participants range 
357   //and the corresponding flag to kTRUE if the cut is used.
358   fNParticipantsMin = low;
359   fNParticipantsMax = high;
360   fNParticipantsFlag = kTRUE;
361 }
362
363 //___________________________________________________________________________
364 void AliEventTagCuts::SetImpactParamRange(Float_t low, Float_t high) {
365   //Sets the impact parameter range 
366   //and the corresponding flag to kTRUE if the cut is used.
367   fImpactParamMin = low;
368   fImpactParamMax = high;
369   fImpactParamFlag = kTRUE;
370 }
371  
372
373 //___________________________________________________________________________
374 void AliEventTagCuts::SetPrimaryVertexFlag(Int_t flag) {
375   //Sets the primary vertex flag cut 
376   //and the corresponding flag to kTRUE if the cut is used.
377   fPrimaryVertexFlag = flag;
378   fPVFlag = kTRUE;
379 }
380
381 //___________________________________________________________________________
382 void AliEventTagCuts::SetZDCNeutron1Range(Float_t low, Float_t high) {
383   //Sets the ZDC's neutron energy range 
384   //and the corresponding flag to kTRUE if the cut is used.
385   fZDCNeutron1EnergyMin = low;
386   fZDCNeutron1EnergyMax = high;
387   fZDCNeutron1EnergyFlag = kTRUE;
388 }
389
390 //___________________________________________________________________________
391 void AliEventTagCuts::SetZDCProton1Range(Float_t low, Float_t high) {
392   //Sets the ZDC's proton energy range 
393   //and the corresponding flag to kTRUE if the cut is used.
394   fZDCProton1EnergyMin = low;
395   fZDCProton1EnergyMax = high;
396   fZDCProton1EnergyFlag = kTRUE;
397 }
398 //___________________________________________________________________________
399 void AliEventTagCuts::SetZDCNeutron2Range(Float_t low, Float_t high) {
400   //Sets the ZDC's neutron energy range 
401   //and the corresponding flag to kTRUE if the cut is used.
402   fZDCNeutron2EnergyMin = low;
403   fZDCNeutron2EnergyMax = high;
404   fZDCNeutron2EnergyFlag = kTRUE;
405 }
406 //___________________________________________________________________________
407 void AliEventTagCuts::SetZDCProton2Range(Float_t low, Float_t high) {
408   //Sets the ZDC's proton energy range 
409   //and the corresponding flag to kTRUE if the cut is used.
410   fZDCProton2EnergyMin = low;
411   fZDCProton2EnergyMax = high;
412   fZDCProton2EnergyFlag = kTRUE;
413 }
414 //___________________________________________________________________________
415 void AliEventTagCuts::SetZDCEMRange(Float_t low, Float_t high) {
416   //Sets the ZDC's e/m energy range 
417   //and the corresponding flag to kTRUE if the cut is used.
418   fZDCEMEnergyMin = low;
419   fZDCEMEnergyMax = high;
420   fZDCEMEnergyFlag = kTRUE;
421 }
422
423 //___________________________________________________________________________
424 void AliEventTagCuts::SetT0VertexZRange(Float_t low, Float_t high) {
425   //Sets the T0's Vz range 
426   //and the corresponding flag to kTRUE if the cut is used.
427   fT0VertexZMin = low;
428   fT0VertexZMax = high;
429   fT0VertexZFlag = kTRUE;
430 }
431
432 //___________________________________________________________________________
433 void AliEventTagCuts::SetPosMultiplicityRange(Int_t low, Int_t high) {
434   //Sets the positive multiplicity range 
435   //and the corresponding flag to kTRUE if the cut is used.
436   fPosMultMin = low;
437   fPosMultMax = high;
438   fPosMultFlag = kTRUE;
439 }
440
441
442 //___________________________________________________________________________
443 void AliEventTagCuts::SetNegMultiplicityRange(Int_t low, Int_t high) {
444   //Sets the negative multiplicity range 
445   //and the corresponding flag to kTRUE if the cut is used.
446   fNegMultMin = low;
447   fNegMultMax = high;
448   fNegMultFlag = kTRUE;
449 }
450
451
452 //___________________________________________________________________________
453 void AliEventTagCuts::SetNeutrMultiplicityRange(Int_t low, Int_t high) {
454   //Sets the neutral particle multiplicity range 
455   //and the corresponding flag to kTRUE if the cut is used.
456   fNeutrMultMin = low;
457   fNeutrMultMax = high;
458   fNeutrMultFlag = kTRUE;
459 }
460
461 //___________________________________________________________________________
462 void AliEventTagCuts::SetNV0sRange(Int_t low, Int_t high) {
463   //Sets the v0s multiplicity range 
464   //and the corresponding flag to kTRUE if the cut is used.
465   fNV0sMin = low;
466   fNV0sMax = high;
467   fNV0sFlag = kTRUE;
468 }
469
470 //___________________________________________________________________________
471 void AliEventTagCuts::SetNCascadesRange(Int_t low, Int_t high) {
472   //Sets the cascades multiplicity range 
473   //and the corresponding flag to kTRUE if the cut is used.
474   fNCascadesMin = low;
475   fNCascadesMax = high;
476   fNCascadesFlag = kTRUE;
477 }
478
479 //___________________________________________________________________________
480 void AliEventTagCuts::SetNKinksRange(Int_t low, Int_t high) {
481   //Sets the kinks multiplicity range 
482   //and the corresponding flag to kTRUE if the cut is used.
483   fNKinksMin = low;
484   fNKinksMax = high;
485   fNKinksFlag = kTRUE;
486 }
487
488 //___________________________________________________________________________
489 void AliEventTagCuts::SetNPMDTracksRange(Int_t low, Int_t high) {
490   //Sets the number of PMD tracks range 
491   //and the corresponding flag to kTRUE if the cut is used.
492   fNPMDTracksMin = low;
493   fNPMDTracksMax = high;
494   fNPMDTracksFlag = kTRUE;
495 }
496
497 //___________________________________________________________________________
498 void AliEventTagCuts::SetNFMDTracksRange(Int_t low, Int_t high) {
499   //Sets the number of FMD tracks range 
500   //and the corresponding flag to kTRUE if the cut is used.
501   fNFMDTracksMin = low;
502   fNFMDTracksMax = high;
503   fNFMDTracksFlag = kTRUE;
504 }
505
506 //___________________________________________________________________________
507 void AliEventTagCuts::SetNPHOSClustersRange(Int_t low, Int_t high) {
508   //Sets the number of PHOS clusters range 
509   //and the corresponding flag to kTRUE if the cut is used.
510   fNPHOSClustersMin = low;
511   fNPHOSClustersMax = high;
512   fNPHOSClustersFlag = kTRUE;
513 }
514
515 //___________________________________________________________________________
516 void AliEventTagCuts::SetNEMCALClustersRange(Int_t low, Int_t high) {
517   //Sets the number of EMCAL clusters range 
518   //and the corresponding flag to kTRUE if the cut is used.
519   fNEMCALClustersMin = low;
520   fNEMCALClustersMax = high;
521   fNEMCALClustersFlag = kTRUE;
522 }
523
524 //___________________________________________________________________________
525 void AliEventTagCuts::SetNJetCandidatesRange(Int_t low, Int_t high) {
526   //Sets the number of jet candidates range 
527   //and the corresponding flag to kTRUE if the cut is used.
528   fNJetCandidatesMin = low;
529   fNJetCandidatesMax = high;
530   fNJetCandidatesFlag = kTRUE;
531 }
532
533 //___________________________________________________________________________
534 void AliEventTagCuts::SetTopJetEnergyMin(Float_t low) {
535   //Sets the lower limit of the maximum jet energy
536   //and the corresponding flag to kTRUE if the cut is used.
537   fTopJetEnergyMin = low; 
538   fTopJetEnergyMinFlag = kTRUE;
539 }
540 //___________________________________________________________________________
541 void AliEventTagCuts::SetTopNeutralEnergyMin(Float_t low) {
542   //Sets the lower limit of the maximum neutral jet energy
543   //and the corresponding flag to kTRUE if the cut is used.
544   fTopNeutralEnergyMin = low; 
545   fTopNeutralEnergyMinFlag = kTRUE;
546 }
547 //___________________________________________________________________________
548 void AliEventTagCuts::SetNHardPhotonsRange(Int_t low, Int_t high) {
549   //Sets the hard photons multiplicity range
550   //and the corresponding flag to kTRUE if the cut is used.
551   fNHardPhotonCandidatesMin = low;
552   fNHardPhotonCandidatesMax = high;
553   fNHardPhotonCandidatesFlag = kTRUE;
554
555
556 //___________________________________________________________________________
557 void AliEventTagCuts::SetNChargedAbove1GeVRange(Int_t low, Int_t high) {
558   //Sets the number of charged above 1GeV range
559   //and the corresponding flag to kTRUE if the cut is used.
560   fNChargedAbove1GeVMin = low;
561   fNChargedAbove1GeVMax = high;
562   fNChargedAbove1GeVFlag = kTRUE;
563 }
564
565 //___________________________________________________________________________
566  void AliEventTagCuts::SetNChargedAbove3GeVRange(Int_t low, Int_t high) {
567   //Sets the number of charged above 3GeV range
568   //and the corresponding flag to kTRUE if the cut is used.
569   fNChargedAbove3GeVMin = low;
570   fNChargedAbove3GeVMax = high;
571   fNChargedAbove3GeVFlag = kTRUE;
572 }
573
574
575 //___________________________________________________________________________
576 void AliEventTagCuts::SetNChargedAbove10GeVRange(Int_t low, Int_t high) {
577   //Sets the number of charged above 10GeV range
578   //and the corresponding flag to kTRUE if the cut is used.
579   fNChargedAbove10GeVMin = low;
580   fNChargedAbove10GeVMax = high;
581   fNChargedAbove10GeVFlag = kTRUE;
582 }
583
584
585 //___________________________________________________________________________
586 void AliEventTagCuts::SetNMuonsAbove1GeVRange(Int_t low, Int_t high) {
587   //Sets the number of muons above 1GeV range
588   //and the corresponding flag to kTRUE if the cut is used.
589   fNMuonsAbove1GeVMin = low;
590   fNMuonsAbove1GeVMax = high;
591   fNMuonsAbove1GeVFlag = kTRUE;
592 }
593
594
595 //___________________________________________________________________________
596 void AliEventTagCuts::SetNMuonsAbove3GeVRange(Int_t low, Int_t high) {
597   //Sets the number of muons above 3GeV range
598   //and the corresponding flag to kTRUE if the cut is used.
599   fNMuonsAbove3GeVMin = low;
600   fNMuonsAbove3GeVMax = high;
601   fNMuonsAbove3GeVFlag = kTRUE;
602
603
604 //___________________________________________________________________________
605 void AliEventTagCuts::SetNMuonsAbove10GeVRange(Int_t low, Int_t high) {
606   //Sets the number of muons above 10GeV range
607   //and the corresponding flag to kTRUE if the cut is used.
608   fNMuonsAbove10GeVMin = low;
609   fNMuonsAbove10GeVMax = high; 
610   fNMuonsAbove10GeVFlag = kTRUE;
611 }
612
613
614 //___________________________________________________________________________
615 void AliEventTagCuts::SetNElectronsAbove1GeVRange(Int_t low, Int_t high) {
616   //Sets the number of electrons above 1GeV range
617   //and the corresponding flag to kTRUE if the cut is used.
618   fNElectronsAbove1GeVMin = low;
619   fNElectronsAbove1GeVMax = high;
620   fNElectronsAbove1GeVFlag = kTRUE;
621 }
622
623 //___________________________________________________________________________
624 void AliEventTagCuts::SetNElectronsAbove3GeVRange(Int_t low, Int_t high) {
625   //Sets the number of electrons above 3GeV range
626   //and the corresponding flag to kTRUE if the cut is used.
627   fNElectronsAbove3GeVMin = low;
628   fNElectronsAbove3GeVMax = high;
629   fNElectronsAbove3GeVFlag = kTRUE;
630 }
631
632 //___________________________________________________________________________
633 void AliEventTagCuts::SetNElectronsAbove10GeVRange(Int_t low, Int_t high) {  
634   //Sets the number of electrons above 10GeV range
635   //and the corresponding flag to kTRUE if the cut is used.
636   fNElectronsAbove10GeVMin = low;
637   fNElectronsAbove10GeVMax = high;
638   fNElectronsAbove10GeVFlag = kTRUE;
639 }
640 //___________________________________________________________________________
641 void AliEventTagCuts::SetNElectronRange(Int_t low, Int_t high) {
642   //Sets the electron multiplicity range
643   //and the corresponding flag to kTRUE if the cut is used.
644   fNElectronsMin = low;
645   fNElectronsMax = high;
646   fNElectronsFlag = kTRUE;
647 }
648 //___________________________________________________________________________
649 void AliEventTagCuts::SetNMuonRange(Int_t low, Int_t high) {
650   //Sets the muon multiplicity range
651   //and the corresponding flag to kTRUE if the cut is used.
652   fNMuonsMin = low;
653   fNMuonsMax = high;
654   fNMuonsFlag = kTRUE;
655 }
656
657 //___________________________________________________________________________
658 void AliEventTagCuts::SetNPionRange(Int_t low, Int_t high) {
659   //Sets the pion multiplicity range
660   //and the corresponding flag to kTRUE if the cut is used.
661   fNPionsMin = low;
662   fNPionsMax = high;
663   fNPionsFlag = kTRUE;
664
665
666 //___________________________________________________________________________
667 void AliEventTagCuts::SetNKaonRange(Int_t low, Int_t high) {
668   //Sets the kaon multiplicity range
669   //and the corresponding flag to kTRUE if the cut is used.
670   fNKaonsMin = low;
671   fNKaonsMax = high;
672   fNKaonsFlag = kTRUE;
673 }
674
675 //___________________________________________________________________________
676 void AliEventTagCuts::SetNProtonRange(Int_t low, Int_t high) {
677   //Sets the proton multiplicity range
678   //and the corresponding flag to kTRUE if the cut is used.
679   fNProtonsMin = low;
680   fNProtonsMax = high;
681   fNProtonsFlag = kTRUE;
682
683
684 //___________________________________________________________________________
685 void AliEventTagCuts::SetNLambdaRange(Int_t low, Int_t high) {
686   //Sets the lambda multiplicity range
687   //and the corresponding flag to kTRUE if the cut is used.
688   fNLambdasMin = low;
689   fNLambdasMax = high;
690   fNLambdasFlag = kTRUE;
691
692 //___________________________________________________________________________
693 void AliEventTagCuts::SetNPhotonRange(Int_t low, Int_t high) {
694   //Sets the photon multiplicity range
695   //and the corresponding flag to kTRUE if the cut is used.
696   fNPhotonsMin = low;
697   fNPhotonsMax = high;
698   fNPhotonFlag = kTRUE;
699
700 //___________________________________________________________________________
701 void AliEventTagCuts::SetNPi0Range(Int_t low, Int_t high) {
702   //Sets the pi0 multiplicity range
703   //and the corresponding flag to kTRUE if the cut is used.
704   fNPi0sMin = low;
705   fNPi0sMax = high; 
706   fNPi0sFlag = kTRUE;
707 }  
708
709 //___________________________________________________________________________
710 void AliEventTagCuts::SetNNeutronRange(Int_t low, Int_t high) {
711   //Sets the neutron multiplicity range
712   //and the corresponding flag to kTRUE if the cut is used.
713   fNNeutronsMin = low;
714   fNNeutronsMax = high; 
715   fNNeutronsFlag = kTRUE;
716 }
717
718 //___________________________________________________________________________
719 void AliEventTagCuts::SetNKaon0Range(Int_t low, Int_t high) {  
720   //Sets the K0s multiplicity range
721   //and the corresponding flag to kTRUE if the cut is used.
722   fNKaon0sMin = low;
723   fNKaon0sMax = high; 
724   fNKaon0sFlag = kTRUE;
725 }
726
727 //___________________________________________________________________________
728 void AliEventTagCuts::SetTotalPRange(Float_t low, Float_t high) {
729   //Sets the total momentum range
730   //and the corresponding flag to kTRUE if the cut is used.
731   fTotalPMin = low;
732   fTotalPMax = high;
733   fTotalPFlag = kTRUE;
734 }
735
736 //___________________________________________________________________________
737 void AliEventTagCuts::SetMeanPtRange(Float_t low, Float_t high) {
738   //Sets the mean Pt range
739   //and the corresponding flag to kTRUE if the cut is used.
740   fMeanPtMin = low;
741   fMeanPtMax = high;
742   fMeanPtFlag = kTRUE;
743 }  
744
745 //___________________________________________________________________________
746 void AliEventTagCuts::SetTopPtMin(Float_t low) {
747   //Sets the lower limit of the max Pt value
748   //and the corresponding flag to kTRUE if the cut is used.
749   fTopPtMin = low; 
750   fTopPtMinFlag = kTRUE;
751 }
752
753 //___________________________________________________________________________
754 void AliEventTagCuts::SetTotalNeutralPRange(Float_t low, Float_t high) {  
755   //Sets the total momentum of neutral particles range
756   //and the corresponding flag to kTRUE if the cut is used.
757   fTotalNeutralPMin =low ;
758   fTotalNeutralPMax = high;  
759   fTotalNeutralPFlag = kTRUE;
760 }
761 //___________________________________________________________________________
762 void AliEventTagCuts::SetMeanNeutralPtPRange(Float_t low, Float_t high) {  
763   //Sets the mean Pt of neutral particles range
764   //and the corresponding flag to kTRUE if the cut is used.
765   fMeanNeutralPtMin = low;
766   fMeanNeutralPtMax = high; 
767   fMeanNeutralPtFlag = kTRUE;
768 }
769 //___________________________________________________________________________
770 void AliEventTagCuts::SetTopNeutralPtMin(Float_t low) {  
771   //Sets the lower limit of the maximum Pt of neutral particles
772   //and the corresponding flag to kTRUE if the cut is used.
773   fTopNeutralPtMin = low; 
774   fTopNeutralPtMinFlag = kTRUE;
775 }
776
777 //___________________________________________________________________________
778 void AliEventTagCuts::SetEventPlaneAngleRange(Float_t low, Float_t high) {
779   //Sets the event plane range
780   //and the corresponding flag to kTRUE if the cut is used.
781   fEventPlaneAngleMin = low;
782   fEventPlaneAngleMax = high; 
783   fEventPlaneAngleFlag = kTRUE;
784 }
785
786 //___________________________________________________________________________
787 void AliEventTagCuts::SetHBTRadiiRange(Float_t low, Float_t high) {
788   //Sets the HBT radii range
789   //and the corresponding flag to kTRUE if the cut is used.
790   fHBTRadiiMin = low;
791   fHBTRadiiMax = high; 
792   fHBTRadiiFlag = kTRUE;
793 }
794
795 //___________________________________________________________________________
796 Bool_t AliEventTagCuts::IsAccepted(AliEventTag *EvTag) const {
797   //Returns true if the event is accepted otherwise false.
798   if(fVzFlag)
799     if((EvTag->GetVertexZ() < fVzMin) || (EvTag->GetVertexZ() > fVzMax))
800       return kFALSE;
801   
802   if(fVyFlag)
803     if((EvTag->GetVertexY() < fVyMin) || (EvTag->GetVertexY() > fVyMax))
804       return kFALSE;
805   
806   if(fVxFlag)
807     if((EvTag->GetVertexX() < fVxMin) || (EvTag->GetVertexX() > fVxMax))
808       return kFALSE;
809   
810   if(fNParticipantsFlag)
811     if((EvTag->GetNumOfParticipants() < fNParticipantsMin) || (EvTag->GetNumOfParticipants() > fNParticipantsMax))
812       return kFALSE; 
813   
814   if(fImpactParamFlag)
815     if((EvTag->GetImpactParameter() < fImpactParamMin) || (EvTag->GetImpactParameter() > fImpactParamMax))
816       return kFALSE; 
817   
818   if(fPVFlag)
819     if((EvTag->GetVertexFlag() != fPrimaryVertexFlag))
820       return kFALSE; 
821   
822   if(fPVzErrorFlag)
823     if((EvTag->GetVertexZError() < fPrimaryVertexZErrorMin) || (EvTag->GetVertexZError() > fPrimaryVertexZErrorMax))
824       return kFALSE; 
825   if(fTriggerMaskFlag)
826     if((EvTag->GetTriggerMask() != fTriggerMask))
827       return kFALSE; 
828   if(fTriggerClusterFlag)
829     if((EvTag->GetTriggerMask() != fTriggerMask))
830       return kFALSE; 
831
832   if(fZDCNeutron1EnergyFlag)
833     if((EvTag->GetZDCNeutron1Energy() < fZDCNeutron1EnergyMin) || (EvTag->GetZDCNeutron1Energy() > fZDCNeutron1EnergyMax))
834       return kFALSE; 
835   
836   if(fZDCProton1EnergyFlag)
837     if((EvTag->GetZDCProton1Energy() < fZDCProton1EnergyMin) || (EvTag->GetZDCProton1Energy() > fZDCProton1EnergyMax))
838       return kFALSE; 
839   
840   if(fZDCNeutron2EnergyFlag)
841     if((EvTag->GetZDCNeutron2Energy() < fZDCNeutron2EnergyMin) || (EvTag->GetZDCNeutron2Energy() > fZDCNeutron2EnergyMax))
842       return kFALSE; 
843   
844   if(fZDCProton2EnergyFlag)
845     if((EvTag->GetZDCProton2Energy() < fZDCProton2EnergyMin) || (EvTag->GetZDCProton2Energy() > fZDCProton2EnergyMax))
846       return kFALSE; 
847   
848   if(fZDCEMEnergyFlag)
849     if((EvTag->GetZDCEMEnergy() < fZDCEMEnergyMin) || (EvTag->GetZDCEMEnergy() > fZDCEMEnergyMax))
850       return kFALSE; 
851   
852   if(fT0VertexZFlag)
853     if((EvTag->GetT0VertexZ() < fT0VertexZMin) || (EvTag->GetT0VertexZ() > fT0VertexZMax))
854       return kFALSE; 
855   
856   if(fMultFlag)
857     if((EvTag->GetNumOfTracks() < fMultMin) || (EvTag->GetNumOfTracks() > fMultMax))
858       return kFALSE; 
859   
860   if(fPosMultFlag)
861     if((EvTag->GetNumOfPosTracks() < fPosMultMin) || (EvTag->GetNumOfPosTracks() > fPosMultMax))
862       return kFALSE; 
863   
864   if(fNegMultFlag)
865     if((EvTag->GetNumOfNegTracks() < fNegMultMin) || (EvTag->GetNumOfNegTracks() > fNegMultMax))
866       return kFALSE; 
867   
868   if(fNeutrMultFlag)
869     if((EvTag->GetNumOfNeutrTracks() < fNeutrMultMin) || (EvTag->GetNumOfNeutrTracks() > fNeutrMultMax))
870       return kFALSE; 
871   
872   if(fNV0sFlag)
873     if((EvTag->GetNumOfV0s() < fNV0sMin) || (EvTag->GetNumOfV0s() > fNV0sMax))
874       return kFALSE; 
875   
876   if(fNCascadesFlag)
877     if((EvTag->GetNumOfCascades() < fNCascadesMin) || (EvTag->GetNumOfCascades() > fNCascadesMax))
878       return kFALSE; 
879   
880   if(fNKinksFlag)
881     if((EvTag->GetNumOfKinks() < fNKinksMin) || (EvTag->GetNumOfKinks() > fNKinksMax))
882       return kFALSE; 
883
884
885   if(fNPMDTracksFlag)
886     if((EvTag->GetNumOfPMDTracks() < fNPMDTracksMin) || (EvTag->GetNumOfPMDTracks() > fNPMDTracksMax))
887       return kFALSE; 
888   if(fNFMDTracksFlag)
889     if((EvTag->GetNumOfFMDTracks() < fNFMDTracksMin) || (EvTag->GetNumOfFMDTracks() > fNFMDTracksMax))
890       return kFALSE; 
891   if(fNPHOSClustersFlag)
892     if((EvTag->GetNumOfPHOSClusters() < fNPHOSClustersMin) || (EvTag->GetNumOfPHOSClusters() > fNPHOSClustersMax))
893       return kFALSE; 
894   if(fNEMCALClustersFlag)
895     if((EvTag->GetNumOfEMCALClusters() < fNEMCALClustersMin) || (EvTag->GetNumOfEMCALClusters() > fNEMCALClustersMax))
896       return kFALSE; 
897   if(fNJetCandidatesFlag)
898     if((EvTag->GetNumOfJetCandidates() < fNJetCandidatesMin) || (EvTag->GetNumOfJetCandidates() > fNJetCandidatesMax))
899       return kFALSE; 
900
901
902   if(fTopJetEnergyMinFlag)
903     if((EvTag->GetMaxJetEnergy() < fTopJetEnergyMin))
904       return kFALSE; 
905   
906   if(fNHardPhotonCandidatesFlag)
907     if((EvTag->GetNumOfHardPhotonsCandidates() < fNHardPhotonCandidatesMin) || (EvTag->GetNumOfHardPhotonsCandidates() > fNHardPhotonCandidatesMax))
908       return kFALSE; 
909   
910   if(fTopNeutralEnergyMinFlag)
911     if((EvTag->GetMaxNeutralEnergy() < fTopNeutralEnergyMin))
912       return kFALSE; 
913   
914   if(fNChargedAbove1GeVFlag)
915     if((EvTag->GetNumOfChargedAbove1GeV() < fNChargedAbove1GeVMin) || (EvTag->GetNumOfChargedAbove1GeV() > fNChargedAbove1GeVMax))
916       return kFALSE; 
917   
918   if(fNChargedAbove3GeVFlag)
919     if((EvTag->GetNumOfChargedAbove3GeV() < fNChargedAbove3GeVMin) || (EvTag->GetNumOfChargedAbove3GeV() > fNChargedAbove3GeVMax))
920       return kFALSE; 
921   
922   if(fNChargedAbove10GeVFlag)
923     if((EvTag->GetNumOfChargedAbove10GeV() < fNChargedAbove10GeVMin) || (EvTag->GetNumOfChargedAbove10GeV() > fNChargedAbove10GeVMax))
924       return kFALSE; 
925   
926   if(fNMuonsAbove1GeVFlag)
927     if((EvTag->GetNumOfMuonsAbove1GeV() < fNMuonsAbove1GeVMin) || (EvTag->GetNumOfMuonsAbove1GeV() > fNMuonsAbove1GeVMax))
928       return kFALSE; 
929   
930   if(fNMuonsAbove3GeVFlag)
931     if((EvTag->GetNumOfMuonsAbove3GeV() < fNMuonsAbove3GeVMin) || (EvTag->GetNumOfMuonsAbove3GeV() > fNMuonsAbove3GeVMax))
932       return kFALSE; 
933   
934   if(fNMuonsAbove10GeVFlag)
935     if((EvTag->GetNumOfMuonsAbove10GeV() < fNMuonsAbove10GeVMin) || (EvTag->GetNumOfMuonsAbove10GeV() > fNMuonsAbove10GeVMax))
936       return kFALSE; 
937   
938   if(fNElectronsAbove1GeVFlag)
939     if((EvTag->GetNumOfElectronsAbove1GeV()  < fNElectronsAbove1GeVMin) || (EvTag->GetNumOfElectronsAbove1GeV()  > fNElectronsAbove1GeVMax))
940       return kFALSE; 
941   
942   if(fNElectronsAbove3GeVFlag)
943     if((EvTag->GetNumOfElectronsAbove3GeV() < fNElectronsAbove3GeVMin) || (EvTag->GetNumOfElectronsAbove3GeV() > fNElectronsAbove3GeVMax))
944       return kFALSE; 
945   
946   if(fNElectronsAbove10GeVFlag)
947     if((EvTag->GetNumOfElectronsAbove10GeV() < fNElectronsAbove10GeVMin) || (EvTag->GetNumOfElectronsAbove10GeV() > fNElectronsAbove10GeVMax))
948       return kFALSE; 
949   
950   if(fNElectronsFlag)
951     if((EvTag->GetNumOfElectrons() < fNElectronsMin) || (EvTag->GetNumOfElectrons() > fNElectronsMax))
952       return kFALSE; 
953   
954   if(fNMuonsFlag)
955     if((EvTag->GetNumOfMuons() < fNMuonsMin) || (EvTag->GetNumOfMuons() > fNMuonsMax))
956       return kFALSE; 
957   
958   if(fNPionsFlag)
959     if((EvTag->GetNumOfPions() < fNPionsMin) || (EvTag->GetNumOfPions() > fNPionsMax))
960       return kFALSE; 
961   
962   if(fNKaonsFlag)
963     if((EvTag->GetNumOfKaons() < fNKaonsMin) || (EvTag->GetNumOfKaons() > fNKaonsMax))
964       return kFALSE; 
965   
966   if(fNProtonsFlag)
967     if((EvTag->GetNumOfProtons() < fNProtonsMin) || (EvTag->GetNumOfProtons() > fNProtonsMax))
968       return kFALSE; 
969   
970   if(fNLambdasFlag)
971     if((EvTag->GetNumOfLambdas() < fNLambdasMin) || (EvTag->GetNumOfLambdas() > fNLambdasMax))
972       return kFALSE; 
973   
974   if(fNPhotonFlag)
975     if((EvTag->GetNumOfPhotons() < fNPhotonsMin) || (EvTag->GetNumOfPhotons() > fNPhotonsMax))
976       return kFALSE; 
977   
978   if(fNPi0sFlag)
979     if((EvTag->GetNumOfPi0s() < fNPi0sMin) || (EvTag->GetNumOfPi0s() > fNPi0sMax))
980       return kFALSE; 
981   
982   if(fNNeutronsFlag)
983     if((EvTag->GetNumOfNeutrons() < fNNeutronsMin) || (EvTag->GetNumOfNeutrons() > fNNeutronsMax))
984       return kFALSE; 
985   
986   if(fNKaon0sFlag)
987     if((EvTag->GetNumOfKaon0s() < fNKaon0sMin) || (EvTag->GetNumOfKaon0s() > fNKaon0sMax))
988       return kFALSE; 
989   
990   if(fTotalPFlag)
991     if((EvTag->GetTotalMomentum() < fTotalPMin) || (EvTag->GetTotalMomentum() > fTotalPMax))
992       return kFALSE; 
993   
994   if(fMeanPtFlag)
995     if((EvTag->GetMeanPt() < fMeanPtMin) || (EvTag->GetMeanPt() > fMeanPtMax))
996       return kFALSE; 
997   
998   if(fTopPtMinFlag)
999     if((EvTag->GetMaxPt() < fTopPtMin))
1000       return kFALSE; 
1001   
1002   if(fTotalNeutralPFlag)
1003     if((EvTag->GetNeutralTotalMomentum() < fTotalNeutralPMin) || (EvTag->GetNeutralTotalMomentum() > fTotalNeutralPMax))
1004       return kFALSE; 
1005   
1006   if(fMeanNeutralPtFlag)
1007     if((EvTag->GetNeutralMeanPt() < fMeanNeutralPtMin) || (EvTag->GetNeutralMeanPt() >fMeanNeutralPtMax ))
1008       return kFALSE; 
1009   
1010   if(fTopNeutralPtMinFlag)
1011     if((EvTag->GetNeutralMaxPt() < fTopNeutralPtMin))
1012       return kFALSE; 
1013   
1014   if(fEventPlaneAngleFlag)
1015     if((EvTag->GetEventPlaneAngle() < fEventPlaneAngleMin) || (EvTag->GetEventPlaneAngle() > fEventPlaneAngleMax))
1016       return kFALSE; 
1017   
1018   if(fHBTRadiiFlag)
1019     if((EvTag->GetHBTRadii() < fHBTRadiiMin) || (EvTag->GetHBTRadii() > fHBTRadiiMax))
1020       return kFALSE; 
1021   
1022   return kTRUE;
1023 }