added InSANEPolradDiffXSec typdef ... not sure if it is correct

	modified:   include/InSANEInelasticRadiativeTail.h
	modified:   include/InSANEPOLRADInternalPolarizedDiffXSec.h
	modified:   include/InSANEPhysics_LinkDef.h
	modified:   main/radiative_corrections/sub_el_tail_1D_perp59.cxx

include/InSANEInelasticRadiativeTail.h

@@ -44,18 +44,14 @@ class InSANEInelasticRadiativeTail : public InSANERadiativeTail {
 
          Double_t sig_rad = 0.0;
 
-
-         std::cout << "test 1" << std::endl;
          // Using the POLRAD IRT
          //sig_rad = fDiffXSec0->EvaluateXSec(x);
 
          // Using the Equiv. Rad. Method
          sig_rad = fRADCOR->ContinuumStragglingStripApprox_InternalEquivRad(Ebeam,Eprime,theta,phi);
-         std::cout << "test 2" << std::endl;
          Double_t sig_rad2 = 0.0;
          if(fAddRegion4) sig_rad2 = fRADCOR->Internal2DEnergyIntegral(Ebeam,Eprime,theta,phi);
          sig_rad += sig_rad2;
-         std::cout << "test 3" << std::endl;
 
          if( IncludeJacobian() ) sig_rad = sig_rad*TMath::Sin(theta);
          if(sig_rad <0.0 || TMath::IsNaN(sig_rad)) sig_rad = 0.0;
@@ -74,8 +70,8 @@ class InSANEInelasticRadiativeTail : public InSANERadiativeTail {
 class InSANEFullInelasticRadiativeTail : public InSANERadiativeTail {
 
    protected:
-      InSANEInelasticRadiativeTail * fDiffXSec;  //Internally radiated tail
-      InSANEPOLRADInternalPolarizedDiffXSec * fDiffXSec0;  //  born cross section used by RADCOR to calculate
+      InSANEInelasticRadiativeTail           * fDiffXSec;  //Internally radiated tail
+      InSANEPOLRADInternalPolarizedDiffXSec  * fDiffXSec0;  //  born cross section used by RADCOR to calculate
     
    public:
       InSANEFullInelasticRadiativeTail();
@@ -90,17 +86,16 @@ class InSANEFullInelasticRadiativeTail : public InSANERadiativeTail {
             return(0.0);
          }
 
-         std::cout << "test 5" << std::endl;
          Double_t Eprime  = x[0];
          Double_t theta   = x[1];
          Double_t phi     = x[2];
          Double_t Ebeam   = GetBeamEnergy();
          fDiffXSec->SetBeamEnergy(Ebeam);
+         fDiffXSec0->SetBeamEnergy(Ebeam);
          Double_t nu      = Ebeam-Eprime;
          Double_t Mtarg   = fPOLRAD->GetTargetMass();  
          Double_t sig_rad = fRADCOR->ContinuumStragglingStripApprox(Ebeam,Eprime,theta,phi);
-         std::cout << "test 6" << std::endl;
-         std::cout << " sig_rad  = " << sig_rad << std::endl;
+         //std::cout << " sig_rad  = " << sig_rad << std::endl;
          if( IncludeJacobian() ) sig_rad = sig_rad*TMath::Sin(theta);
          if(sig_rad <0.0 || TMath::IsNaN(sig_rad)) sig_rad = 0.0;
          return sig_rad;

include/InSANEPOLRADInternalPolarizedDiffXSec.h

@@ -81,6 +81,8 @@ class InSANEPOLRADInternalPolarizedDiffXSec: public InSANEInclusiveDiffXSec {
 
       ClassDef(InSANEPOLRADInternalPolarizedDiffXSec,3)
 };
+
+typedef InSANEPOLRADInternalPolarizedDiffXSec InSANEPolradDiffXSec;
 //____________________________________________________________________________________
 
 /** Born Cross Section for elastic scattering.

include/InSANEPhysics_LinkDef.h

@@ -160,6 +160,7 @@
 #pragma link off function InSANEPOLRADUltraRelativistic::AdaptiveSimpsonAux(double (*)(double &), double,double,double,double,double,double,double,int); 
 
 #pragma link C++ class InSANEPOLRADInternalPolarizedDiffXSec+;
+#pragma link C++ typedef InSANEPolradDiffXSec;
 #pragma link C++ class InSANEPOLRADElasticDiffXSec+;
 #pragma link C++ class InSANEPOLRADBornDiffXSec+;
 #pragma link C++ class InSANEPOLRADRadiatedDiffXSec+;

main/radiative_corrections/sub_el_tail_1D_perp59.cxx

@@ -3,6 +3,13 @@
  */
 Int_t sub_el_tail_1D_perp59(Int_t paraRunGroup = 23,Int_t aNumber=23){
 
+   const TString weh("sub_el_tail_1D_perp59()");
+   if( gROOT->LoadMacro("asym/sane_data_bins.cxx") )  {
+      Error(weh, "Failed loading asym/sane_data_bins.cxx in compiled mode.");
+      return -1;
+   }
+   sane_data_bins();
+
    Double_t theta_target = 80.0*degree;
    Double_t fBeamEnergy  = 5.89; 
 
@@ -41,7 +48,7 @@ Int_t sub_el_tail_1D_perp59(Int_t paraRunGroup = 23,Int_t aNumber=23){
       for(Int_t i=0; i<= xBinmax; i++){
          bin   = fA_x->GetBin(i);
          std::cout << " -------------------------------------" << std::endl;
-         std::cout << " Bin " << bin << "/" << xBinmax << std::endl; 
+         std::cout << " Bin " << bin << "/" << xBinmax <<  " of Asymmetry " << jj << "/" << fAsymmetries->GetEntries() << std::endl; 
 
          //fA_x->GetBinXYZ(bin,binx,biny,binz);
          //x     = fA_x->GetBinCenter(bin);