/[PAMELA software]/DarthVader/ToFLevel2/src/ToFLevel2.cpp
ViewVC logotype

Diff of /DarthVader/ToFLevel2/src/ToFLevel2.cpp

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

revision 1.42 by pam-fi, Fri Oct 17 07:54:17 2014 UTC revision 1.48 by pam-fi, Wed May 20 10:57:02 2015 UTC
# Line 147  void ToFLevel2::Delete(Option_t *t){ //E Line 147  void ToFLevel2::Delete(Option_t *t){ //E
147    //    //
148  }; //ELENA  }; //ELENA
149    
150    /**
151     * Retrieves the itrk-th tof track stored in the array
152     * @param itrk Array index (min 0, max ToFLevel2::ntrk())
153     *
154     */
155  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t itrk){  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t itrk){
156    //        //    
157    if(itrk >= ntrk()){    if(itrk >= ntrk()){
# Line 163  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t Line 168  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t
168    
169  /**  /**
170   * Retrieves the tof track matching the seqno-th tracker stored track.   * Retrieves the tof track matching the seqno-th tracker stored track.
171     * @param seqno Track sequential number
172     * (seqno = -1 for standalone tof track, seqno=0-TrkLevel2::ntrk() for tof tracks associated to a tracker track)
173   *   *
174   */   */
175  ToFTrkVar *ToFLevel2::GetToFStoredTrack(int seqno){  ToFTrkVar *ToFLevel2::GetToFStoredTrack(int seqno){
# Line 252  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit) Line 259  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
259         false) return true;         false) return true;
260      else return false;      else return false;
261  };  };
262    
263  /**  /**
264   * Method to get the number of hit paddles on a ToF plane.   * Strict method to get the number of hit paddles on a ToF plane.
265     * The method uses "HitPaddle" which checks if there is a TDC signal
266     * from both PMTs.
267   * @param plane Plane ID (11, 12, 21, 22, 31, 32) or Plane index (0,1,2,3,4,5).   * @param plane Plane ID (11, 12, 21, 22, 31, 32) or Plane index (0,1,2,3,4,5).
268   */   */
269  Int_t ToFLevel2::GetNHitPaddles(Int_t plane){  Int_t ToFLevel2::GetNHitPaddles(Int_t plane){
# Line 263  Int_t ToFLevel2::GetNHitPaddles(Int_t pl Line 273  Int_t ToFLevel2::GetNHitPaddles(Int_t pl
273  };  };
274    
275  /**  /**
276   * Method to get the number of hit paddles on a ToF plane.   * Optional method to get the number of hit paddles on a ToF plane.
277     * The method does NOT check if there is a signal from both PMTs, it only
278     * checks if there is some PMT signal in a paddle
279   * @param plane Plane ID (11, 12, 21, 22, 31, 32) or Plane index (0,1,2,3,4,5).   * @param plane Plane ID (11, 12, 21, 22, 31, 32) or Plane index (0,1,2,3,4,5).
280   */   */
281  Int_t ToFLevel2::GetTrueNHitPaddles(Int_t plane){  Int_t ToFLevel2::GetTrueNHitPaddles(Int_t plane){
# Line 283  Int_t ToFLevel2::GetTrueNHitPaddles(Int_ Line 295  Int_t ToFLevel2::GetTrueNHitPaddles(Int_
295      return npad;      return npad;
296  };  };
297    
298    //new, wm Feb 15
299  //wm Nov 08  //wm Nov 08
300  //gf Apr 07  //gf Apr 07
301  /**  /**
302   * Method to get the mean dEdx from a ToF layer - ATTENTION:   * Method to get the mean dEdx from a ToF layer
303   * It will sum up the dEdx of all the paddles, but since by definition   * By definition there should be PMTs with dEdx values only in one paddle of a layer
304   * only the paddle hitted by the track gets a dEdx value and the other   * (the paddle hitted by the track), this method looks for the hitted paddle
305   * paddles are set to zero, the output is just the dEdx of the hitted   * and gives the mean dEdx of that paddle as the output
306   * paddle in each layer!   * The method was modified for the "ToF-standalone" part in february 2015
307   * The "adcfl" option is not very useful (an artificial dEdx is per   * The "adcfl" option is not very useful (an artificial dEdx is per
308   * definition= 1 mip and not a real measurement), anyway left in the code   * definition= 1 mip and not a real measurement), anyway left in the code
309   * @param notrack Track Number   * @param notrack Track Number (array index, ranging from 0 to ntrk())
310   * @param plane Plane index (0,1,2,3,4,5)   * @param plane Plane index (0,1,2,3,4,5)
311   * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )   * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
312   */   */
313  Float_t ToFLevel2::GetdEdx(ToFTrkVar *trk, Int_t plane, Int_t adcfl){  Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){
314    //  printf("fiffi\n");      ToFTrkVar *trk = GetToFTrkVar(notrack);
315    Float_t dedx = 0.;  //    ToFTrkVar *trk = GetToFStoredTrack(seqno);//Elena 2015
316    Float_t PadEdx =0.;    return this->GetdEdx(trk, plane, adcfl);
317    Int_t SatWarning;  }
   Int_t pad=-1;  
   //  
   if(!trk) return 0; //ELENA  
   //  
   if ( trk->trkseqno == -1 ){ //standalone, only paddles along the track, or about...  
     //    printf("ciccio\n");  
     Float_t xleft=0;  
     Float_t xright=0;  
     Float_t yleft=0;  
     Float_t yright=0;  
     Float_t xtof_temp[6]={100.,100.,100.,100.,100.,100.};  
     Float_t ytof_temp[6]={100.,100.,100.,100.,100.,100.};  
     ToFTrkVar *t_tof = trk;  
       
     xtof_temp[0]=t_tof->xtofpos[0];  
     ytof_temp[0]=t_tof->ytofpos[0];  
     xtof_temp[1]=t_tof->xtofpos[0];  
     ytof_temp[1]=t_tof->ytofpos[0];  
   
     xtof_temp[2]=t_tof->xtofpos[1];  
     ytof_temp[2]=t_tof->ytofpos[1];  
     xtof_temp[3]=t_tof->xtofpos[1];  
     ytof_temp[3]=t_tof->ytofpos[1];  
   
     xtof_temp[4]=t_tof->xtofpos[2];  
     ytof_temp[4]=t_tof->ytofpos[2];  
     xtof_temp[5]=t_tof->xtofpos[2];  
     ytof_temp[5]=t_tof->ytofpos[2];  
   
     if(t_tof->xtofpos[0]<100. && t_tof->ytofpos[0]<100.){  
       xtof_temp[1]=t_tof->xtofpos[0];  
       ytof_temp[0]=t_tof->ytofpos[0];  
     }else if(t_tof->xtofpos[0]>=100. && t_tof->ytofpos[0]<100.){  
       ytof_temp[0]=t_tof->ytofpos[0];  
       this->GetPaddleGeometry(0,(Int_t)log2(this->tof_j_flag[0]),xleft, xright, yleft, yright);  
       xtof_temp[1]=xleft+2.55;  
     }else if(t_tof->ytofpos[0]>=100. && t_tof->xtofpos[0]<100.){  
       xtof_temp[1]=t_tof->xtofpos[0];  
       this->GetPaddleGeometry(1,(Int_t)log2(this->tof_j_flag[1]),xleft, xright, yleft, yright);  
       ytof_temp[0]=yleft+2.75;  
     }  
       
     if(t_tof->xtofpos[1]<100. && t_tof->ytofpos[1]<100.){  
       xtof_temp[2]=t_tof->xtofpos[1];  
       ytof_temp[3]=t_tof->ytofpos[1];  
     }else if(t_tof->xtofpos[1]>=100. && t_tof->ytofpos[1]<100.){  
       ytof_temp[3]=t_tof->ytofpos[1];  
       this->GetPaddleGeometry(3,(Int_t)log2(this->tof_j_flag[3]),xleft, xright, yleft, yright);  
       xtof_temp[2]=xleft+4.5;  
     }else if(t_tof->ytofpos[1]>=100. && t_tof->xtofpos[1]<100.){  
       xtof_temp[2]=t_tof->xtofpos[1];  
       this->GetPaddleGeometry(2,(Int_t)log2(this->tof_j_flag[2]),xleft, xright, yleft, yright);  
       ytof_temp[3]=yleft+3.75;  
     }  
       
     if(t_tof->xtofpos[2]<100. && t_tof->ytofpos[2]<100.){  
       xtof_temp[5]=t_tof->xtofpos[2];  
       ytof_temp[4]=t_tof->ytofpos[2];  
     }else if(t_tof->xtofpos[2]>=100. && t_tof->ytofpos[2]<100.){  
       ytof_temp[4]=t_tof->ytofpos[2];  
       this->GetPaddleGeometry(4,(Int_t)log2(this->tof_j_flag[4]),xleft, xright, yleft, yright);  
       xtof_temp[5]=xleft+3;  
     }else if(t_tof->ytofpos[2]>=100. && t_tof->xtofpos[2]<100.){  
       xtof_temp[5]=t_tof->xtofpos[2];  
       this->GetPaddleGeometry(5,(Int_t)log2(this->tof_j_flag[5]),xleft, xright, yleft, yright);  
       ytof_temp[4]=yleft+2.5;  
     }  
   
     if ( (xtof_temp[0])> 17.85 && fabs(xtof_temp[0])<17.85+4. ) xtof_temp[0] =17.84;  
     if ( (xtof_temp[1])> 17.85 && fabs(xtof_temp[1])<17.85+4. ) xtof_temp[1] =17.84;  
     if ( (ytof_temp[0])> 13.75 && fabs(ytof_temp[0])<13.75+4. ) ytof_temp[0] =13.74;  
     if ( (ytof_temp[1])> 13.75 && fabs(ytof_temp[1])<13.75+4. ) ytof_temp[1] =13.74;  
     if ( (xtof_temp[0])< -17.85 && fabs(xtof_temp[0])>-(17.85+4.) ) xtof_temp[0] =-17.84;  
     if ( (xtof_temp[1])< -17.85 && fabs(xtof_temp[1])>-(17.85+4.) ) xtof_temp[1] =-17.84;  
     if ( (ytof_temp[0])< -13.75 && fabs(ytof_temp[0])>-(13.75+4.) ) ytof_temp[0] =-13.74;  
     if ( (ytof_temp[1])< -13.75 && fabs(ytof_temp[1])>-(13.75+4.) ) ytof_temp[1] =-13.74;  
   
     if ( (xtof_temp[2])> 4.5 && (xtof_temp[2])<4.5+4. ) xtof_temp[2] =4.4;  
     if ( (xtof_temp[3])> 4.45 && (xtof_temp[3])<4.5+4. ) xtof_temp[3] =4.4;  
     if ( (ytof_temp[2])> 3.75 && (ytof_temp[2])<3.75+4. ) ytof_temp[2] =3.74;  
     if ( (ytof_temp[3])> 3.75 && (ytof_temp[3])<3.75+4. ) ytof_temp[3] =3.74;  
     if ( (xtof_temp[2])< -4.5 && (xtof_temp[2])>-(4.5+4.) ) xtof_temp[2] =-4.4;  
     if ( (xtof_temp[3])< -4.45 && (xtof_temp[3])>-(4.5+4.) ) xtof_temp[3] =-4.4;  
     if ( (ytof_temp[2])< -3.75 && (ytof_temp[2])>-(3.75+4.) ) ytof_temp[2] =-3.74;  
     if ( (ytof_temp[3])< -3.75 && (ytof_temp[3])>-(3.75+4.) ) ytof_temp[3] =-3.74;  
   
     if ( (xtof_temp[4])> 6. && (xtof_temp[4])<6.+4. ) xtof_temp[2] =5.9;  
     if ( (xtof_temp[5])> 6. && (xtof_temp[5])<6.+4. ) xtof_temp[3] =5.9;  
     if ( (ytof_temp[4])> 5. && (ytof_temp[4])<5.+4. ) ytof_temp[2] =4.9;  
     if ( (ytof_temp[5])> 5. && (ytof_temp[5])<5.+4. ) ytof_temp[3] =4.9;  
     if ( (xtof_temp[4])< -6. && (xtof_temp[4])>-(6.+4.) ) xtof_temp[2] =-5.9;  
     if ( (xtof_temp[5])< -6. && (xtof_temp[5])>-(6.+4.) ) xtof_temp[3] =-5.9;  
     if ( (ytof_temp[4])< -5. && (ytof_temp[4])>-(5.+4.) ) ytof_temp[2] =-4.9;  
     if ( (ytof_temp[5])< -5. && (ytof_temp[5])>-(5.+4.) ) ytof_temp[3] =-4.9;  
   
   
     for (Int_t ii=0; ii<GetNPaddle(plane); ii++){  
       Int_t paddleid=ii;  
       pad = GetPaddleid(plane,paddleid);  
       Int_t IpaddleT=-1;  
       IpaddleT=this->GetPaddleIdOfTrack(xtof_temp[plane],ytof_temp[plane], plane,0.0);  
 //      printf("SSS %i %i %f %f  %f %f \n",IpaddleT,paddleid,xtof_temp[plane],ytof_temp[plane],t_tof->xtofpos[plane],t_tof->ytofpos[plane]);  
       //      if ( IpaddleT == paddleid || IpaddleT-1 == paddleid || IpaddleT+1 == paddleid ){  
       if ( IpaddleT == paddleid || GetTrueNHitPaddles(plane) == 1 ){  
         //IpaddleT-1 == paddleid || IpaddleT+1 == paddleid ){  
         GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);  
         dedx += PadEdx;  
       }  
     };  
   } else {  
     for (Int_t ii=0; ii<GetNPaddle(plane); ii++){  
       Int_t paddleid=ii;  
       pad = GetPaddleid(plane,paddleid);  
       GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);  
       dedx += PadEdx;  
 //      printf("TTT %i %i %f\n",paddleid,plane,PadEdx);  
     };  
   }  
   //  
   return(dedx);  
 };  
   
