--- calo/flight/CaloBragg/src/CaloBragg.cpp	2011/06/15 12:28:45	1.11
+++ calo/flight/CaloBragg/src/CaloBragg.cpp	2011/06/17 11:56:42	1.12
@@ -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");