/[PAMELA software]/DarthVader/CalorimeterLevel2/src/CaloProcessing.cpp
ViewVC logotype

Diff of /DarthVader/CalorimeterLevel2/src/CaloProcessing.cpp

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

revision 1.12 by mocchiut, Tue Jan 23 11:52:25 2007 UTC revision 1.13 by mocchiut, Tue Feb 13 16:55:08 2007 UTC
# Line 67  CaloProcessing::CaloProcessing(){ Line 67  CaloProcessing::CaloProcessing(){
67    calopar1 = true;    calopar1 = true;
68    calopar2 = true;    calopar2 = true;
69    calopar3 = true;    calopar3 = true;
70      crosst = true;
71    ftcalopar1 = 0;    ftcalopar1 = 0;
72    ttcalopar1 = 0;    ttcalopar1 = 0;
73    ftcalopar2 = 0;    ftcalopar2 = 0;
# Line 75  CaloProcessing::CaloProcessing(){ Line 76  CaloProcessing::CaloProcessing(){
76    ttcalopar3 = 0;    ttcalopar3 = 0;
77  }  }
78    
79    void CaloProcessing::SetCrossTalk(Bool_t ct){
80      crosst = ct;
81    };
82    
83  /**  /**
84   * Initialize CaloProcessing object   * Initialize CaloProcessing object
85  **/  **/
# Line 575  Int_t CaloProcessing::Calibrate(Int_t ei Line 580  Int_t CaloProcessing::Calibrate(Int_t ei
580                };                };
581              };              };
582            };            };
583            if (ck == 1){            if ( crosst ){
584              if (ip[i]%2 == 0) {              if (ck == 1){
585                ipre = ip[i] + 1;                if (ip[i]%2 == 0) {
586              } else {                  ipre = ip[i] + 1;
587                ipre = ip[i] - 1;                } else {
588              };                  ipre = ip[i] - 1;
589              for (Int_t j = ipre*16 ; j < (ipre+1)*16 ; j++){                };
590                clevel1->estrip[j][m][l] += (qpre[ipre] - qpre[ip[i]]) * 0.00478;                for (Int_t j = ipre*16 ; j < (ipre+1)*16 ; j++){
591                    clevel1->estrip[j][m][l] += (qpre[ipre] - qpre[ip[i]]) * 0.00478;
592                  };
593              };              };
594            };              if (ck == 2){
595            if (ck == 2){                for (Int_t j = i*32 ; j < (i+1)*32 ; j++){
596              for (Int_t j = i*32 ; j < (i+1)*32 ; j++){                  ipre = j/16 + 1;
597                ipre = j/16 + 1;                  clevel1->estrip[j][m][l] +=  qpre[ipre] * 0.00478;
598                clevel1->estrip[j][m][l] +=  qpre[ipre] * 0.00478;                };
599              };              };
600            };            };
601          };          };
# Line 607  Int_t CaloProcessing::Calibrate(Int_t ei Line 614  Int_t CaloProcessing::Calibrate(Int_t ei
614            jjj++;            jjj++;
615            j4++;            j4++;
616            if ( j < 96 ) ene[j] = clevel1->estrip[j][m][l];            if ( j < 96 ) ene[j] = clevel1->estrip[j][m][l];
617            if ( jj >= 0 && jj < 96 ){            if ( crosst ){
618              if ( jj != 0 && jj != 32 && jj != 64 ) ene[jj-1] += -clevel1->estrip[jj][m][l] * 0.01581;              if ( jj >= 0 && jj < 96 ){
619              if ( jj != 31 && jj != 63 && jj != 95 ) ene[jj+1] += -clevel1->estrip[jj][m][l] * 0.01581;                                  if ( jj != 0 && jj != 32 && jj != 64 ) ene[jj-1] += -clevel1->estrip[jj][m][l] * 0.01581;
620            };                if ( jj != 31 && jj != 63 && jj != 95 ) ene[jj+1] += -clevel1->estrip[jj][m][l] * 0.01581;                        
621            if ( jjj >= 0 && jjj < 96 ){              };
622              if ( jjj != 0 && jjj != 32 && jjj != 64 ) clevel1->estrip[jjj-1][m][l] +=  -ene[jjj] * 0.01581;              if ( jjj >= 0 && jjj < 96 ){
623              if ( jjj != 31 && jjj != 63 && jjj != 95 ) clevel1->estrip[jjj+1][m][l] +=  -ene[jjj] * 0.01581;                if ( jjj != 0 && jjj != 32 && jjj != 64 ) clevel1->estrip[jjj-1][m][l] +=  -ene[jjj] * 0.01581;
624                  if ( jjj != 31 && jjj != 63 && jjj != 95 ) clevel1->estrip[jjj+1][m][l] +=  -ene[jjj] * 0.01581;
625                };
626            };            };
627            if ( j4 >= 0 && j4 < 96 ){            if ( j4 >= 0 && j4 < 96 ){
628              //              //
Legend:
Removed from v.1.12  
changed lines
  Added in v.1.13
  ViewVC Help
Powered by ViewVC 1.1.23