/[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.36 by mocchiut, Thu Apr 12 12:27:10 2012 UTC revision 1.46 by pam-fi, Tue May 19 10:49:14 2015 UTC
# Line 113  void ToFLevel2::Set(){//ELENA Line 113  void ToFLevel2::Set(){//ELENA
113      if(!PMT)PMT = new TClonesArray("ToFPMT",12); //ELENA      if(!PMT)PMT = new TClonesArray("ToFPMT",12); //ELENA
114      if(!ToFTrk)ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA      if(!ToFTrk)ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA
115  }//ELENA  }//ELENA
116    //--------------------------------------
117    //
118    //
119    //--------------------------------------
120    void ToFLevel2::SetTrackArray(TClonesArray *track){//ELENA
121        if(track && strcmp(track->GetClass()->GetName(),"ToFTrkVar")==0){
122            if(ToFTrk)ToFTrk->Clear("C");
123            ToFTrk = track;
124        }
125    }
126    
127  void ToFLevel2::Clear(Option_t *t){  void ToFLevel2::Clear(Option_t *t){
128    //    //
# Line 151  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t Line 161  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t
161    return toftrack;    return toftrack;
162  }  }
163    
164    /**
165     * Retrieves the tof track matching the seqno-th tracker stored track.
166     *
167     */
168    ToFTrkVar *ToFLevel2::GetToFStoredTrack(int seqno){
169    
170      if( ntrk()==0 ){
171        printf("ToFLevel2::GetToFStoredTrack(int) : requested tracker SeqNo %i but no ToFrimeter tracks are stored\n",seqno);
172        return NULL;
173      };
174      
175      ToFTrkVar *c = 0;
176      Int_t it_tof=0;
177        
178      do {
179        c = GetToFTrkVar(it_tof);
180        it_tof++;
181      } while( c && seqno != c->trkseqno && it_tof < ntrk());      
182      
183      if(!c || seqno != c->trkseqno){
184        c = 0;
185        if(seqno!=-1 ) printf("ToFLevel2::GetToFStoredTrack(int) : requested tracker SeqNo %i does not match ToFrimeter stored tracks\n",seqno);
186      };
187      return c;
188        
189    }
190    
191    
192  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit){  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit){
193    //        //    
194    if(ihit >= npmt()){    if(ihit >= npmt()){
# Line 214  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit) Line 252  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
252         false) return true;         false) return true;
253      else return false;      else return false;
254  };  };
255    
256  /**  /**
257   * 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.
258     * The method uses "HitPaddle" which checks if there is a TDC signal
259     * from both PMTs.
260   * @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).
261   */   */
262  Int_t ToFLevel2::GetNHitPaddles(Int_t plane){  Int_t ToFLevel2::GetNHitPaddles(Int_t plane){
# Line 225  Int_t ToFLevel2::GetNHitPaddles(Int_t pl Line 266  Int_t ToFLevel2::GetNHitPaddles(Int_t pl
266  };  };
267    
268  /**  /**
269   * 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.
270     * The method does NOT check if there is a signal from both PMTs, it only
271     * checks if there is some PMT signal in a paddle
272   * @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).
273   */   */
274  Int_t ToFLevel2::GetTrueNHitPaddles(Int_t plane){  Int_t ToFLevel2::GetTrueNHitPaddles(Int_t plane){
# Line 245  Int_t ToFLevel2::GetTrueNHitPaddles(Int_ Line 288  Int_t ToFLevel2::GetTrueNHitPaddles(Int_
288      return npad;      return npad;
289  };  };
290    
291    //new, wm Feb 15
292  //wm Nov 08  //wm Nov 08
293  //gf Apr 07  //gf Apr 07
294  /**  /**
295   * Method to get the mean dEdx from a ToF layer - ATTENTION:   * Method to get the mean dEdx from a ToF layer
296   * 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
297   * 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
298   * 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
299   * paddle in each layer!   * The method was modified for the "ToF-standalone" part in february 2015
300   * The "adcfl" option is not very useful (an artificial dEdx is per   * The "adcfl" option is not very useful (an artificial dEdx is per
301   * 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
302   * @param notrack Track Number   * @param notrack Track Number
# Line 260  Int_t ToFLevel2::GetTrueNHitPaddles(Int_ Line 304  Int_t ToFLevel2::GetTrueNHitPaddles(Int_
304   * @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; )
305   */   */
306  Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){  Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){
307    //  printf("fiffi\n");  //  ToFTrkVar *trk = GetToFTrkVar(notrack);
308        ToFTrkVar *trk = GetToFStoredTrack(notrack);//Elena 2015
309      return this->GetdEdx(trk, plane, adcfl);
310    }
311    
312    //new, wm Feb 15
313    //wm Nov 08
314    //gf Apr 07
315    /**
316     * Method to get the mean dEdx from a ToF layer
317     * By definition there should be PMTs with dEdx values only in one paddle of a layer
318     * (the paddle hitted by the track), this method looks for the hitted paddle
319     * and gives the mean dEdx of that paddle as the output
320     * The method was modified for the "ToF-standalone" part in february 2015
321     * The "adcfl" option is not very useful (an artificial dEdx is per
322     * definition= 1 mip and not a real measurement), anyway left in the code
323     * @param trk Pointer to TofTrkVar object
324     * @param plane Plane index (0,1,2,3,4,5)
325     * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
326     */
327    Float_t ToFLevel2::GetdEdx(ToFTrkVar *trk, Int_t plane, Int_t adcfl){
328      
329    Float_t dedx = 0.;    Float_t dedx = 0.;
330    Float_t PadEdx =0.;    Float_t PadEdx =0.;
331    Int_t SatWarning;    Int_t SatWarning;
332    Int_t pad=-1;    Int_t pad=-1;
333    //    //
334    ToFTrkVar *trk = GetToFTrkVar(notrack);    if(!trk) cout << "ToFLevel2::GetdEdx(...) ---> NULL ToFTrkVar obj "<<endl;
335    if(!trk) return 0; //ELENA    if(!trk) return 0; //ELENA
336    //    //
337    if ( trk->trkseqno == -1 ){ //standalone, only paddles along the track, or about...    // ToF standalone part
338      //    printf("ciccio\n");    //
339      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;  
340            
341      xtof_temp[0]=t_tof->xtofpos[0];      //    ToFTrkVar *t_tof = trk;
342      ytof_temp[0]=t_tof->ytofpos[0];      
343      xtof_temp[1]=t_tof->xtofpos[0];      // Find the hitted paddle  (two good TDC values) using the tof_j_flag (from tofl2com.for)
344      ytof_temp[1]=t_tof->ytofpos[0];      
345        Int_t Ipaddle=-1;
346      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
347      ytof_temp[2]=t_tof->ytofpos[1];      if (this->tof_j_flag[plane] > 0)  Ipaddle = (Int_t)log2(this->tof_j_flag[plane]) ;
348      xtof_temp[3]=t_tof->xtofpos[1];      
349      ytof_temp[3]=t_tof->ytofpos[1];      Ipaddle =  (Int_t)log2(this->tof_j_flag[plane]) ;
   
     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;  
     }  
350            
351      if(t_tof->xtofpos[1]<100. && t_tof->ytofpos[1]<100.){      // Get the dEdx of this paddle using "GetdEdxPaddle"
352        xtof_temp[2]=t_tof->xtofpos[1];      if (Ipaddle>-1) {
353        ytof_temp[3]=t_tof->ytofpos[1];        Int_t pad = GetPaddleid(plane,Ipaddle);
354      }else if(t_tof->xtofpos[1]>=100. && t_tof->ytofpos[1]<100.){        GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
355        ytof_temp[3]=t_tof->ytofpos[1];        dedx = PadEdx;
       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;  
356      }      }
357            
358      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
359        xtof_temp[5]=t_tof->xtofpos[2];      // PMT entries in the PMT-class (found with "GetTrueNHitPaddles", use the dEdx of this paddle
360        ytof_temp[4]=t_tof->ytofpos[2];      
361      }else if(t_tof->xtofpos[2]>=100. && t_tof->ytofpos[2]<100.){      if ((Ipaddle<0) && (GetTrueNHitPaddles(plane)==1)) {
362        ytof_temp[4]=t_tof->ytofpos[2];        // find the paddle by looping over the paddles in each layer
363        this->GetPaddleGeometry(4,(Int_t)log2(this->tof_j_flag[4]),xleft, xright, yleft, yright);        // since GetTrueNHitPaddles==1 this is OK
364        xtof_temp[5]=xleft+3;        for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
365      }else if(t_tof->ytofpos[2]>=100. && t_tof->xtofpos[2]<100.){          Int_t paddleid=ii;
366        xtof_temp[5]=t_tof->xtofpos[2];          Int_t pad = GetPaddleid(plane,paddleid);
367        this->GetPaddleGeometry(5,(Int_t)log2(this->tof_j_flag[5]),xleft, xright, yleft, yright);          GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
368        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;  
369        }        }
370      };      }
371    } else {    } else {
372        // track dependent dEdx: simple, there will be only one paddle hitted in    each layer
373        // so just loop over the paddles in each layer
374      for (Int_t ii=0; ii<GetNPaddle(plane); ii++){      for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
375        Int_t paddleid=ii;        Int_t paddleid=ii;
376        pad = GetPaddleid(plane,paddleid);        pad = GetPaddleid(plane,paddleid);
377        GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);        GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
378        dedx += PadEdx;        dedx += PadEdx;
379  //      printf("TTT %i %i %f\n",paddleid,plane,PadEdx);      }
     };  
