Methods added to distinguish muons from charm and beauty
authorauras <auras@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 30 Sep 2012 22:01:27 +0000 (22:01 +0000)
committerauras <auras@f7af4fe6-9843-0410-8265-dc069ae4e863>
Sun, 30 Sep 2012 22:01:27 +0000 (22:01 +0000)
MFT/AliMuonForwardTrack.cxx
MFT/AliMuonForwardTrack.h
MFT/runReconstruction.C

index 16c3e9b..b12130c 100644 (file)
@@ -421,3 +421,106 @@ Double_t AliMuonForwardTrack::GetOffset(Double_t x, Double_t y, Double_t z) {
 
 //====================================================================================================================================================
 
+Bool_t AliMuonForwardTrack::IsFromResonance() {
+
+  Bool_t result = kFALSE;
+
+  if ( GetParentPDGCode(0) ==    113 ||
+       GetParentPDGCode(0) ==    221 ||
+       GetParentPDGCode(0) ==    223 ||
+       GetParentPDGCode(0) ==    331 ||
+       GetParentPDGCode(0) ==    333 ||
+       GetParentPDGCode(0) ==    443 ||
+       GetParentPDGCode(0) == 100443 ||
+       GetParentPDGCode(0) ==    553 ||
+       GetParentPDGCode(0) == 100553 ) result = kTRUE;
+  
+  if (result) AliDebug(1, Form("Muon comes from a resonance %d", GetParentPDGCode(0)));
+  
+  return result; 
+  
+}
+
+//====================================================================================================================================================
+
+Bool_t AliMuonForwardTrack::IsFromCharm() {
+
+  Bool_t result = kFALSE;
+
+  if ( GetParentPDGCode(0) ==   411 ||
+       GetParentPDGCode(0) ==   421 ||
+       GetParentPDGCode(0) == 10411 ||
+       GetParentPDGCode(0) == 10421 ||
+       GetParentPDGCode(0) ==   413 ||
+       GetParentPDGCode(0) ==   423 ||
+       GetParentPDGCode(0) == 10413 ||
+       GetParentPDGCode(0) == 10423 ||
+       GetParentPDGCode(0) == 20413 ||
+       GetParentPDGCode(0) == 20423 ||
+       GetParentPDGCode(0) ==   415 ||
+       GetParentPDGCode(0) ==   425 ||
+       GetParentPDGCode(0) ==   431 ||
+       GetParentPDGCode(0) == 10431 ||
+       GetParentPDGCode(0) ==   433 ||
+       GetParentPDGCode(0) == 10433 ||
+       GetParentPDGCode(0) == 20433 ||
+       GetParentPDGCode(0) ==   435 ) result = kTRUE;
+  
+  if (result) AliDebug(1, Form("Muon comes from a charmed hadron %d", GetParentPDGCode(0)));
+  
+  return result; 
+  
+}
+
+//====================================================================================================================================================
+
+Bool_t AliMuonForwardTrack::IsFromBeauty() {
+
+  Bool_t result = kFALSE;
+
+  if ( GetParentPDGCode(0) ==   511 ||
+       GetParentPDGCode(0) ==   521 ||
+       GetParentPDGCode(0) == 10511 ||
+       GetParentPDGCode(0) == 10521 ||
+       GetParentPDGCode(0) ==   513 ||
+       GetParentPDGCode(0) ==   523 ||
+       GetParentPDGCode(0) == 10513 ||
+       GetParentPDGCode(0) == 10523 ||
+       GetParentPDGCode(0) == 20513 ||
+       GetParentPDGCode(0) == 20523 ||
+       GetParentPDGCode(0) ==   515 ||
+       GetParentPDGCode(0) ==   525 ||
+       GetParentPDGCode(0) ==   531 ||
+       GetParentPDGCode(0) == 10531 ||
+       GetParentPDGCode(0) ==   533 ||
+       GetParentPDGCode(0) == 10533 ||
+       GetParentPDGCode(0) == 20533 ||
+       GetParentPDGCode(0) ==   535 ||
+       GetParentPDGCode(0) ==   541 ||
+       GetParentPDGCode(0) == 10541 ||
+       GetParentPDGCode(0) ==   543 ||
+       GetParentPDGCode(0) == 10543 ||
+       GetParentPDGCode(0) == 20543 ||
+       GetParentPDGCode(0) ==   545 ) result = kTRUE;
+  
+  if (result) AliDebug(1, Form("Muon comes from a beauty hadron %d", GetParentPDGCode(0)));
+  
+  return result; 
+  
+}
+
+//====================================================================================================================================================
+
+Bool_t AliMuonForwardTrack::IsFromBackground() {
+
+  Bool_t result = kFALSE;
+
+  if (!IsFromResonance() && !IsFromCharm() && !IsFromBeauty()) result = kTRUE;
+
+  if (result) AliDebug(1, Form("Muon comes from a background source %d", GetParentPDGCode(0)));
+
+  return result;
+
+}
+
+//====================================================================================================================================================
index d434406..3bdc152 100644 (file)
@@ -83,6 +83,11 @@ public:
   void SetTrackMCId(Int_t id) { fTrackMCId = id; }
   Int_t GetTrackMCId() { return fTrackMCId; }
   
+  Bool_t IsFromResonance();
+  Bool_t IsFromCharm();
+  Bool_t IsFromBeauty();
+  Bool_t IsFromBackground();
+
 protected:
 
   static const Int_t fNMaxPlanes = AliMFTConstants::fNMaxPlanes;        // max number of MFT planes
index c005af7..7bf7f3d 100644 (file)
@@ -39,8 +39,8 @@ void runReconstruction(Int_t seed, const Char_t *recOptions) {
   reco->SetRunReconstruction("MUON MFT");\r
   reco->SetRunLocalReconstruction("MUON MFT");\r
   reco->SetOption("MUON MFT",recOptions);\r
-  reco->SetRunQA("DetectorList:ActionList");\r
-  reco->SetQAWriteExpert(AliQAv1::kMUON);\r
+  //  reco->SetRunQA("DetectorList:ActionList");\r
+  //  reco->SetQAWriteExpert(AliQAv1::kMUON);\r
 \r
   reco->SetWriteESDfriend(kFALSE);\r
   reco->SetStopOnError(kFALSE);\r