/[PAMELA software]/PamCut/CaloCuts/CaloNotIntCut/CaloNotIntCut.cpp
ViewVC logotype

Diff of /PamCut/CaloCuts/CaloNotIntCut/CaloNotIntCut.cpp

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

revision 1.1 by pam-fi, Wed May 27 13:30:09 2009 UTC revision 1.4 by pam-fi, Tue Nov 24 14:42:31 2009 UTC
# Line 10  Line 10 
10    
11  int CaloNotIntCut::Check(PamLevel2 *event) {  int CaloNotIntCut::Check(PamLevel2 *event) {
12    
13    if ( (!(_xCaloAxis) || !(_yCaloAxis)) && !(_trkTrack) ) { // check that one of the two possible initializations has been correctly done    if ( (!(_xCaloAxis) || !(_yCaloAxis)) && !(_pamTrack) ) { // check that one of the two possible initializations has been correctly done
14      return 0;      return 0;
15    }    }
16        
17    if (!(_trkTrack)) { // CALO track    if (!(_pamTrack)) { // CALOAXIS track
18        
19      // energy released along the event axis in the calo / total energy release in the calo      // energy released along the event axis in the calo / total energy release in the calo
20            
21      float qratio = ( _xCaloAxis->GetQaxis() + _yCaloAxis->GetQaxis() ) / event->GetCaloLevel2()->qtot;      Double_t qtot = event->GetCaloLevel2()->qtot;
22        if ( qtot>0 ) {
23        
24          Double_t qratio = ( _xCaloAxis->GetQaxis() + _yCaloAxis->GetQaxis() ) / qtot;
25    
26      if (! (qratio > _qRatioMin) ) {        if (! (qratio > _qRatioMin) ) {
27        return 1;          return 1;
28      }        }
29        
30        }
31    }    }
32        
33    if (_trkTrack) { // TRK track    if (_pamTrack) { // PAMELA track
34      
35     return 2; // temporary      Double_t qtot = event->GetCaloLevel2()->qtot;
36      // Implement here the check using the tracker track.      if ( qtot>0 ) {
37    
38          Double_t qratio = ( (*_pamTrack)->GetCaloTrack()->qtrack ) / qtot;
39    
40          if (! (qratio > _qRatioMin) ) {  
41            return 2;
42          }
43          
44        }
45    
46    }    }
47    
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.23