Changeset 1352 for trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc

Timestamp:

06/10/02 08:49:28 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc (modified) (6 diffs)

trunk/WuerzburgSoft/Thomas/mphys/MPairProduction.cc

@@ -32 +32 @@
 #include <TF1.h>
 #include <TMath.h>
+#include <TRandom.h>
 
 #include "TList.h"
@@ -70 +71 @@
 
 // --------------------------------------------------------------------------
-Bool_t MPairProduction::Process(MParticle *gamma, MParticle *phot, TList *list)
+Bool_t MPairProduction::Process(MParticle *gamma, MParticle *phot, const Double_t theta, TList *list)
 {
     //
@@ -78 +79 @@
     const Double_t E0     = 511e-6;              // [GeV]
 
-    const Double_t theta  = phot->GetAngle();    // [2pi]
+    //const Double_t theta  = phot->GetAngle();    // [2pi]
     const Double_t Ep     = phot->GetEnergy();   // [GeV]
     const Double_t Eg     = gamma->GetEnergy();  // [GeV]
@@ -120 +121 @@
     const Double_t dE = E*Bcosd;
 
-    MElectron *p[2];
-    p[0] = new MElectron(E+dE, gamma->GetZ());
-    p[1] = new MElectron(E-dE, gamma->GetZ());
+    MElectron &p0 = *new MElectron(E+dE, gamma->GetZ());
+    MElectron &p1 = *new MElectron(E-dE, gamma->GetZ());
+    p0 = *gamma; // Set Position and direction
+    p1 = *gamma; // Set Position and direction
 
     const Double_t E1 = E0/(E+dE);
@@ -130 +132 @@
     const Double_t beta2 = sqrt(1.-E2*E2);
 
-    p[0]->SetBeta(beta1);
-    p[1]->SetBeta(beta2);
+    p0.SetBeta(beta1);
+    p1.SetBeta(beta2);
 
     const Double_t Bscp = Bsind*cphi;
@@ -140 +142 @@
     const Double_t tan2 = (spg*(Bcosd-1) + Bscp)/((cpg*(Bcosd-1) - Bscp));
 
-    p[0]->SetAngle(atan(tan1));
-    p[1]->SetAngle(atan(tan2));
+    static TRandom rand(0);
+    Double_t rnd = rand.Uniform(TMath::Pi() * 2);
+    p0.SetNewDirection(atan(tan1), rnd);
+    p1.SetNewDirection(atan(tan2), rnd+TMath::Pi());
 
-    list->Add(p[0]);
-    list->Add(p[1]);
+    list->Add(&p0);
+    list->Add(&p1);
 
     return kTRUE;