380    }    }
381    //    //
382    return(dedx);    return(dedx);
383  };  }
384    
385  /**  /**
386   * 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
# Line 460  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt Line 455  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt
455   */   */
456  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){
457    //    //
458    short tof[4][24] = {    static const short tof[4][24] = {
459      {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},
460      {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},
461      {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 496  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_ Line 491  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_
491   */   */
492  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){
493    //    //
494    short tof[4][24] = {    static const short tof[4][24] = {
495      {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},
496      {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},
497      {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 536  void ToFLevel2::GetPMTIndex(Int_t ind, I Line 531  void ToFLevel2::GetPMTIndex(Int_t ind, I
531   */   */
532  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){
533    
534    //  ToFTrkVar *trk = GetToFTrkVar(notrack);
535        ToFTrkVar *trk = GetToFStoredTrack(notrack); //Elena 2015
536      this->GetdEdxPaddle(trk, paddleid, adcfl, PadEdx, SatWarning);
537      
538    };
539    
540    //
541    //  wm Nov 08 revision - saturation values included
542    /// gf Apr 07
543    /**
544     * Method to get the dEdx from a given ToF paddle.
545     * If two PMTs are good, the mean dEdx of both PMTs is taken, otherwise
546     * just the dEdx of the "good" PMT. If both PMTs are above saturation => dEdx=1000
547     * @param notrack Track Number
548     * @param Paddle index (0,1,...,23).
549     * @param adcflag in the paddle (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
550     * @param PadEdx dEdx from a given ToF paddle
551     * @param SatWarning 1 if the PMT ios near saturation region (adcraw ~3000)
552     */
553    void ToFLevel2::GetdEdxPaddle(ToFTrkVar *trk, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){
554    
555    /*    /*
556      Float_t  PMTsat[48] = {      Float_t  PMTsat[48] = {
557      3162.14, 3165.48, 3153.85, 3085.73, 3089.65, 3107.64, 3097.52, 3078.37,      3162.14, 3165.48, 3153.85, 3085.73, 3089.65, 3107.64, 3097.52, 3078.37,
# Line 566  void ToFLevel2::GetdEdxPaddle(Int_t notr Line 582  void ToFLevel2::GetdEdxPaddle(Int_t notr
582    Int_t pmt_id = -1;    Int_t pmt_id = -1;
583    Float_t adcraw[48];    Float_t adcraw[48];
584    //    //
585    ToFTrkVar *trk = GetToFTrkVar(notrack);    if(!trk)cout << "ToFLevel2::GetdEdxPaddle(...) ---> NULL ToFTrkVar obj "<<endl;
586    if(!trk) return; //ELENA    if(!trk) return; //ELENA
587    //    //
588    
# Line 622  void ToFLevel2::GetdEdxPaddle(Int_t notr Line 638  void ToFLevel2::GetdEdxPaddle(Int_t notr
638    if(dEdx[pmtleft]<1000 && dEdx[pmtright]==1000) PadEdx = dEdx[pmtleft];    if(dEdx[pmtleft]<1000 && dEdx[pmtright]==1000) PadEdx = dEdx[pmtleft];
639        
640  };  };
 //  
   
641    
642  // gf Apr 07  // gf Apr 07
643    
# Line 639  TString ToFLevel2::GetPMTName(Int_t ind, Line 653  TString ToFLevel2::GetPMTName(Int_t ind,
653        
654    TString pmtname = " ";    TString pmtname = " ";
655        
656    TString photoS[48] = {    static const TString photoS[48] = {
657      "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",
658      "S11_4B",      "S11_4B",
659      "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 1099  Int_t ToFLevel2::GetNPaddle(Int_t plane) Line 1113  Int_t ToFLevel2::GetNPaddle(Int_t plane)
1113   * @param cut on chi2   * @param cut on chi2
1114   */   */
1115    
 Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){  
1116    
1117  //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;  Float_t ToFTrkVar::CalcBeta( Float_t resmax, Float_t qualitycut, Float_t chi2cut){
1118    
1119    
1120    Float_t bxx = 100.;    Float_t bxx = 100.;
1121    //    //
1122    ToFTrkVar *trk = GetToFTrkVar(notrack);    ToFTrkVar *trk = this;
   if(!trk) return 0; //ELENA  
1123    
1124    
1125    Float_t chi2,xhelp,beta_mean;    Float_t chi2,xhelp,beta_mean;
# Line 1133  Float_t  w_il[6]; Line 1146  Float_t  w_il[6];
1146    for (Int_t i=0; i<trk->npmttdc; i++){    for (Int_t i=0; i<trk->npmttdc; i++){
1147      //      //
1148      pmt_id = (trk->pmttdc).At(i);      pmt_id = (trk->pmttdc).At(i);
1149      pmt_plane = GetPlaneIndex(pmt_id);      pmt_plane = ToFLevel2::GetPlaneIndex(pmt_id);
1150      tdcfl = (trk->tdcflag).At(i);      tdcfl = (trk->tdcflag).At(i);
1151      if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag      if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag
1152                                       };                                       };
# Line 1220  Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2, Line 1233  Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2,
1233    //    //
1234    return(bxx);    return(bxx);
1235  };  };
1236    ////////////////////////////////////////////////////
1237    ////////////////////////////////////////////////////
1238    /**
1239     * See ToFTrkVar::CalcBeta(Float_t,Float_t, Float_t).
1240     */
1241    Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){
1242    
1243    //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;
1244    
1245      ToFTrkVar *trk = GetToFTrkVar(notrack);
1246      if(!trk) return 0; //ELENA
1247    
1248      return trk->CalcBeta(resmax,qualitycut,chi2cut);
1249    
1250    };
1251    
1252    
1253  ////////////////////////////////////////////////////  ////////////////////////////////////////////////////
# Line 1309  Int_t ToFLevel2::Process(TrkLevel2 *trk, Line 1337  Int_t ToFLevel2::Process(TrkLevel2 *trk,
1337  //   if ( !dbc->IsConnected() ) return 1;  //   if ( !dbc->IsConnected() ) return 1;
1338  //   stringstream myquery;  //   stringstream myquery;
1339  //   myquery.str("");  //   myquery.str("");
1340  //   myquery << "SET time_zone='+0:00'";  //   myquery << "SET time_zone='+0:00';";
1341  //   dbc->Query(myquery.str().c_str());  //   dbc->Query(myquery.str().c_str());
1342    //   delete dbc->Query("SET sql_mode = 'NO_UNSIGNED_SUBTRACTION';");
1343  //   GL_PARAM *glparam = new GL_PARAM();  //   GL_PARAM *glparam = new GL_PARAM();
1344  //   glparam->Query_GL_PARAM(1,1,dbc); // parameters stored in DB in GL_PRAM table  //   glparam->Query_GL_PARAM(1,1,dbc); // parameters stored in DB in GL_PRAM table
1345  //   trk->LoadField(glparam->PATH+glparam->NAME);  //   trk->LoadField(glparam->PATH+glparam->NAME);
# Line 1784  void ToFdEdx::Init(pamela::tof::TofEvent Line 1813  void ToFdEdx::Init(pamela::tof::TofEvent
1813      for (Int_t hh=0; hh<12;hh++){      for (Int_t hh=0; hh<12;hh++){
1814        //          tofinput_.tdc[hh][gg]=tofEvent->tdc[gg][hh];                  //          tofinput_.tdc[hh][gg]=tofEvent->tdc[gg][hh];          
1815        int mm = tf.GetPMTid(gg,hh);                int mm = tf.GetPMTid(gg,hh);        
1816        adc[mm]=tofl0->adc[gg][hh];        adc[mm]= (0xFFF & tofl0->adc[gg][hh]); // EM, exclude warning bits
1817      };            };      
1818    };    };
1819        
# Line 1899  void ToFdEdx::Process(UInt_t atime, Floa Line 1928  void ToFdEdx::Process(UInt_t atime, Floa
1928      //    printf(" e qui? \n");      //    printf(" e qui? \n");
1929      //---------------------------------------------------- Z reconstruction      //---------------------------------------------------- Z reconstruction
1930    
1931      double adcHe, adcnorm, adclin, dEdx, Zeta;      double adcHe, adcnorm, adclin, dEdx;//, Zeta; // EM GCC4.7
1932    
1933      adcHe=-2;      adcHe=-2;
1934      adcnorm=-2;      adcnorm=-2;
1935      adclin=-2;      adclin=-2;
1936      dEdx=-2;      dEdx=-2;
1937      Zeta=-2;      //    Zeta=-2;//EM GCC4.7
1938      Double_t correction = 1.;      Double_t correction = 1.;
1939    
1940      if(Aconn==1 && (ii==0 || ii==20 || ii==22 || ii==24)){      if(Aconn==1 && (ii==0 || ii==20 || ii==22 || ii==24)){
Legend:
Removed from v.1.36  
changed lines
  Added in v.1.46
  ViewVC Help
Powered by ViewVC 1.1.23