change print out info

[u/mrichter/AliRoot.git] / LHC / AliLhcProcessIBS.cxx

diff --git a/LHC/AliLhcProcessIBS.cxx b/LHC/AliLhcProcessIBS.cxx
index fbb0994bd0f4ceed8ed2fd9b415b400b0af36a59..a5263d22cf3a03602eb7d1129415bf8fd1c823b4 100644 (file)
--- a/LHC/AliLhcProcessIBS.cxx
+++ b/LHC/AliLhcProcessIBS.cxx
@@ -23,17 +23,18 @@
 // andreas.morsch@cern.ch
 //
 
+#include <TCanvas.h>
+#include <TF1.h>
+#include <TGraph.h>
+#include <TH1F.h>
+#include <TMath.h>
+#include <TMultiGraph.h>
+
 #include "AliLhcProcessIBS.h"
 #include "AliLHC.h"
 #include "AliLhcIRegion.h"
 #include "AliLhcBeam.h"
 
-#include <TF1.h>
-#include <TMath.h>
-#include <TCanvas.h>
-#include <TGraph.h>
-#include <TMultiGraph.h>
-
 ClassImp(AliLhcProcessIBS)
 
   Double_t func(Double_t *x, Double_t *par);
@@ -113,9 +114,9 @@ void AliLhcProcessIBS::Evolve(Float_t dt)
      // impact parameter
      Float_t d    = sige*avd/sigx;
      // 
-     Float_t at = sige*TMath::Sqrt(1.-d*d)/(gamma*ssigx);
-     Float_t bt = sige*TMath::Sqrt(1.-d*d)/(gamma*ssigy);
-     Float_t ct = sige*TMath::Sqrt(1.-d*d)*TMath::Sqrt(4.*sigy/fR[i]);
+     Float_t at = sige*TMath::Sqrt((1.-d)*(1.+d))/(gamma*ssigx);
+     Float_t bt = sige*TMath::Sqrt((1.-d)*(1.+d))/(gamma*ssigy);
+     Float_t ct = sige*TMath::Sqrt((1.-d)*(1.+d))*TMath::Sqrt(4.*sigy/fR[i]);
      //
      Double_t par[3];
      par[0] = at;
@@ -126,8 +127,9 @@ void AliLhcProcessIBS::Evolve(Float_t dt)
      Float_t f = (Float_t) 8.0*TMath::Pi()*
         fct->Integral(0., 1., par, 1.e-5);
      //
-     fTaux =  1./(asd*f*(d*d-at*at/2.));
-     fTaue =  1./(asd*f*(1.-d*d));
+     const Double_t osq2=1./TMath::Sqrt(2.);
+     fTaux =  1./(asd*f*(d-osq2*at)*(d+osq2*at)); 
+     fTaue =  1./(asd*f*((1.-d)*(1.+d)));
      //     printf("\n taux, taue %f %f", taux, taue);
      //     Float_t tauy = -2./at*at/asd/f;
      fBeam[i]->IncreaseEmittance(dt/fTaux, dt/fTaue);