Modified PrimaryGeneratorAction to read both PID and energy from text dump file

run1.mac

@@ -10,6 +10,6 @@
 #/process/list
 #/tracking/verbose 1
 
-/run/beamOn 4
+/run/beamOn 10
 
 #/run/particle/dumpList

src/PrimaryGeneratorAction.cc

@@ -320,8 +320,7 @@ void PrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent)
 
 
   // FASERnu text read:
-  float tmp1, tmp2, tmp3, tmp4;
-  int i1, i2, i3;
+  float tmp1, tmp2;
   bool start = false;
   string line;
   while ( getline(*source,line) ) {
@@ -342,11 +341,15 @@ void PrimaryGeneratorAction::GeneratePrimaries(G4Event* anEvent)
       start = true;
     }*/
     else if(start) {
-      sscanf(line.data(),"%f",&tmp1);
-      fParticleGun->SetParticleDefinition(fElectron); 
+      sscanf(line.data(),"%f %f",&tmp1,&tmp2);
+      if(tmp1==11) fParticleGun->SetParticleDefinition(fElectron);
+      else if(tmp1==-11) fParticleGun->SetParticleDefinition(fPositron);
+      else if(tmp1==13) fParticleGun->SetParticleDefinition(fMuonM);
+      else fParticleGun->SetParticleDefinition(fMuonP);
+
       //printf("Xin: %f\n",tmp1);
 
-      fParticleGun->SetParticleEnergy(tmp1*GeV);
+      fParticleGun->SetParticleEnergy(tmp2*GeV);
       fParticleGun->SetParticleMomentumDirection(G4ThreeVector(0,0,1));
       fParticleGun->SetParticlePosition(G4ThreeVector(0,0,-worldZHalfLength));
       fParticleGun->GeneratePrimaryVertex(anEvent);

src/StackingAction.cc

@@ -47,29 +47,29 @@ int flag = 0;
 G4ClassificationOfNewTrack
 StackingAction::ClassifyNewTrack(const G4Track* track)
 {
-//  std::cout << "Calling ClassifyNewTrack" << std::endl;
-  G4cout << "flag = " << flag << G4endl;
-  G4cout << "Tracks: Total: " << stackManager->GetNTotalTrack() << " Urgent: " << stackManager->GetNUrgentTrack() << " Waiting: " << stackManager->GetNWaitingTrack() << " Postponed: " << stackManager->GetNPostponedTrack() << G4endl;
-  G4cout << "TrackID = " << track->GetTrackID() << " ParentID = " << track->GetParentID() << " Particle = " << track->GetParticleDefinition()->GetParticleName() << G4endl;
+
+//  G4cout << "flag = " << flag << G4endl;
+//  G4cout << "Tracks: Total: " << stackManager->GetNTotalTrack() << " Urgent: " << stackManager->GetNUrgentTrack() << " Waiting: " << stackManager->GetNWaitingTrack() << " Postponed: " << stackManager->GetNPostponedTrack() << G4endl;
+//  G4cout << "TrackID = " << track->GetTrackID() << " ParentID = " << track->GetParentID() << " Particle = " << track->GetParticleDefinition()->GetParticleName() << G4endl;
   //keep primary particle
   if (flag == 0) return fPostpone;
 
   //kill secondary neutrino
-//  if (track->GetDefinition() == G4NeutrinoE::NeutrinoE()) return fKill;
+
   else return fUrgent;
 }
 void StackingAction::NewStage()
 {
-G4cout << "Calling NewStage()" << G4endl;
-  G4cout << "Tracks: Total: " << stackManager->GetNTotalTrack() << " Urgent: " << stackManager->GetNUrgentTrack() << " Waiting: " << stackManager->GetNWaitingTrack() << " Postponed: " << stackManager->GetNPostponedTrack() << G4endl;
+//G4cout << "Calling NewStage()" << G4endl;
+//  G4cout << "Tracks: Total: " << stackManager->GetNTotalTrack() << " Urgent: " << stackManager->GetNUrgentTrack() << " Waiting: " << stackManager->GetNWaitingTrack() << " Postponed: " << stackManager->GetNPostponedTrack() << G4endl;
 } 
 void StackingAction::PrepareNewEvent()
 {
 G4cout << "flag = " << flag << G4endl;
-G4cout << "Calling PrepareNewEvent()" << G4endl;
-  G4cout << "Tracks: Total: " << stackManager->GetNTotalTrack() << " Urgent: " << stackManager->GetNUrgentTrack() << " Waiting: " << stackManager->GetNWaitingTrack() << " Postponed: " << stackManager->GetNPostponedTrack() << G4endl;
+//G4cout << "Calling PrepareNewEvent()" << G4endl;
+//  G4cout << "Tracks: Total: " << stackManager->GetNTotalTrack() << " Urgent: " << stackManager->GetNUrgentTrack() << " Waiting: " << stackManager->GetNWaitingTrack() << " Postponed: " << stackManager->GetNPostponedTrack() << G4endl;
 if(stackManager->GetNPostponedTrack()>0) flag = 1;
 else flag = 0;
-G4cout << "flag = " << flag << G4endl;
+//G4cout << "flag = " << flag << G4endl;
 }
 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......