Changeset 1134 in svn for trunk/examples/DelphesCMSFWLite.cpp

Timestamp:

Jun 5, 2013, 11:43:57 PM (11 years ago)

Author:

Pavel Demin

Message:

fix particles status selection in DelphesCMSFWLite.cpp

File:

• trunk/examples/DelphesCMSFWLite.cpp (modified) (1 diff)

trunk/examples/DelphesCMSFWLite.cpp

@@ -60 +60 @@
     const reco::GenParticle &particle = (*handleParticle)[i];
 
-     pid = particle.pdgId();
-     status = particle.status();
-     px = particle.px(), py = particle.py(), pz = particle.pz(), e = particle.energy();
-     x = particle.vx(), y = particle.vy(), z = particle.vz();
-
-    if(status == 1 || status == 2)
-    {
-      candidate = factory->NewCandidate();
-
-      candidate->PID = pid;
-      pdgCode = TMath::Abs(candidate->PID);
-
-      candidate->Status = status;
-
-      pdgParticle = pdg->GetParticle(pid);
-      candidate->Charge = pdgParticle ? Int_t(pdgParticle->Charge()/3.0) : -999;
-      candidate->Mass = pdgParticle ? pdgParticle->Mass() : -999.9;
-
-      candidate->Momentum.SetPxPyPzE(px, py, pz, e);
-
-      candidate->Position.SetXYZT(x, y, z, 0.0);
-
-      allParticleOutputArray->Add(candidate);
-
-      if(!pdgParticle) return;
-
-      if(status == 1)
-      {
-        stableParticleOutputArray->Add(candidate);
-      }
-      else if(pdgCode <= 5 || pdgCode == 21 || pdgCode == 15)
-      {
-        partonOutputArray->Add(candidate);
-      }
+    pid = particle.pdgId();
+    status = particle.status();
+    px = particle.px(); py = particle.py(); pz = particle.pz(); e = particle.energy();
+    x = particle.vx(); y = particle.vy(); z = particle.vz();
+
+    candidate = factory->NewCandidate();
+
+    candidate->PID = pid;
+    pdgCode = TMath::Abs(candidate->PID);
+
+    candidate->Status = status;
+
+    pdgParticle = pdg->GetParticle(pid);
+    candidate->Charge = pdgParticle ? Int_t(pdgParticle->Charge()/3.0) : -999;
+    candidate->Mass = pdgParticle ? pdgParticle->Mass() : -999.9;
+
+    candidate->Momentum.SetPxPyPzE(px, py, pz, e);
+
+    candidate->Position.SetXYZT(x, y, z, 0.0);
+
+    allParticleOutputArray->Add(candidate);
+
+    if(!pdgParticle) return;
+
+    if(status == 1)
+    {
+      stableParticleOutputArray->Add(candidate);
+    }
+    else if(pdgCode <= 5 || pdgCode == 21 || pdgCode == 15)
+    {
+      partonOutputArray->Add(candidate);
     }
   }