Changeset a09b75f in git for modules/LLPFilter.cc

Timestamp:

Dec 9, 2021, 10:38:13 PM (3 years ago)

Author:

christinaw97 <christina.wang@…>

Children:

c61b5ce

Parents:

29b722a

Message:

cleaned modules

File:

• modules/LLPFilter.cc (modified) (4 diffs)

modules/LLPFilter.cc

@@ -151 +151 @@
     if (fDaughterNumber > 0)
     {
-      if (candidate->D2-candidate->D1 != fDaughterNumber) continue;//require 3 daughters
+      if (candidate->D2-candidate->D1 != fDaughterNumber) continue;//require at least fDaughterNumber daughters
 
     }
@@ -160 +160 @@
     candidate->Eem = 0.0;
     candidate->Ehad = 0.0;
-
     fItParticleInputArray->Reset();
 
@@ -166 +165 @@
     {
 
-      // check mother is LLP
-      //check ID is 123456 or neutrinos or charged leptons
       daughterPdg = daughter->PID;
       if (daughter->Status != 1)continue;
-
       if (daughter->IsPU)continue;
-
       if (abs(daughterPdg)==12 || abs(daughterPdg)==14 || abs(daughterPdg)==16 || abs(daughterPdg)==13)continue; // ignore neutrinos and muons
-
-
       if (abs(daughterPdg) > 1000000) continue;//ignore BSM particles
+      
       const TLorentzVector &daughterProdPosition = daughter->Position;
       const TLorentzVector &daughterMomentum = daughter->Momentum;
-
       const TLorentzVector distance = daughterProdPosition - candidateDecayPosition;
 
-      if (sqrt(pow(distance.X(), 2) + pow(distance.Y(), 2) + pow(distance.Z(), 2))>100) continue; //if vertices are close, then matched, this 1 m only works for those decay in muon system, would be too large for tracker related
+      // look for mother until find LLP or reach the top of the tree
       tempCandidate = daughter;
       while(tempCandidate->M1 != -1 && tempCandidate->M1 != index)
       {
         tempCandidate  = static_cast<Candidate *>(fParticleInputArray->At(tempCandidate->M1));
-        // if (tempCandidate->PID == pdgCode)break;
-
       }
       if (tempCandidate->M1 == -1) continue;
@@ -194 +185 @@
       if (abs(daughterPdg)==11 || abs(daughterPdg)==22 || abs(daughterPdg)==111)candidate->Eem += daughterMomentum.E();
       else candidate->Ehad += daughterMomentum.E();
-
-
     }