Changeset 84a097e in git for readers/DelphesPythia8.cpp

Timestamp:

Dec 2, 2019, 6:55:22 PM (5 years ago)

Author:

Michele Selvaggi <michele.selvaggi@…>

Branches:

Timing

Children:

2ad823e, 6fc566b

Parents:

c614dd7

Message:

first iteration of adding timing in Delphes

File:

• readers/DelphesPythia8.cpp (modified) (1 diff)

readers/DelphesPythia8.cpp

@@ -122 +122 @@
 
     pdgParticle = pdg->GetParticle(pid);
-    candidate->Charge = pdgParticle ? Int_t(pdgParticle->Charge() / 3.0) : -999;
+
+    if(pdgParticle)
+    {
+      candidate->Charge = int(pdgParticle->Charge()/3.0);
+    }
+    else if(abs(pid) == 1000612) candidate->Charge = TMath::Sign(1, pid);
+    else if(abs(pid) == 1000622) candidate->Charge = 0;
+    else if(abs(pid) == 1000632) candidate->Charge = TMath::Sign(1, pid);
+    else if(abs(pid) == 1000642) candidate->Charge = 0;
+    else if(abs(pid) == 1000652) candidate->Charge = TMath::Sign(1, pid);
+    else if(abs(pid) == 1006113) candidate->Charge = 0;
+    else if(abs(pid) == 1006211) candidate->Charge = TMath::Sign(1, pid);
+    else if(abs(pid) == 1006213) candidate->Charge = TMath::Sign(1, pid);
+    else if(abs(pid) == 1006223) candidate->Charge = TMath::Sign(2, pid);
+    else if(abs(pid) == 1006311) candidate->Charge = 0;
+    else if(abs(pid) == 1006313) candidate->Charge = 0;
+    else if(abs(pid) == 1006321) candidate->Charge = TMath::Sign(1, pid);
+    else if(abs(pid) == 1006323) candidate->Charge = TMath::Sign(1, pid);
+    else if(abs(pid) == 1006333) candidate->Charge = 0;
+    else {
+      candidate->Charge = -999;
+      cout << "[WARNING] Unrecognised particle. PID: " << pid << endl;
+    }
+
     candidate->Mass = mass;