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

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

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

revision 1.39 by pam-fi, Wed Aug 22 07:03:45 2007 UTC revision 1.43 by pam-fi, Tue Jan 22 08:55:07 2008 UTC
# Line 52  TrkTrack::TrkTrack(){ Line 52  TrkTrack::TrkTrack(){
52          ayv[ip]    = 0;          ayv[ip]    = 0;
53          dedx_x[ip] = 0;          dedx_x[ip] = 0;
54          dedx_y[ip] = 0;          dedx_y[ip] = 0;
55            multmaxx[ip] = 0;
56            multmaxy[ip] = 0;
57            seedx[ip]  = 0;  
58            seedy[ip]  = 0;
59            xpu[ip]    = 0;  
60            ypu[ip]    = 0;  
61    
62      };      };
63  //    clx = 0;  
64  //    cly = 0;  //     TrkParams::SetTrackingMode();
65  //    clx = new TRefArray(6,0); //forse causa memory leak???  //     TrkParams::SetPrecisionFactor();
66  //    cly = new TRefArray(6,0); //forse causa memory leak???  //     TrkParams::SetStepMin();
67  //    clx = TRefArray(6,0);      TrkParams::SetMiniDefault();
 //    cly = TRefArray(6,0);  
   
     TrkParams::SetTrackingMode();  
     TrkParams::SetPrecisionFactor();  
     TrkParams::SetStepMin();  
68      TrkParams::SetPFA();      TrkParams::SetPFA();
69    
70  };  };
# Line 96  TrkTrack::TrkTrack(const TrkTrack& t){ Line 98  TrkTrack::TrkTrack(const TrkTrack& t){
98          ayv[ip]    = t.ayv[ip];          ayv[ip]    = t.ayv[ip];
99          dedx_x[ip] = t.dedx_x[ip];          dedx_x[ip] = t.dedx_x[ip];
100          dedx_y[ip] = t.dedx_y[ip];          dedx_y[ip] = t.dedx_y[ip];
101            multmaxx[ip] = t.multmaxx[ip];
102            multmaxy[ip] = t.multmaxy[ip];
103            seedx[ip]    = t.seedx[ip];  
104            seedy[ip]    = t.seedy[ip];
105            xpu[ip]      = t.xpu[ip];  
106            ypu[ip]      = t.ypu[ip];  
107      };      };
108  //    clx = 0;  
109  //    cly = 0;  //     TrkParams::SetTrackingMode();
110  //    if(t.clx)clx = new TRefArray(*(t.clx));  //     TrkParams::SetPrecisionFactor();
111  //    if(t.cly)cly = new TRefArray(*(t.cly));  //     TrkParams::SetStepMin();  
112  //    clx = TRefArray(t.clx);      TrkParams::SetMiniDefault();
 //    cly = TRefArray(t.cly);  
   
     TrkParams::SetTrackingMode();  
     TrkParams::SetPrecisionFactor();  
     TrkParams::SetStepMin();    
113      TrkParams::SetPFA();      TrkParams::SetPFA();
114    
115  };  };
# Line 141  void TrkTrack::Copy(TrkTrack& t){ Line 144  void TrkTrack::Copy(TrkTrack& t){
144          t.ayv[ip]    = ayv[ip];          t.ayv[ip]    = ayv[ip];
145          t.dedx_x[ip] = dedx_x[ip];          t.dedx_x[ip] = dedx_x[ip];
146          t.dedx_y[ip] = dedx_y[ip];          t.dedx_y[ip] = dedx_y[ip];
147            t.multmaxx[ip] = multmaxx[ip];
148            t.multmaxy[ip] = multmaxy[ip];
149            t.seedx[ip]    = seedx[ip];  
150            t.seedy[ip]    = seedy[ip];
151            t.xpu[ip]      = xpu[ip];  
152            t.ypu[ip]      = ypu[ip];  
153                            
154      };      };
155    
 //    t.clx = TRefArray(clx);  
 //    t.cly = TRefArray(cly);  
