]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - STEER/AliEventTagCuts.cxx
POI's and RP's for LeeYang Zeroes eventplane
[u/mrichter/AliRoot.git] / STEER / AliEventTagCuts.cxx
index c446a9a2ec2c7b291e74dc2631c6ca3e3b2cfbce..fb666827fdaf01b723f1d0771846c5ba50ffe58d 100644 (file)
@@ -117,6 +117,8 @@ AliEventTagCuts::AliEventTagCuts() :
   fNElectronsAbove10GeVFlag(kFALSE),
   fNElectronsMin(-1), fNElectronsMax(100000),
   fNElectronsFlag(kFALSE),
+  fNFWMuonsMin(-1), fNFWMuonsMax(100000),
+  fNFWMuonsFlag(kFALSE),
   fNMuonsMin(-1), fNMuonsMax(100000),
   fNMuonsFlag(kFALSE),
   fNPionsMin(-1), fNPionsMax(100000),
@@ -209,6 +211,7 @@ void AliEventTagCuts::Reset() {
   fNElectronsAbove3GeVFlag = kFALSE;
   fNElectronsAbove10GeVFlag = kFALSE;
   fNElectronsFlag = kFALSE;
+  fNFWMuonsFlag = kFALSE;
   fNMuonsFlag = kFALSE;
   fNPionsFlag = kFALSE;
   fNKaonsFlag = kFALSE;
@@ -271,6 +274,7 @@ void AliEventTagCuts::Reset() {
   fNElectronsAbove3GeVMin = -1; fNElectronsAbove3GeVMax = 100000;
   fNElectronsAbove10GeVMin = -1; fNElectronsAbove10GeVMax = 100000;
   fNElectronsMin = -1; fNElectronsMax = 100000;
+  fNFWMuonsMin = -1; fNFWMuonsMax = 100000;
   fNMuonsMin = -1; fNMuonsMax = 100000;
   fNPionsMin = -1; fNPionsMax = 100000;
   fNKaonsMin = -1; fNKaonsMax = 100000;
@@ -645,6 +649,15 @@ void AliEventTagCuts::SetNElectronRange(Int_t low, Int_t high) {
   fNElectronsMax = high;
   fNElectronsFlag = kTRUE;
 }
+//___________________________________________________________________________
+void AliEventTagCuts::SetNFWMuonRange(Int_t low, Int_t high) {
+  //Sets the forward muon multiplicity range
+  //and the corresponding flag to kTRUE if the cut is used.
+  fNFWMuonsMin = low;
+  fNFWMuonsMax = high;
+  fNFWMuonsFlag = kTRUE;
+}
+
 //___________________________________________________________________________
 void AliEventTagCuts::SetNMuonRange(Int_t low, Int_t high) {
   //Sets the muon multiplicity range
@@ -951,6 +964,10 @@ Bool_t AliEventTagCuts::IsAccepted(AliEventTag *EvTag) const {
     if((EvTag->GetNumOfElectrons() < fNElectronsMin) || (EvTag->GetNumOfElectrons() > fNElectronsMax))
       return kFALSE; 
   
+  if(fNFWMuonsFlag)
+    if((EvTag->GetNumOfFWMuons() < fNFWMuonsMin) || (EvTag->GetNumOfFWMuons() > fNFWMuonsMax))
+      return kFALSE; 
+  
   if(fNMuonsFlag)
     if((EvTag->GetNumOfMuons() < fNMuonsMin) || (EvTag->GetNumOfMuons() > fNMuonsMax))
       return kFALSE;