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
  36
  37
  38   fVxMin(-1000.0), fVxMax(1000.0),
  39   fVxFlag(kFALSE),
  40   fVyMin(-1000.0), fVyMax(1000.0),
  41   fVyFlag(kFALSE),
  42   fVzMin(-1000.0), fVzMax(1000.0),
  43   fVzFlag(kFALSE),
  44   fParticipantsMin(-1), fParticipantMax(10000),
  45   fParticipantsFlag(kFALSE),
  46   fImpactParamMin(-1.0), fImpactParamMax(1000.0),
  47   fImpactParamFlag(kFALSE),
  48   fPrimaryVertexFlag(1),
  49   fPVFlag(kFALSE),
  50
  51   fPrimaryVertexZErrorMin(-10000.), fPrimaryVertexZErrorMax(10000.),
  52   fPVzErrorFlag(kFALSE),
  53   fTriggerMask(0),
  54   fTriggerMaskFlag(kFALSE),
  55   fTriggerCluster(0),
  56   fTriggerClusterFlag(kFALSE),
  57
  58   fZDCNeutron1EnergyMin(-1.0), fZDCNeutron1EnergyMax(100000.0),
  59   fZDCNeutron1EnergyFlag(kFALSE),
  60   fZDCProton1EnergyMin(-1.0), fZDCProton1EnergyMax(100000.0),
  61   fZDCProton1EnergyFlag(kFALSE),
  62   fZDCNeutron2EnergyMin(-1.0), fZDCNeutron2EnergyMax(100000.0),
  63   fZDCNeutron2EnergyFlag(kFALSE),
  64   fZDCProton2EnergyMin(-1.0), fZDCProton2EnergyMax(100000.0),
  65   fZDCProton2EnergyFlag(kFALSE),
  66   fZDCEMEnergyMin(-1.0), fZDCEMEnergyMax(100000.0),
  67   fZDCEMEnergyFlag(kFALSE),
  68   fT0VertexZMin(-10000.0), fT0VertexZMax(10000.0),
  69   fT0VertexZFlag(kFALSE),
  70   fMultMin(0), fMultMax(100000),
  71   fMultFlag(kFALSE),
  72   fMultPosMin(-1), fMultPosMax(100000),
  73   fMultPosFlag(kFALSE),
  74   fMultNegMin(-1), fMultNegMax(100000),
  75   fMultNegFlag(kFALSE),
  76   fMultNeutrMin(-1), fMultNeutrMax(100000),
  77   fMultNeutrFlag(kFALSE),
  78   fV0sMin(-1), fV0sMax(1000000),
  79   fV0sFlag(kFALSE),
  80   fCascadesMin(-1), fCascadesMax(100000),
  81   fCascadesFlag(kFALSE),
  82   fkinksMin(-1), fkinksMax(1000000),
  83   fkinksFlag(kFALSE),
  84
  85   fPMDTracksMin(-1), fPMDTracksMax(100000),
  86   fPMDTracksFlag(kFALSE),
  87   fFMDTracksMin(-1), fFMDTracksMax(100000),
  88   fFMDTracksFlag(kFALSE),
  89   fPHOSClustersMin(-1), fPHOSClustersMax(100000),
  90   fPHOSClustersFlag(kFALSE),
  91   fEMCALClustersMin(-1), fEMCALClustersMax(100000),
  92   fEMCALClustersFlag(kFALSE),
  93   fJetCandidatesMin(-1), fJetCandidatesMax(100000),
  94   fJetCandidatesFlag(kFALSE),
  95
  96   fMaxJetEnergy(-1.0),
  97   fMaxJetEnergyFlag(kFALSE),
  98   fNHardPhotonsCandidatesMin(-1), fNHardPhotonsCandidatesMax(100000),
  99   fNHardPhotonsCandidatesFlag(kFALSE),
 100   fMaxNeutralEnergy(-1.0),
 101   fMaxNeutralFlag(kFALSE),
 102   fChargedAbove1GeVMin(-1), fChargedAbove1GeVMax(100000),
 103   fChargedAbove1GeVFlag(kFALSE),
 104   fChargedAbove3GeVMin(-1), fChargedAbove3GeVMax(100000),
 105   fChargedAbove3GeVFlag(kFALSE),
 106   fChargedAbove10GeVMin(-1), fChargedAbove10GeVMax(100000),
 107   fChargedAbove10GeVFlag(kFALSE),
 108   fMuonsAbove1GeVMin(-1), fMuonsAbove1GeVMax(100000),
 109   fMuonsAbove1GeVFlag(kFALSE),
 110   fMuonsAbove3GeVMin(-1), fMuonsAbove3GeVMax(100000),
 111   fMuonsAbove3GeVFlag(kFALSE),
 112   fMuonsAbove10GeVMin(-1), fMuonsAbove10GeVMax(100000),
 113   fMuonsAbove10GeVFlag(kFALSE),
 114   fElectronsAbove1GeVMin(-1), fElectronsAbove1GeVMax(100000),
 115   fElectronsAbove1GeVFlag(kFALSE),
 116   fElectronsAbove3GeVMin(-1), fElectronsAbove3GeVMax(100000),
 117   fElectronsAbove3GeVFlag(kFALSE),
 118   fElectronsAbove10GeVMin(-1), fElectronsAbove10GeVMax(100000),
 119   fElectronsAbove10GeVFlag(kFALSE),
 120   fElectronsMin(-1), fElectronsMax(100000),
 121   fElectronsFlag(kFALSE),
 122   fMuonsMin(-1), fMuonsMax(100000),
 123   fMuonsFlag(kFALSE),
 124   fPionsMin(-1), fPionsMax(100000),
 125   fPionsFlag(kFALSE),
 126   fKaonsMin(-1), fKaonsMax(100000),
 127   fKaonsFlag(kFALSE),
 128   fProtonsMin(-1), fProtonsMax(100000),
 129   fProtonsFlag(kFALSE),
 130   fLambdasMin(-1), fLambdasMax(100000),
 131   fLambdasFlag(kFALSE),
 132   fPhotonsMin(-1), fPhotonsMax(100000),
 133   fPhotonFlag(kFALSE),
 134   fPi0sMin(-1), fPi0sMax(100000),
 135   fPi0sFlag(kFALSE),
 136   fNeutronsMin(-1), fNeutronsMax(100000),
 137   fNeutronsFlag(kFALSE),
 138   fKaon0sMin(-1), fKaon0sMax(100000),
 139   fKaon0sFlag(kFALSE),
 140   fTotalPMin(-1.0), fTotalPMax(1000000.0),
 141   fTotalPFlag(kFALSE),
 142   fMeanPtMin(-1.0), fMeanPtMax(100000.0),
 143   fMeanPtFlag(kFALSE),
 144   fMaxPt(-1.0),
 145   fMaxPtFlag(kFALSE),
 146   fTotalNeutralPMin(-1.0), fTotalNeutralPMax(1000000.0),
 147   fTotalNeutralPFlag(kFALSE),
 148   fMeanNeutralPtMin(-1.0), fMeanNeutralPtMax(1000000.0),
 149   fMeanNeutralPtFlag(kFALSE),
 150   fMaxNeutralPt(-1.0),
 151   fMaxNeutralPtFlag(kFALSE),
 152   fEventPlaneAngleMin(-10000000.0), fEventPlaneAngleMax(10000000.0),
 153   fEventPlaneAngleFlag(kFALSE),
 154   fHBTRadiiMin(-1.0), fHBTRadiiMax(100000.0),
 155   fHBTRadiiFlag(kFALSE)
 156 {
 157   //Default constructor which calls the Reset method.
 158   Reset();
 159 }
 160
 161 //___________________________________________________________________________
 162 AliEventTagCuts::~AliEventTagCuts() {
 163   //Defaut destructor.
 164 }
 165
 166 //___________________________________________________________________________
 167 void AliEventTagCuts::Reset() {
 168   //Sets dummy values to every private member.
 169   fVxFlag = kFALSE;
 170   fVyFlag = kFALSE;
 171   fVzFlag = kFALSE;
 172   fParticipantsFlag = kFALSE;
 173   fImpactParamFlag = kFALSE;
 174   fPVFlag = kFALSE;
 175
 176   fPVzErrorFlag = kFALSE;
 177   fTriggerMaskFlag = kFALSE;
 178   fTriggerClusterFlag = kFALSE;
 179
 180   fZDCNeutron1EnergyFlag = kFALSE;
 181   fZDCProton1EnergyFlag = kFALSE;
 182   fZDCNeutron2EnergyFlag = kFALSE;
 183   fZDCProton2EnergyFlag = kFALSE;
 184   fZDCEMEnergyFlag = kFALSE;
 185   fT0VertexZFlag = kFALSE;
 186   fMultFlag = kFALSE;
 187   fMultPosFlag = kFALSE;
 188   fMultNegFlag = kFALSE;
 189   fMultNeutrFlag = kFALSE;
 190   fV0sFlag = kFALSE;
 191   fCascadesFlag = kFALSE;
 192   fkinksFlag = kFALSE;
 193
 194   fPMDTracksFlag = kFALSE;
 195   fFMDTracksFlag = kFALSE;
 196   fPHOSClustersFlag = kFALSE;
 197   fEMCALClustersFlag = kFALSE;
 198   fJetCandidatesFlag = kFALSE;
 199
 200   fMaxJetEnergyFlag = kFALSE;
 201   fNHardPhotonsCandidatesFlag = kFALSE;
 202   fMaxNeutralFlag = kFALSE;
 203   fChargedAbove1GeVFlag = kFALSE;
 204   fChargedAbove3GeVFlag = kFALSE;
 205   fChargedAbove10GeVFlag = kFALSE;
 206   fMuonsAbove1GeVFlag = kFALSE;
 207   fMuonsAbove3GeVFlag = kFALSE;
 208   fMuonsAbove10GeVFlag = kFALSE;
 209   fElectronsAbove1GeVFlag = kFALSE;
 210   fElectronsAbove3GeVFlag = kFALSE;
 211   fElectronsAbove10GeVFlag = kFALSE;
 212   fElectronsFlag = kFALSE;
 213   fMuonsFlag = kFALSE;
 214   fPionsFlag = kFALSE;
 215   fKaonsFlag = kFALSE;
 216   fProtonsFlag = kFALSE;
 217   fLambdasFlag = kFALSE;
 218   fPhotonFlag = kFALSE;
 219   fPi0sFlag = kFALSE;
 220   fNeutronsFlag = kFALSE;
 221   fKaon0sFlag = kFALSE;
 222   fTotalPFlag = kFALSE;
 223   fMeanPtFlag = kFALSE;
 224   fMaxPtFlag = kFALSE;
 225   fTotalNeutralPFlag = kFALSE;
 226   fMeanNeutralPtFlag = kFALSE;
 227   fMaxNeutralPtFlag = kFALSE;
 228   fEventPlaneAngleFlag = kFALSE;
 229   fHBTRadiiFlag = kFALSE;
 230
 231   fVxMin = -1000.0; fVxMax = 1000.0;
 232   fVyMin = -1000.0; fVyMax = 1000.0;
 233   fVzMin = -1000.0; fVzMax = 1000.0;
 234   fParticipantsMin = -1; fParticipantMax = 10000;
 235   fImpactParamMin = -1.0; fImpactParamMax = 1000.0;
 236   fPrimaryVertexFlag = 1;
 237
 238   fPrimaryVertexZErrorMin = -10000.; fPrimaryVertexZErrorMax = 10000.;
 239   fTriggerMask = 0;
 240   fTriggerCluster = 0;
 241
 242   fZDCNeutron1EnergyMin = -1.0; fZDCNeutron1EnergyMax = 100000.0;
 243   fZDCProton1EnergyMin = -1.0; fZDCProton1EnergyMax = 100000.0;
 244   fZDCNeutron2EnergyMin = -1.0; fZDCNeutron2EnergyMax = 100000.0;
 245   fZDCProton2EnergyMin = -1.0; fZDCProton2EnergyMax = 100000.0;
 246   fZDCEMEnergyMin = -1.0; fZDCEMEnergyMax = 100000.0;
 247   fT0VertexZMin = -10000.0; fT0VertexZMax = 10000.0;
 248   fMultMin = 0; fMultMax = 100000;
 249   fMultPosMin = -1; fMultPosMax = 100000;
 250   fMultNegMin = -1; fMultNegMax = 100000;
 251   fMultNeutrMin = -1; fMultNeutrMax = 100000;
 252   fV0sMin = -1; fV0sMax = 1000000;
 253   fCascadesMin = -1; fCascadesMax = 100000;
 254   fkinksMin = -1; fkinksMax = 1000000;
 255
 256   fPMDTracksMin = -1, fPMDTracksMax = 100000;
 257   fFMDTracksMin = -1, fFMDTracksMax = 100000;
 258   fPHOSClustersMin = -1, fPHOSClustersMax = 100000;
 259   fEMCALClustersMin = -1, fEMCALClustersMax = 100000;
 260   fJetCandidatesMin = -1, fJetCandidatesMax = 100000;
 261
 262   fMaxJetEnergy = -1.0;
 263   fNHardPhotonsCandidatesMin = -1; fNHardPhotonsCandidatesMax = 100000;
 264   fMaxNeutralEnergy = -1.0;
 265   fChargedAbove1GeVMin = -1; fChargedAbove1GeVMax = 100000;
 266   fChargedAbove3GeVMin = -1; fChargedAbove3GeVMax = 100000;
 267   fChargedAbove10GeVMin = -1; fChargedAbove10GeVMax = 100000;
 268   fMuonsAbove1GeVMin = -1; fMuonsAbove1GeVMax = 100000;
 269   fMuonsAbove3GeVMin = -1; fMuonsAbove3GeVMax = 100000;
 270   fMuonsAbove10GeVMin = -1; fMuonsAbove10GeVMax = 100000;
 271   fElectronsAbove1GeVMin = -1; fElectronsAbove1GeVMax = 100000;
 272   fElectronsAbove3GeVMin = -1; fElectronsAbove3GeVMax = 100000;
 273   fElectronsAbove10GeVMin = -1; fElectronsAbove10GeVMax = 100000;
 274   fElectronsMin = -1; fElectronsMax = 100000;
 275   fMuonsMin = -1; fMuonsMax = 100000;
 276   fPionsMin = -1; fPionsMax = 100000;
 277   fKaonsMin = -1; fKaonsMax = 100000;
 278   fProtonsMin = -1; fProtonsMax = 100000;
 279   fLambdasMin = -1; fLambdasMax = 100000;
 280   fPhotonsMin = -1; fPhotonsMax = 100000;
 281   fPi0sMin = -1; fPi0sMax = 100000;
 282   fNeutronsMin = -1; fNeutronsMax = 100000;
 283   fKaon0sMin = -1; fKaon0sMax = 100000;
 284   fTotalPMin = -1.0; fTotalPMax = 1000000.0;
 285   fMeanPtMin = -1.0; fMeanPtMax = 100000.0;
 286   fMaxPt = -1.0; fTotalNeutralPMin = -1.0;
 287   fTotalNeutralPMax = 1000000.0;
 288   fMeanNeutralPtMin = -1.0; fMeanNeutralPtMax = 1000000.0;
 289   fMaxNeutralPt = -1.0;
 290   fEventPlaneAngleMin = -10000000.0; fEventPlaneAngleMax = 10000000.0;
 291   fHBTRadiiMin = -1.0; fHBTRadiiMax = 100000.0;
 292 }
 293
 294 //___________________________________________________________________________
 295 void AliEventTagCuts::SetPrimaryVertexXRange(Float_t low, Float_t high) {
 296   //Sets the primary vertex x range
 297   //and the corresponding flag to kTRUE if the cut is used.
 298   fVxMin = low;
 299   fVxMax = high;
 300   fVxFlag = kTRUE;
 301 }
 302
 303 //___________________________________________________________________________
 304 void AliEventTagCuts::SetPrimaryVertexYRange(Float_t low, Float_t high) {
 305   //Sets the primary vertex y range
 306   //and the corresponding flag to kTRUE if the cut is used.
 307   fVyMin = low;
 308   fVyMax = high;
 309   fVyFlag = kTRUE;
 310 }
 311
 312 //___________________________________________________________________________
 313 void AliEventTagCuts::SetPrimaryVertexZRange(Float_t low, Float_t high) {
 314   //Sets the primary vertex z range
 315   //and the corresponding flag to kTRUE if the cut is used.
 316   fVzMin = low;
 317   fVzMax = high;
 318   fVzFlag = kTRUE;
 319 }
 320
 321 //___________________________________________________________________________
 322 void AliEventTagCuts::SetPrimaryVertexZErrorRange(Float_t low, Float_t high) {
 323   //Sets the primary vertex z error range
 324   //and the corresponding flag to kTRUE if the cut is used.
 325   fPrimaryVertexZErrorMin = low;
 326   fPrimaryVertexZErrorMax = high;
 327   fPVzErrorFlag = kTRUE;
 328 }
 329
 330 //___________________________________________________________________________
 331 void AliEventTagCuts::SetTriggerMask(ULong64_t trmask) {
 332   //Sets the trigger mask
 333   //and the corresponding flag to kTRUE if the cut is used.
 334   fTriggerMask = trmask;
 335   fTriggerMaskFlag = kTRUE;
 336 }
 337
 338 //___________________________________________________________________________
 339 void AliEventTagCuts::SetTriggerCluster(UChar_t trcluster) {
 340   //Sets the trigger cluster
 341   //and the corresponding flag to kTRUE if the cut is used.
 342   fTriggerCluster = trcluster;
 343   fTriggerClusterFlag = kTRUE;
 344 }
 345
 346 //___________________________________________________________________________
 347 void AliEventTagCuts::SetMultiplicityRange(Int_t low, Int_t high) {
 348   //Sets the primary multiplicity range
 349   //and the corresponding flag to kTRUE if the cut is used.
 350   fMultMin = low;
 351   fMultMax = high;
 352   fMultFlag = kTRUE;
 353 }
 354
 355 //___________________________________________________________________________
 356 void AliEventTagCuts::SetNParticipantsRange(Int_t low, Int_t high) {
 357   //Sets the number of participants range
 358   //and the corresponding flag to kTRUE if the cut is used.
 359   fParticipantsMin = low;
 360   fParticipantMax = high;
 361   fParticipantsFlag = kTRUE;
 362 }
 363
 364 //___________________________________________________________________________
 365 void AliEventTagCuts::SetImpactParamRange(Float_t low, Float_t high) {
 366   //Sets the impact parameter range
 367   //and the corresponding flag to kTRUE if the cut is used.
 368   fImpactParamMin = low;
 369   fImpactParamMax = high;
 370   fImpactParamFlag = kTRUE;
 371 }
 372
 373
 374 //___________________________________________________________________________
 375 void AliEventTagCuts::SetPrimaryVertexFlag(Int_t flag) {
 376   //Sets the primary vertex flag cut
 377   //and the corresponding flag to kTRUE if the cut is used.
 378   fPrimaryVertexFlag = flag;
 379   fPVFlag = kTRUE;
 380 }
 381
 382 //___________________________________________________________________________
 383 void AliEventTagCuts::SetZDCNeutron1Range(Float_t low, Float_t high) {
 384   //Sets the ZDC's neutron energy range
 385   //and the corresponding flag to kTRUE if the cut is used.
 386   fZDCNeutron1EnergyMin = low;
 387   fZDCNeutron1EnergyMax = high;
 388   fZDCNeutron1EnergyFlag = kTRUE;
 389 }
 390
 391 //___________________________________________________________________________
 392 void AliEventTagCuts::SetZDCProton1Range(Float_t low, Float_t high) {
 393   //Sets the ZDC's proton energy range
 394   //and the corresponding flag to kTRUE if the cut is used.
 395   fZDCProton1EnergyMin = low;
 396   fZDCProton1EnergyMax = high;
 397   fZDCProton1EnergyFlag = kTRUE;
 398 }
 399 //___________________________________________________________________________
 400 void AliEventTagCuts::SetZDCNeutron2Range(Float_t low, Float_t high) {
 401   //Sets the ZDC's neutron energy range
 402   //and the corresponding flag to kTRUE if the cut is used.
 403   fZDCNeutron2EnergyMin = low;
 404   fZDCNeutron2EnergyMax = high;
 405   fZDCNeutron2EnergyFlag = kTRUE;
 406 }
 407 //___________________________________________________________________________
 408 void AliEventTagCuts::SetZDCProton2Range(Float_t low, Float_t high) {
 409   //Sets the ZDC's proton energy range
 410   //and the corresponding flag to kTRUE if the cut is used.
 411   fZDCProton2EnergyMin = low;
 412   fZDCProton2EnergyMax = high;
 413   fZDCProton2EnergyFlag = kTRUE;
 414 }
 415 //___________________________________________________________________________
 416 void AliEventTagCuts::SetZDCEMRange(Float_t low, Float_t high) {
 417   //Sets the ZDC's e/m energy range
 418   //and the corresponding flag to kTRUE if the cut is used.
 419   fZDCEMEnergyMin = low;
 420   fZDCEMEnergyMax = high;
 421   fZDCEMEnergyFlag = kTRUE;
 422 }
 423
 424 //___________________________________________________________________________
 425 void AliEventTagCuts::SetT0VertexZRange(Float_t low, Float_t high) {
 426   //Sets the T0's Vz range
 427   //and the corresponding flag to kTRUE if the cut is used.
 428   fT0VertexZMin = low;
 429   fT0VertexZMax = high;
 430   fT0VertexZFlag = kTRUE;
 431 }
 432
 433 //___________________________________________________________________________
 434 void AliEventTagCuts::SetPosMultiplicityRange(Int_t low, Int_t high) {
 435   //Sets the positive multiplicity range
 436   //and the corresponding flag to kTRUE if the cut is used.
 437   fMultPosMin = low;
 438   fMultPosMax = high;
 439   fMultPosFlag = kTRUE;
 440 }
 441
 442
 443 //___________________________________________________________________________
 444 void AliEventTagCuts::SetNegMultiplicityRange(Int_t low, Int_t high) {
 445   //Sets the negative multiplicity range
 446   //and the corresponding flag to kTRUE if the cut is used.
 447   fMultNegMin = low;
 448   fMultNegMax = high;
 449   fMultNegFlag = kTRUE;
 450 }
 451
 452
 453 //___________________________________________________________________________
 454 void AliEventTagCuts::SetNeutrMultiplicityRange(Int_t low, Int_t high) {
 455   //Sets the neutral particle multiplicity range
 456   //and the corresponding flag to kTRUE if the cut is used.
 457   fMultNeutrMin = low;
 458   fMultNeutrMax = high;
 459   fMultNeutrFlag = kTRUE;
 460 }
 461
 462 //___________________________________________________________________________
 463 void AliEventTagCuts::SetNV0sRange(Int_t low, Int_t high) {
 464   //Sets the v0s multiplicity range
 465   //and the corresponding flag to kTRUE if the cut is used.
 466   fV0sMin = low;
 467   fV0sMax = high;
 468   fV0sFlag = kTRUE;
 469 }
 470
 471 //___________________________________________________________________________
 472 void AliEventTagCuts::SetNCascadesRange(Int_t low, Int_t high) {
 473   //Sets the cascades multiplicity range
 474   //and the corresponding flag to kTRUE if the cut is used.
 475   fCascadesMin = low;
 476   fCascadesMax = high;
 477   fCascadesFlag = kTRUE;
 478 }
 479
 480 //___________________________________________________________________________
 481 void AliEventTagCuts::SetNKinksRange(Int_t low, Int_t high) {
 482   //Sets the kinks multiplicity range
 483   //and the corresponding flag to kTRUE if the cut is used.
 484   fkinksMin = low;
 485   fkinksMax = high;
 486   fkinksFlag = kTRUE;
 487 }
 488
 489 //___________________________________________________________________________
 490 void AliEventTagCuts::SetNPMDTracksRange(Int_t low, Int_t high) {
 491   //Sets the number of PMD tracks range
 492   //and the corresponding flag to kTRUE if the cut is used.
 493   fPMDTracksMin = low;
 494   fPMDTracksMax = high;
 495   fPMDTracksFlag = kTRUE;
 496 }
 497
 498 //___________________________________________________________________________
 499 void AliEventTagCuts::SetNFMDTracksRange(Int_t low, Int_t high) {
 500   //Sets the number of FMD tracks range
 501   //and the corresponding flag to kTRUE if the cut is used.
 502   fFMDTracksMin = low;
 503   fFMDTracksMax = high;
 504   fFMDTracksFlag = kTRUE;
 505 }
 506
 507 //___________________________________________________________________________
 508 void AliEventTagCuts::SetNPHOSClustersRange(Int_t low, Int_t high) {
 509   //Sets the number of PHOS clusters range
 510   //and the corresponding flag to kTRUE if the cut is used.
 511   fPHOSClustersMin = low;
 512   fPHOSClustersMax = high;
 513   fPHOSClustersFlag = kTRUE;
 514 }
 515
 516 //___________________________________________________________________________
 517 void AliEventTagCuts::SetNEMCALClustersRange(Int_t low, Int_t high) {
 518   //Sets the number of EMCAL clusters range
 519   //and the corresponding flag to kTRUE if the cut is used.
 520   fEMCALClustersMin = low;
 521   fEMCALClustersMax = high;
 522   fEMCALClustersFlag = kTRUE;
 523 }
 524
 525 //___________________________________________________________________________
 526 void AliEventTagCuts::SetNJetCandidatesRange(Int_t low, Int_t high) {
 527   //Sets the number of jet candidates range
 528   //and the corresponding flag to kTRUE if the cut is used.
 529   fJetCandidatesMin = low;
 530   fJetCandidatesMax = high;
 531   fJetCandidatesFlag = kTRUE;
 532 }
 533
 534 //___________________________________________________________________________
 535 void AliEventTagCuts::SetTopJetEnergyMin(Float_t low) {
 536   //Sets the lower limit of the maximum jet energy
 537   //and the corresponding flag to kTRUE if the cut is used.
 538   fMaxJetEnergy = low;
 539   fMaxJetEnergyFlag = kTRUE;
 540 }
 541 //___________________________________________________________________________
 542 void AliEventTagCuts::SetTopNeutralEnergyMin(Float_t low) {
 543   //Sets the lower limit of the maximum neutral jet energy
 544   //and the corresponding flag to kTRUE if the cut is used.
 545   fMaxNeutralEnergy = low;
 546   fMaxNeutralFlag = kTRUE;
 547 }
 548 //___________________________________________________________________________
 549 void AliEventTagCuts::SetNHardPhotonsRange(Int_t low, Int_t high) {
 550   //Sets the hard photons multiplicity range
 551   //and the corresponding flag to kTRUE if the cut is used.
 552   fNHardPhotonsCandidatesMin = low;
 553   fNHardPhotonsCandidatesMax = high;
 554   fNHardPhotonsCandidatesFlag = kTRUE;
 555 }
 556
 557 //___________________________________________________________________________
 558 void AliEventTagCuts::SetNChargedAbove1GeVRange(Int_t low, Int_t high) {
 559   //Sets the number of charged above 1GeV range
 560   //and the corresponding flag to kTRUE if the cut is used.
 561   fChargedAbove1GeVMin = low;
 562   fChargedAbove1GeVMax = high;
 563   fChargedAbove1GeVFlag = kTRUE;
 564 }
 565
 566 //___________________________________________________________________________
 567  void AliEventTagCuts::SetNChargedAbove3GeVRange(Int_t low, Int_t high) {
 568   //Sets the number of charged above 3GeV range
 569   //and the corresponding flag to kTRUE if the cut is used.
 570   fChargedAbove3GeVMin = low;
 571   fChargedAbove3GeVMax = high;
 572   fChargedAbove3GeVFlag = kTRUE;
 573 }
 574
 575
 576 //___________________________________________________________________________
 577 void AliEventTagCuts::SetNChargedAbove10GeVRange(Int_t low, Int_t high) {
 578   //Sets the number of charged above 10GeV range
 579   //and the corresponding flag to kTRUE if the cut is used.
 580   fChargedAbove10GeVMin = low;
 581   fChargedAbove10GeVMax = high;
 582   fChargedAbove10GeVFlag = kTRUE;
 583 }
 584
 585
 586 //___________________________________________________________________________
 587 void AliEventTagCuts::SetNMuonsAbove1GeVRange(Int_t low, Int_t high) {
 588   //Sets the number of muons above 1GeV range
 589   //and the corresponding flag to kTRUE if the cut is used.
 590   fMuonsAbove1GeVMin = low;
 591   fMuonsAbove1GeVMax = high;
 592   fMuonsAbove1GeVFlag = kTRUE;
 593 }
 594
 595
 596 //___________________________________________________________________________
 597 void AliEventTagCuts::SetNMuonsAbove3GeVRange(Int_t low, Int_t high) {
 598   //Sets the number of muons above 3GeV range
 599   //and the corresponding flag to kTRUE if the cut is used.
 600   fMuonsAbove3GeVMin = low;
 601   fMuonsAbove3GeVMax = high;
 602   fMuonsAbove3GeVFlag = kTRUE;
 603 }
 604
 605 //___________________________________________________________________________
 606 void AliEventTagCuts::SetNMuonsAbove10GeVRange(Int_t low, Int_t high) {
 607   //Sets the number of muons above 10GeV range
 608   //and the corresponding flag to kTRUE if the cut is used.
 609   fMuonsAbove10GeVMin = low;
 610   fMuonsAbove10GeVMax = high;
 611   fMuonsAbove10GeVFlag = kTRUE;
 612 }
 613
 614
 615 //___________________________________________________________________________
 616 void AliEventTagCuts::SetNElectronsAbove1GeVRange(Int_t low, Int_t high) {
 617   //Sets the number of electrons above 1GeV range
 618   //and the corresponding flag to kTRUE if the cut is used.
 619   fElectronsAbove1GeVMin = low;
 620   fElectronsAbove1GeVMax = high;
 621   fElectronsAbove1GeVFlag = kTRUE;
 622 }
 623
 624 //___________________________________________________________________________
 625 void AliEventTagCuts::SetNElectronsAbove3GeVRange(Int_t low, Int_t high) {
 626   //Sets the number of electrons above 3GeV range
 627   //and the corresponding flag to kTRUE if the cut is used.
 628   fElectronsAbove3GeVMin = low;
 629   fElectronsAbove3GeVMax = high;
 630   fElectronsAbove3GeVFlag = kTRUE;
 631 }
 632
 633 //___________________________________________________________________________
 634 void AliEventTagCuts::SetNElectronsAbove10GeVRange(Int_t low, Int_t high) {
 635   //Sets the number of electrons above 10GeV range
 636   //and the corresponding flag to kTRUE if the cut is used.
 637   fElectronsAbove10GeVMin = low;
 638   fElectronsAbove10GeVMax = high;
 639   fElectronsAbove10GeVFlag = kTRUE;
 640 }
 641 //___________________________________________________________________________
 642 void AliEventTagCuts::SetNElectronRange(Int_t low, Int_t high) {
 643   //Sets the electron multiplicity range
 644   //and the corresponding flag to kTRUE if the cut is used.
 645   fElectronsMin = low;
 646   fElectronsMax = high;
 647   fElectronsFlag = kTRUE;
 648 }
 649 //___________________________________________________________________________
 650 void AliEventTagCuts::SetNMuonRange(Int_t low, Int_t high) {
 651   //Sets the muon multiplicity range
 652   //and the corresponding flag to kTRUE if the cut is used.
 653   fMuonsMin = low;
 654   fMuonsMax = high;
 655   fMuonsFlag = kTRUE;
 656 }
 657
 658 //___________________________________________________________________________
 659 void AliEventTagCuts::SetNPionRange(Int_t low, Int_t high) {
 660   //Sets the pion multiplicity range
 661   //and the corresponding flag to kTRUE if the cut is used.
 662   fPionsMin = low;
 663   fPionsMax = high;
 664   fPionsFlag = kTRUE;
 665 }
 666
 667 //___________________________________________________________________________
 668 void AliEventTagCuts::SetNKaonRange(Int_t low, Int_t high) {
 669   //Sets the kaon multiplicity range
 670   //and the corresponding flag to kTRUE if the cut is used.
 671   fKaonsMin = low;
 672   fKaonsMax = high;
 673   fKaonsFlag = kTRUE;
 674 }
 675
 676 //___________________________________________________________________________
 677 void AliEventTagCuts::SetNProtonRange(Int_t low, Int_t high) {
 678   //Sets the proton multiplicity range
 679   //and the corresponding flag to kTRUE if the cut is used.
 680   fProtonsMin = low;
 681   fProtonsMax = high;
 682   fProtonsFlag = kTRUE;
 683 }
 684
 685 //___________________________________________________________________________
 686 void AliEventTagCuts::SetNLambdaRange(Int_t low, Int_t high) {
 687   //Sets the lambda multiplicity range
 688   //and the corresponding flag to kTRUE if the cut is used.
 689   fLambdasMin = low;
 690   fLambdasMax = high;
 691   fLambdasFlag = kTRUE;
 692 }
 693 //___________________________________________________________________________
 694 void AliEventTagCuts::SetNPhotonRange(Int_t low, Int_t high) {
 695   //Sets the photon multiplicity range
 696   //and the corresponding flag to kTRUE if the cut is used.
 697   fPhotonsMin = low;
 698   fPhotonsMax = high;
 699   fPhotonFlag = kTRUE;
 700 }
 701 //___________________________________________________________________________
 702 void AliEventTagCuts::SetNPi0Range(Int_t low, Int_t high) {
 703   //Sets the pi0 multiplicity range
 704   //and the corresponding flag to kTRUE if the cut is used.
 705   fPi0sMin = low;
 706   fPi0sMax = high;
 707   fPi0sFlag = kTRUE;
 708 }
 709
 710 //___________________________________________________________________________
 711 void AliEventTagCuts::SetNNeutronRange(Int_t low, Int_t high) {
 712   //Sets the neutron multiplicity range
 713   //and the corresponding flag to kTRUE if the cut is used.
 714   fNeutronsMin = low;
 715   fNeutronsMax = high;
 716   fNeutronsFlag = kTRUE;
 717 }
 718
 719 //___________________________________________________________________________
 720 void AliEventTagCuts::SetNKaon0Range(Int_t low, Int_t high) {
 721   //Sets the K0s multiplicity range
 722   //and the corresponding flag to kTRUE if the cut is used.
 723   fKaon0sMin = low;
 724   fKaon0sMax = high;
 725   fKaon0sFlag = kTRUE;
 726 }
 727
 728 //___________________________________________________________________________
 729 void AliEventTagCuts::SetTotalPRange(Float_t low, Float_t high) {
 730   //Sets the total momentum range
 731   //and the corresponding flag to kTRUE if the cut is used.
 732   fTotalPMin = low;
 733   fTotalPMax = high;
 734   fTotalPFlag = kTRUE;
 735 }
 736
 737 //___________________________________________________________________________
 738 void AliEventTagCuts::SetMeanPtRange(Float_t low, Float_t high) {
 739   //Sets the mean Pt range
 740   //and the corresponding flag to kTRUE if the cut is used.
 741   fMeanPtMin = low;
 742   fMeanPtMax = high;
 743   fMeanPtFlag = kTRUE;
 744 }
 745
 746 //___________________________________________________________________________
 747 void AliEventTagCuts::SetTopPtMin(Float_t low) {
 748   //Sets the lower limit of the max Pt value
 749   //and the corresponding flag to kTRUE if the cut is used.
 750   fMaxPt = low;
 751   fMaxPtFlag = kTRUE;
 752 }
 753
 754 //___________________________________________________________________________
 755 void AliEventTagCuts::SetTotalNeutralPRange(Float_t low, Float_t high) {
 756   //Sets the total momentum of neutral particles range
 757   //and the corresponding flag to kTRUE if the cut is used.
 758   fTotalNeutralPMin =low ;
 759   fTotalNeutralPMax = high;
 760   fTotalNeutralPFlag = kTRUE;
 761 }
 762 //___________________________________________________________________________
 763 void AliEventTagCuts::SetMeanNeutralPtPRange(Float_t low, Float_t high) {
 764   //Sets the mean Pt of neutral particles range
 765   //and the corresponding flag to kTRUE if the cut is used.
 766   fMeanNeutralPtMin = low;
 767   fMeanNeutralPtMax = high;
 768   fMeanNeutralPtFlag = kTRUE;
 769 }
 770 //___________________________________________________________________________
 771 void AliEventTagCuts::SetTopNeutralPtMin(Float_t low) {
 772   //Sets the lower limit of the maximum Pt of neutral particles
 773   //and the corresponding flag to kTRUE if the cut is used.
 774   fMaxNeutralPt = low;
 775   fMaxNeutralPtFlag = kTRUE;
 776 }
 777
 778 //___________________________________________________________________________
 779 void AliEventTagCuts::SetEvPlaneAngleRange(Float_t low, Float_t high) {
 780   //Sets the event plane range
 781   //and the corresponding flag to kTRUE if the cut is used.
 782   fEventPlaneAngleMin = low;
 783   fEventPlaneAngleMax = high;
 784   fEventPlaneAngleFlag = kTRUE;
 785 }
 786
 787 //___________________________________________________________________________
 788 void AliEventTagCuts::SetHBTRadiiRange(Float_t low, Float_t high) {
 789   //Sets the HBT radii range
 790   //and the corresponding flag to kTRUE if the cut is used.
 791   fHBTRadiiMin = low;
 792   fHBTRadiiMax = high;
 793   fHBTRadiiFlag = kTRUE;
 794 }
 795
 796 //___________________________________________________________________________
 797 Bool_t AliEventTagCuts::IsAccepted(AliEventTag *EvTag) const {
 798   //Returns true if the event is accepted otherwise false.
 799   if(fVzFlag)
 800     if((EvTag->GetVertexZ() < fVzMin) || (EvTag->GetVertexZ() > fVzMax))
 801       return kFALSE;
 802
 803   if(fVyFlag)
 804     if((EvTag->GetVertexY() < fVyMin) || (EvTag->GetVertexY() > fVyMax))
 805       return kFALSE;
 806
 807   if(fVxFlag)
 808     if((EvTag->GetVertexX() < fVxMin) || (EvTag->GetVertexX() > fVxMax))
 809       return kFALSE;
 810
 811   if(fParticipantsFlag)
 812     if((EvTag->GetNumOfParticipants() < fParticipantsMin) || (EvTag->GetNumOfParticipants() > fParticipantMax))
 813       return kFALSE;
 814
 815   if(fImpactParamFlag)
 816     if((EvTag->GetImpactParameter() < fImpactParamMin) || (EvTag->GetImpactParameter() > fImpactParamMax))
 817       return kFALSE;
 818
 819   if(fPVFlag)
 820     if((EvTag->GetVertexFlag() != fPrimaryVertexFlag))
 821       return kFALSE;
 822
 823   if(fPVzErrorFlag)
 824     if((EvTag->GetVertexZError() < fPrimaryVertexZErrorMin) || (EvTag->GetVertexZError() > fPrimaryVertexZErrorMax))
 825       return kFALSE;
 826   if(fTriggerMaskFlag)
 827     if((EvTag->GetTriggerMask() != fTriggerMask))
 828       return kFALSE;
 829   if(fTriggerClusterFlag)
 830     if((EvTag->GetTriggerMask() != fTriggerMask))
 831       return kFALSE;
 832
 833   if(fZDCNeutron1EnergyFlag)
 834     if((EvTag->GetZDCNeutron1Energy() < fZDCNeutron1EnergyMin) || (EvTag->GetZDCNeutron1Energy() > fZDCNeutron1EnergyMax))
 835       return kFALSE;
 836
 837   if(fZDCProton1EnergyFlag)
 838     if((EvTag->GetZDCProton1Energy() < fZDCProton1EnergyMin) || (EvTag->GetZDCProton1Energy() > fZDCProton1EnergyMax))
 839       return kFALSE;
 840
 841   if(fZDCNeutron2EnergyFlag)
 842     if((EvTag->GetZDCNeutron2Energy() < fZDCNeutron2EnergyMin) || (EvTag->GetZDCNeutron2Energy() > fZDCNeutron2EnergyMax))
 843       return kFALSE;
 844
 845   if(fZDCProton2EnergyFlag)
 846     if((EvTag->GetZDCProton2Energy() < fZDCProton2EnergyMin) || (EvTag->GetZDCProton2Energy() > fZDCProton2EnergyMax))
 847       return kFALSE;
 848
 849   if(fZDCEMEnergyFlag)
 850     if((EvTag->GetZDCEMEnergy() < fZDCEMEnergyMin) || (EvTag->GetZDCEMEnergy() > fZDCEMEnergyMax))
 851       return kFALSE;
 852
 853   if(fT0VertexZFlag)
 854     if((EvTag->GetT0VertexZ() < fT0VertexZMin) || (EvTag->GetT0VertexZ() > fT0VertexZMax))
 855       return kFALSE;
 856
 857   if(fMultFlag)
 858     if((EvTag->GetNumOfTracks() < fMultMin) || (EvTag->GetNumOfTracks() > fMultMax))
 859       return kFALSE;
 860   if(fMultPosFlag)
 861     if((EvTag->GetNumOfPosTracks() < fMultPosMin) || (EvTag->GetNumOfPosTracks() > fMultPosMax))
 862       return kFALSE;
 863
 864   if(fMultNegFlag)
 865     if((EvTag->GetNumOfNegTracks() < fMultNegMin) || (EvTag->GetNumOfNegTracks() > fMultNegMax))
 866       return kFALSE;
 867
 868   if(fMultNeutrFlag)
 869     if((EvTag->GetNumOfNeutrTracks() < fMultNeutrMin) || (EvTag->GetNumOfNeutrTracks() > fMultNeutrMax))
 870       return kFALSE;
 871
 872   if(fV0sFlag)
 873     if((EvTag->GetNumOfV0s() < fV0sMin) || (EvTag->GetNumOfV0s() > fV0sMax))
 874       return kFALSE;
 875
 876   if(fCascadesFlag)
 877     if((EvTag->GetNumOfCascades() < fCascadesMin) || (EvTag->GetNumOfCascades() > fCascadesMax))
 878       return kFALSE;
 879
 880   if(fkinksFlag)
 881     if((EvTag->GetNumOfKinks() < fkinksMin) || (EvTag->GetNumOfKinks() > fkinksMax))
 882       return kFALSE;
 883
 884
 885   if(fPMDTracksFlag)
 886     if((EvTag->GetNumOfPMDTracks() < fPMDTracksMin) || (EvTag->GetNumOfPMDTracks() > fPMDTracksMax))
 887       return kFALSE;
 888   if(fFMDTracksFlag)
 889     if((EvTag->GetNumOfFMDTracks() < fFMDTracksMin) || (EvTag->GetNumOfFMDTracks() > fFMDTracksMax))
 890       return kFALSE;
 891   if(fPHOSClustersFlag)
 892     if((EvTag->GetNumOfPHOSClusters() < fPHOSClustersMin) || (EvTag->GetNumOfPHOSClusters() > fPHOSClustersMax))
 893       return kFALSE;
 894   if(fEMCALClustersFlag)
 895     if((EvTag->GetNumOfEMCALClusters() < fEMCALClustersMin) || (EvTag->GetNumOfEMCALClusters() > fEMCALClustersMax))
 896       return kFALSE;
 897   if(fJetCandidatesFlag)
 898     if((EvTag->GetNumOfJetCandidates() < fJetCandidatesMin) || (EvTag->GetNumOfJetCandidates() > fJetCandidatesMax))
 899       return kFALSE;
 900
 901
 902   if(fMaxJetEnergyFlag)
 903     if((EvTag->GetMaxJetEnergy() < fMaxJetEnergy))
 904       return kFALSE;
 905
 906   if(fNHardPhotonsCandidatesFlag)
 907     if((EvTag->GetNumOfHardPhotonsCandidates() < fNHardPhotonsCandidatesMin) || (EvTag->GetNumOfHardPhotonsCandidates() > fNHardPhotonsCandidatesMax))
 908       return kFALSE;
 909
 910   if(fMaxNeutralFlag)
 911     if((EvTag->GetMaxNeutralEnergy() < fMaxNeutralEnergy))
 912       return kFALSE;
 913
 914   if(fChargedAbove1GeVFlag)
 915     if((EvTag->GetNumOfChargedAbove1GeV() < fChargedAbove1GeVMin) || (EvTag->GetNumOfChargedAbove1GeV() > fChargedAbove1GeVMax))
 916       return kFALSE;
 917
 918   if(fChargedAbove3GeVFlag)
 919     if((EvTag->GetNumOfChargedAbove3GeV() < fChargedAbove3GeVMin) || (EvTag->GetNumOfChargedAbove3GeV() > fChargedAbove3GeVMax))
 920       return kFALSE;
 921
 922   if(fChargedAbove10GeVFlag)
 923     if((EvTag->GetNumOfChargedAbove10GeV() < fChargedAbove10GeVMin) || (EvTag->GetNumOfChargedAbove10GeV() > fChargedAbove10GeVMax))
 924       return kFALSE;
 925
 926   if(fMuonsAbove1GeVFlag)
 927     if((EvTag->GetNumOfMuonsAbove1GeV() < fMuonsAbove1GeVMin) || (EvTag->GetNumOfMuonsAbove1GeV() > fMuonsAbove1GeVMax))
 928       return kFALSE;
 929
 930   if(fMuonsAbove3GeVFlag)
 931     if((EvTag->GetNumOfMuonsAbove3GeV() < fMuonsAbove3GeVMin) || (EvTag->GetNumOfMuonsAbove3GeV() > fMuonsAbove3GeVMax))
 932       return kFALSE;
 933
 934   if(fMuonsAbove10GeVFlag)
 935     if((EvTag->GetNumOfMuonsAbove10GeV() < fMuonsAbove10GeVMin) || (EvTag->GetNumOfMuonsAbove10GeV() > fMuonsAbove10GeVMax))
 936       return kFALSE;
 937
 938   if(fElectronsAbove1GeVFlag)
 939     if((EvTag->GetNumOfElectronsAbove1GeV()  < fElectronsAbove1GeVMin) || (EvTag->GetNumOfElectronsAbove1GeV()  > fElectronsAbove1GeVMax))
 940       return kFALSE;
 941
 942   if(fElectronsAbove3GeVFlag)
 943     if((EvTag->GetNumOfElectronsAbove3GeV() < fElectronsAbove3GeVMin) || (EvTag->GetNumOfElectronsAbove3GeV() > fElectronsAbove3GeVMax))
 944       return kFALSE;
 945
 946   if(fElectronsAbove10GeVFlag)
 947     if((EvTag->GetNumOfElectronsAbove10GeV() < fElectronsAbove10GeVMin) || (EvTag->GetNumOfElectronsAbove10GeV() > fElectronsAbove10GeVMax))
 948       return kFALSE;
 949
 950   if(fElectronsFlag)
 951     if((EvTag->GetNumOfElectrons() < fElectronsMin) || (EvTag->GetNumOfElectrons() > fElectronsMax))
 952       return kFALSE;
 953
 954   if(fMuonsFlag)
 955     if((EvTag->GetNumOfMuons() < fMuonsMin) || (EvTag->GetNumOfMuons() > fMuonsMax))
 956       return kFALSE;
 957
 958   if(fPionsFlag)
 959     if((EvTag->GetNumOfPions() < fPionsMin) || (EvTag->GetNumOfPions() > fPionsMax))
 960       return kFALSE;
 961
 962   if(fKaonsFlag)
 963     if((EvTag->GetNumOfKaons() < fKaonsMin) || (EvTag->GetNumOfKaons() > fKaonsMax))
 964       return kFALSE;
 965
 966   if(fProtonsFlag)
 967     if((EvTag->GetNumOfProtons() < fProtonsMin) || (EvTag->GetNumOfProtons() > fProtonsMax))
 968       return kFALSE;
 969
 970   if(fLambdasFlag)
 971     if((EvTag->GetNumOfLambdas() < fLambdasMin) || (EvTag->GetNumOfLambdas() > fLambdasMax))
 972       return kFALSE;
 973
 974   if(fPhotonFlag)
 975     if((EvTag->GetNumOfPhotons() < fPhotonsMin) || (EvTag->GetNumOfPhotons() > fPhotonsMax))
 976       return kFALSE;
 977
 978   if(fPi0sFlag)
 979     if((EvTag->GetNumOfPi0s() < fPi0sMin) || (EvTag->GetNumOfPi0s() > fPi0sMax))
 980       return kFALSE;
 981
 982   if(fNeutronsFlag)
 983     if((EvTag->GetNumOfNeutrons() < fNeutronsMin) || (EvTag->GetNumOfNeutrons() > fNeutronsMax))
 984       return kFALSE;
 985
 986   if(fKaon0sFlag)
 987     if((EvTag->GetNumOfKaon0s() < fKaon0sMin) || (EvTag->GetNumOfKaon0s() > fKaon0sMax))
 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(fMaxPtFlag)
 999     if((EvTag->GetMaxPt() < fMaxPt))
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(fMaxNeutralPtFlag)
1011     if((EvTag->GetNeutralMaxPt() < fMaxNeutralPt))
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 }