Update to trunk of hepforge

[u/mrichter/AliRoot.git] / STARLIGHT / starlight / src / incoherentVMCrossSection.cpp

diff --git a/STARLIGHT/starlight/src/incoherentVMCrossSection.cpp b/STARLIGHT/starlight/src/incoherentVMCrossSection.cpp
index 32d75c8fb99c9fe264f329b3e9df1c7b340941d0..8c0c57dfc54f1d3bac1b90b7aea1ae656c3a23cb 100644 (file)
--- a/STARLIGHT/starlight/src/incoherentVMCrossSection.cpp
+++ b/STARLIGHT/starlight/src/incoherentVMCrossSection.cpp
@@ -69,9 +69,11 @@ incoherentVMCrossSection::crossSectionCalculation(const double)  // _bwnormsave
        // double Av,Wgp,cs,cvma;
        double W,dY;
        double y1,y2,y12,ega1,ega2,ega12;
+        // double y1lab,y2lab,y12lab;
        // double t,tmin,tmax;
        double csgA1,csgA2,csgA12,int_r,dR,rate;
         double Wgp,csVN,csVA; 
+        double Wgpm; 
        double tmp;
        // double ax,bx;
        double Eth;
@@ -94,6 +96,7 @@ incoherentVMCrossSection::crossSectionCalculation(const double)  // _bwnormsave
   
        for(J=0;J<=(NY-1);J++){
     
+               // This is the fdefault
                y1  = _narrowYmin + double(J)*dY;
                y2  = _narrowYmin + double(J+1)*dY;
                y12 = 0.5*(y1+y2);
@@ -101,6 +104,28 @@ incoherentVMCrossSection::crossSectionCalculation(const double)  // _bwnormsave
                ega1  = 0.5*W*exp(y1);
                ega2  = 0.5*W*exp(y2);
                ega12 = 0.5*W*exp(y12);
+
+                // This is for checking things in the lab frame 
+                // y1lab  = _narrowYmin + double(J)*dY;
+                // y2lab  = _narrowYmin + double(J+1)*dY;
+                // y12lab = 0.5*(y1lab+y2lab);
+
+                // p+Pb
+                // y1 = y1lab + 0.465;
+                // y2 = y2lab + 0.465;
+                // y12 = y12lab + 0.465; 
+                // ega1  = 0.5*W*exp(y1);
+                // ega2  = 0.5*W*exp(y2);
+                // ega12 = 0.5*W*exp(y12);
+
+                // Pb+p
+                // y1 = y1lab - 0.465;
+                // y2 = y2lab - 0.465;
+                // y12 = y12lab - 0.465; 
+                // ega1  = 0.5*W*exp(-y1);
+                // ega2  = 0.5*W*exp(-y2);
+                // ega12 = 0.5*W*exp(-y12);
+
     
                if(ega1 < Eth)   
                        continue;
@@ -120,6 +145,7 @@ incoherentVMCrossSection::crossSectionCalculation(const double)  // _bwnormsave
                // Middle point 
                 Wgp = sqrt(2.*ega12*(_Ep+sqrt(_Ep*_Ep-starlightConstants::protonMass*starlightConstants::protonMass))
                                  +starlightConstants::protonMass*starlightConstants::protonMass);
+                Wgpm = Wgp; 
                 csVN = sigma_N(Wgp);            
                 csVA = sigma_A(csVN); 
                 csgA12 = (csVA/csVN)*sigmagp(Wgp); 
@@ -143,6 +169,7 @@ incoherentVMCrossSection::crossSectionCalculation(const double)  // _bwnormsave
                 dR = dR*(dY/6.); 
 
                // cout<<" y: "<<y12<<" egamma: "<<ega12<<" flux: "<<photonFlux(ega12)<<" sigma_gA: "<<10000000.*csgA12<<" dsig/dy (microb): "<<10000.*dR/dY<<endl;
+                // cout<<" y: "<<y12lab<<" egamma: "<<ega12<<" flux: "<<ega12*photonFlux(ega12)<<" W: "<<Wgpm<<" Wflux: "<<Wgpm*photonFlux(ega12)<<" sigma_gA (nb): "<<10000000.*csgA12<<" dsig/dy (microb): "<<10000.*ega12*photonFlux(ega12)*csgA12<<endl;
 
                // The 2 accounts for the 2 beams    
                if(getbbs().beam1().A()==getbbs().beam2().A()){