--- calo/flight/CaloNuclei/src/CaloNuclei.cpp 2007/09/28 10:05:00 1.8 +++ calo/flight/CaloNuclei/src/CaloNuclei.cpp 2008/03/12 14:19:21 1.12 @@ -29,11 +29,14 @@ void CaloNuclei::Clear(){ // + UN = 0; tr = 0; sntr = 0; interplane = 0; preq = 0.; postq = 0.; + stdedx1 = 0.; + ethr = 0.; dedx1 = 0.; dedx3 = 0.; qpremean = 0.; @@ -106,6 +109,7 @@ if ( debug ) printf(" Processing event at OBT %u PKT %u time %u \n",OBT,PKT,atime); // Clear(); + // if ( debug ) printf(" Always calculate stdedx1 \n"); // // Always calculate stdedx1 @@ -138,7 +142,7 @@ // find energy released along the strip of maximum on the first plane and on the two neighbour strips // if ( indx > 0 ){ - Int_t mindx = (Int_t)TMath::LocMax(indx,stfpl); + Int_t mindx = (Int_t)TMath::LocMax(indx,vfpl); for (Int_t ii=0; ii=0 && stfpl[ii] == stfpl[mindx-1] ) stdedx1 += vfpl[ii]; @@ -171,8 +175,8 @@ }; } else { if ( ntr >= 0 ){ - printf(" ERROR: you asked not to use a track but you are looking for track number %i !\n",ntr); - printf(" ERROR: CaloNuclei variables not completely filled \n"); + if ( debug ) printf(" ERROR: you asked not to use a track but you are looking for track number %i !\n",ntr); + if ( debug ) printf(" ERROR: CaloNuclei variables not completely filled \n"); return; }; }; @@ -500,6 +504,7 @@ }; // qpremeanN /= (Float_t)RN; + UN = RN; // if ( debug ) printf(" charge is %f \n",sqrt(qpremean)); //