/[PAMELA software]/calo/flight/CaloNuclei/src/CaloNuclei.cpp
ViewVC logotype

Diff of /calo/flight/CaloNuclei/src/CaloNuclei.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.17 by marcelli, Thu Dec 18 09:43:46 2008 UTC revision 1.18 by mocchiut, Thu Dec 18 10:31:00 2008 UTC
# Line 501  void CaloNuclei::Process(Int_t ntr){ Line 501  void CaloNuclei::Process(Int_t ntr){
501      Int_t uplim2 = interplane-1;      Int_t uplim2 = interplane-1;
502      //      //
503      while ( l < uplim && ind < interplane ){      while ( l < uplim && ind < interplane ){
504        qm = TMath::KOrdStat(interplane,qme,ind,work);        qm = TMath::KOrdStat((Long64_t)interplane,qme,(Long64_t)ind,work);
505        if ( qm >= qmt ){        if ( qm >= qmt ){
506          if ( l < 3 ){          if ( l < 3 ){
507            qpremean += qm;            qpremean += qm;
# Line 518  void CaloNuclei::Process(Int_t ntr){ Line 518  void CaloNuclei::Process(Int_t ntr){
518      l = 0;      l = 0;
519      RN = 0;      RN = 0;
520      while ( l < uplim && ind < interplane ){      while ( l < uplim && ind < interplane ){
521        qm2 = TMath::KOrdStat(interplane,qme2,ind,work);        qm2 = TMath::KOrdStat((Long64_t)interplane,qme2,(Long64_t)ind,work);
522        if ( qm2 >= qmt ){                if ( qm2 >= qmt ){        
523          if ( l < N ){          if ( l < N ){
524            qpremeanN += qm2;            qpremeanN += qm2;
# Line 539  void CaloNuclei::Process(Int_t ntr){ Line 539  void CaloNuclei::Process(Int_t ntr){
539      RN = 0;      RN = 0;
540      S2=0;      S2=0;
541      while ( l < uplim2 && ind<interplane){      while ( l < uplim2 && ind<interplane){
542        qm2 = TMath::KOrdStat(interplane,qme2,ind,work);        qm2 = TMath::KOrdStat((Long64_t)interplane,qme2,(Long64_t)ind,work);
543        if ( qm2 < qmt ) S2++;        if ( qm2 < qmt ) S2++;
544        ind++;        ind++;
545      }      }
# Line 548  void CaloNuclei::Process(Int_t ntr){ Line 548  void CaloNuclei::Process(Int_t ntr){
548      l = 0;      l = 0;
549      RN = 0;      RN = 0;
550      while ( l < uplim2 && ind < interplane ){      while ( l < uplim2 && ind < interplane ){
551        qm2 = TMath::KOrdStat(interplane,qme2,ind,work);        qm2 = TMath::KOrdStat((Long64_t)interplane,qme2,(Long64_t)ind,work);
552        if ( qm2 >= qmt ){                if ( qm2 >= qmt ){        
553          if ( l < (interplane - 1 - S2)){            if ( l < (interplane - 1 - S2)){  
554            qNmin1_w += qm2;            qNmin1_w += qm2;
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.18
  ViewVC Help
Powered by ViewVC 1.1.23