Charm and Beauty definitions updated
authorauras <auras@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 8 Feb 2013 18:13:53 +0000 (18:13 +0000)
committerauras <auras@f7af4fe6-9843-0410-8265-dc069ae4e863>
Fri, 8 Feb 2013 18:13:53 +0000 (18:13 +0000)
MFT/AliMuonForwardTrack.cxx
MFT/AliMuonForwardTrack.h

index 04881e1..11d2e90 100644 (file)
@@ -497,7 +497,7 @@ Bool_t AliMuonForwardTrack::IsFromResonance() {
 
 //====================================================================================================================================================
 
-Bool_t AliMuonForwardTrack::IsFromCharm() {
+Bool_t AliMuonForwardTrack::IsDirectCharm() {
 
   Bool_t result = kFALSE;
 
@@ -511,7 +511,7 @@ Bool_t AliMuonForwardTrack::IsFromCharm() {
 
 //====================================================================================================================================================
 
-Bool_t AliMuonForwardTrack::IsFromBeauty() {
+Bool_t AliMuonForwardTrack::IsDirectBeauty() {
 
   Bool_t result = kFALSE;
 
@@ -525,6 +525,20 @@ Bool_t AliMuonForwardTrack::IsFromBeauty() {
 
 //====================================================================================================================================================
 
+Bool_t AliMuonForwardTrack::IsChainBeauty() {
+
+  Bool_t result = kFALSE;
+
+  if (IsPDGCharm(GetParentPDGCode(0)) && IsPDGBeauty(GetParentPDGCode(1))) result = kTRUE;
+  
+  if (result) AliDebug(1, Form("Muon comes from a charmed hadron %d which comes from a beauty hadron %d", GetParentPDGCode(0), GetParentPDGCode(1)));
+  
+  return result; 
+  
+}
+
+//====================================================================================================================================================
+
 Bool_t AliMuonForwardTrack::IsPDGCharm(Int_t pdg) {
 
   Bool_t result = kFALSE;
@@ -593,7 +607,7 @@ Bool_t AliMuonForwardTrack::IsFromBackground() {
 
   Bool_t result = kFALSE;
 
-  if (!IsFromResonance() && !IsFromCharm() && !IsFromBeauty()) result = kTRUE;
+  if (!IsFromResonance() && !IsDirectCharm() && !IsDirectBeauty() && !IsChainBeauty()) result = kTRUE;
 
   if (result) AliDebug(1, Form("Muon comes from a background source %d", GetParentPDGCode(0)));
 
index 2d1c523..5706b08 100644 (file)
@@ -98,8 +98,11 @@ public:
   Int_t GetTrackMCId() { return fTrackMCId; }
   
   Bool_t IsFromResonance();
-  Bool_t IsFromCharm();
-  Bool_t IsFromBeauty();
+  Bool_t IsDirectCharm();
+  Bool_t IsDirectBeauty();
+  Bool_t IsChainBeauty();
+  Bool_t IsFromCharm()  { return IsDirectCharm(); }
+  Bool_t IsFromBeauty() { return IsDirectBeauty() || IsChainBeauty(); }
   Bool_t IsPDGCharm(Int_t pdg);
   Bool_t IsPDGBeauty(Int_t pdg);
   Bool_t IsFromBackground();