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

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

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

revision 1.2 by pam-ts, Wed Jun 4 07:57:03 2014 UTC revision 1.6 by mocchiut, Thu Jul 31 10:14:46 2014 UTC
# Line 24  ExtTrkingAlg::ExtTrkingAlg(Int_t id){ Line 24  ExtTrkingAlg::ExtTrkingAlg(Int_t id){
24    
25    _whichAlg = id;    _whichAlg = id;
26    
   SetDebug(false);  
27    
28      SetDebug(false);
29        
30    if(id == 0){    if(id == 0){
31        if(_debug){
32      cout << "ExtTrkingAlg::ExtTrkingAlg("<<id<<")"<<endl;        cout << "ExtTrkingAlg::ExtTrkingAlg("<<id<<")"<<endl;
33      cout << "Creating array of TrkTrack objects "<<endl;  //      cout << "Creating array of TrkTrack objects "<<endl;
34      cout << "WARNING!!! tracking not accurate!! bug not fixed yet... "<<endl;        cout << "WARNING!!! tracking not accurate!! bug not fixed yet... "<<endl;
35            }
36      _trkArray = new TClonesArray("TrkTrack");      _trkArray = new TClonesArray("TrkTrack");
37      //      //
38      _sel_nClstrMAX  = 20. ; // selection parameter: maximum number of cluster      _sel_nClstrMAX  = 20. ; // selection parameter: maximum number of cluster
# Line 55  ExtTrkingAlg::ExtTrkingAlg(Int_t id){ Line 55  ExtTrkingAlg::ExtTrkingAlg(Int_t id){
55    }else if ( (id >= 100 && id < 144) ||    }else if ( (id >= 100 && id < 144) ||
56               (id >= 200 && id < 244) ||               (id >= 200 && id < 244) ||
57               false){ //add n calo planes               false){ //add n calo planes
58         if(_debug){
59      cout << "ExtTrkingAlg::ExtTrkingAlg("<<id<<")"<<endl;        cout << "ExtTrkingAlg::ExtTrkingAlg("<<id<<")"<<endl;
60      cout << "Creating array of ExtTrack objects "<<endl;  //       cout << "Creating array of ExtTrack objects "<<endl;
61         }
62      _trkArray = new TClonesArray("ExtTrack");      _trkArray = new TClonesArray("ExtTrack");
63      //      //
64      _alg_nClFixX   = 3.;     // algorythm parameter: n.hits required on X view      _alg_nClFixX   = 3.;     // algorythm parameter: n.hits required on X view
# Line 70  ExtTrkingAlg::ExtTrkingAlg(Int_t id){ Line 71  ExtTrkingAlg::ExtTrkingAlg(Int_t id){
71      if((int)(id/200)==1)_alg_nViewCal = id - 200;  // algorythm parameter: n.calorimeter planes included      if((int)(id/200)==1)_alg_nViewCal = id - 200;  // algorythm parameter: n.calorimeter planes included
72    
73      //      //
74      _sel_nClstrMAX  = 10. ; // selection parameter: maximum number of cluster      _sel_nClstrMAX  = 15.;//10. ; // selection parameter: maximum number of cluster
75      _sel_nPlaneXMIN = (_alg_nViewCal ? 2 : 3);   // selection parameter: minimum number of hit x-views      _sel_nPlaneXMIN = (_alg_nViewCal ? 2 : 3);   // selection parameter: minimum number of hit x-views
76      _sel_nPlaneYMIN = 2 - (int)(_alg_nViewCal+1)/2 ;   // selection parameter: minimum number of hit y-views      _sel_nPlaneYMIN = 2 - (int)(_alg_nViewCal+1)/2 ;   // selection parameter: minimum number of hit y-views
77    
# Line 95  ExtTrkingAlg::ExtTrkingAlg(Int_t id){ Line 96  ExtTrkingAlg::ExtTrkingAlg(Int_t id){
96      // _extTrack->SetZ(4,ZTRK5);      // _extTrack->SetZ(4,ZTRK5);
97      // _extTrack->SetZ(5,ZTRK6);      // _extTrack->SetZ(5,ZTRK6);
98            
99    
100    
101      _caloStripRoto.clear();      _caloStripRoto.clear();
102      for(int iv=0; iv<_alg_nViewCal; iv++){ //loop over calorimater tracking planes      for(int iv=0; iv<_alg_nViewCal; iv++){ //loop over calorimater tracking planes
103        for(int is=0; is<9; is++){ //loop over calorimeter sensors        for(int is=0; is<9; is++){ //loop over calorimeter sensors
# Line 106  ExtTrkingAlg::ExtTrkingAlg(Int_t id){ Line 109  ExtTrkingAlg::ExtTrkingAlg(Int_t id){
109      //------------------------------------------      //------------------------------------------
110      // read alignment parameters      // read alignment parameters
111      //------------------------------------------      //------------------------------------------
112      cout << "Reading calorimeter alignment parameters"<<endl;      //    printf("qui debug : %i \n",_debug);
113        if(_debug) cout << "Reading calorimeter alignment parameters"<<endl;
114      const char *pamca = gSystem->Getenv("PAM_CALIB");      const char *pamca = gSystem->Getenv("PAM_CALIB");
115      TString filep = "/trk-param/align_param_calo-0/CaloAlignParams.txt";      TString filep = "/trk-param/align_param_calo-0/CaloAlignParams.txt";
116      filep.Prepend(pamca);      filep.Prepend(pamca);
117    
118        if(_debug)cout << " --->  "<<filep.Data()<<endl;
119      std::fstream fs;      std::fstream fs;
120      fs.open (filep.Data(), std::fstream::in );      fs.open (filep.Data(), std::fstream::in );
121      do{  
122        if(!fs.is_open()){
123            
124            cout <<" WARNING!!!!!!!!!!!! error opening file "<<filep.Data()<<endl;
125            cout <<" No alignment parameters loaded "<<endl;
126        }
127    //    do{
128        while(fs.is_open() && !fs.eof()){
129    
130          int view,plane,sensor;          int view,plane,sensor;
131          float par[5],err[5];          float par[5],err[5];
132          fs>>view >> plane>> sensor;          fs>>view >> plane>> sensor;
# Line 123  ExtTrkingAlg::ExtTrkingAlg(Int_t id){ Line 137  ExtTrkingAlg::ExtTrkingAlg(Int_t id){
137    
138    
139              int index = (plane*2+1-view)*9 + sensor;              int index = (plane*2+1-view)*9 + sensor;
140              cout<<index<<" >> ";              if(_debug){
141              cout<<view << " "<<plane<<" "<<sensor<<"          ";                cout<<index<<" >> ";
142              for(int i=0; i<5; i++)cout<<setw(10)<<par[i]<<" "<<setw(10)<<err[i]<<" ";                cout<<view << " "<<plane<<" "<<sensor<<"          ";
143              cout<<endl;                for(int i=0; i<5; i++)cout<<setw(10)<<par[i]<<" "<<setw(10)<<err[i]<<" ";
144                  cout<<endl;
145                }
146    
147              _caloStripRoto[index].SetAligParams__fPitch(par[0]);              _caloStripRoto[index].SetAligParams__fPitch(par[0]);
148              _caloStripRoto[index].SetAligParams__shift(par+1);              _caloStripRoto[index].SetAligParams__shift(par+1);
# Line 135  ExtTrkingAlg::ExtTrkingAlg(Int_t id){ Line 151  ExtTrkingAlg::ExtTrkingAlg(Int_t id){
151  //          _caloStripRoto[index].SetView(view);  //          _caloStripRoto[index].SetView(view);
152  //          _caloStripRoto[index].SetPlane(plane);  //          _caloStripRoto[index].SetPlane(plane);
153    
154          }          };
155      }while(!fs.eof());  //    }while(!fs.eof());
156        };
157      fs.close();      fs.close();
158    
159      // -----------------------------------      // -----------------------------------
# Line 156  ExtTrkingAlg::ExtTrkingAlg(Int_t id){ Line 173  ExtTrkingAlg::ExtTrkingAlg(Int_t id){
173    
174      for(int ip =0 ; ip<_extTrack->nplanes; ip++ )_extTrack->SetZ(ip,_zMech[ip]);      for(int ip =0 ; ip<_extTrack->nplanes; ip++ )_extTrack->SetZ(ip,_zMech[ip]);
175    
176      cout <<" Extended-track Z-coordinates: "<<endl;      if ( _debug ){
177      for(int ip =0 ; ip<_extTrack->nplanes; ip++ )cout <<  _extTrack->zm[ip]<<endl;        cout <<" Extended-track Z-coordinates: "<<endl;
178              for(int ip =0 ; ip<_extTrack->nplanes; ip++ )cout <<  _extTrack->zm[ip]<<endl;
179        }
180    }else{    }else{
181            if ( _debug ){
182      cout << "ExtTrkingAlg(Int_t id) -- algorythm id="<<id<<" not valid "<<endl;        cout << "ExtTrkingAlg(Int_t id) -- algorythm id="<<id<<" not valid "<<endl;
183      cout << "--> Track array not created "<<endl;        cout << "--> Track array not created "<<endl;
184        }
185      _sel_nClstrMAX  = 0.;      _sel_nClstrMAX  = 0.;
186      _sel_nPlaneXMIN = 0.;        _sel_nPlaneXMIN = 0.;  
187      _sel_nPlaneYMIN = 0.;        _sel_nPlaneYMIN = 0.;  
# Line 277  TClonesArray* ExtTrkingAlg::GetTrackArra Line 296  TClonesArray* ExtTrkingAlg::GetTrackArra
296   */   */
297  Bool_t ExtTrkingAlg::CheckEvent(){  Bool_t ExtTrkingAlg::CheckEvent(){
298    
299    if(!_trk_l2)return true;    if(!_trk_l1)return false;
300      TrkLevel1 *trkl1 = _trk_l1;//event->GetTrkLevel1();
   TrkLevel2 *trkl2 = _trk_l2;//event->GetTrkLevel2();  
301    int nClstrMAX  = _sel_nClstrMAX; //maximum number of cluster    int nClstrMAX  = _sel_nClstrMAX; //maximum number of cluster
302    int nPlaneXMIN = _sel_nPlaneXMIN;    /////////////////////////////////////////////////////////////////////
303    int nPlaneYMIN = _sel_nPlaneYMIN;    /// dump selection
304      /////////////////////////////////////////////////////////////////////
305      if(_debug){
306        cout << " n.clusters "<<trkl1->nclstr()<<endl;      
307      }  
308      if(trkl1->nclstr() > nClstrMAX) return false;
309    /////////////////////////////////////////////////////////////////////    /////////////////////////////////////////////////////////////////////
310    /// count number of hit planes    /// count number of hit planes
311    /////////////////////////////////////////////////////////////////////    /////////////////////////////////////////////////////////////////////
312      int nPlaneXMIN = _sel_nPlaneXMIN;
313      int nPlaneYMIN = _sel_nPlaneYMIN;
314    int nHitX[] = {0,0,0,0,0,0};    int nHitX[] = {0,0,0,0,0,0};
315    int nHitY[] = {0,0,0,0,0,0};    int nHitY[] = {0,0,0,0,0,0};
316    for(int ix=0; ix<trkl2->nclsx(); ix++)nHitX[trkl2->GetSingletX(ix)->plane-1]++;        for(int ic=0; ic<trkl1->nclstr(); ic++){      
317    for(int iy=0; iy<trkl2->nclsy(); iy++)nHitY[trkl2->GetSingletY(iy)->plane-1]++;        if      (trkl1->GetCluster(ic)->IsX() ) nHitX[trkl1->GetCluster(ic)->GetPlane()-1]++;    
318          else if (trkl1->GetCluster(ic)->IsY() ) nHitY[trkl1->GetCluster(ic)->GetPlane()-1]++;  
319      }
320    /////////////////////////////////////////////////////////////////////    /////////////////////////////////////////////////////////////////////
321    /// set minimum condition 3x+2y    /// set minimum condition 3x+2y
322    /////////////////////////////////////////////////////////////////////    /////////////////////////////////////////////////////////////////////
# Line 297  Bool_t ExtTrkingAlg::CheckEvent(){ Line 324  Bool_t ExtTrkingAlg::CheckEvent(){
324    for(int ip=0; ip<6; ip++)if(nHitX[ip])nPlaneX++;          for(int ip=0; ip<6; ip++)if(nHitX[ip])nPlaneX++;      
325    int nPlaneY=0;    int nPlaneY=0;
326    for(int ip=0; ip<6; ip++)if(nHitY[ip])nPlaneY++;          for(int ip=0; ip<6; ip++)if(nHitY[ip])nPlaneY++;      
327    
328    /////////////////////////////////////////////////////////////////////    /////////////////////////////////////////////////////////////////////
329    /// dump selection    /// dump selection
330    /////////////////////////////////////////////////////////////////////    /////////////////////////////////////////////////////////////////////
331    if(_debug){    if(_debug){
332      cout << " n.tracks "<<trkl2->GetNTracks()<<endl;        cout << " n.x-hit planes  "<<nPlaneX<<endl;
333      cout << " n.x-singles "<<trkl2->nclsx()<<endl;        cout << " n.y-hit planes  "<<nPlaneY<<endl;
334      cout << " n.y-singles "<<trkl2->nclsy()<<endl;    }  
335      cout << " n.x-planes with singles "<<nPlaneX<<endl;    if( nPlaneX<nPlaneXMIN ) return false;
336      cout << " n.y-planes with singles "<<nPlaneY<<endl;    if( nPlaneY<nPlaneYMIN ) return false;
337    }  
338    
339    //   if(!_trk_l2)return true;
340    //   TrkLevel2 *trkl2 = _trk_l2;//event->GetTrkLevel2();
341    //   /////////////////////////////////////////////////////////////////////
342    //   /// dump selection
343    //   /////////////////////////////////////////////////////////////////////
344    //   if(_debug){
345    //     cout << " n.tracks "<<trkl2->GetNTracks()<<endl;
346    //     cout << " n.x-singles "<<trkl2->nclsx()<<endl;
347    //     cout << " n.y-singles "<<trkl2->nclsy()<<endl;
348    //     cout << " n.y-singles "<<trkl2->nclsy()<<endl;
349    //   }
350    //   if(trkl2->nclsy()+trkl2->nclsx() > nClstrMAX) return false;
351    //   /////////////////////////////////////////////////////////////////////
352    //   /// count number of hit planes
353    //   /////////////////////////////////////////////////////////////////////
354    //   int nPlaneXMIN = _sel_nPlaneXMIN;
355    //   int nPlaneYMIN = _sel_nPlaneYMIN;
356    //   int nHitX[] = {0,0,0,0,0,0};
357    //   int nHitY[] = {0,0,0,0,0,0};
358    //   for(int ix=0; ix<trkl2->nclsx(); ix++)nHitX[trkl2->GetSingletX(ix)->plane-1]++;    
359    //   for(int iy=0; iy<trkl2->nclsy(); iy++)nHitY[trkl2->GetSingletY(iy)->plane-1]++;
360    //   /////////////////////////////////////////////////////////////////////
361    //   /// set minimum condition 3x+2y
362    //   /////////////////////////////////////////////////////////////////////
363    //   int nPlaneX=0;
364    //   for(int ip=0; ip<6; ip++)if(nHitX[ip])nPlaneX++;  
365    //   int nPlaneY=0;
366    //   for(int ip=0; ip<6; ip++)if(nHitY[ip])nPlaneY++;  
367    //   /////////////////////////////////////////////////////////////////////
368    //   /// dump selection
369    //   /////////////////////////////////////////////////////////////////////
370    //   if(_debug){
371    //     cout << " n.x-planes with singles "<<nPlaneX<<endl;
372    //     cout << " n.y-planes with singles "<<nPlaneY<<endl;
373    //   }
374    ////////////////////////////////////////////////////////////////////////    ////////////////////////////////////////////////////////////////////////
375    //    //
376    // check condition for retracking    // check condition for retracking
377    //    //
378    ////////////////////////////////////////////////////////////////////////    ////////////////////////////////////////////////////////////////////////
379    if(trkl2->GetNTracks()==0 && nPlaneX<nPlaneXMIN) return false;  //   if(trkl2->GetNTracks()==0 && nPlaneX<nPlaneXMIN) return false;
380    if(trkl2->GetNTracks()==0 && nPlaneY<nPlaneYMIN) return false;  //   if(trkl2->GetNTracks()==0 && nPlaneY<nPlaneYMIN) return false;
   if(trkl2->nclsy()+trkl2->nclsx() > nClstrMAX) return false;  
381    
382    return true;    return true;
383    
# Line 479  void ExtTrkingAlg::ProcessEvent0(Bool_t Line 542  void ExtTrkingAlg::ProcessEvent0(Bool_t
542    }//end iterations    }//end iterations
543        
544    
545    if( trackCandidates[mapIndex].size() > nTrackMAX ){    if( trackCandidates[mapIndex].size() > (uint)nTrackMAX ){ // EM,  warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
546      if(_debug)cout << "n.candidates "<<trackCandidates[mapIndex].size()<<" > "<<nTrackMAX<<endl;      if(_debug)cout << "n.candidates "<<trackCandidates[mapIndex].size()<<" > "<<nTrackMAX<<endl;
547      return;//to many candidates      return;//to many candidates
548    }    }
# Line 491  void ExtTrkingAlg::ProcessEvent0(Bool_t Line 554  void ExtTrkingAlg::ProcessEvent0(Bool_t
554        cout << " >> ";        cout << " >> ";
555        cout << " X: ";        cout << " X: ";
556        for(int ip=0; ip<6; ip++)        for(int ip=0; ip<6; ip++)
557          if(trackCand.GetClusterX_ID(ip)>=0)cout << trackCand.GetClusterX_ID(ip)<<":"<<trackCand.GetSensor(ip)<<"|";          if(trackCand.GetClusterX_ID(ip)>=0)cout << trackCand.GetClusterX_ID(ip)<<":"<<trackCand.GetSensor(ip)<<"|";
558        cout << " Y: ";        cout << " Y: ";
559        for(int ip=0; ip<6; ip++)        for(int ip=0; ip<6; ip++)
560          if(trackCand.GetClusterY_ID(ip)>=0)cout << trackCand.GetClusterY_ID(ip)<<":"<<trackCand.GetSensor(ip)<<"|";              if(trackCand.GetClusterY_ID(ip)>=0)cout << trackCand.GetClusterY_ID(ip)<<":"<<trackCand.GetSensor(ip)<<"|";    
561        cout << endl;        cout << endl;
562      }      }
563    }    }
# Line 999  void ExtTrkingAlg::FillClusterMap(multim Line 1062  void ExtTrkingAlg::FillClusterMap(multim
1062      Float_t mip = l1->DecodeEstrip(ih,view,plane,strip);      Float_t mip = l1->DecodeEstrip(ih,view,plane,strip);
1063            
1064      // if(strip == 0 || strip == 95)cout <<" strip "<<strip<<" OK "<<endl;      // if(strip == 0 || strip == 95)cout <<" strip "<<strip<<" OK "<<endl;
1065      if(strip < 0 || strip > 95)cout <<" strip "<<strip<<" AHI AHI "<<endl;      if((strip < 0 || strip > 95) && _debug)cout <<" strip "<<strip<<" AHI AHI "<<endl;
1066            
1067      if(view<0)continue;      if(view<0)continue;
1068      if(plane<0)continue;      if(plane<0)continue;
# Line 1293  bool ExtTrkingAlg::EvaluateClusterPositi Line 1356  bool ExtTrkingAlg::EvaluateClusterPositi
1356    //  int plane = 6 + cluster.view;            //0...21    //  int plane = 6 + cluster.view;            //0...21
1357    int plane = (int)(cluster.view / 2 );    //0...21    int plane = (int)(cluster.view / 2 );    //0...21
1358    int ladder = (int)(cluster.coordPU / 32);//0...2    int ladder = (int)(cluster.coordPU / 32);//0...2
1359    int view = cluster.view;//0...43    //  int view = cluster.view;//0...43    // EM, unused variable
1360    bool isX = (cluster.view)%2;    bool isX = (cluster.view)%2;
1361    bool isY = !isX;    bool isY = !isX;
1362        
# Line 1356  bool ExtTrkingAlg::EvaluateClusterPositi Line 1419  bool ExtTrkingAlg::EvaluateClusterPositi
1419  //       cout << " extst.Set("<<(int)isY<<", "<<plane<<", "<<cluster.coordPU<<", "<<sensor<<") "<<endl;  //       cout << " extst.Set("<<(int)isY<<", "<<plane<<", "<<cluster.coordPU<<", "<<sensor<<") "<<endl;
1420  //       extst.Set((int)isY, plane, cluster.coordPU, sensor);  //       extst.Set((int)isY, plane, cluster.coordPU, sensor);
1421    
1422        if( sisensor != extst.GetSiSensor() )cout << " extst.GetSiSensor() "<<extst.GetSiSensor()<<" != "<<sisensor<<endl;        if( sisensor != extst.GetSiSensor() && _debug )cout << " extst.GetSiSensor() "<<extst.GetSiSensor()<<" != "<<sisensor<<endl;
1423        if( plane != extst.GetPlane() )cout << " extst.GetPlane() "<<extst.GetPlane()<<" != "<<plane<<endl;        if( plane != extst.GetPlane() && _debug )cout << " extst.GetPlane() "<<extst.GetPlane()<<" != "<<plane<<endl;
1424        if( (int)isY != extst.GetView() )cout << " extst.GetView() "<<extst.GetView()<<" != "<<isY<<endl;        if( (int)isY != extst.GetView() && _debug )cout << " extst.GetView() "<<extst.GetView()<<" != "<<isY<<endl;
1425  //      cout << " extst.SetStrip("<<cluster.coordPU<<") "<<endl;  //      cout << " extst.SetStrip("<<cluster.coordPU<<") "<<endl;
1426        extst.SetStrip(cluster.coordPU);        extst.SetStrip(cluster.coordPU);
1427                
# Line 1442  bool ExtTrkingAlg::EvaluateClusterPositi Line 1505  bool ExtTrkingAlg::EvaluateClusterPositi
1505            
1506            int view = cluster.view; //0... 43            int view = cluster.view; //0... 43
1507    
1508            int nW = (int)((view + 1)/2);            //      int nW = (int)((view + 1)/2); // EM, unused variable
1509            float dW = 0.74*factor;//X0            float dW = 0.74*factor;//X0
1510            float dW_cm = 0.26*factor;//cm            float dW_cm = 0.26*factor;//cm
1511            float dSi_cm = 0.; //cm            float dSi_cm = 0.; //cm
# Line 1462  bool ExtTrkingAlg::EvaluateClusterPositi Line 1525  bool ExtTrkingAlg::EvaluateClusterPositi
1525            // if(isY)          VMS[nW] = rms*rms;            // if(isY)          VMS[nW] = rms*rms;
1526            // if(_debug)cout <<endl<<view<<" nW "<<nW<<" res(MS) "<<rms<<" factor "<<factor;            // if(_debug)cout <<endl<<view<<" nW "<<nW<<" res(MS) "<<rms<<" factor "<<factor;
1527                
1528            if(view>=3)cout << " --> NB error computation not accurate ";            if(view>=3 && _debug)cout << " --> NB error computation not accurate ";
1529            res_ms = rms;            res_ms = rms;
1530    
1531        }        }
# Line 1472  bool ExtTrkingAlg::EvaluateClusterPositi Line 1535  bool ExtTrkingAlg::EvaluateClusterPositi
1535    
1536    }    }
1537    
1538      return true; // EM, warning: no return statement in function returning non-void [-Wreturn-type]
1539  }  }
1540    
1541    
# Line 1762  void ExtTrkingAlg::ProcessEvent1(Bool_t Line 1825  void ExtTrkingAlg::ProcessEvent1(Bool_t
1825    
1826            }            }
1827    
1828            if( trackCandidates[mapIndex].size() > nTrackMAX ){            if( trackCandidates[mapIndex].size() > (uint)nTrackMAX ){ // EM, compilation warning comparison between signed and unsigned
1829              if(_debug)cout << "n.candidates "<<trackCandidates[mapIndex].size()<<" > "<<nTrackMAX<<endl;              if(_debug)cout << "n.candidates "<<trackCandidates[mapIndex].size()<<" > "<<nTrackMAX<<endl;
1830              return;//to many candidates              return;//to many candidates
1831            }            }
# Line 2242  void ExtTrkingAlg::ProcessEvent1(Bool_t Line 2305  void ExtTrkingAlg::ProcessEvent1(Bool_t
2305        // and        // and
2306        // evaluated other track info        // evaluated other track info
2307        // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~          // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  
2308        float VMS[22];VMS[0]=0.;        //      float VMS[22];VMS[0]=0.; // EM, 2245:13: warning: variable ?VMS? set but not used [-Wunused-but-set-variable]
2309        for(int ip=0; ip<NEXTPLANES; ip++){        for(int ip=0; ip<NEXTPLANES; ip++){
2310    
2311          if(!t_track.XGood(ip)&&!t_track.YGood(ip))continue;          if(!t_track.XGood(ip)&&!t_track.YGood(ip))continue;
# Line 2823  void ExtTrkingAlg::ProcessEvent2(Bool_t Line 2886  void ExtTrkingAlg::ProcessEvent2(Bool_t
2886    
2887                        }                        }
2888    
2889                        if( trackCandidates[mapIndex+1].size() > nTrackMAX ){                        if( trackCandidates[mapIndex+1].size() > (uint)nTrackMAX ){ //EM, 2826:50: warning: comparison between signed and unsigned integer expressions [-Wsign-compare]
2890                            if(_debug)cout << "n.candidates "<<trackCandidates[mapIndex+1].size()<<" > "<<nTrackMAX<<endl;                            if(_debug)cout << "n.candidates "<<trackCandidates[mapIndex+1].size()<<" > "<<nTrackMAX<<endl;
2891                            return;//to many candidates                            return;//to many candidates
2892                        }                        }
# Line 3292  void ExtTrkingAlg::ProcessEvent2(Bool_t Line 3355  void ExtTrkingAlg::ProcessEvent2(Bool_t
3355                // and                // and
3356                // evaluated other track info                // evaluated other track info
3357                // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~                  // ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  
3358                float VMS[22];VMS[0]=0.;                //float VMS[22];VMS[0]=0.;  //EM, 3295:14: warning: variable ?VMS? set but not used [-Wunused-but-set-variable]
3359                for(int ip=0; ip<NEXTPLANES; ip++){                for(int ip=0; ip<NEXTPLANES; ip++){
3360    
3361                    if(!t_track.XGood(ip)&&!t_track.YGood(ip))continue;                    if(!t_track.XGood(ip)&&!t_track.YGood(ip))continue;
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.23