318    
319    //new, wm Feb 15
320  //wm Nov 08  //wm Nov 08
321  //gf Apr 07  //gf Apr 07
322  /**  /**
323   * Method to get the mean dEdx from a ToF layer - ATTENTION:   * Method to get the mean dEdx from a ToF layer
324   * It will sum up the dEdx of all the paddles, but since by definition   * By definition there should be PMTs with dEdx values only in one paddle of a layer
325   * only the paddle hitted by the track gets a dEdx value and the other   * (the paddle hitted by the track), this method looks for the hitted paddle
326   * paddles are set to zero, the output is just the dEdx of the hitted   * and gives the mean dEdx of that paddle as the output
327   * paddle in each layer!   * The method was modified for the "ToF-standalone" part in february 2015
328   * The "adcfl" option is not very useful (an artificial dEdx is per   * The "adcfl" option is not very useful (an artificial dEdx is per
329   * definition= 1 mip and not a real measurement), anyway left in the code   * definition= 1 mip and not a real measurement), anyway left in the code
330   * @param notrack Track Number   * @param trk Pointer to TofTrkVar object
331   * @param plane Plane index (0,1,2,3,4,5)   * @param plane Plane index (0,1,2,3,4,5)
332   * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )   * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
333   */   */
334  Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){  Float_t ToFLevel2::GetdEdx(ToFTrkVar *trk, Int_t plane, Int_t adcfl){
335    //  printf("fiffi\n");    
336    Float_t dedx = 0.;    Float_t dedx = 0.;
337    Float_t PadEdx =0.;    Float_t PadEdx =0.;
338    Int_t SatWarning;    Int_t SatWarning;
339    Int_t pad=-1;    Int_t pad=-1;
340    //    //
341    ToFTrkVar *trk = GetToFTrkVar(notrack);    if(!trk) cout << "ToFLevel2::GetdEdx(...) ---> NULL ToFTrkVar obj "<<endl;
342    if(!trk) return 0; //ELENA    if(!trk) return 0; //ELENA
343    //    //
344    if ( trk->trkseqno == -1 ){ //standalone, only paddles along the track, or about...    // ToF standalone part
345      //    printf("ciccio\n");    //
346      Float_t xleft=0;    if ( trk->trkseqno == -1 ){
     Float_t xright=0;  
     Float_t yleft=0;  
     Float_t yright=0;  
     Float_t xtof_temp[6]={100.,100.,100.,100.,100.,100.};  
     Float_t ytof_temp[6]={100.,100.,100.,100.,100.,100.};  
     ToFTrkVar *t_tof = trk;  
347            
348      xtof_temp[0]=t_tof->xtofpos[0];      //    ToFTrkVar *t_tof = trk;
349      ytof_temp[0]=t_tof->ytofpos[0];      
350      xtof_temp[1]=t_tof->xtofpos[0];      // Find the hitted paddle  (two good TDC values) using the tof_j_flag (from tofl2com.for)
351      ytof_temp[1]=t_tof->ytofpos[0];      
352        Int_t Ipaddle=-1;
353      xtof_temp[2]=t_tof->xtofpos[1];      // if tof_j_flag == 0: no paddle was hitted. Otherwise decode tof_j_flag to get the paddle
354      ytof_temp[2]=t_tof->ytofpos[1];      if (this->tof_j_flag[plane] > 0)  Ipaddle = (Int_t)log2(this->tof_j_flag[plane]) ;
     xtof_temp[3]=t_tof->xtofpos[1];  
     ytof_temp[3]=t_tof->ytofpos[1];  
   
     xtof_temp[4]=t_tof->xtofpos[2];  
     ytof_temp[4]=t_tof->ytofpos[2];  
     xtof_temp[5]=t_tof->xtofpos[2];  
     ytof_temp[5]=t_tof->ytofpos[2];  
   
     if(t_tof->xtofpos[0]<100. && t_tof->ytofpos[0]<100.){  
       xtof_temp[1]=t_tof->xtofpos[0];  
       ytof_temp[0]=t_tof->ytofpos[0];  
     }else if(t_tof->xtofpos[0]>=100. && t_tof->ytofpos[0]<100.){  
       ytof_temp[0]=t_tof->ytofpos[0];  
       this->GetPaddleGeometry(0,(Int_t)log2(this->tof_j_flag[0]),xleft, xright, yleft, yright);  
       xtof_temp[1]=xleft+2.55;  
     }else if(t_tof->ytofpos[0]>=100. && t_tof->xtofpos[0]<100.){  
       xtof_temp[1]=t_tof->xtofpos[0];  
       this->GetPaddleGeometry(1,(Int_t)log2(this->tof_j_flag[1]),xleft, xright, yleft, yright);  
       ytof_temp[0]=yleft+2.75;  
     }  
355            
356      if(t_tof->xtofpos[1]<100. && t_tof->ytofpos[1]<100.){      Ipaddle =  (Int_t)log2(this->tof_j_flag[plane]) ;
357        xtof_temp[2]=t_tof->xtofpos[1];      
358        ytof_temp[3]=t_tof->ytofpos[1];      // Get the dEdx of this paddle using "GetdEdxPaddle"
359      }else if(t_tof->xtofpos[1]>=100. && t_tof->ytofpos[1]<100.){      if (Ipaddle>-1) {
360        ytof_temp[3]=t_tof->ytofpos[1];        Int_t pad = GetPaddleid(plane,Ipaddle);
361        this->GetPaddleGeometry(3,(Int_t)log2(this->tof_j_flag[3]),xleft, xright, yleft, yright);        GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
362        xtof_temp[2]=xleft+4.5;        dedx = PadEdx;
     }else if(t_tof->ytofpos[1]>=100. && t_tof->xtofpos[1]<100.){  
       xtof_temp[2]=t_tof->xtofpos[1];  
       this->GetPaddleGeometry(2,(Int_t)log2(this->tof_j_flag[2]),xleft, xright, yleft, yright);  
       ytof_temp[3]=yleft+3.75;  
363      }      }
364            
365      if(t_tof->xtofpos[2]<100. && t_tof->ytofpos[2]<100.){      // If there was no correct hitted paddle, but there was one (and only one) paddle with some
366        xtof_temp[5]=t_tof->xtofpos[2];      // PMT entries in the PMT-class (found with "GetTrueNHitPaddles", use the dEdx of this paddle
367        ytof_temp[4]=t_tof->ytofpos[2];      
368      }else if(t_tof->xtofpos[2]>=100. && t_tof->ytofpos[2]<100.){      if ((Ipaddle<0) && (GetTrueNHitPaddles(plane)==1)) {
369        ytof_temp[4]=t_tof->ytofpos[2];        // find the paddle by looping over the paddles in each layer
370        this->GetPaddleGeometry(4,(Int_t)log2(this->tof_j_flag[4]),xleft, xright, yleft, yright);        // since GetTrueNHitPaddles==1 this is OK
371        xtof_temp[5]=xleft+3;        for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
372      }else if(t_tof->ytofpos[2]>=100. && t_tof->xtofpos[2]<100.){          Int_t paddleid=ii;
373        xtof_temp[5]=t_tof->xtofpos[2];          Int_t pad = GetPaddleid(plane,paddleid);
374        this->GetPaddleGeometry(5,(Int_t)log2(this->tof_j_flag[5]),xleft, xright, yleft, yright);          GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
375        ytof_temp[4]=yleft+2.5;          dedx += PadEdx;
     }  
   
     if ( (xtof_temp[0])> 17.85 && fabs(xtof_temp[0])<17.85+4. ) xtof_temp[0] =17.84;  
     if ( (xtof_temp[1])> 17.85 && fabs(xtof_temp[1])<17.85+4. ) xtof_temp[1] =17.84;  
     if ( (ytof_temp[0])> 13.75 && fabs(ytof_temp[0])<13.75+4. ) ytof_temp[0] =13.74;  
     if ( (ytof_temp[1])> 13.75 && fabs(ytof_temp[1])<13.75+4. ) ytof_temp[1] =13.74;  
     if ( (xtof_temp[0])< -17.85 && fabs(xtof_temp[0])>-(17.85+4.) ) xtof_temp[0] =-17.84;  
     if ( (xtof_temp[1])< -17.85 && fabs(xtof_temp[1])>-(17.85+4.) ) xtof_temp[1] =-17.84;  
     if ( (ytof_temp[0])< -13.75 && fabs(ytof_temp[0])>-(13.75+4.) ) ytof_temp[0] =-13.74;  
     if ( (ytof_temp[1])< -13.75 && fabs(ytof_temp[1])>-(13.75+4.) ) ytof_temp[1] =-13.74;  
   
     if ( (xtof_temp[2])> 4.5 && (xtof_temp[2])<4.5+4. ) xtof_temp[2] =4.4;  
     if ( (xtof_temp[3])> 4.45 && (xtof_temp[3])<4.5+4. ) xtof_temp[3] =4.4;  
     if ( (ytof_temp[2])> 3.75 && (ytof_temp[2])<3.75+4. ) ytof_temp[2] =3.74;  
     if ( (ytof_temp[3])> 3.75 && (ytof_temp[3])<3.75+4. ) ytof_temp[3] =3.74;  
     if ( (xtof_temp[2])< -4.5 && (xtof_temp[2])>-(4.5+4.) ) xtof_temp[2] =-4.4;  
     if ( (xtof_temp[3])< -4.45 && (xtof_temp[3])>-(4.5+4.) ) xtof_temp[3] =-4.4;  
     if ( (ytof_temp[2])< -3.75 && (ytof_temp[2])>-(3.75+4.) ) ytof_temp[2] =-3.74;  
     if ( (ytof_temp[3])< -3.75 && (ytof_temp[3])>-(3.75+4.) ) ytof_temp[3] =-3.74;  
   
     if ( (xtof_temp[4])> 6. && (xtof_temp[4])<6.+4. ) xtof_temp[2] =5.9;  
     if ( (xtof_temp[5])> 6. && (xtof_temp[5])<6.+4. ) xtof_temp[3] =5.9;  
     if ( (ytof_temp[4])> 5. && (ytof_temp[4])<5.+4. ) ytof_temp[2] =4.9;  
     if ( (ytof_temp[5])> 5. && (ytof_temp[5])<5.+4. ) ytof_temp[3] =4.9;  
     if ( (xtof_temp[4])< -6. && (xtof_temp[4])>-(6.+4.) ) xtof_temp[2] =-5.9;  
     if ( (xtof_temp[5])< -6. && (xtof_temp[5])>-(6.+4.) ) xtof_temp[3] =-5.9;  
     if ( (ytof_temp[4])< -5. && (ytof_temp[4])>-(5.+4.) ) ytof_temp[2] =-4.9;  
     if ( (ytof_temp[5])< -5. && (ytof_temp[5])>-(5.+4.) ) ytof_temp[3] =-4.9;  
   
   
     for (Int_t ii=0; ii<GetNPaddle(plane); ii++){  
       Int_t paddleid=ii;  
       pad = GetPaddleid(plane,paddleid);  
       Int_t IpaddleT=-1;  
       IpaddleT=this->GetPaddleIdOfTrack(xtof_temp[plane],ytof_temp[plane], plane,0.0);  
 //      printf("SSS %i %i %f %f  %f %f \n",IpaddleT,paddleid,xtof_temp[plane],ytof_temp[plane],t_tof->xtofpos[plane],t_tof->ytofpos[plane]);  
       //      if ( IpaddleT == paddleid || IpaddleT-1 == paddleid || IpaddleT+1 == paddleid ){  
       if ( IpaddleT == paddleid || GetTrueNHitPaddles(plane) == 1 ){  
         //IpaddleT-1 == paddleid || IpaddleT+1 == paddleid ){  
         GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);  
         dedx += PadEdx;  
376        }        }
377      };      }
378    } else {    } else {
379        // track dependent dEdx: simple, there will be only one paddle hitted in    each layer
380        // so just loop over the paddles in each layer
381      for (Int_t ii=0; ii<GetNPaddle(plane); ii++){      for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
382        Int_t paddleid=ii;        Int_t paddleid=ii;
383        pad = GetPaddleid(plane,paddleid);        pad = GetPaddleid(plane,paddleid);
384        GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);        GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
385        dedx += PadEdx;        dedx += PadEdx;
386  //      printf("TTT %i %i %f\n",paddleid,plane,PadEdx);      }
     };  
387    }    }
388    //    //
389    return(dedx);    return(dedx);
390  };  }
391    
392  /**  /**
393   * Method to fill the ADC_C 4x12 matrix with the dEdx values and the TDC 4x12 matrix   * Method to fill the ADC_C 4x12 matrix with the dEdx values and the TDC 4x12 matrix
394   * with the time-walk corrected TDC values.   * with the time-walk corrected TDC values.
395   * @param notrack Track Number   * @param notrack Track Number (arry index, ranging from 0 to ntrk())
396   * @param adc  ADC_C matrix with dEdx values   * @param adc  ADC_C matrix with dEdx values
397   * @param tdc  TDC matrix   * @param tdc  TDC matrix
398   */   */
# Line 640  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt Line 462  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt
462   */   */
463  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){
464    //    //
465    short tof[4][24] = {    static const short tof[4][24] = {
466      {4, 4,  4,  4,  1,  1, 2, 2,  3,  3, 3, 3,  3,  3, 1, 1,  1,  1, 2, 3,  3, 3, 3,  4},      {4, 4,  4,  4,  1,  1, 2, 2,  3,  3, 3, 3,  3,  3, 1, 1,  1,  1, 2, 3,  3, 3, 3,  4},
467      {1, 3,  5,  7, 10, 12, 2, 4,  2,  4, 6, 8, 10, 12, 1, 5,  3,  9, 7, 9, 11, 1, 5,  9},      {1, 3,  5,  7, 10, 12, 2, 4,  2,  4, 6, 8, 10, 12, 1, 5,  3,  9, 7, 9, 11, 1, 5,  9},
468      {2, 2,  2,  2,  1,  1, 1, 1,  4,  4, 4, 4,  4,  4, 2, 1,  2,  1, 2, 2,  2, 3, 3,  4},      {2, 2,  2,  2,  1,  1, 1, 1,  4,  4, 4, 4,  4,  4, 2, 1,  2,  1, 2, 2,  2, 3, 3,  4},
# Line 676  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_ Line 498  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_
498   */   */
499  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){
500    //    //
501    short tof[4][24] = {    static const short tof[4][24] = {
502      {4, 4,  4,  4,  1,  1, 2, 2,  3,  3, 3, 3,  3,  3, 1, 1,  1,  1, 2, 3,  3, 3, 3,  4},      {4, 4,  4,  4,  1,  1, 2, 2,  3,  3, 3, 3,  3,  3, 1, 1,  1,  1, 2, 3,  3, 3, 3,  4},
503      {1, 3,  5,  7, 10, 12, 2, 4,  2,  4, 6, 8, 10, 12, 1, 5,  3,  9, 7, 9, 11, 1, 5,  9},      {1, 3,  5,  7, 10, 12, 2, 4,  2,  4, 6, 8, 10, 12, 1, 5,  3,  9, 7, 9, 11, 1, 5,  9},
504      {2, 2,  2,  2,  1,  1, 1, 1,  4,  4, 4, 4,  4,  4, 2, 1,  2,  1, 2, 2,  2, 3, 3,  4},      {2, 2,  2,  2,  1,  1, 1, 1,  4,  4, 4, 4,  4,  4, 2, 1,  2,  1, 2, 2,  2, 3, 3,  4},
# Line 708  void ToFLevel2::GetPMTIndex(Int_t ind, I Line 530  void ToFLevel2::GetPMTIndex(Int_t ind, I
530   * Method to get the dEdx from a given ToF paddle.   * Method to get the dEdx from a given ToF paddle.
531   * If two PMTs are good, the mean dEdx of both PMTs is taken, otherwise   * If two PMTs are good, the mean dEdx of both PMTs is taken, otherwise
532   * just the dEdx of the "good" PMT. If both PMTs are above saturation => dEdx=1000   * just the dEdx of the "good" PMT. If both PMTs are above saturation => dEdx=1000
533   * @param notrack Track Number   * @param notrack Track Number (-1 for standalone info)
534   * @param Paddle index (0,1,...,23).   * @param Paddle index (0,1,...,23).
535   * @param adcflag in the paddle (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )   * @param adcflag in the paddle (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
536   * @param PadEdx dEdx from a given ToF paddle   * @param PadEdx dEdx from a given ToF paddle
# Line 716  void ToFLevel2::GetPMTIndex(Int_t ind, I Line 538  void ToFLevel2::GetPMTIndex(Int_t ind, I
538   */   */
539  void ToFLevel2::GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){  void ToFLevel2::GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){
540    
   /*  
     Float_t  PMTsat[48] = {  
     3162.14, 3165.48, 3153.85, 3085.73, 3089.65, 3107.64, 3097.52, 3078.37,  
     3130.05, 3087.07, 3112.22, 3102.92, 3080.58, 3092.55, 3087.94, 3125.03,  
     3094.09, 3143.16, 3125.51, 3181.27, 3092.09, 3124.98, 3069.3, 3095.53,  
     3097.11, 3133.53, 3114.73, 3113.01, 3091.19, 3097.99, 3033.84, 3134.98,  
     3081.37, 3111.04, 3066.77, 3108.17, 3133, 3111.06, 3052.52, 3140.66,  
     3106.33, 3094.85, 3150.85, 3118.8, 3096.24, 3118.47,3111.36, 3117.11 } ;  
   */  
   
   // new values from Napoli dec 2008  
   Float_t  PMTsat[48] = {  
     3176.35,3178.19,3167.38,3099.73,3117.00,3126.29,3111.44,3092.27,  
     3146.48,3094.41,3132.13,3115.37,3099.32,3110.97,3111.80,3143.14,  
     3106.72,3153.44,3136.00,3188.96,3104.73,3140.45,3073.18,3106.62,  
     3112.48,3146.92,3127.24,3136.52,3109.59,3112.89,3045.15,3147.26,  
     3095.92,3121.05,3083.25,3123.62,3150.92,3125.30,3067.60,3160.18,  
     3119.36,3108.92,3164.77,3133.64,3111.47,3131.98,3128.87,3135.56 };  
   
   for (Int_t i=0; i<48;i++) PMTsat[i] = PMTsat[i] - 5.;  // safety margin  
   
   
   PadEdx = 0.;  
   //  SatWarning = 1000;  
   SatWarning = 0;   // 0=good, increase for each bad PMT  
   
   Float_t dEdx[48] = {0};  
   Int_t pmt_id = -1;  
   Float_t adcraw[48];  
   //  
541    ToFTrkVar *trk = GetToFTrkVar(notrack);    ToFTrkVar *trk = GetToFTrkVar(notrack);
542    if(!trk) return; //ELENA  //    ToFTrkVar *trk = GetToFStoredTrack(seqno); //Elena 2015
543    //    this->GetdEdxPaddle(trk, paddleid, adcfl, PadEdx, SatWarning);
   
   Int_t pmtleft=-1;  
   Int_t pmtright=-1;  
   GetPaddlePMT(paddleid, pmtleft, pmtright);  
   
   adcraw[pmtleft] = 4095;  
   adcraw[pmtright] = 4095;  
   
     
   for (Int_t jj=0; jj<npmt(); jj++){  
       
     ToFPMT *pmt = GetToFPMT(jj);  
     if(!pmt)break; //ELENA  
       
     pmt_id = pmt->pmt_id;  
     if(pmt_id==pmtleft){  
       adcraw[pmtleft] = pmt->adc;  
     }  
       
     if(pmt_id==pmtright){  
       adcraw[pmtright] = pmt->adc;  
     }  
   }  
   
     
   for (Int_t i=0; i<trk->npmtadc; i++){  
   
     if((trk->adcflag).At(i)==0 || adcfl==100){  
       if((trk->pmtadc).At(i) == pmtleft)dEdx[pmtleft] = (trk->dedx).At(i);  
       if((trk->pmtadc).At(i) == pmtright)dEdx[pmtright] = (trk->dedx).At(i);  
     }else{  
       if((trk->pmtadc).At(i) == pmtleft)dEdx[pmtleft] = 0.;  
       if((trk->pmtadc).At(i) == pmtright)dEdx[pmtright] = 0.;  
     }  
   }  
   
   
   //  if( adcraw[pmtleft] >3000 || adcraw[pmtright] >3000)SatWarning=1;  //old version  
   
   // Increase SatWarning Counter for each PMT>Sat  
   if( adcraw[pmtleft] > PMTsat[pmtleft])SatWarning++;    
   if( adcraw[pmtright] > PMTsat[pmtright])SatWarning++;  
   
   // if ADC  > sat set dEdx=1000  
   if( adcraw[pmtleft] > PMTsat[pmtleft]) dEdx[pmtleft] = 1000.;  
   if( adcraw[pmtright] > PMTsat[pmtright]) dEdx[pmtright] = 1000. ;  
   
   // if two PMT are good, take mean dEdx, otherwise only the good dEdx  
   if(dEdx[pmtleft]<1000 && dEdx[pmtright]<1000) PadEdx = (dEdx[pmtleft]+dEdx[pmtright])*0.5;  
   if(dEdx[pmtleft]==1000 && dEdx[pmtright]<1000) PadEdx = dEdx[pmtright];    
   if(dEdx[pmtleft]<1000 && dEdx[pmtright]==1000) PadEdx = dEdx[pmtleft];  
544        
545  };  };
546    
# Line 848  void ToFLevel2::GetdEdxPaddle(ToFTrkVar Line 589  void ToFLevel2::GetdEdxPaddle(ToFTrkVar
589    Int_t pmt_id = -1;    Int_t pmt_id = -1;
590    Float_t adcraw[48];    Float_t adcraw[48];
591    //    //
592      if(!trk)cout << "ToFLevel2::GetdEdxPaddle(...) ---> NULL ToFTrkVar obj "<<endl;
593    if(!trk) return; //ELENA    if(!trk) return; //ELENA
594    //    //
595    
# Line 918  TString ToFLevel2::GetPMTName(Int_t ind, Line 660  TString ToFLevel2::GetPMTName(Int_t ind,
660        
661    TString pmtname = " ";    TString pmtname = " ";
662        
663    TString photoS[48] = {    static const TString photoS[48] = {
664      "S11_1A", "S11_1B", "S11_2A", "S11_2B", "S11_3A", "S11_3B", "S11_4A",      "S11_1A", "S11_1B", "S11_2A", "S11_2B", "S11_3A", "S11_3B", "S11_4A",
665      "S11_4B",      "S11_4B",
666      "S11_5A", "S11_5B", "S11_6A", "S11_6B", "S11_7A", "S11_7B", "S11_8A",      "S11_5A", "S11_5B", "S11_6A", "S11_6B", "S11_7A", "S11_7B", "S11_8A",
# Line 1378  Int_t ToFLevel2::GetNPaddle(Int_t plane) Line 1120  Int_t ToFLevel2::GetNPaddle(Int_t plane)
1120   * @param cut on chi2   * @param cut on chi2
1121   */   */
1122    
 Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){  
1123    
1124  //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;  Float_t ToFTrkVar::CalcBeta( Float_t resmax, Float_t qualitycut, Float_t chi2cut){
1125    
1126    
1127    Float_t bxx = 100.;    Float_t bxx = 100.;
1128    //    //
1129    ToFTrkVar *trk = GetToFTrkVar(notrack);    ToFTrkVar *trk = this;
   if(!trk) return 0; //ELENA  
1130    
1131    
1132    Float_t chi2,xhelp,beta_mean;    Float_t chi2,xhelp,beta_mean;
# Line 1412  Float_t  w_il[6]; Line 1153  Float_t  w_il[6];
1153    for (Int_t i=0; i<trk->npmttdc; i++){    for (Int_t i=0; i<trk->npmttdc; i++){
1154      //      //
1155      pmt_id = (trk->pmttdc).At(i);      pmt_id = (trk->pmttdc).At(i);
1156      pmt_plane = GetPlaneIndex(pmt_id);      pmt_plane = ToFLevel2::GetPlaneIndex(pmt_id);
1157      tdcfl = (trk->tdcflag).At(i);      tdcfl = (trk->tdcflag).At(i);
1158      if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag      if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag
1159                                       };                                       };
# Line 1499  Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2, Line 1240  Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2,
1240    //    //
1241    return(bxx);    return(bxx);
1242  };  };
1243    ////////////////////////////////////////////////////
1244    ////////////////////////////////////////////////////
1245    /**
1246     * See ToFTrkVar::CalcBeta(Float_t,Float_t, Float_t).
1247     * @param notrack Track Number (arry index, ranging from 0 to ntrk())
1248     */
1249    Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){
1250    
1251    //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;
1252    
1253      ToFTrkVar *trk = GetToFTrkVar(notrack);
1254    //    ToFTrkVar *trk = GetToFStoredTrack(seqno);//Elena Apr 2015
1255      if(!trk) return 0; //ELENA
1256    
1257      return trk->CalcBeta(resmax,qualitycut,chi2cut);
1258    
1259    };
1260    
1261    
1262  ////////////////////////////////////////////////////  ////////////////////////////////////////////////////
Legend:
Removed from v.1.42  
changed lines
  Added in v.1.48
  ViewVC Help
Powered by ViewVC 1.1.23