Pluto track mother index fix

sim/transport/generators/CbmPlutoGenerator.cxx

@@ -9,27 +9,17 @@
 #include "CbmPlutoGenerator.h"
 
 #include "CbmFileUtils.h"
-
 #include "FairPrimaryGenerator.h"  // for FairPrimaryGenerator
-#include <Logger.h>
-
-#include "TArchiveFile.h"    // for TArchiveFile
+#include "PDataBase.h"             // for PDataBase
+#include "PParticle.h"             // for PParticle
+#include "RtypesCore.h"
 #include "TChain.h"          // for TChain
 #include "TClonesArray.h"    // for TClonesArray
 #include "TDatabasePDG.h"    // for TDatabasePDG
-#include "TFile.h"           // for TFile
 #include "TLorentzVector.h"  // for TLorentzVector
-#include "TTree.h"           // for TTree
 #include "TVector3.h"        // for TVector3
 
-#include <iosfwd>  // for ostream
-
-#include "PDataBase.h"    // for PDataBase
-#include "PParticle.h"    // for PParticle
-#include "PStaticData.h"  // for PStaticData
-
-//#include <stddef.h>                     // for NULL
-#include <iostream>  // for operator<<, basic_ostream, etc
+#include <Logger.h>
 
 #include <sys/stat.h>
 
@@ -108,6 +98,9 @@ Bool_t CbmPlutoGenerator::ReadEvent(FairPrimaryGenerator* primGen)
   // Get number of particle in TClonesrray
   Int_t nParts = fParticles->GetEntriesFast();
 
+  // define a dummy value
+  Int_t dummyPdg = 0;
+
   // Loop over particles in TClonesArray
   for (Int_t iPart = 0; iPart < nParts; iPart++) {
     PParticle* part = (PParticle*) fParticles->At(iPart);
@@ -118,8 +111,8 @@ Bool_t CbmPlutoGenerator::ReadEvent(FairPrimaryGenerator* primGen)
 
     // Check if particle type is known to database
     if (!found) {
-      LOG(warn) << "CbmPlutoGenerator: Unknown type " << part->ID() << ", skipping particle.";
-      continue;
+      LOG(warn) << "CbmPlutoGenerator: Unknown type " << part->ID() << ", setting PDG code to " << dummyPdg << ".";
+      pdgType = &dummyPdg;
     }
     LOG(info) << iPart << " Particle (geant " << part->ID() << " PDG " << *pdgType << " -> "
               << dataBase->GetParticle(*pdgType)->GetName();
@@ -153,7 +146,7 @@ Bool_t CbmPlutoGenerator::ReadEvent(FairPrimaryGenerator* primGen)
     Double_t vz     = vertex.z();
 
     Bool_t wanttracking = kTRUE;
-    if (idx > -1) wanttracking = kFALSE;  // only tracking for decay products
+    if (!found || idx > -1) wanttracking = kFALSE;  // only tracking for decay products that are known
     Int_t parent = parIdx;
     LOG(info) << "Add particle with parent at index " << parIdx << " and do tracking " << wanttracking;
 

sim/transport/generators/CbmPlutoGenerator.h

@@ -21,16 +21,13 @@
 #define FAIR_PLUTOGENERATOR_H
 
 #include "FairGenerator.h"  // for FairGenerator
-
+#include "PStaticData.h"    // for PStaticData
 #include "Rtypes.h"         // for Char_t, etc
 #include "TClonesArray.h"   // for TClonesArray
 
 #include <string>
 #include <vector>
 
-#include "PParticle.h"    // for PParticle
-#include "PStaticData.h"  // for PStaticData
-
 class FairPrimaryGenerator;
 
 class TChain;