/[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.13 by mocchiut, Mon Jan 22 14:09:32 2007 UTC revision 1.19 by mocchiut, Mon Mar 3 09:51:02 2008 UTC
# Line 1  Line 1 
1    /**
2     * \file ToFLevel2.cpp
3     * \author Gianfranca DeRosa, Wolfgang Menn
4     */
5    
6  #include <TObject.h>  #include <TObject.h>
7  #include <ToFLevel2.h>  #include <ToFLevel2.h>
8  #include <iostream>  #include <iostream>
# Line 10  ToFPMT::ToFPMT(){ Line 15  ToFPMT::ToFPMT(){
15    pmt_id = 0;    pmt_id = 0;
16    adc = 0.;    adc = 0.;
17    tdc_tw = 0.;    tdc_tw = 0.;
18      tdc = 0.;
19  }  }
20    
21  ToFPMT::ToFPMT(const ToFPMT &t){  ToFPMT::ToFPMT(const ToFPMT &t){
22    pmt_id = t.pmt_id;    pmt_id = t.pmt_id;
23    adc = t.adc;    adc = t.adc;
24    tdc_tw = t.tdc_tw;    tdc_tw = t.tdc_tw;
25      tdc = t.tdc;
26  }  }
27    
28  void ToFPMT::Clear(){  void ToFPMT::Clear(Option_t *t){
29    pmt_id = 0;    pmt_id = 0;
30    adc = 0.;    adc = 0.;
31    tdc_tw = 0.;    tdc_tw = 0.;
32      tdc = 0.;
33  }  }
34    
35    
# Line 40  ToFTrkVar::ToFTrkVar() { Line 48  ToFTrkVar::ToFTrkVar() {
48    memset(beta,  0, 13*sizeof(Float_t));    memset(beta,  0, 13*sizeof(Float_t));
49    memset(xtofpos,  0, 3*sizeof(Float_t));    memset(xtofpos,  0, 3*sizeof(Float_t));
50    memset(ytofpos,  0, 3*sizeof(Float_t));    memset(ytofpos,  0, 3*sizeof(Float_t));
51      memset(xtr_tof,  0, 6*sizeof(Float_t));
52      memset(ytr_tof,  0, 6*sizeof(Float_t));
53    //    //
54  };  };
55    
56  void ToFTrkVar::Clear() {  void ToFTrkVar::Clear(Option_t *t) {
57    trkseqno = 0;    trkseqno = 0;
58    npmttdc = 0;    npmttdc = 0;
59    npmtadc = 0;    npmtadc = 0;
# Line 56  void ToFTrkVar::Clear() { Line 66  void ToFTrkVar::Clear() {
66    memset(beta,  0, 13*sizeof(Float_t));    memset(beta,  0, 13*sizeof(Float_t));
67    memset(xtofpos,  0, 3*sizeof(Float_t));    memset(xtofpos,  0, 3*sizeof(Float_t));
68    memset(ytofpos,  0, 3*sizeof(Float_t));    memset(ytofpos,  0, 3*sizeof(Float_t));
69      memset(xtr_tof,  0, 6*sizeof(Float_t));
70      memset(ytr_tof,  0, 6*sizeof(Float_t));
71    //    //
72  };  };
73    
# Line 74  ToFTrkVar::ToFTrkVar(const ToFTrkVar &t) Line 86  ToFTrkVar::ToFTrkVar(const ToFTrkVar &t)
86    memcpy(beta,t.beta,sizeof(beta));    memcpy(beta,t.beta,sizeof(beta));
87    memcpy(xtofpos,t.xtofpos,sizeof(xtofpos));    memcpy(xtofpos,t.xtofpos,sizeof(xtofpos));
88    memcpy(ytofpos,t.ytofpos,sizeof(ytofpos));    memcpy(ytofpos,t.ytofpos,sizeof(ytofpos));
89      memcpy(xtr_tof,t.xtr_tof,sizeof(xtr_tof));
90      memcpy(ytr_tof,t.ytr_tof,sizeof(ytr_tof));
91    //    //
92  };  };
93    
# Line 93  void ToFLevel2::Set(){//ELENA Line 107  void ToFLevel2::Set(){//ELENA
107      if(!ToFTrk)ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA      if(!ToFTrk)ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA
108  }//ELENA  }//ELENA
109    
110  void ToFLevel2::Clear(){  void ToFLevel2::Clear(Option_t *t){
111    //    //
112    if(ToFTrk)ToFTrk->Delete(); //ELENA    if(ToFTrk)ToFTrk->Delete(); //ELENA
113    if(PMT)PMT->Delete(); //ELENA    if(PMT)PMT->Delete(); //ELENA
# Line 102  void ToFLevel2::Clear(){ Line 116  void ToFLevel2::Clear(){
116    //    //
117  };  };
118    
119  void ToFLevel2::Delete(){ //ELENA  void ToFLevel2::Delete(Option_t *t){ //ELENA
120    //    //
121    if(ToFTrk){    if(ToFTrk){
122        ToFTrk->Delete(); //ELENA        ToFTrk->Delete(); //ELENA
# Line 148  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit) Line 162  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
162  //--------------------------------------  //--------------------------------------
163  /**  /**
164   * Method to get the plane ID (11 12 21 22 31 32) from the plane index (0 1 2 3 4 5)   * Method to get the plane ID (11 12 21 22 31 32) from the plane index (0 1 2 3 4 5)
165     * @param Plane index (0,1,2,3,4,5).
166   */   */
167    Int_t  ToFLevel2::GetToFPlaneID(Int_t ip){    Int_t  ToFLevel2::GetToFPlaneID(Int_t ip){
168        if(ip>=0 && ip<6)return 10*((int)(ip/2+1.1))+(ip%2)+1;        if(ip>=0 && ip<6)return 10*((int)(ip/2+1.1))+(ip%2)+1;
# Line 155  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit) Line 170  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
170    };    };
171  /**  /**
172   * Method to get the plane index (0 1 2 3 4 5) from the plane ID (11 12 21 22 31 32)   * Method to get the plane index (0 1 2 3 4 5) from the plane ID (11 12 21 22 31 32)
173     * @param plane Plane ID (11, 12, 21, 22, 31, 32)
174   */   */
175    Int_t  ToFLevel2::GetToFPlaneIndex(Int_t plane_id){    Int_t  ToFLevel2::GetToFPlaneIndex(Int_t plane_id){
176        if(        if(
# Line 169  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit) Line 185  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
185    };    };
186  /**  /**
187   * Method to know if a given ToF paddle was hit, that is there is a TDC signal   * Method to know if a given ToF paddle was hit, that is there is a TDC signal
188   * from both PMTs   * from both PMTs. The method uses the "tof_j_flag" variable.
189   * @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).
190   * @param paddle_id Paddle ID.   * @param paddle_id Paddle ID.
191   * @return 1 if the paddle was hit.   * @return 1 if the paddle was hit.
# Line 201  Int_t ToFLevel2::GetNHitPaddles(Int_t pl Line 217  Int_t ToFLevel2::GetNHitPaddles(Int_t pl
217  };  };
218    
219    
220  Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane){  //gf Apr 07
221    /**
222     * Method to get the mean dEdx from a given ToF plane. This current version
223     * is just summing up all PMT signals, which will not give proper results,
224     *  and needs a revision.
225     * @param notrack Track Number
226     * @param plane Plane index (0,1,2,3,4,5)
227     * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
228     */
229    Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){
230    
231    Float_t dedx = 0.;    Float_t dedx = 0.;
232    Int_t ip = 0;    Float_t PadEdx =0.;
233    Int_t pmt_id = 0;    Int_t SatWarning;
234    Int_t pl = 0;    Int_t pad=-1;
   if ( plane >= 6 ){  
     ip = GetToFPlaneIndex(plane);  
   } else {  
     ip = plane;  
   };  
235    //    //
236    ToFTrkVar *trk = GetToFTrkVar(notrack);    ToFTrkVar *trk = GetToFTrkVar(notrack);
237    if(!trk) return 0; //ELENA    if(!trk) return 0; //ELENA
238    //    //
239    for (Int_t i=0; i<trk->npmtadc; i++){    for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
240      //      Int_t paddleid=ii;
241      pmt_id = (trk->pmtadc).At(i);      pad = GetPaddleid(plane,paddleid);
242      //      GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);
243      pl = GetPlaneIndex(pmt_id);      dedx += PadEdx;
     //  
     if ( pl == ip ) dedx += (trk->dedx).At(i);  
     //  
244    };    };
245    //    //
246    return(dedx);    return(dedx);
247  };  };
248    
249    /**
250     * Method to fill the ADC_C 4x12 matrix with the dEdx values and the TDC 4x12 matrix
251     * with the time-walk corrected TDC values.
252     * @param notrack Track Number
253     * @param adc  ADC_C matrix with dEdx values
254     * @param tdc  TDC matrix
255     */
256  void ToFLevel2::GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]){  void ToFLevel2::GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]){
257    //    //
258    for (Int_t aa=0; aa<4;aa++){    for (Int_t aa=0; aa<4;aa++){
# Line 270  void ToFLevel2::GetMatrix(Int_t notrack, Line 293  void ToFLevel2::GetMatrix(Int_t notrack,
293  };  };
294    
295    
296    /**
297     * Method to get the plane index (0 - 5) for the PMT_ID as input
298     * @param pmt_id  PMT_ID (0 - 47)
299     */
300  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt_id){  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt_id){
301    TString pmtname = GetPMTName(pmt_id);    TString pmtname = GetPMTName(pmt_id);
302    pmtname.Resize(3);    pmtname.Resize(3);
# Line 284  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt Line 310  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt
310  };  };
311    
312    
313    /**
314     * Method to get the PMT_ID if the index (4,12) is given. We have 4 channels on
315     * each of the 12 half-boards, this method decodes which PMT is cables to which
316     * channel.
317     * @param hh Channel
318     * @param kk HalfBoard
319     */
320  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){
321    //    //
322    short tof[4][24] = {    short tof[4][24] = {
# Line 312  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_ Line 345  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_
345    return ind;    return ind;
346  };  };
347    
 TString ToFLevel2::GetPMTName(Int_t ind){  
     
   TString pmtname = " ";  
     
   TString photoS[48] = {  
     "S11_1A", "S11_1B", "S11_2A", "S11_2B", "S11_3A", "S11_3B", "S11_4A", "S11_4B",  
     "S11_5A", "S11_5B", "S11_6A", "S11_6B", "S11_7A", "S11_7B", "S11_8A", "S11_8B",  
     "S12_1A", "S12_1B", "S12_2A", "S12_2B", "S12_3A", "S12_3B", "S12_4A", "S12_4B", "S12_5A",  "S12_5B", "S12_6A", "S12_6B",  
     "S21_1A", "S21_1B", "S21_2A", "S21_2B",  
     "S22_1A", "S22_1B", "S22_2A", "S22_2B",  
     "S31_1A", "S31_1B", "S31_2A", "S31_2B", "S31_3A", "S31_3B",  
     "S32_1A", "S32_1B", "S32_2A", "S32_2B", "S32_3A", "S32_3B"  
   };  
     
   
   pmtname = photoS[ind].Data();  
   
   return pmtname;  
 };  
   
348    
349    /**
350     * Method to get the PMT index if the PMT ID is given. This method is the
351     * "reverse" of method "GetPMTid"
352     * @param ind  PMT_ID (0 - 47)
353     * @param hb   HalfBoard
354     * @param ch   Channel
355     */
356  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){
357    //    //
358    short tof[4][24] = {    short tof[4][24] = {
# Line 359  void ToFLevel2::GetPMTIndex(Int_t ind, I Line 379  void ToFLevel2::GetPMTIndex(Int_t ind, I
379    return;    return;
380  };  };
381    
382    
383    
384    
385    /// gf Apr 07
386    
387    /**
388     * Method to get the dEdx from a given ToF paddle.
389     * @param notrack Track Number
390     * @param Paddle index (0,1,...,23).
391     * @param adcflag in the paddle (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
392     * @param PadEdx dEdx from a given ToF paddle
393     * @param SatWarning 1 if the PMT ios near saturation region (adcraw ~3000)
394     */
395    void ToFLevel2::GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){
396    
397      PadEdx = 0.;
398      SatWarning = 1000;
399    
400      Float_t dEdx[48] = {0};
401      Int_t pmt_id = -1;
402      Float_t adcraw[48];
403      //
404      ToFTrkVar *trk = GetToFTrkVar(notrack);
405      if(!trk) return; //ELENA
406      //
407    
408      Int_t pmtleft=-1;
409      Int_t pmtright=-1;
410      GetPaddlePMT(paddleid, pmtleft, pmtright);
411    
412      adcraw[pmtleft] = 4095;
413      adcraw[pmtright] = 4095;
414    
415      
416      for (Int_t jj=0; jj<npmt(); jj++){
417        
418        ToFPMT *pmt = GetToFPMT(jj);
419        if(!pmt)break; //ELENA
420        
421        pmt_id = pmt->pmt_id;
422        if(pmt_id==pmtleft){
423          adcraw[pmtleft] = pmt->adc;
424        }
425        
426        if(pmt_id==pmtright){
427          adcraw[pmtright] = pmt->adc;
428        }
429      }
430      
431      for (Int_t i=0; i<trk->npmtadc; i++){
432    
433        if((trk->adcflag).At(i)==0 || adcfl==100){
434          if((trk->pmtadc).At(i) == pmtleft)dEdx[pmtleft] = (trk->dedx).At(i);
435          if((trk->pmtadc).At(i) == pmtright)dEdx[pmtright] = (trk->dedx).At(i);
436        }else{
437          if((trk->pmtadc).At(i) == pmtleft)dEdx[pmtleft] = 0.;
438          if((trk->pmtadc).At(i) == pmtright)dEdx[pmtright] = 0.;
439        }
440      }
441    
442      if( adcraw[pmtleft] >3000 || adcraw[pmtright] >3000)SatWarning=1;
443        
444      if(dEdx[pmtleft]!=0 && dEdx[pmtright]!=0){
445        PadEdx = (dEdx[pmtleft]+dEdx[pmtright])*0.5;
446      }
447      if(dEdx[pmtleft]==0 && dEdx[pmtright]!=0){
448        PadEdx = dEdx[pmtright];
449      }
450      if(dEdx[pmtleft]!=0 && dEdx[pmtright]==0){
451        PadEdx = dEdx[pmtleft];
452      }
453      
454      return;
455    };
456    //
457    
458    
459    // gf Apr 07
460    
461    /**
462     * Method to get the PMT name (like "S11_1A") if the PMT_ID is given.
463     * Indexes of corresponding  plane, paddle and  pmt are also given as output.
464     * @param ind  PMT_ID (0 - 47)
465     * @param iplane plane index (0 - 5)
466     * @param ipaddle paddle index (relative to the plane)
467     * @param ipmt pmt index (0(A), 1(B))
468     */
469    TString ToFLevel2::GetPMTName(Int_t ind, Int_t &iplane, Int_t &ipaddle,Int_t &ipmt){
470      
471      TString pmtname = " ";
472      
473      TString photoS[48] = {
474        "S11_1A", "S11_1B", "S11_2A", "S11_2B", "S11_3A", "S11_3B", "S11_4A",
475        "S11_4B",
476        "S11_5A", "S11_5B", "S11_6A", "S11_6B", "S11_7A", "S11_7B", "S11_8A",
477        "S11_8B",
478        "S12_1A", "S12_1B", "S12_2A", "S12_2B", "S12_3A", "S12_3B", "S12_4A",
479        "S12_4B", "S12_5A",  "S12_5B", "S12_6A", "S12_6B",
480        "S21_1A", "S21_1B", "S21_2A", "S21_2B",
481        "S22_1A", "S22_1B", "S22_2A", "S22_2B",
482        "S31_1A", "S31_1B", "S31_2A", "S31_2B", "S31_3A", "S31_3B",
483        "S32_1A", "S32_1B", "S32_2A", "S32_2B", "S32_3A", "S32_3B"
484      };
485      
486      
487      pmtname = photoS[ind].Data();
488      
489      TString ss = pmtname(1,2);
490      iplane  = (int)(atoi(ss.Data())/10)*2-3+atoi(ss.Data())%10;
491      ss = pmtname(4);
492      ipaddle = atoi(ss.Data())-1 ;
493      if( pmtname.Contains("A") )ipmt=0;
494      if( pmtname.Contains("B") )ipmt=1;
495      
496      return pmtname;
497    };
498    /**
499     * Method to get the PMT name (like "S11_1A") if the PMT_ID is given
500     * @param ind  PMT_ID (0 - 47)
501     */
502    TString ToFLevel2::GetPMTName(Int_t ind){
503    
504      Int_t iplane  = -1;
505      Int_t ipaddle = -1;
506      Int_t ipmt    = -1;
507      return GetPMTName(ind,iplane,ipaddle,ipmt);
508      
509    };
510    
511    
512    // gf Apr 07
513    Int_t ToFLevel2::GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane){
514    
515      Double_t xt,yt,xl,xh,yl,yh;
516      
517      Float_t tof11_x[8] = {-17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85};
518      Float_t tof12_y[6] = { -13.75,-8.25,-2.75,2.75,8.25,13.75};
519      Float_t tof21_y[2] = { 3.75,-3.75};
520      Float_t tof22_x[2] = { -4.5,4.5};
521      Float_t tof31_x[3] = { -6.0,0.,6.0};
522      Float_t tof32_y[3] = { -5.0,0.0,5.0};
523      
524      //  S11 8 paddles  33.0 x 5.1 cm
525      //  S12 6 paddles  40.8 x 5.5 cm
526      //  S21 2 paddles  18.0 x 7.5 cm
527      //  S22 2 paddles  15.0 x 9.0 cm
528      //  S31 3 paddles  15.0 x 6.0 cm
529      //  S32 3 paddles  18.0 x 5.0 cm
530      
531      Int_t paddleidoftrack=-1;
532      //
533      
534      //--- S11 ------
535      
536      if(plane==0){
537        xt = xtr;
538        yt = ytr;
539        paddleidoftrack=-1;
540        yl = -33.0/2. ;
541        yh =  33.0/2. ;
542        if ((yt>yl)&&(yt<yh)) {
543          for (Int_t i1=0; i1<8;i1++){
544            xl = tof11_x[i1] - (5.1-0.4)/2. ;
545            xh = tof11_x[i1] + (5.1-0.4)/2. ;
546            if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;
547          }
548        }
549      }
550      //      cout<<"S11  "<<paddleidoftrack[0]<<"\n";
551      
552      //--- S12 -------
553      if(plane==1){
554        xt = xtr;
555        yt = ytr;
556        paddleidoftrack=-1;
557        xl = -40.8/2. ;
558        xh =  40.8/2. ;
559        
560        if ((xt>xl)&&(xt<xh)) {
561          for (Int_t i1=0; i1<6;i1++){
562            yl = tof12_y[i1] - (5.5-0.4)/2. ;
563            yh = tof12_y[i1] + (5.5-0.4)/2. ;
564            if ((yt>yl)&&(yt<yh))  paddleidoftrack=i1;
565          }
566        }
567      }
568      
569      //--- S21 ------
570    
571      if(plane==2){
572        xt = xtr;
573        yt = ytr;
574        paddleidoftrack=-1;
575        xl = -18./2. ;
576        xh =  18./2. ;
577        
578        if ((xt>xl)&&(xt<xh)) {
579          for (Int_t i1=0; i1<2;i1++){
580            yl = tof21_y[i1] - (7.5-0.4)/2. ;
581            yh = tof21_y[i1] + (7.5-0.4)/2. ;
582            if ((yt>yl)&&(yt<yh))  paddleidoftrack=i1;
583          }
584        }
585      }
586      
587      //--- S22 ------
588      if(plane==3){
589        xt = xtr;
590        yt = ytr;
591        paddleidoftrack=-1;
592        yl = -15./2. ;
593        yh =  15./2. ;
594        
595        if ((yt>yl)&&(yt<yh)) {
596          for (Int_t i1=0; i1<2;i1++){
597            xl = tof22_x[i1] - (9.0-0.4)/2. ;
598            xh = tof22_x[i1] + (9.0-0.4)/2. ;
599            if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;
600          }
601        }
602      }  
603      
604      //--- S31 ------
605      if(plane==4){
606        xt = xtr;
607        yt = ytr;
608        paddleidoftrack=-1;
609        yl = -15.0/2. ;
610        yh =  15.0/2. ;
611        
612        if ((yt>yl)&&(yt<yh)) {
613          for (Int_t i1=0; i1<3;i1++){
614            xl = tof31_x[i1] - (6.0-0.4)/2. ;
615            xh = tof31_x[i1] + (6.0-0.4)/2. ;
616            if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;
617          }
618        }
619      }  
620      
621      //---  S32 ------
622      if(plane==5){
623        xt = xtr;
624        yt = ytr;
625        paddleidoftrack=-1;
626        xl = -18.0/2. ;
627        xh =  18.0/2. ;
628        
629        if ((xt>xl)&&(xt<xh)) {
630          for (Int_t i1=0; i1<3;i1++){
631            yl = tof32_y[i1] - (5.0-0.4)/2. ;
632            yh = tof32_y[i1] + (5.0-0.4)/2. ;
633            if ((yt>yl)&&(yt<yh)) paddleidoftrack=i1;
634          }
635        }
636      }
637      
638      return paddleidoftrack;
639    
640    }  
641    
642    //
643    
644    // gf Apr 07
645    
646    void ToFLevel2::GetPMTPaddle(Int_t pmt_id, Int_t &plane, Int_t &paddle){
647      
648      plane = GetPlaneIndex(pmt_id);
649    
650      if(plane == 0){
651        if(pmt_id==0 || pmt_id==1)paddle=0;
652        if(pmt_id==2 || pmt_id==3)paddle=1;
653        if(pmt_id==4 || pmt_id==5)paddle=2;
654        if(pmt_id==6 || pmt_id==7)paddle=3;
655        if(pmt_id==8 || pmt_id==9)paddle=4;
656        if(pmt_id==10 || pmt_id==11)paddle=5;
657        if(pmt_id==12 || pmt_id==13)paddle=6;
658        if(pmt_id==14 || pmt_id==15)paddle=7;
659      }
660      
661      if(plane == 1){
662        if(pmt_id==16 || pmt_id==17)paddle=0;
663        if(pmt_id==18 || pmt_id==19)paddle=1;
664        if(pmt_id==20 || pmt_id==21)paddle=2;
665        if(pmt_id==22 || pmt_id==23)paddle=3;
666        if(pmt_id==24 || pmt_id==25)paddle=4;
667        if(pmt_id==26 || pmt_id==27)paddle=5;
668      }
669      
670      if(plane == 2){
671        if(pmt_id==28 || pmt_id==29)paddle=0;
672        if(pmt_id==30 || pmt_id==31)paddle=1;
673      }
674      
675      if(plane == 3){
676        if(pmt_id==32 || pmt_id==33)paddle=0;
677        if(pmt_id==34 || pmt_id==35)paddle=1;
678      }
679      
680      if(plane == 4){
681        if(pmt_id==36 || pmt_id==37)paddle=0;
682        if(pmt_id==38 || pmt_id==39)paddle=1;
683        if(pmt_id==40 || pmt_id==41)paddle=2;
684      }
685      
686      if(plane == 5){
687        if(pmt_id==42 || pmt_id==43)paddle=0;
688        if(pmt_id==44 || pmt_id==45)paddle=1;
689        if(pmt_id==46 || pmt_id==47)paddle=2;
690      }
691      return;
692    }
693    
694    //
695    
696    // gf Apr 07
697    
698    void ToFLevel2::GetPaddlePMT(Int_t paddle, Int_t &pmtleft, Int_t &pmtright){
699    
700      if(paddle==0){
701        pmtleft=0;
702        pmtright=1;
703      }
704    
705      if(paddle==1){
706        pmtleft=2;
707        pmtright=3;
708      }
709    
710      if(paddle==2){
711        pmtleft=4;
712        pmtright=5;
713      }
714    
715      if(paddle==3){
716        pmtleft=6;
717        pmtright=7;
718      }
719    
720      if(paddle==4){
721        pmtleft=8;
722        pmtright=9;
723      }
724    
725      if(paddle==5){
726        pmtleft=10;
727        pmtright=11;
728      }
729    
730      if(paddle==6){
731        pmtleft=12;
732        pmtright=13;
733      }
734    
735      if(paddle==7){
736        pmtleft=14;
737        pmtright=15;
738      }
739    
740      if(paddle==8){
741        pmtleft=16;
742        pmtright=17;
743      }
744    
745      if(paddle==9){
746        pmtleft=18;
747        pmtright=19;
748      }
749    
750      if(paddle==10){
751        pmtleft=20;
752        pmtright=21;
753      }
754    
755      if(paddle==11){
756        pmtleft=22;
757        pmtright=23;
758      }
759    
760      if(paddle==12){
761        pmtleft=24;
762        pmtright=25;
763      }
764    
765      if(paddle==13){
766        pmtleft=26;
767        pmtright=27;
768      }
769    
770      if(paddle==14){
771        pmtleft=28;
772        pmtright=29;
773      }
774    
775      if(paddle==15){
776        pmtleft=30;
777        pmtright=31;
778      }
779    
780      if(paddle==16){
781        pmtleft=32;
782        pmtright=33;
783      }
784    
785      if(paddle==17){
786        pmtleft=34;
787        pmtright=35;
788      }
789    
790      if(paddle==18){
791        pmtleft=36;
792        pmtright=37;
793      }
794    
795      if(paddle==19){
796        pmtleft=38;
797        pmtright=39;
798      }
799    
800      if(paddle==20){
801        pmtleft=40;
802        pmtright=41;
803      }
804    
805      if(paddle==21){
806        pmtleft=42;
807        pmtright=43;
808      }
809    
810      if(paddle==22){
811        pmtleft=44;
812        pmtright=45;
813      }
814    
815      if(paddle==23){
816        pmtleft=46;
817        pmtright=47;
818      }
819      
820      return;
821    }
822    
823    //
824    
825    
826    
827    // // gf Apr 07
828    
829    void ToFLevel2::GetPaddleGeometry(Int_t plane, Int_t paddle, Float_t &xleft, Float_t &xright, Float_t &yleft, Float_t &yright){
830      
831      Int_t i1;
832    
833      Float_t tof11_x[8] = {-17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85};
834      Float_t tof12_y[6] = { -13.75,-8.25,-2.75,2.75,8.25,13.75};
835      Float_t tof21_y[2] = { 3.75,-3.75};
836      Float_t tof22_x[2] = { -4.5,4.5};
837      Float_t tof31_x[3] = { -6.0,0.,6.0};
838      Float_t tof32_y[3] = { -5.0,0.0,5.0};
839            
840      //  S11 8 paddles  33.0 x 5.1 cm
841      //  S12 6 paddles  40.8 x 5.5 cm
842      //  S21 2 paddles  18.0 x 7.5 cm
843      //  S22 2 paddles  15.0 x 9.0 cm
844      //  S31 3 paddles  15.0 x 6.0 cm
845      //  S32 3 paddles  18.0 x 5.0 cm
846    
847      if(plane==0)
848        {
849          for (i1=0; i1<8;i1++){
850            if(i1 == paddle){
851              xleft = tof11_x[i1] - 5.1/2.;
852              xright = tof11_x[i1] + 5.1/2.;
853              yleft = -33.0/2.;
854              yright = 33.0/2.;
855            }
856          }
857        }
858      
859      if(plane==1)
860        {
861          for (i1=0; i1<6;i1++){
862            if(i1 == paddle){
863              xleft = -40.8/2.;
864              xright = 40.8/2.;
865              yleft = tof12_y[i1] - 5.5/2.;
866              yright = tof12_y[i1] + 5.5/2.;
867            }
868          }
869        }
870    
871      if(plane==2)
872        {
873          for (i1=0; i1<2;i1++){
874            if(i1 == paddle){
875              xleft =  -18./2.;
876              xright = 18./2.;
877              yleft = tof21_y[i1] - 7.5/2.;
878              yright = tof21_y[i1] + 7.5/2.;
879            }
880          }
881        }
882      
883      if(plane==3)
884        {
885          for (i1=0; i1<2;i1++){
886            if(i1 == paddle){
887              xleft = tof22_x[i1] - 9.0/2.;
888              xright = tof22_x[i1] + 9.0/2.;
889              yleft = -15./2.;
890              yright = 15./2.;
891            }
892          }
893        }
894    
895    
896      if(plane==4)
897        {
898          for (i1=0; i1<3;i1++){
899            if(i1 == paddle){
900              xleft = tof31_x[i1] - 6.0/2.;
901              xright = tof31_x[i1] + 6.0/2.;
902              yleft = -15./2.;
903              yright = 15./2.;
904            }
905          }
906        }
907    
908      if(plane==5)
909        {
910          for (i1=0; i1<3;i1++){
911            if(i1 == paddle){
912              xleft = -18.0/2.;
913              xright = 18.0/2.;
914              yleft = tof32_y[i1] - 5.0/2.;
915              yright = tof32_y[i1] + 5.0/2.;
916            }
917          }
918        }
919      return;
920    }
921    
922    // gf Apr 07
923    /**
924     * Method to get the paddle index (0,...23) if the plane ID and the paddle id in the plane is given.
925     * This method is the
926     * "reverse" of method "GetPaddlePlane"
927     * @param plane    (0 - 5)
928     * @param paddle   (plane=0, paddle = 0,...5)
929     * @param padid    (0 - 23)
930     */
931    Int_t ToFLevel2::GetPaddleid(Int_t plane, Int_t paddle)
932    {
933    
934      Int_t padid=-1;
935      Int_t pads11=8;
936      Int_t pads12=6;
937      Int_t pads21=2;
938      Int_t pads22=2;
939      Int_t pads31=3;
940      //  Int_t pads32=3;
941    
942    
943      if(plane == 0){
944        padid=paddle;
945      }
946    
947      if(plane == 1){
948        padid=pads11+paddle;
949      }
950    
951      if(plane == 2){
952        padid=pads11+pads12+paddle;
953      }
954    
955      if(plane == 3){
956        padid=pads11+pads12+pads21+paddle;
957      }
958    
959      if(plane == 4){
960        padid=pads11+pads12+pads21+pads22+paddle;
961      }
962    
963      if(plane == 5){
964        padid=pads11+pads12+pads21+pads22+pads31+paddle;
965      }
966    
967      return padid;
968    
969    }
970    
971    
972    // gf Apr 07
973    /**
974     * Method to get the plane ID and the paddle id in the plane if the paddle index (0,...23) is given.
975     * This method is the
976     * "reverse" of method "GetPaddleid"
977     * @param pad      (0 - 23)
978     * @param plane    (0 - 5)
979     * @param paddle   (plane=0, paddle = 0,...5)
980     */
981    void ToFLevel2::GetPaddlePlane(Int_t pad, Int_t &plane, Int_t &paddle)
982    {
983    
984      Int_t pads11=8;
985      Int_t pads12=6;
986      Int_t pads21=2;
987      Int_t pads22=2;
988      Int_t pads31=3;
989      // Int_t pads32=3;
990    
991      if(pad<8){
992        plane=0;
993        paddle=pad;
994        return;
995      }
996    
997      if(7<pad<14){
998        plane=1;
999        paddle=pad-pads11;
1000        return;
1001      }
1002      
1003      if(13<pad<16){
1004        plane=2;
1005        paddle=pad-pads11-pads12;
1006        return;
1007      }
1008    
1009      if(15<pad<18){
1010        plane=3;
1011        paddle=pad-pads11-pads12-pads21;
1012        return;
1013      }
1014    
1015      if(17<pad<21){
1016        plane=4;
1017        paddle=pad-pads11-pads12-pads21-pads22;
1018        return;
1019      }
1020    
1021      if(20<pad<24){
1022        plane=5;
1023        paddle=pad-pads11-pads12-pads21-pads22-pads31;
1024        return;
1025      }  
1026    
1027    }
1028    
1029    
1030    Int_t ToFLevel2::GetNPaddle(Int_t plane){
1031    
1032      Int_t npaddle=-1;
1033    
1034      Int_t pads11=8;
1035      Int_t pads12=6;
1036      Int_t pads21=2;
1037      Int_t pads22=2;
1038      Int_t pads31=3;
1039      Int_t pads32=3;
1040    
1041      if(plane==0)npaddle=pads11;
1042      if(plane==1)npaddle=pads12;
1043      if(plane==2)npaddle=pads21;
1044      if(plane==3)npaddle=pads22;
1045      if(plane==4)npaddle=pads31;
1046      if(plane==5)npaddle=pads32;
1047    
1048      return npaddle;
1049    
1050    }
1051    
1052    
1053    
1054    /// wm feb 08
1055    
1056    /**
1057     * Method to calculate Beta from the 12 single measurements
1058     * we check the individual weights for artificial TDC values, then calculate
1059     * am mean beta for the first time. In a second step we loop again through
1060     * the single measurements, checking for the residual from the mean
1061     * The cut on the residual reject measurements > "x"-sigma. A chi2 value is
1062     * calculated, furthermore a "quality" value by adding the weights which
1063     * are finally used. If all measurements are taken, "quality" will be = 22.47.
1064     * A chi2 cut around 3-4 and a quality-cut > 20 is needed for clean beta
1065     * measurements like antiprotons etc.
1066     * The Level2 output is derived in the fortran routines using: 10.,10.,20.
1067     * @param notrack Track Number
1068     * @param cut on residual: difference between single measurement and mean
1069     * @param cut on "quality"
1070     * @param cut on chi2
1071     */
1072    
1073    Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){
1074    
1075    //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;
1076    
1077      Float_t bxx = 100.;
1078      //
1079      ToFTrkVar *trk = GetToFTrkVar(notrack);
1080      if(!trk) return 0; //ELENA
1081    
1082    
1083      Float_t chi2,xhelp,beta_mean;
1084      Float_t w_i[12],quality,sw,sxw,res,betachi,beta_mean_inv;
1085      Float_t b[12],tdcfl;
1086      Int_t  pmt_id,pmt_plane;
1087    
1088      for (Int_t i=0; i<12; i++){
1089        b[i] = trk->beta[i];
1090                                  }
1091          
1092    
1093    //========================================================================
1094    //---  Find out ToF layers with artificial TDC values & fill vector    ---
1095    //========================================================================
1096    
1097    Float_t  w_il[6];
1098    
1099         for (Int_t jj=0; jj<6;jj++) {
1100             w_il[jj] = 1000.;
1101                                     }
1102    
1103    
1104      for (Int_t i=0; i<trk->npmttdc; i++){
1105        //
1106        pmt_id = (trk->pmttdc).At(i);
1107        pmt_plane = GetPlaneIndex(pmt_id);
1108        tdcfl = (trk->tdcflag).At(i);
1109        if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag
1110                                         };
1111      
1112    //========================================================================
1113    //---  Set weights for the 12 measurements using information for top and bottom:
1114    //---  if no measurements: weight = set to very high value=> not used
1115    //---  top or bottom artificial: weight*sqrt(2)
1116    //---  top and bottom artificial: weight*sqrt(2)*sqrt(2)
1117    //========================================================================
1118    
1119    Int_t itop[12] = {0,0,1,1,2,2,3,3,0,0,1,1};
1120    Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2,3};
1121    
1122         xhelp= 1E09;
1123      
1124         for (Int_t jj=0; jj<12;jj++) {
1125         if (jj<4)           xhelp = 0.11;    // S1-S3
1126         if ((jj>3)&&(jj<8)) xhelp = 0.18;    // S2-S3
1127         if (jj>7)           xhelp = 0.28;    // S1-S2
1128         if ((w_il[itop[jj]] == 1000.) && (w_il[ibot[jj]] == 1000.)) xhelp = 1E09;
1129         if ((w_il[itop[jj]] == 1) || (w_il[ibot[jj]] == 1.)) xhelp = xhelp*1.414 ;
1130         if ((w_il[itop[jj]] == 1) && (w_il[ibot[jj]] == 1.)) xhelp = xhelp*2. ;
1131    
1132         w_i[jj] = 1./xhelp;
1133                                      }
1134    
1135    
1136    //========================================================================
1137    //--- Calculate mean beta for the first time -----------------------------
1138    //--- We are using "1/beta" since its error is gaussian ------------------
1139    //========================================================================
1140    
1141          Int_t icount=0;
1142          sw=0.;
1143          sxw=0.;
1144          beta_mean=100.;
1145    
1146              for (Int_t jj=0; jj<12;jj++){
1147            if ((fabs(1./b[jj])>0.1)&&(fabs(1./b[jj])<15.))
1148             {
1149                icount= icount+1;
1150                sxw=sxw + (1./b[jj])*w_i[jj]*w_i[jj] ;
1151                sw =sw + w_i[jj]*w_i[jj] ;
1152    
1153             }
1154             }
1155    
1156          if (icount>0) beta_mean=1./(sxw/sw);
1157          beta_mean_inv = 1./beta_mean;
1158    
1159    //========================================================================
1160    //--- Calculate beta for the second time, use residuals of the single
1161    //--- measurements to get a chi2 value
1162    //========================================================================
1163    
1164          icount=0;
1165          sw=0.;
1166          sxw=0.;
1167          betachi = 100.;
1168          chi2 = 0.;
1169          quality=0.;
1170    
1171    
1172              for (Int_t jj=0; jj<12;jj++){
1173           if ((fabs(1./b[jj])>0.1)&&(fabs(1./b[jj])<15.)&&(w_i[jj]>0.01)) {
1174                res = beta_mean_inv - (1./b[jj]) ;
1175                if (fabs(res*w_i[jj])<resmax)          {;
1176                chi2 = chi2 + pow((res*w_i[jj]),2) ;
1177                icount= icount+1;
1178                sxw=sxw + (1./b[jj])*w_i[jj]*w_i[jj] ;
1179                sw =sw + w_i[jj]*w_i[jj] ;
1180                                                   }
1181                                                                            }
1182                                          }
1183          quality = sqrt(sw) ;
1184    
1185          if (icount==0) chi2 = 1000.;
1186          if (icount>0) chi2 = chi2/(icount) ;
1187          if (icount>0) betachi=1./(sxw/sw);
1188    
1189       bxx = 100.;
1190       if ((chi2 < chi2cut)&&(quality>qualitycut)) bxx = betachi;
1191      //
1192      return(bxx);
1193    };
1194    
1195    
1196    ////////////////////////////////////////////////////
1197    ////////////////////////////////////////////////////
1198    
1199    
1200  /**  /**
1201   * Fills a struct cToFLevel2 with values from a ToFLevel2 object (to put data into a F77 common).   * Fills a struct cToFLevel2 with values from a ToFLevel2 object (to put data into a F77 common).
1202   */   */
# Line 389  void ToFLevel2::GetLevel2Struct(cToFLeve Line 1227  void ToFLevel2::GetLevel2Struct(cToFLeve
1227                l2->xtofpos[i][j]=((ToFTrkVar*)ToFTrk->At(j))->xtofpos[i];                l2->xtofpos[i][j]=((ToFTrkVar*)ToFTrk->At(j))->xtofpos[i];
1228                l2->ytofpos[i][j]=((ToFTrkVar*)ToFTrk->At(j))->ytofpos[i];                l2->ytofpos[i][j]=((ToFTrkVar*)ToFTrk->At(j))->ytofpos[i];
1229            }            }
1230              for(Int_t i=0;i<6;i++){
1231                  l2->xtr_tof[i][j]=((ToFTrkVar*)ToFTrk->At(j))->xtr_tof[i];
1232                  l2->ytr_tof[i][j]=((ToFTrkVar*)ToFTrk->At(j))->ytr_tof[i];
1233              }
1234        }        }
1235    } //ELENA    } //ELENA
1236            
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.19
  ViewVC Help
Powered by ViewVC 1.1.23