diff --git a/NPLib/CATS/TCATSPhysics.cxx b/NPLib/CATS/TCATSPhysics.cxx
index a229129ef3da7eff096fcbd4ffaa14e803fd025d..af36a2a654b513623a5073464d6f8ee238f73865 100644
--- a/NPLib/CATS/TCATSPhysics.cxx
+++ b/NPLib/CATS/TCATSPhysics.cxx
@@ -47,6 +47,7 @@ ClassImp(TCATSPhysics)
     m_PreTreatedData    = new TCATSData ;
     m_EventPhysics 			= this			    ;
     m_NumberOfCATS      = 0             ;
+    m_Spectra           = NULL          ;
   }
 
 ///////////////////////////////////////////////////////////////////////////
diff --git a/NPLib/Charissa/TCharissaPhysics.cxx b/NPLib/Charissa/TCharissaPhysics.cxx
index 46d3f3c143805270db0029f6dae5789344908961..cd5ec2f2e5a62f44241ce9c1d52bb1f5a70ef7c7 100755
--- a/NPLib/Charissa/TCharissaPhysics.cxx
+++ b/NPLib/Charissa/TCharissaPhysics.cxx
@@ -149,7 +149,7 @@ void TCharissaPhysics::BuildPhysicalEvent(){
       double tLayer2_Si_Y_T= -1000;
       double tLayer2_Si_E= -1000;
       double tLayer2_Si_T = -1000;
-
+      
       for(unsigned int ll = 0 ; ll < Layer2_couple.size() ; ++ll){	
         int Layer2_N = m_PreTreatedData->GetCharissaLayer2StripXEDetectorNbr(Layer2_couple[ll].X()) ;
         // Same detector
@@ -195,18 +195,18 @@ void TCharissaPhysics::BuildPhysicalEvent(){
 
       double CsIE = -1000;
       int CsIN = -1000;
+      
       // Look for an associate CsI E
       for(unsigned int t = 0 ; t < m_CsIEMult ; ++t ){
-        if(m_PreTreatedData->GetCharissaLayer2StripXEDetectorNbr( Layer1_couple[i].X() ) == m_PreTreatedData->GetCharissaCsIEDetectorNbr(t)){
+        if(m_PreTreatedData->GetCharissaLayer1StripXEDetectorNbr( Layer1_couple[i].X() ) == m_PreTreatedData->GetCharissaCsIEDetectorNbr(t)){
           CsIE = m_PreTreatedData->GetCharissaCsIEEnergy(t);
           CsIN = m_PreTreatedData->GetCharissaCsIECristalNbr(t);
         }
       }
-
       double CsIT = -1000;
       // Look for an associate CsI T
       for(unsigned int t = 0 ; t < m_CsITMult ; ++t ){
-          if(m_PreTreatedData->GetCharissaLayer2StripXEDetectorNbr( Layer1_couple[i].X() ) == m_PreTreatedData->GetCharissaCsITDetectorNbr(t))
+          if(m_PreTreatedData->GetCharissaLayer1StripXEDetectorNbr( Layer1_couple[i].X() ) == m_PreTreatedData->GetCharissaCsITDetectorNbr(t))
             CsIT = m_PreTreatedData->GetCharissaCsITTime(t);
       }
       
diff --git a/NPLib/Tiara/TTiaraBarrelPhysics.cxx b/NPLib/Tiara/TTiaraBarrelPhysics.cxx
index d75250b7082508c57c97933666ffe6cd1adeb66d..07b4d2f56a1868ea2b23314ba94f2231a108eb58 100644
--- a/NPLib/Tiara/TTiaraBarrelPhysics.cxx
+++ b/NPLib/Tiara/TTiaraBarrelPhysics.cxx
@@ -57,6 +57,7 @@ ClassImp(TTiaraBarrelPhysics)
     m_Back_E_Threshold = 0.40 ;
     m_Maximum_FrontBack_Difference = 0.4;
     m_OuterBack_E_Threshold =0;
+    m_Spectra = NULL ;
   }
 
 ///////////////////////////////////////////////////////////////////////////
diff --git a/NPLib/Tiara/TTiaraHyballPhysics.cxx b/NPLib/Tiara/TTiaraHyballPhysics.cxx
index fd790053b6cb71abe7caa8a4ab4bc8936bfb4c4a..b948affd8513bad7a318fc681c0a703df2259c92 100644
--- a/NPLib/Tiara/TTiaraHyballPhysics.cxx
+++ b/NPLib/Tiara/TTiaraHyballPhysics.cxx
@@ -65,6 +65,7 @@ ClassImp(TTiaraHyballPhysics)
 
     m_Take_E_Ring=false;
     m_Take_T_Sector=true;
+    m_Spectra = NULL;
   }
 
 ///////////////////////////////////////////////////////////////////////////