--- calo/flight/CaloBragg/src/CaloBragg.cpp	2011/06/15 10:04:59	1.10
+++ calo/flight/CaloBragg/src/CaloBragg.cpp	2011/06/17 11:56:42	1.12
@@ -46,7 +46,7 @@
   lpetot = 0.;
   lppskip = 0.;
   memset(calorimetro,0,44*2*sizeof(Float_t));
-  memset(spessore,0,3*sizeof(Float_t));
+  memset(spessore,0,4*sizeof(Float_t));
   memset(estremi,0,2*2*sizeof(Float_t));
   Integrale=0.;
  
@@ -367,17 +367,58 @@
   return integ;
 }
 
+
 void CaloBragg::Draw(){
 
   Process();
 
+  this->Draw(0.,0.);
+
+}
+
+void CaloBragg::Draw(Int_t Z, Float_t enetot){
+
   //  Float_t dEpianimean[44];
   Float_t dEpianiloop[44];
   Float_t Depth[44];
   //  Int_t tz=(Int_t)qtz;
-  Int_t tz1=(Int_t)lpz;
-  //  Enetrack(&tz, &qtetot, &estremi[0][0],&estremi[1][0], dEpianimean);//calcola rilascio energetico sui piani da media troncata
-  Enetrack(&tz1, &lpetot, &estremi[0][0],&estremi[1][0], dEpianiloop);//calcola rilascio energetico sui piani da loop
+  Int_t tz1= Z; 
+  Float_t enet = enetot;
+  //  Float_t enet = lpetot;
+  
+  if ( Z > 0. && enetot > 0. ){
+    estremi[0][0] = 0;
+    estremi[1][0] = 43;
+
+
+    Float_t    ytgx = 0.;
+    Float_t    ytgy = 0.;
+    
+    //lunghezza effettiva di silicio attraversata (mm)
+    Float_t    SiCross = sqrt(SQ(ySi) + SQ(ytgx) + SQ(ytgy));
+    
+    spessore[0] = (SiCross/10.) * rhoSi; //spessore silicio in g/cm2
+    
+    /*tungsteno*/
+    
+    //rapporto tra rilasci energetici nei due materiali
+    Float_t WCross = sqrt((yW*yW) + (ytgx*ytgx) + (ytgy*ytgy));//mm* rapporto lunghezze rad
+    //gcm2W = WCross/10. * rhoW;
+    
+    //       (g/cm2W)/(g/cm2Si)
+    spessore[3] =  (WCross/10.) * rhoW;
+    Float_t a=(WCross/SiCross)*(rhoW/rhoSi)*(1.145/1.664);  //(gcm2W)/(SiCross/10. * rhoSi)* (1.145/1.664);
+    spessore[1] =  a;
+    //riscala mip allo spessore attraversato
+    spessore[2] = MIP*(SiCross/ySi);   
+
+  } else {
+    tz1=(Int_t)lpz;
+    enet = lpetot;
+    //  Enetrack(&tz, &qtetot, &estremi[0][0],&estremi[1][0], dEpianimean);//calcola rilascio energetico sui piani da media troncata
+  
+  }
+  Enetrack(&tz1, &enet, &estremi[0][0],&estremi[1][0], dEpianiloop);//calcola rilascio energetico sui piani da loop
 
   Float_t sp= spessore[0]*spessore[1];
   for(Int_t i=0;i<44;i++)Depth[i]=i*sp;
@@ -430,7 +471,10 @@
   //  tc->cd(2);
   tc->cd();
   //
-  for(Int_t i=0;i<=estremi[1][0];i++)th3->Fill(Depth[i],dEpianiloop[i]);
+  for(Int_t i=0;i<=estremi[1][0];i++){
+    th3->Fill(Depth[i],dEpianiloop[i]);
+    //    printf(" i %i Depth %f depianiloop %f \n",i,Depth[i],dEpianiloop[i]);
+  }
   th3->Draw();
   th2->Draw("same");
 
@@ -552,15 +596,11 @@
   //gcm2W = WCross/10. * rhoW;
 
   //       (g/cm2W)/(g/cm2Si)
-  if ( usenewBB ){
-    spessore[1] =  (WCross/10.) * rhoW;
-  } else {
-
-    a=(WCross/SiCross)*(rhoW/rhoSi)*(1.145/1.664);  //(gcm2W)/(SiCross/10. * rhoSi)* (1.145/1.664);
-    spessore[1] =  a;
-    //riscala mip allo spessore attraversato
-    spessore[2] = MIP*(SiCross/ySi);   
-  };
+  spessore[3] =  (WCross/10.) * rhoW;
+  a=(WCross/SiCross)*(rhoW/rhoSi)*(1.145/1.664);  //(gcm2W)/(SiCross/10. * rhoSi)* (1.145/1.664);
+  spessore[1] =  a;
+  //riscala mip allo spessore attraversato
+  spessore[2] = MIP*(SiCross/ySi);   
 };//end conversione
 
 
@@ -585,9 +625,9 @@
   Float_t lg =0.;
   Float_t Energia=0.;
   Float_t C=0.;
-  Float_t INo = II;
+  Float_t INo = ISi;
 
-  if ( !usenewBB ) INo = ISi;
+  if ( usenewBB ) INo = II;
 
   eta = (*gam)*(*Bet);
 
@@ -693,11 +733,12 @@
       Float_t sp = 0.;
       Float_t II = ISi;
       if ( usenewBB ){
-        sp = spessore[1];        
+        sp = spessore[3];        
         II = IW;
       } else {
         sp = spessore[0]*spessore[1]; //((gcm2Si)*(WinSi))//spessore attraversato  in g/cm2
       }
+      //      printf(" sp %f II %f \n",sp,II);
       ELOSS(&sp, Z, &Ezero, &dE,II);
       if((Ezero-dE) <= Massa){//se l'energia depositata e' maggiore dell'energia della perticella stop
 	return;