AliAODEvent::GetHeader now return AliVHeader
[u/mrichter/AliRoot.git] / PWG / muon / AliUtilityMuonAncestor.cxx
index 968a10a..eea875b 100644 (file)
-/**************************************************************************\r
- * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *\r
- *                                                                        *\r
- * Author: The ALICE Off-line Project.                                    *\r
- * Contributors are mentioned in the code where appropriate.              *\r
- *                                                                        *\r
- * Permission to use, copy, modify and distribute this software and its   *\r
- * documentation strictly for non-commercial purposes is hereby granted   *\r
- * without fee, provided that the above copyright notice appears in all   *\r
- * copies and that both the copyright notice and this permission notice   *\r
- * appear in the supporting documentation. The authors make no claims     *\r
- * about the suitability of this software for any purpose. It is          *\r
- * provided "as is" without express or implied warranty.                  *\r
- **************************************************************************/\r
-\r
-/* $Id: AliUtilityMuonAncestor.cxx 47782 2011-02-24 18:37:31Z martinez $ */\r
-\r
-//-----------------------------------------------------------------------------\r
-/// \class AliUtilityMuonAncestor\r
-/// Static utilities to get the muon ancestor in MC\r
-///\r
-/// \author Diego Stocco\r
-//-----------------------------------------------------------------------------\r
-\r
-#include "AliUtilityMuonAncestor.h"\r
-\r
-// ROOT includes\r
-#include "TDatabasePDG.h"\r
-#include "TParticlePDG.h"\r
-#include "TPDGCode.h"\r
-#include "TMCProcess.h"\r
-#include "TMath.h"\r
-\r
-// STEER includes\r
-#include "AliMCEvent.h"\r
-#include "AliVParticle.h"\r
-#include "AliMCParticle.h"\r
-#include "AliAODMCParticle.h"\r
-#include "AliLog.h"\r
-\r
-// PWGmuon includes\r
-#include "AliAnalysisMuonUtility.h"\r
-\r
-/// \cond CLASSIMP\r
-ClassImp(AliUtilityMuonAncestor) // Class implementation in ROOT context\r
-/// \endcond\r
-\r
-//_________________________________________________________\r
-AliUtilityMuonAncestor::AliUtilityMuonAncestor() :\r
-TObject(),\r
-fPx(0.),\r
-fPy(0.),\r
-fPz(0.),\r
-fMask(0),\r
-fAncestor(-999)\r
-{\r
-  /// Default constructor\r
-}\r
-\r
-\r
-//_________________________________________________________\r
-AliUtilityMuonAncestor::~AliUtilityMuonAncestor()\r
-{\r
-  /// Default destructor\r
-}\r
-\r
-//_________________________________________________________\r
-AliUtilityMuonAncestor::AliUtilityMuonAncestor(const AliUtilityMuonAncestor& obj) :\r
-TObject(obj),\r
-fPx(obj.fPx),\r
-fPy(obj.fPy),\r
-fPz(obj.fPz),\r
-fMask(obj.fMask),\r
-fAncestor(obj.fAncestor)\r
-{\r
-  /// Copy constructor\r
-}\r
-\r
-//_________________________________________________________\r
-AliUtilityMuonAncestor& AliUtilityMuonAncestor::operator=(const AliUtilityMuonAncestor& obj)\r
-{\r
-  /// Copy operator\r
-  if ( this != &obj ) {\r
-    TObject::operator=(obj);\r
-    fPx = obj.fPx;\r
-    fPy = obj.fPy;\r
-    fPz = obj.fPz;\r
-    fMask = obj.fMask;\r
-    fAncestor = obj.fAncestor;\r
-  }\r
-  return *this;\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::BuildAncestor ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Build ancestor\r
-  \r
-  // If track is the same as the one in memory, do not re-build the track ancestor\r
-  if ( track->Px() == fPx && track->Py() == fPy && track->Pz() == fPz ) return kTRUE;\r
-  fPx = track->Px();\r
-  fPy = track->Py();\r
-  fPz = track->Pz();\r
-  fMask = 0;\r
-  fAncestor = -999;\r
-  if ( ! mcEvent ) return kFALSE;\r
-  \r
-  AliVParticle* mcParticle = 0x0;\r
-  \r
-  if ( AliAnalysisMuonUtility::IsMCTrack(track) ) mcParticle = track;\r
-  else {\r
-    Int_t trackLabel = track->GetLabel();\r
-    if ( trackLabel < 0 ) return kFALSE;\r
-    mcParticle = mcEvent->GetTrack(trackLabel);\r
-  }\r
-  \r
-  // Track is MC (or matches a MC track)\r
-  SETBIT(fMask,kIsID);\r
-  \r
-  Int_t recoPdg = mcParticle->PdgCode();\r
-  \r
-  // Track is not a muon\r
-  if ( TMath::Abs(recoPdg) == 13 ) SETBIT(fMask,kIsMuon);\r
-  \r
-  Int_t imother = AliAnalysisMuonUtility::GetMotherIndex(mcParticle);\r
-\r
-  while ( imother >= 0 ) {\r
-    AliVParticle* part = mcEvent->GetTrack(imother);\r
-    \r
-    Int_t absPdg = TMath::Abs(part->PdgCode());\r
-    \r
-    // This is a quark\r
-    if ( absPdg < 10 ) return kTRUE;\r
-    \r
-    fAncestor = imother;\r
-    Bool_t isPrimary = AliAnalysisMuonUtility::IsPrimary(part, mcEvent);\r
-    \r
-    if ( isPrimary ) {\r
-      Int_t mpdg = absPdg%100000;\r
-      if ( mpdg >= 100 && mpdg < 10000 ) {\r
-        Int_t flv  = Int_t ( mpdg / TMath::Power(10, Int_t(TMath::Log10(mpdg) )));\r
-        if ( flv < 4 ) SETBIT(fMask,kHasLightParent);\r
-        else if ( flv >= 6 ) continue;\r
-        else {\r
-          TParticlePDG* partPdg = TDatabasePDG::Instance()->GetParticle(part->PdgCode());\r
-          if ( partPdg && ! partPdg->AntiParticle() ) SETBIT(fMask,kHasQuarkoniumParent);\r
-          else if ( flv == 4 ) SETBIT(fMask,kHasCharmParent);\r
-          else SETBIT(fMask,kHasBeautyParent);\r
-        }\r
-      } // absPdg within 100 and 10000\r
-    } // is primary\r
-    else {\r
-      UInt_t mcProcess = AliAnalysisMuonUtility::GetMCProcess(part);\r
-      if ( mcProcess == kPHadronic ||\r
-           ( mcProcess == 0 && part->Zv() < -90. ) ) {\r
-        // The MC process is not well computed in the AODs of old MC productions\r
-        // In this case, declare the particle as "secondary" if it is produced inside the front absorber\r
-        SETBIT(fMask,kIsSecondary);\r
-//        return kTRUE;\r
-      }\r
-    } // is secondary\r
-    \r
-    imother = AliAnalysisMuonUtility::GetMotherIndex(part);\r
-    \r
-  } // loop on mothers\r
-  return kTRUE;\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::CheckAncestor ( AliVParticle* track, const AliMCEvent* mcEvent, Int_t ancestorPdg, Bool_t matchAbsPdg )\r
-{\r
-  /// Check ancestor\r
-  Int_t pdgCode = GetAncestorPdg(track, mcEvent);\r
-  if ( matchAbsPdg ) {\r
-    pdgCode = TMath::Abs(pdgCode);\r
-    ancestorPdg = TMath::Abs(ancestorPdg);\r
-  }\r
-  return ( pdgCode == ancestorPdg );\r
-}\r
-\r
-//_________________________________________________________\r
-Int_t AliUtilityMuonAncestor::GetAncestor ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Return ancestor (compute it if necessary)\r
-  BuildAncestor(track,mcEvent);\r
-  return fAncestor;\r
-}\r
-\r
-//_________________________________________________________\r
-Int_t AliUtilityMuonAncestor::GetAncestorPdg ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Return ancestor Pdg\r
-  if ( BuildAncestor(track,mcEvent) ) {\r
-    if ( fAncestor >= 0 ) return mcEvent->GetTrack(fAncestor)->PdgCode();\r
-  }\r
-  return 0;\r
-}\r
-\r
-//_________________________________________________________\r
-Long64_t AliUtilityMuonAncestor::GetMask ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Return mask\r
-  BuildAncestor(track,mcEvent);\r
-  return fMask;\r
-}\r
-\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsBeautyMu ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Muon from beauty decays\r
-  Long64_t mask = GetMask(track,mcEvent);\r
-  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kHasBeautyParent) & ! TESTBIT(mask,kHasLightParent) );\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsBJpsiMu ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Muon B->J/psi decays\r
-  if ( IsBeautyMu(track,mcEvent) ) {\r
-    Int_t imother = AliAnalysisMuonUtility::GetMotherIndex(track);\r
-    if ( imother >= 0 ) return ( mcEvent->GetTrack(imother)->PdgCode() == 443 );\r
-  }\r
-  return kFALSE;\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsCharmMu ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Muon from charm decays\r
-  Long64_t mask = GetMask(track,mcEvent);\r
-  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kHasCharmParent) && ! TESTBIT(mask,kHasBeautyParent) && ! TESTBIT(mask,kHasLightParent) );\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsDecayMu ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Muon from light hadron decays\r
-  Long64_t mask = GetMask(track,mcEvent);\r
-  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kHasLightParent) && ! TESTBIT(mask,kIsSecondary) );\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsHadron ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Reconstructed hadron\r
-  Long64_t mask = GetMask(track,mcEvent);\r
-  return ( TESTBIT(mask,kIsID) && ! TESTBIT(mask,kIsMuon) );\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsMuon ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Track is muon\r
-  Long64_t mask = GetMask(track,mcEvent);\r
-  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) );\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsQuarkoniumMu ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Mu from quarkonium decay\r
-  Long64_t mask = GetMask(track,mcEvent);\r
-  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kHasQuarkoniumParent) );\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsSecondaryMu ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Muon from secondary decays in absorber\r
-  Long64_t mask = GetMask(track,mcEvent);\r
-  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kIsSecondary) );\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsUnidentified ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Unidentified muon\r
-  Long64_t mask = GetMask(track,mcEvent);\r
-  return ( ! TESTBIT(mask,kIsID) );\r
-}\r
-\r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsWBosonMu ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Muon from W boson decays\r
-  return ( IsMuon(track,mcEvent) && CheckAncestor(track,mcEvent,24) );\r
-}\r
-          \r
-//_________________________________________________________\r
-Bool_t AliUtilityMuonAncestor::IsZBosonMu ( AliVParticle* track, const AliMCEvent* mcEvent )\r
-{\r
-  /// Muon from Z boson decays\r
-  return ( IsMuon(track,mcEvent) && CheckAncestor(track,mcEvent,kZ0) );\r
-}\r
+/**************************************************************************
+ * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
+ *                                                                        *
+ * Author: The ALICE Off-line Project.                                    *
+ * Contributors are mentioned in the code where appropriate.              *
+ *                                                                        *
+ * Permission to use, copy, modify and distribute this software and its   *
+ * documentation strictly for non-commercial purposes is hereby granted   *
+ * without fee, provided that the above copyright notice appears in all   *
+ * copies and that both the copyright notice and this permission notice   *
+ * appear in the supporting documentation. The authors make no claims     *
+ * about the suitability of this software for any purpose. It is          *
+ * provided "as is" without express or implied warranty.                  *
+ **************************************************************************/
+
+/* $Id: AliUtilityMuonAncestor.cxx 47782 2011-02-24 18:37:31Z martinez $ */
+
+//-----------------------------------------------------------------------------
+/// \class AliUtilityMuonAncestor
+/// Static utilities to get the muon ancestor in MC
+///
+/// \author Diego Stocco
+//-----------------------------------------------------------------------------
+
+#include "AliUtilityMuonAncestor.h"
+
+// ROOT includes
+#include "TDatabasePDG.h"
+#include "TParticlePDG.h"
+#include "TPDGCode.h"
+#include "TMCProcess.h"
+#include "TMath.h"
+
+// STEER includes
+#include "AliMCEvent.h"
+#include "AliVParticle.h"
+#include "AliMCParticle.h"
+#include "AliAODMCParticle.h"
+#include "AliLog.h"
+
+// PWGmuon includes
+#include "AliAnalysisMuonUtility.h"
+
+/// \cond CLASSIMP
+ClassImp(AliUtilityMuonAncestor) // Class implementation in ROOT context
+/// \endcond
+
+//_________________________________________________________
+AliUtilityMuonAncestor::AliUtilityMuonAncestor() :
+TObject(),
+fPx(0.),
+fPy(0.),
+fPz(0.),
+fMask(0),
+fAncestor(-999)
+{
+  /// Default constructor
+}
+
+
+//_________________________________________________________
+AliUtilityMuonAncestor::~AliUtilityMuonAncestor()
+{
+  /// Default destructor
+}
+
+//_________________________________________________________
+AliUtilityMuonAncestor::AliUtilityMuonAncestor(const AliUtilityMuonAncestor& obj) :
+TObject(obj),
+fPx(obj.fPx),
+fPy(obj.fPy),
+fPz(obj.fPz),
+fMask(obj.fMask),
+fAncestor(obj.fAncestor)
+{
+  /// Copy constructor
+}
+
+//_________________________________________________________
+AliUtilityMuonAncestor& AliUtilityMuonAncestor::operator=(const AliUtilityMuonAncestor& obj)
+{
+  /// Copy operator
+  if ( this != &obj ) {
+    TObject::operator=(obj);
+    fPx = obj.fPx;
+    fPy = obj.fPy;
+    fPz = obj.fPz;
+    fMask = obj.fMask;
+    fAncestor = obj.fAncestor;
+  }
+  return *this;
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::BuildAncestor ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Build ancestor
+  
+  // If track is the same as the one in memory, do not re-build the track ancestor
+  if ( track->Px() == fPx && track->Py() == fPy && track->Pz() == fPz ) return kTRUE;
+  fPx = track->Px();
+  fPy = track->Py();
+  fPz = track->Pz();
+  fMask = 0;
+  fAncestor = -999;
+  if ( ! mcEvent ) return kFALSE;
+  
+  AliVParticle* mcParticle = 0x0;
+  
+  if ( AliAnalysisMuonUtility::IsMCTrack(track) ) mcParticle = track;
+  else {
+    Int_t trackLabel = track->GetLabel();
+    if ( trackLabel < 0 ) return kFALSE;
+    mcParticle = mcEvent->GetTrack(trackLabel);
+  }
+  
+  // Track is MC (or matches a MC track)
+  SETBIT(fMask,kIsID);
+  
+  Int_t recoPdg = mcParticle->PdgCode();
+  
+  // Track is not a muon
+  if ( TMath::Abs(recoPdg) == 13 ) SETBIT(fMask,kIsMuon);
+  
+  Int_t imother = AliAnalysisMuonUtility::GetMotherIndex(mcParticle);
+
+  while ( imother >= 0 ) {
+    AliVParticle* part = mcEvent->GetTrack(imother);
+    
+    Int_t absPdg = TMath::Abs(part->PdgCode());
+    
+    // This is a quark
+    if ( absPdg < 10 ) return kTRUE;
+    
+    fAncestor = imother;
+    Bool_t isPrimary = AliAnalysisMuonUtility::IsPrimary(part, mcEvent);
+    
+    if ( isPrimary ) {
+      Int_t mpdg = absPdg%100000;
+      if ( mpdg >= 100 && mpdg < 10000 ) {
+        Int_t flv  = Int_t ( mpdg / TMath::Power(10, Int_t(TMath::Log10(mpdg) )));
+        if ( flv < 4 ) SETBIT(fMask,kHasLightParent);
+        else if ( flv >= 6 ) continue;
+        else {
+          TParticlePDG* partPdg = TDatabasePDG::Instance()->GetParticle(part->PdgCode());
+          if ( partPdg && ! partPdg->AntiParticle() ) SETBIT(fMask,kHasQuarkoniumParent);
+          else if ( flv == 4 ) SETBIT(fMask,kHasCharmParent);
+          else SETBIT(fMask,kHasBeautyParent);
+        }
+      } // absPdg within 100 and 10000
+    } // is primary
+    else {
+      UInt_t mcProcess = AliAnalysisMuonUtility::GetMCProcess(part);
+      if ( mcProcess == kPHadronic ||
+           ( mcProcess == 0 && part->Zv() < -90. ) ) {
+        // The MC process is not well computed in the AODs of old MC productions
+        // In this case, declare the particle as "secondary" if it is produced inside the front absorber
+        SETBIT(fMask,kIsSecondary);
+//        return kTRUE;
+      }
+    } // is secondary
+    
+    imother = AliAnalysisMuonUtility::GetMotherIndex(part);
+    
+  } // loop on mothers
+  return kTRUE;
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::CheckAncestor ( AliVParticle* track, const AliMCEvent* mcEvent, Int_t ancestorPdg, Bool_t matchAbsPdg )
+{
+  /// Check ancestor
+  Int_t pdgCode = GetAncestorPdg(track, mcEvent);
+  if ( matchAbsPdg ) {
+    pdgCode = TMath::Abs(pdgCode);
+    ancestorPdg = TMath::Abs(ancestorPdg);
+  }
+  return ( pdgCode == ancestorPdg );
+}
+
+//_________________________________________________________
+Int_t AliUtilityMuonAncestor::GetAncestor ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Return ancestor (compute it if necessary)
+  BuildAncestor(track,mcEvent);
+  return fAncestor;
+}
+
+//_________________________________________________________
+Int_t AliUtilityMuonAncestor::GetAncestorPdg ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Return ancestor Pdg
+  if ( BuildAncestor(track,mcEvent) ) {
+    if ( fAncestor >= 0 ) return mcEvent->GetTrack(fAncestor)->PdgCode();
+  }
+  return 0;
+}
+
+//_________________________________________________________
+Long64_t AliUtilityMuonAncestor::GetMask ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Return mask
+  BuildAncestor(track,mcEvent);
+  return fMask;
+}
+
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsBeautyMu ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Muon from beauty decays
+  Long64_t mask = GetMask(track,mcEvent);
+  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kHasBeautyParent) & ! TESTBIT(mask,kHasLightParent) );
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsBJpsiMu ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Muon B->J/psi decays
+  if ( IsBeautyMu(track,mcEvent) ) {
+    Int_t imother = AliAnalysisMuonUtility::GetMotherIndex(track);
+    if ( imother >= 0 ) return ( mcEvent->GetTrack(imother)->PdgCode() == 443 );
+  }
+  return kFALSE;
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsCharmMu ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Muon from charm decays
+  Long64_t mask = GetMask(track,mcEvent);
+  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kHasCharmParent) && ! TESTBIT(mask,kHasBeautyParent) && ! TESTBIT(mask,kHasLightParent) );
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsDecayMu ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Muon from light hadron decays
+  Long64_t mask = GetMask(track,mcEvent);
+  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kHasLightParent) && ! TESTBIT(mask,kIsSecondary) );
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsHadron ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Reconstructed hadron
+  Long64_t mask = GetMask(track,mcEvent);
+  return ( TESTBIT(mask,kIsID) && ! TESTBIT(mask,kIsMuon) );
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsMuon ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Track is muon
+  Long64_t mask = GetMask(track,mcEvent);
+  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) );
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsQuarkoniumMu ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Mu from quarkonium decay
+  Long64_t mask = GetMask(track,mcEvent);
+  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kHasQuarkoniumParent) );
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsSecondaryMu ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Muon from secondary decays in absorber
+  Long64_t mask = GetMask(track,mcEvent);
+  return ( TESTBIT(mask,kIsID) && TESTBIT(mask,kIsMuon) && TESTBIT(mask,kIsSecondary) );
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsUnidentified ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Unidentified muon
+  Long64_t mask = GetMask(track,mcEvent);
+  return ( ! TESTBIT(mask,kIsID) );
+}
+
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsWBosonMu ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Muon from W boson decays
+  return ( IsMuon(track,mcEvent) && CheckAncestor(track,mcEvent,24) );
+}
+          
+//_________________________________________________________
+Bool_t AliUtilityMuonAncestor::IsZBosonMu ( AliVParticle* track, const AliMCEvent* mcEvent )
+{
+  /// Muon from Z boson decays
+  return ( IsMuon(track,mcEvent) && CheckAncestor(track,mcEvent,kZ0) );
+}