/[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.19 by mocchiut, Mon Dec 14 14:42:53 2009 UTC
# Line 28  CaloNuclei::CaloNuclei(PamLevel2 *l2p){ Line 28  CaloNuclei::CaloNuclei(PamLevel2 *l2p){
28    debug = false;    debug = false;
29    // debug = true;    // debug = true;
30    usetrack = true;    usetrack = true;
31      usepl18x = false;
32    //    //
33  };  };
34    
# Line 146  void CaloNuclei::Process(Int_t ntr){ Line 147  void CaloNuclei::Process(Int_t ntr){
147      //      //
148      mip = L2->GetCaloLevel1()->DecodeEstrip(i,view,plane,strip);      mip = L2->GetCaloLevel1()->DecodeEstrip(i,view,plane,strip);
149      //      //
150        if ( !usepl18x && view==0 && plane==18 ) mip = 0.;
151        //
152        //
153      // put in vfpl vector the energy release on the first plane      // put in vfpl vector the energy release on the first plane
154      //      //
155      if ( strip != -1 && view == 1 && plane == 0 ) {      if ( strip != -1 && view == 1 && plane == 0 ) {
# Line 214  void CaloNuclei::Process(Int_t ntr){ Line 218  void CaloNuclei::Process(Int_t ntr){
218      //      //
219      mip = L2->GetCaloLevel1()->DecodeEstrip(i,view,plane,strip);      mip = L2->GetCaloLevel1()->DecodeEstrip(i,view,plane,strip);
220      //      //
221        if ( !usepl18x && view==0 && plane==18 ) mip = 0.;
222        //
223      if ( ntr >= 0 ){      if ( ntr >= 0 ){
224        //        //
225        if ( strip != -1 &&        if ( strip != -1 &&
# Line 283  void CaloNuclei::Process(Int_t ntr){ Line 289  void CaloNuclei::Process(Int_t ntr){
289      //      //
290      mip = L2->GetCaloLevel1()->DecodeEstrip(ii,view,plane,strip);          mip = L2->GetCaloLevel1()->DecodeEstrip(ii,view,plane,strip);    
291      //      //
292        if ( !usepl18x && view==0 && plane==18 ) mip = 0.;
293        //
294        //
295      if ( ntr >= 0 ){      if ( ntr >= 0 ){
296        if ( strip != -1 && mip > ethr && !wmulthit[view] && !wgap[view] &&        if ( strip != -1 && mip > ethr && !wmulthit[view] && !wgap[view] &&
297             ( strip == (track->tibar[plane][view]-1) || strip == (track->tibar[plane][view]-2) || strip == (track->tibar[plane][view]) )             ( strip == (track->tibar[plane][view]-1) || strip == (track->tibar[plane][view]-2) || strip == (track->tibar[plane][view]) )
# Line 411  void CaloNuclei::Process(Int_t ntr){ Line 420  void CaloNuclei::Process(Int_t ntr){
420        //        //
421        mip = L2->GetCaloLevel1()->DecodeEstrip(ii,view,plane,strip);            mip = L2->GetCaloLevel1()->DecodeEstrip(ii,view,plane,strip);    
422        //        //
423          if ( !usepl18x && view==0 && plane==18 ) mip = 0.;
424          //
425          //
426        if ( strip != -1 ){        if ( strip != -1 ){
427          if ( view == 0 ){          if ( view == 0 ){
428            ipl = (1 + plane) * 2;            ipl = (1 + plane) * 2;
# Line 501  void CaloNuclei::Process(Int_t ntr){ Line 513  void CaloNuclei::Process(Int_t ntr){
513      Int_t uplim2 = interplane-1;      Int_t uplim2 = interplane-1;
514      //      //
515      while ( l < uplim && ind < interplane ){      while ( l < uplim && ind < interplane ){
516        qm = TMath::KOrdStat(interplane,qme,ind,work);        qm = TMath::KOrdStat((Long64_t)interplane,qme,(Long64_t)ind,work);
517        if ( qm >= qmt ){        if ( qm >= qmt ){
518          if ( l < 3 ){          if ( l < 3 ){
519            qpremean += qm;            qpremean += qm;
# Line 518  void CaloNuclei::Process(Int_t ntr){ Line 530  void CaloNuclei::Process(Int_t ntr){
530      l = 0;      l = 0;
531      RN = 0;      RN = 0;
532      while ( l < uplim && ind < interplane ){      while ( l < uplim && ind < interplane ){
533        qm2 = TMath::KOrdStat(interplane,qme2,ind,work);        qm2 = TMath::KOrdStat((Long64_t)interplane,qme2,(Long64_t)ind,work);
534        if ( qm2 >= qmt ){                if ( qm2 >= qmt ){        
535          if ( l < N ){          if ( l < N ){
536            qpremeanN += qm2;            qpremeanN += qm2;
# Line 539  void CaloNuclei::Process(Int_t ntr){ Line 551  void CaloNuclei::Process(Int_t ntr){
551      RN = 0;      RN = 0;
552      S2=0;      S2=0;
553      while ( l < uplim2 && ind<interplane){      while ( l < uplim2 && ind<interplane){
554        qm2 = TMath::KOrdStat(interplane,qme2,ind,work);        qm2 = TMath::KOrdStat((Long64_t)interplane,qme2,(Long64_t)ind,work);
555        if ( qm2 < qmt ) S2++;        if ( qm2 < qmt ) S2++;
556        ind++;        ind++;
557      }      }
# Line 548  void CaloNuclei::Process(Int_t ntr){ Line 560  void CaloNuclei::Process(Int_t ntr){
560      l = 0;      l = 0;
561      RN = 0;      RN = 0;
562      while ( l < uplim2 && ind < interplane ){      while ( l < uplim2 && ind < interplane ){
563        qm2 = TMath::KOrdStat(interplane,qme2,ind,work);        qm2 = TMath::KOrdStat((Long64_t)interplane,qme2,(Long64_t)ind,work);
564        if ( qm2 >= qmt ){                if ( qm2 >= qmt ){        
565          if ( l < (interplane - 1 - S2)){            if ( l < (interplane - 1 - S2)){  
566            qNmin1_w += qm2;            qNmin1_w += qm2;
Legend:
Removed from v.1.17  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.23