/[PAMELA software]/DarthVader/TrackerLevel2/src/TrkLevel1.cpp
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/TrkLevel1.cpp

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

revision 1.28 by pam-ts, Wed Jun 4 07:57:04 2014 UTC revision 1.29 by mocchiut, Thu Jul 3 10:38:00 2014 UTC
# Line 908  int TrkLevel1::ProcessEvent(){ Line 908  int TrkLevel1::ProcessEvent(){
908  void TrkLevel1::Set(TrkLevel1 *trkl1, float mipCut, float fCut){  void TrkLevel1::Set(TrkLevel1 *trkl1, float mipCut, float fCut){
909    
910    
   
911            
912      if(!trkl1)return;      if(!trkl1)return;
913    
# Line 930  void TrkLevel1::Set(TrkLevel1 *trkl1, fl Line 929  void TrkLevel1::Set(TrkLevel1 *trkl1, fl
929      Cluster = new TClonesArray("TrkCluster");      Cluster = new TClonesArray("TrkCluster");
930      TClonesArray &t = *Cluster;      TClonesArray &t = *Cluster;
931    
932      int isel=0;      if ( trkl1->GetClusters() ){ // EM: crash nuclei alg when  DSP status: -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
933      for(int icl=0 ; icl< trkl1->GetClusters()->GetEntries(); icl++){        int isel=0;
934          for(int icl=0 ; icl< trkl1->GetClusters()->GetEntries(); icl++){
935          TrkCluster *cl = trkl1->GetCluster(icl);          TrkCluster *cl = trkl1->GetCluster(icl);
936            
937          float mip = TrkParams::GetMIP(cl->GetLadder()-1,cl->view-1);          float mip = TrkParams::GetMIP(cl->GetLadder()-1,cl->view-1);
938          float smip = cl->GetSignal()/(mip>0.?mip:1.);          float smip = cl->GetSignal()/(mip>0.?mip:1.);
939          float smax =  cl->clsignal[cl->indmax]/(mip>0.?mip:1.);          float smax =  cl->clsignal[cl->indmax]/(mip>0.?mip:1.);
# Line 941  void TrkLevel1::Set(TrkLevel1 *trkl1, fl Line 941  void TrkLevel1::Set(TrkLevel1 *trkl1, fl
941          if(smip<mipCut)continue;          if(smip<mipCut)continue;
942          if(smax<0.5*mipCut)continue;          if(smax<0.5*mipCut)continue;
943                    
944            
945            
946          new(t[isel]) TrkCluster(*cl); // <<< store cluster          new(t[isel]) TrkCluster(*cl); // <<< store cluster
947          isel++;          isel++;
948          }
949      }      }
950    
951    
Legend:
Removed from v.1.28  
changed lines
  Added in v.1.29
  ViewVC Help
Powered by ViewVC 1.1.23