156            
157  };  };
158  //--------------------------------------  //--------------------------------------
# Line 492  void TrkTrack::Dump(){ Line 499  void TrkTrack::Dump(){
499      cout << endl << "           "; for(int i=0; i<5; i++)cout << coval[4][i]<<" ";      cout << endl << "           "; for(int i=0; i<5; i++)cout << coval[4][i]<<" ";
500      cout << endl << "dedx_x   : "; for(int i=0; i<6; i++)cout << dedx_x[i] << " ";      cout << endl << "dedx_x   : "; for(int i=0; i<6; i++)cout << dedx_x[i] << " ";
501      cout << endl << "dedx_y   : "; for(int i=0; i<6; i++)cout << dedx_y[i] << " ";      cout << endl << "dedx_y   : "; for(int i=0; i<6; i++)cout << dedx_y[i] << " ";
502        cout << endl << "maxs x   : "; for(int i=0; i<6; i++)cout << GetClusterX_MaxStrip(i) << " ";
503        cout << endl << "maxs y   : "; for(int i=0; i<6; i++)cout << GetClusterY_MaxStrip(i) << " ";
504        cout << endl << "mult x   : "; for(int i=0; i<6; i++)cout << GetClusterX_Multiplicity(i) << " ";
505        cout << endl << "mult y   : "; for(int i=0; i<6; i++)cout << GetClusterY_Multiplicity(i) << " ";
506        cout << endl << "seed x   : "; for(int i=0; i<6; i++)cout << GetClusterX_Seed(i) << " ";
507        cout << endl << "seed y   : "; for(int i=0; i<6; i++)cout << GetClusterY_Seed(i) << " ";
508        cout << endl << "xpu      : "; for(int i=0; i<6; i++)cout << xpu[i] << " ";
509        cout << endl << "ypu      : "; for(int i=0; i<6; i++)cout << ypu[i] << " ";
510    
511      cout << endl;      cout << endl;
512  }  }
513  /**  /**
# Line 607  void TrkTrack::LoadField(TString path){ Line 623  void TrkTrack::LoadField(TString path){
623  //     path_.error   = 0;  //     path_.error   = 0;
624  //     readb_();  //     readb_();
625    
626      TrkParams::SetTrackingMode();  //     TrkParams::SetTrackingMode();
627      TrkParams::SetPrecisionFactor();  //     TrkParams::SetPrecisionFactor();
628      TrkParams::SetStepMin();  //     TrkParams::SetStepMin();
629        TrkParams::SetMiniDefault();
630    
631      TrkParams::Set(path,1);      TrkParams::Set(path,1);
632      TrkParams::Load(1);      TrkParams::Load(1);
# Line 710  void TrkTrack::SetFromMiniStruct(cMini2t Line 727  void TrkTrack::SetFromMiniStruct(cMini2t
727   */   */
728  Bool_t TrkTrack::EvaluateClusterPositions(){  Bool_t TrkTrack::EvaluateClusterPositions(){
729            
730  //     cout << "void TrkTrack::GetClusterPositions() "<<endl;  //     cout << "void TrkTrack::GetClusterositions() "<<endl;
731    
732      TrkParams::Load( );      TrkParams::Load( );
733      if( !TrkParams::IsLoaded() )return false;      if( !TrkParams::IsLoaded() )return false;
# Line 758  Bool_t TrkTrack::EvaluateClusterPosition Line 775  Bool_t TrkTrack::EvaluateClusterPosition
775   * @see EvaluateClusterPositions()   * @see EvaluateClusterPositions()
776   *   *
777   * The fitting procedure can be varied by changing the tracking mode,   * The fitting procedure can be varied by changing the tracking mode,
778   * the fit-precision factor and the minimum number of step.   * the fit-precision factor, the minimum number of step, etc.
779   * @see SetTrackingMode(int)   * @see SetTrackingMode(int)
780   * @see SetPrecisionFactor(double)   * @see SetPrecisionFactor(double)
781   * @see SetStepMin(int)   * @see SetStepMin(int)
782     * @see SetDeltaB(int,double)
783   */   */
784  void TrkTrack::Fit(double pfixed, int& fail, int iprint, int froml1){  void TrkTrack::Fit(double pfixed, int& fail, int iprint, int froml1){
785    
# Line 864  void TrkTrack::SetStepMin(int istepmin){ Line 882  void TrkTrack::SetStepMin(int istepmin){
882      track_.istepmin = istepmin;      track_.istepmin = istepmin;
883  }  }
884  /**  /**
885     * Set deltaB parameters (id=0,1). By default they are set to zero.
886     */
887    void TrkTrack::SetDeltaB(int id, double db){
888        if(id!=0 && id!=1)cout << "void TrkTrack::SetDeltaB(int id,double db) -- wrong input parameters: "<<id<<" "<<db<<endl;
889        TrkParams::SetDeltaB(id,db);
890    }
891    
892    /**
893   * Returns 1 if the track is inside the magnet cavity   * Returns 1 if the track is inside the magnet cavity
894   * Set the minimum number of steps for tracking precision   * Set the minimum number of steps for tracking precision
895   */   */
896  Bool_t TrkTrack::IsInsideCavity(){  Bool_t TrkTrack::IsInsideCavity(){
897      float xmagntop, ymagntop, xmagnbottom, ymagnbottom;      float xmagntop, ymagntop, xmagnbottom, ymagnbottom;
898      xmagntop = xv[0] + (ZMAGNHIGH-zv[0])*tan(cos(-1.0)*axv[0]/180.);      xmagntop = xv[0] + (ZMAGNHIGH-zv[0])*tan(acos(-1.0)*axv[0]/180.);
899      ymagntop = yv[0] + (ZMAGNHIGH-zv[0])*tan(cos(-1.0)*ayv[0]/180.);      ymagntop = yv[0] + (ZMAGNHIGH-zv[0])*tan(acos(-1.0)*ayv[0]/180.);
900      xmagnbottom = xv[5] + (ZMAGNLOW-zv[5])*tan(cos(-1.0)*axv[5]/180.);      xmagnbottom = xv[5] + (ZMAGNLOW-zv[5])*tan(acos(-1.0)*axv[5]/180.);
901      ymagnbottom = yv[5] + (ZMAGNLOW-zv[5])*tan(cos(-1.0)*ayv[5]/180.);      ymagnbottom = yv[5] + (ZMAGNLOW-zv[5])*tan(acos(-1.0)*ayv[5]/180.);
902      if( xmagntop>XMAGNLOW && xmagntop<XMAGNHIGH &&      if( xmagntop>XMAGNLOW && xmagntop<XMAGNHIGH &&
903          ymagntop>YMAGNLOW && ymagntop<YMAGNHIGH &&          ymagntop>YMAGNLOW && ymagntop<YMAGNHIGH &&
904          xmagnbottom>XMAGNLOW && xmagnbottom<XMAGNHIGH &&          xmagnbottom>XMAGNLOW && xmagnbottom<XMAGNHIGH &&
# Line 895  Int_t TrkTrack::GetClusterX_ID(int ip){ Line 921  Int_t TrkTrack::GetClusterX_ID(int ip){
921  Int_t TrkTrack::GetClusterY_ID(int ip){  Int_t TrkTrack::GetClusterY_ID(int ip){
922      return ((Int_t)fabs(ygood[ip]))%10000000-1;      return ((Int_t)fabs(ygood[ip]))%10000000-1;
923  };  };
924    
925  /**  /**
926   * Method to retrieve the ladder (0-4, increasing x) traversed by the track on this plane.   * Method to retrieve the ladder (0-4, increasing x) traversed by the track on this plane.
927   * If no ladder is traversed (dead area) the metod retuns -1.   * If no ladder is traversed (dead area) the metod retuns -1.
# Line 943  void TrkTrack::SetYGood(int ip, int clid Line 970  void TrkTrack::SetYGood(int ip, int clid
970      if(bad)ygood[ip]=-ygood[ip];      if(bad)ygood[ip]=-ygood[ip];
971  };  };
972    
973    /**
974     * \brief Average X
975     * Average value of <xv>, evaluated from the first to the last hit x view.
976     */
977    Float_t TrkTrack::GetXav(){
978    
979        int first_plane = -1;
980        int last_plane  = -1;
981        for(Int_t ip=0; ip<6; ip++){
982            if( XGood(ip) && first_plane == -1 )first_plane = ip;
983            if( XGood(ip) && first_plane != -1 )last_plane = ip;
984        }
985        if( first_plane == -1 || last_plane == -1){
986            return -100;
987        }
988        if( last_plane-first_plane+1 ==0 )return -100;
989    
990        Float_t av = 0;    
991        for(int ip=first_plane; ip<=last_plane; ip++)av+=xv[ip];
992        
993        return (av/(last_plane-first_plane+1));    
994    }
995    /**
996     * \brief Average Y
997     * Average value of <yv>, evaluated from the first to the last hit x view.
998     */
999    Float_t TrkTrack::GetYav(){
1000    
1001        int first_plane = -1;
1002        int last_plane  = -1;
1003        for(Int_t ip=0; ip<6; ip++){
1004            if( XGood(ip) && first_plane == -1 )first_plane = ip;
1005            if( XGood(ip) && first_plane != -1 )last_plane = ip;
1006        }
1007        if( first_plane == -1 || last_plane == -1){
1008            return -100;
1009        }
1010        if( last_plane-first_plane+1 ==0 )return -100;
1011    
1012        Float_t av = 0;    
1013        for(int ip=first_plane; ip<=last_plane; ip++)av+=yv[ip];
1014    
1015        return (av/(last_plane-first_plane+1));    
1016    }
1017    /**
1018     * \brief Average Z
1019     * Average value of <zv>, evaluated from the first to the last hit x view.
1020     */
1021    Float_t TrkTrack::GetZav(){
1022    
1023        int first_plane = -1;
1024        int last_plane  = -1;
1025        for(Int_t ip=0; ip<6; ip++){
1026            if( XGood(ip) && first_plane == -1 )first_plane = ip;
1027            if( XGood(ip) && first_plane != -1 )last_plane = ip;
1028        }
1029        if( first_plane == -1 || last_plane == -1){
1030            return -100;
1031        }
1032        if( last_plane-first_plane+1 ==0 )return -100;
1033    
1034        Float_t av = 0;    
1035        for(int ip=first_plane; ip<=last_plane; ip++)av+=zv[ip];
1036        
1037        return (av/(last_plane-first_plane+1));    
1038    }
1039    
1040    /**
1041     * \brief Number of column traversed
1042     */
1043    Int_t TrkTrack::GetNColumns(){
1044        int sensors[] = {0,0,0,0,0,0};
1045        for(int ip=0; ip<6; ip++){
1046            int sensorid = GetLadder(ip)+3*GetSensor(ip);    
1047            if(XGood(ip)||YGood(ip))
1048                if(sensorid>=0 && sensorid<6)sensors[sensorid]=1;
1049        }
1050        int nsensors=0;
1051        for(int is=0; is<6; is++)nsensors += sensors[is];
1052        return nsensors;
1053    };
1054    /**
1055     * \brief Give the maximum energy release
1056     */
1057    Float_t TrkTrack::GetDEDX_max(int ip, int iv){
1058        Float_t max=0;
1059        int pfrom = 0;
1060        int pto   = 6;
1061        int vfrom = 0;
1062        int vto   = 2;
1063        if(ip>=0&&ip<6){
1064            pfrom = ip;
1065            pto   = ip+1;
1066        }
1067        if(iv>=0&&iv<2){
1068            vfrom = iv;
1069            vto   = iv+1;
1070        }
1071        for(int i=pfrom; i<pto; i++)
1072            for(int j=0; j<vto; j++)
1073                if(GetDEDX(i,j)>max)max=GetDEDX(i,j);
1074    
1075        return max;
1076    
1077    };
1078    
1079    /**
1080     * \brief Give the minimum energy release
1081     */
1082    Float_t TrkTrack::GetDEDX_min(int ip, int iv){
1083        Float_t min=100000000;
1084        int pfrom = 0;
1085        int pto   = 6;
1086        int vfrom = 0;
1087        int vto   = 2;
1088        if(ip>=0&&ip<6){
1089            pfrom = ip;
1090            pto   = ip+1;
1091        }
1092        if(iv>=0&&iv<2){
1093            vfrom = iv;
1094            vto   = iv+1;
1095        }
1096        for(int i=pfrom; i<pto; i++)
1097            for(int j=0; j<vto; j++)
1098                if(GetDEDX(i,j)<min)min=GetDEDX(i,j);
1099    
1100        return min;
1101    
1102    };
1103    
1104    /**
1105     * \brief Give the maximum spatial residual release
1106     */
1107    Float_t TrkTrack::GetResidual_max(int ip, int iv){
1108        Float_t max=0;
1109        int pfrom = 0;
1110        int pto   = 6;
1111        int vfrom = 0;
1112        int vto   = 2;
1113        if(ip>=0&&ip<6){
1114            pfrom = ip;
1115            pto   = ip+1;
1116        }
1117        if(iv>=0&&iv<2){
1118            vfrom = iv;
1119            vto   = iv+1;
1120        }
1121        for(int i=pfrom; i<pto; i++){
1122            for(int j=0; j<vto; j++){
1123                if(j==0 && XGood(i) && fabs(xm[i]-xv[i])>fabs(max))max=xv[i]-xm[i];
1124                if(j==1 && YGood(i) && fabs(ym[i]-yv[i])>fabs(max))max=yv[i]-ym[i];
1125            }
1126        }
1127        return max;
1128    
1129    };
1130    
1131    
1132    /**
1133     * \brief Give the maximum multiplicity on the x view
1134     */
1135    Int_t TrkTrack::GetClusterX_Multiplicity_max(){
1136        int max=0;
1137        for(int ip=0; ip<6; ip++)
1138            if(GetClusterX_Multiplicity(ip)>max)max=GetClusterX_Multiplicity(ip);
1139        return max;
1140    };
1141    /**
1142     * \brief Give the minimum multiplicity on the x view
1143     */
1144    Int_t TrkTrack::GetClusterX_Multiplicity_min(){
1145        int min=50;
1146        for(int ip=0; ip<6; ip++)
1147            if(GetClusterX_Multiplicity(ip)<min)min=GetClusterX_Multiplicity(ip);
1148        return min;
1149    };
1150    /**
1151     * \brief Give the maximum multiplicity on the x view
1152     */
1153    Int_t TrkTrack::GetClusterY_Multiplicity_max(){
1154        int max=0;
1155        for(int ip=0; ip<6; ip++)
1156            if(GetClusterY_Multiplicity(ip)>max)max=GetClusterY_Multiplicity(ip);
1157        return max;
1158    };
1159    /**
1160     * \brief Give the minimum multiplicity on the x view
1161     */
1162    Int_t TrkTrack::GetClusterY_Multiplicity_min(){
1163        int min=50;
1164        for(int ip=0; ip<6; ip++)
1165            if(GetClusterY_Multiplicity(ip)<min)min=GetClusterY_Multiplicity(ip);
1166        return min;
1167    };
1168    
1169    /**
1170     * \brief Give the minimum seed on the x view
1171     */
1172    Float_t TrkTrack::GetClusterX_Seed_min(){
1173        Float_t min=100000;
1174        for(int ip=0; ip<6; ip++)
1175            if(XGood(ip) && GetClusterX_Seed(ip)<min)min=GetClusterX_Seed(ip);
1176        return min;
1177    };
1178    /**
1179     * \brief Give the minimum seed on the x view
1180     */
1181    Float_t TrkTrack::GetClusterY_Seed_min(){
1182        Float_t min=100000;
1183        for(int ip=0; ip<6; ip++)
1184            if(YGood(ip) && GetClusterY_Seed(ip)<min)min=GetClusterY_Seed(ip);
1185        return min;
1186    };
1187    
1188    
1189  //--------------------------------------  //--------------------------------------
1190  //  //
1191  //  //
# Line 1078  void TrkLevel2::Dump(){ Line 1321  void TrkLevel2::Dump(){
1321                    
1322          //          //
1323      cout << endl << endl << "=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-";      cout << endl << endl << "=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-";
1324      cout << endl << "good     : "; for(int i=0; i<12; i++) cout << good[i]<<" ";      cout << endl << "good     : "; for(int i=0; i<12; i++) cout << hex <<" 0x"<< good[i]<<dec;
1325      cout << endl << "ntrk()   : " << this->ntrk() ;      cout << endl << "ntrk()   : " << ntrk() ;
1326      cout << endl << "nclsx()  : " << this->nclsx();      cout << endl << "nclsx()  : " << nclsx();
1327      cout << endl << "nclsy()  : " << this->nclsy();      cout << endl << "nclsy()  : " << nclsy();
1328      if(Track){      if(Track){
1329          TClonesArray &t  = *Track;          TClonesArray &t  = *Track;
1330          for(int i=0; i<ntrk(); i++)     ((TrkTrack *)t[i])->Dump();          for(int i=0; i<ntrk(); i++)     ((TrkTrack *)t[i])->Dump();
1331      }            }      
1332      if(SingletX){  //     if(SingletX){
1333          TClonesArray &sx = *SingletX;  //      TClonesArray &sx = *SingletX;
1334          for(int i=0; i<nclsx(); i++) ((TrkSinglet *)sx[i])->Dump();  //      for(int i=0; i<nclsx(); i++) ((TrkSinglet *)sx[i])->Dump();
1335      }  //     }
1336      if(SingletY){  //     if(SingletY){
1337          TClonesArray &sy = *SingletY;  //      TClonesArray &sy = *SingletY;
1338          for(int i=0; i<nclsy(); i++) ((TrkSinglet *)sy[i])->Dump();  //      for(int i=0; i<nclsy(); i++) ((TrkSinglet *)sy[i])->Dump();
1339      }  //     }
1340        cout << endl;
1341  }  }
1342  /**  /**
1343   * \brief Dump processing status   * \brief Dump processing status
# Line 1237  void TrkLevel2::SetFromLevel2Struct(cTrk Line 1481  void TrkLevel2::SetFromLevel2Struct(cTrk
1481              t_track->ayv[ip]    = l2->ayv_nt[i][ip];              t_track->ayv[ip]    = l2->ayv_nt[i][ip];
1482              t_track->dedx_x[ip] = l2->dedx_x[i][ip];              t_track->dedx_x[ip] = l2->dedx_x[i][ip];
1483              t_track->dedx_y[ip] = l2->dedx_y[i][ip];              t_track->dedx_y[ip] = l2->dedx_y[i][ip];
1484                t_track->multmaxx[ip] = l2->multmaxx[i][ip];
1485                t_track->multmaxy[ip] = l2->multmaxy[i][ip];
1486                t_track->seedx[ip]  = l2->seedx[i][ip];  
1487                t_track->seedy[ip]  = l2->seedy[i][ip];
1488                t_track->xpu[ip]    = l2->xpu[i][ip];  
1489                t_track->ypu[ip]    = l2->ypu[i][ip];  
1490              //-----------------------------------------------------              //-----------------------------------------------------
1491              //-----------------------------------------------------              //-----------------------------------------------------
1492              //-----------------------------------------------------              //-----------------------------------------------------
# Line 1451  TRefArray *TrkLevel2::GetTracks_NFitSort Line 1701  TRefArray *TrkLevel2::GetTracks_NFitSort
1701  TrkTrack *TrkLevel2::GetStoredTrack(int is){  TrkTrack *TrkLevel2::GetStoredTrack(int is){
1702    
1703      if(is >= this->ntrk()){      if(is >= this->ntrk()){
1704          cout << "** TrkLevel2 ** Track "<< is << "doen not exits! " << endl;          cout << "TrkTrack *TrkLevel2::GetStoredTrack(int) >> Track "<< is << "doen not exits! " << endl;
1705          cout << "                Stored tracks ntrk() = "<< this->ntrk() << endl;          cout << "Stored tracks ntrk() = "<< this->ntrk() << endl;
1706          return 0;          return 0;
1707      }      }
1708      if(!Track){      if(!Track){
# Line 1473  TrkTrack *TrkLevel2::GetStoredTrack(int Line 1723  TrkTrack *TrkLevel2::GetStoredTrack(int
1723  TrkSinglet *TrkLevel2::GetSingletX(int is){  TrkSinglet *TrkLevel2::GetSingletX(int is){
1724    
1725          if(is >= this->nclsx()){          if(is >= this->nclsx()){
1726                  cout << "** TrkLevel2 ** Singlet "<< is << "doen not exits! " << endl;                  cout << "TrkSinglet *TrkLevel2::GetSingletX(int) >> Singlet "<< is << "doen not exits! " << endl;
1727                  cout << "                Stored x-singlets nclsx() = "<< this->nclsx() << endl;                  cout << "Stored x-singlets nclsx() = "<< this->nclsx() << endl;
1728                  return 0;                  return 0;
1729          }          }
1730          if(!SingletX)return 0;          if(!SingletX)return 0;
# Line 1493  TrkSinglet *TrkLevel2::GetSingletX(int i Line 1743  TrkSinglet *TrkLevel2::GetSingletX(int i
1743  TrkSinglet *TrkLevel2::GetSingletY(int is){  TrkSinglet *TrkLevel2::GetSingletY(int is){
1744    
1745          if(is >= this->nclsy()){          if(is >= this->nclsy()){
1746                  cout << "** TrkLevel2 ** Singlet "<< is << "doen not exits! " << endl;                  cout << "TrkSinglet *TrkLevel2::GetSingletY(int) >> Singlet "<< is << "doen not exits! " << endl;
1747                  cout << "                Stored y-singlets nclsy() = "<< this->nclsx() << endl;                  cout << "Stored y-singlets nclsx() = "<< this->nclsx() << endl;
1748                  return 0;                  return 0;
1749          }          }
1750          if(!SingletY)return 0;          if(!SingletY)return 0;
# Line 1514  TrkSinglet *TrkLevel2::GetSingletY(int i Line 1764  TrkSinglet *TrkLevel2::GetSingletY(int i
1764  TrkTrack *TrkLevel2::GetTrack(int it){  TrkTrack *TrkLevel2::GetTrack(int it){
1765            
1766          if(it >= this->GetNTracks()){          if(it >= this->GetNTracks()){
1767                  cout << "** TrkLevel2 ** Track "<< it << "does not exits! " << endl;                  cout << "TrkTrack *TrkLevel2::GetTrack(int) >> Track "<< it << "does not exits! " << endl;
1768                  cout << "                Physical tracks GetNTracks() = "<< this->ntrk() << endl;                  cout << "Physical tracks GetNTracks() = "<< this->ntrk() << endl;
1769                  return 0;                  return 0;
1770          }          }
1771                    
# Line 1552  Int_t TrkLevel2::GetNTracks(){ Line 1802  Int_t TrkLevel2::GetNTracks(){
1802  TrkTrack *TrkLevel2::GetTrackImage(int it){  TrkTrack *TrkLevel2::GetTrackImage(int it){
1803    
1804      if(it >= this->GetNTracks()){      if(it >= this->GetNTracks()){
1805          cout << "** TrkLevel2 ** Track "<< it << "does not exits! " << endl;          cout << "TrkTrack *TrkLevel2::GetTrackImage(int) >> Track "<< it << "does not exits! " << endl;
1806          cout << "                Physical tracks GetNTracks() = "<< this->ntrk() << endl;          cout << "Physical tracks GetNTracks() = "<< this->ntrk() << endl;
1807          return 0;          return 0;
1808      }      }
1809                    
# Line 1562  TrkTrack *TrkLevel2::GetTrackImage(int i Line 1812  TrkTrack *TrkLevel2::GetTrackImage(int i
1812      TrkTrack *track = (TrkTrack*)sorted->At(it);      TrkTrack *track = (TrkTrack*)sorted->At(it);
1813                    
1814      if(!track->HasImage()){      if(!track->HasImage()){
1815          cout << "** TrkLevel2 ** Track "<< it << "does not have image! " << endl;          cout << "TrkTrack *TrkLevel2::GetTrackImage(int) >> Track "<< it << "does not have image! " << endl;
1816          return 0;          return 0;
1817      }      }
1818      if(!Track)return 0;      if(!Track)return 0;
# Line 1589  void TrkLevel2::LoadField(TString path){ Line 1839  void TrkLevel2::LoadField(TString path){
1839  //     path_.error   = 0;  //     path_.error   = 0;
1840  //     readb_();  //     readb_();
1841    
1842      TrkParams::SetTrackingMode();  //     TrkParams::SetTrackingMode();
1843      TrkParams::SetPrecisionFactor();  //     TrkParams::SetPrecisionFactor();
1844      TrkParams::SetStepMin();  //     TrkParams::SetStepMin();
1845        TrkParams::SetMiniDefault();
1846    
1847      TrkParams::Set(path,1);      TrkParams::Set(path,1);
1848      TrkParams::Load(1);      TrkParams::Load(1);
Legend:
Removed from v.1.39  
changed lines
  Added in v.1.43
  ViewVC Help
Powered by ViewVC 1.1.23