In decays with AliDecayer put long-lived particles undecayed on the stack.
authormorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Jul 2001 09:10:23 +0000 (09:10 +0000)
committermorsch <morsch@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 19 Jul 2001 09:10:23 +0000 (09:10 +0000)
TGeant3/AliGeant3.cxx

index 547f58b..2c8bc63 100644 (file)
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.17  2001/06/15 09:31:23  morsch
+In gudcay: write only first generation decay products to stack to respect the possibility of secondary, tertiary, ... vertices during tracking.
+
 Revision 1.16  2001/05/16 14:57:23  alibrary
 New files for folders and Stack
 
@@ -434,8 +437,8 @@ void gudcay()
 //
     
     TGeant3* geant3=(TGeant3*) gMC;
-  // set decay table
-  gMC->Decayer()->ForceDecay();
+    // set decay table
+    gMC->Decayer()->ForceDecay();
 
 // Initialize 4-momentum vector    
     Int_t ipart = geant3->Gckine()->ipart;
@@ -459,20 +462,39 @@ void gudcay()
     if (np <=1) return;
 
     TParticle *  iparticle = (TParticle *) particles->At(0);
-    Int_t ichF = iparticle->GetFirstDaughter();
-    Int_t ichL = iparticle->GetLastDaughter();    
-
-    for (Int_t i=ichF-1; i < ichL; i++) 
+    Int_t ipF = 0, ipL = 0 ;
+    Int_t i,j;
+
+// Array to flag deselected particles
+    Int_t* pFlag = new Int_t[np];
+    for (i=0; i<np; i++) pFlag[i]=0;
+// Particle loop
+    for (i=1; i < np; i++) 
     {
        iparticle = (TParticle *) particles->At(i);
+       ipF = iparticle->GetFirstDaughter();
+       ipL = iparticle->GetLastDaughter();     
        Int_t kf = iparticle->GetPdgCode();
-//     Int_t ks = iparticle->GetStatusCode();
-//     printf("\n %d %d %d %d %d %d", i, np, kf, ks, ichL, ichF);
-// Final state particles only
-//     if (ks < 1 || ks >  10) continue;
+       Int_t ks = iparticle->GetStatusCode();
+//
+// Deselect daughters of deselected particles
+// and jump skip the current particle
+       if (pFlag[i] == 1) {
+           if (ipF >= 0) for (j=ipF-1; j<ipL; j++) pFlag[j]=1;
+           continue;
+       } // deselected ??
+// Particles with long life-time are put on the stack for further tracking
+// Decay products are deselected
+//     
+       if (ks != 1) { 
+           Double_t lifeTime = gMC->Decayer()->GetLifetime(kf);
+           if (lifeTime > (Double_t) 1.e-15) {
+               if (ipF >= 0) for (j=ipF-1; j<ipL; j++) pFlag[j]=1;
+           } else{
+               continue;
+           }
+       } // ks==1 ?
 // Skip neutrinos
-
-
        if (kf==12 || kf ==-12) continue;
        if (kf==14 || kf ==-14) continue;
        if (kf==16 || kf ==-16) continue;
@@ -480,7 +502,7 @@ void gudcay()
        Int_t index=geant3->Gcking()->ngkine;
 // Put particle on geant stack
 // momentum vector
-
+       
        (geant3->Gcking()->gkin[index][0]) = iparticle->Px();
        (geant3->Gcking()->gkin[index][1]) = iparticle->Py();
        (geant3->Gcking()->gkin[index][2]) = iparticle->Pz();
@@ -495,10 +517,11 @@ void gudcay()
        (geant3->Gckin3()->gpos[index][2]) = geant3->Gctrak()->vect[2];
 // time of flight offset (mm)
        (geant3->Gcking()->tofd[index])    = 0.;
-//     (geant3->Gcking()->tofd[index])    = iparticle->T()/(10*kSpeedOfLight);
 // increase stack counter
        (geant3->Gcking()->ngkine)=index+1;
     }
+//
+    if (pFlag) delete[] pFlag;
 }
 
 //______________________________________________________________________