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

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

Parent Directory | Revision Log | View Patch Patch

-revision 1.13 by mocchiut,
Mon Jan 22 14:09:32 2007 UTC
+revision 1.23 by pamelats,
Thu Dec  4 11:41:08 2008 UTC
 Line 1
+ /**
+  * \file ToFLevel2.cpp
+  * \author Gianfranca DeRosa, Wolfgang Menn
+  *
+  * WM dec 2008: Description of "GetdEdx" changed
+  * WM dec 2008: "GetdEdxPaddle" modified: Now includes saturation limit
+  *              PMTs higher than the saturation limit are not used for dEdx
+  */
  #include <TObject.h>
  #include <ToFLevel2.h>
  #include <iostream>
-Line 10 
 ToFPMT::ToFPMT(){
+Line 19 
 ToFPMT::ToFPMT(){
    pmt_id = 0;
    adc = 0.;
    tdc_tw = 0.;
+   tdc = 0.;
  }
  ToFPMT::ToFPMT(const ToFPMT &t){
    pmt_id = t.pmt_id;
    adc = t.adc;
    tdc_tw = t.tdc_tw;
+   tdc = t.tdc;
  }
- void ToFPMT::Clear(){
+ void ToFPMT::Clear(Option_t *t){
    pmt_id = 0;
    adc = 0.;
    tdc_tw = 0.;
+   tdc = 0.;
  }
-Line 40 
 ToFTrkVar::ToFTrkVar() {
+Line 52 
 ToFTrkVar::ToFTrkVar() {
    memset(beta,  0, 13*sizeof(Float_t));
    memset(xtofpos,  0, 3*sizeof(Float_t));
    memset(ytofpos,  0, 3*sizeof(Float_t));
+   memset(xtr_tof,  0, 6*sizeof(Float_t));
+   memset(ytr_tof,  0, 6*sizeof(Float_t));
    //
  };
- void ToFTrkVar::Clear() {
+ void ToFTrkVar::Clear(Option_t *t) {
    trkseqno = 0;
    npmttdc = 0;
    npmtadc = 0;
-Line 56 
 void ToFTrkVar::Clear() {
+Line 70 
 void ToFTrkVar::Clear() {
    memset(beta,  0, 13*sizeof(Float_t));
    memset(xtofpos,  0, 3*sizeof(Float_t));
    memset(ytofpos,  0, 3*sizeof(Float_t));
+   memset(xtr_tof,  0, 6*sizeof(Float_t));
+   memset(ytr_tof,  0, 6*sizeof(Float_t));
    //
  };
-Line 74 
 ToFTrkVar::ToFTrkVar(const ToFTrkVar &t)
+Line 90 
 ToFTrkVar::ToFTrkVar(const ToFTrkVar &t)
    memcpy(beta,t.beta,sizeof(beta));
    memcpy(xtofpos,t.xtofpos,sizeof(xtofpos));
    memcpy(ytofpos,t.ytofpos,sizeof(ytofpos));
+   memcpy(xtr_tof,t.xtr_tof,sizeof(xtr_tof));
+   memcpy(ytr_tof,t.ytr_tof,sizeof(ytr_tof));
    //
  };
-Line 93 
 void ToFLevel2::Set(){//ELENA
+Line 111 
 void ToFLevel2::Set(){//ELENA
      if(!ToFTrk)ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA
  }//ELENA
- void ToFLevel2::Clear(){
+ void ToFLevel2::Clear(Option_t *t){
    //
    if(ToFTrk)ToFTrk->Delete(); //ELENA
    if(PMT)PMT->Delete(); //ELENA
-Line 102 
 void ToFLevel2::Clear(){
+Line 120 
 void ToFLevel2::Clear(){
    //
  };
- void ToFLevel2::Delete(){ //ELENA
+ void ToFLevel2::Delete(Option_t *t){ //ELENA
    //
    if(ToFTrk){
        ToFTrk->Delete(); //ELENA
-Line 148 
 ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
+Line 166 
 ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
  //--------------------------------------
  /**
   * Method to get the plane ID (11 12 21 22 31 32) from the plane index (0 1 2 3 4 5)
+  * @param Plane index (0,1,2,3,4,5).
   */
    Int_t  ToFLevel2::GetToFPlaneID(Int_t ip){
        if(ip>=0 && ip<6)return 10*((int)(ip/2+1.1))+(ip%2)+1;
-Line 155 
 ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
+Line 174 
 ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
    };
  /**
   * Method to get the plane index (0 1 2 3 4 5) from the plane ID (11 12 21 22 31 32)
+  * @param plane Plane ID (11, 12, 21, 22, 31, 32)
   */
    Int_t  ToFLevel2::GetToFPlaneIndex(Int_t plane_id){
        if(
-Line 169 
 ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
+Line 189 
 ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
    };
  /**
   * Method to know if a given ToF paddle was hit, that is there is a TDC signal
-  * from both PMTs
+  * from both PMTs. The method uses the "tof_j_flag" variable.
   * @param plane Plane ID (11, 12, 21, 22, 31, 32) or Plane index (0,1,2,3,4,5).
   * @param paddle_id Paddle ID.
   * @return 1 if the paddle was hit.
-Line 200 
 Int_t ToFLevel2::GetNHitPaddles(Int_t pl
+Line 220 
 Int_t ToFLevel2::GetNHitPaddles(Int_t pl
      return npad;
  };
+ //wm Nov 08
- Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane){
+ //gf Apr 07
+ /**
+  * Method to get the mean dEdx from a ToF layer - ATTENTION:
+  * It will sum up the dEdx of all the paddles, but since by definition
+  * only the paddle hitted by the track gets a dEdx value and the other
+  * paddles are set to zero, the output is just the dEdx of the hitted
+  * paddle in each layer!
+  * The "adcfl" option is not very useful (an artificial dEdx is per
+  * definition= 1 mip and not a real measurement), anyway left in the code
+  * @param notrack Track Number
+  * @param plane Plane index (0,1,2,3,4,5)
+  * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
+  */
+ Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){
    Float_t dedx = 0.;
-   Int_t ip = 0;
+   Float_t PadEdx =0.;
-   Int_t pmt_id = 0;
+   Int_t SatWarning;
-   Int_t pl = 0;
+   Int_t pad=-1;
-   if ( plane >= 6 ){
-     ip = GetToFPlaneIndex(plane);
-   } else {
-     ip = plane;
-   };
    //
    ToFTrkVar *trk = GetToFTrkVar(notrack);
    if(!trk) return 0; //ELENA
    //
-   for (Int_t i=0; i<trk->npmtadc; i++){
+   for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
-     //
+     Int_t paddleid=ii;
-     pmt_id = (trk->pmtadc).At(i);
+     pad = GetPaddleid(plane,paddleid);
-     //
+     GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);
-     pl = GetPlaneIndex(pmt_id);
+     dedx += PadEdx;
-     //
-     if ( pl == ip ) dedx += (trk->dedx).At(i);
-     //
    };
    //
    return(dedx);
  };
+ /**
+  * Method to fill the ADC_C 4x12 matrix with the dEdx values and the TDC 4x12 matrix
+  * with the time-walk corrected TDC values.
+  * @param notrack Track Number
+  * @param adc  ADC_C matrix with dEdx values
+  * @param tdc  TDC matrix
+  */
  void ToFLevel2::GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]){
    //
    for (Int_t aa=0; aa<4;aa++){
-Line 270 
 void ToFLevel2::GetMatrix(Int_t notrack,
+Line 301 
 void ToFLevel2::GetMatrix(Int_t notrack,
  };
+ /**
+  * Method to get the plane index (0 - 5) for the PMT_ID as input
+  * @param pmt_id  PMT_ID (0 - 47)
+  */
  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt_id){
    TString pmtname = GetPMTName(pmt_id);
    pmtname.Resize(3);
-Line 284 
 Int_t ToFLevel2::GetPlaneIndex(Int_t pmt
+Line 318 
 Int_t ToFLevel2::GetPlaneIndex(Int_t pmt
  };
+ /**
+  * Method to get the PMT_ID if the index (4,12) is given. We have 4 channels on
+  * each of the 12 half-boards, this method decodes which PMT is cables to which
+  * channel.
+  * @param hh Channel
+  * @param kk HalfBoard
+  */
  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){
    //
    short tof[4][24] = {
-Line 312 
 Int_t ToFLevel2::GetPMTid(Int_t hh, Int_
+Line 353 
 Int_t ToFLevel2::GetPMTid(Int_t hh, Int_
    return ind;
  };
- 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;
- };
+ /**
+  * Method to get the PMT index if the PMT ID is given. This method is the
+  * "reverse" of method "GetPMTid"
+  * @param ind  PMT_ID (0 - 47)
+  * @param hb   HalfBoard
+  * @param ch   Channel
+  */
  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){
    //
    short tof[4][24] = {
-Line 359 
 void ToFLevel2::GetPMTIndex(Int_t ind, I
+Line 387 
 void ToFLevel2::GetPMTIndex(Int_t ind, I
    return;
  };
+ //  wm Nov 08 revision - saturation values included
+ /// gf Apr 07
+ /**
+  * Method to get the dEdx from a given ToF paddle.
+  * If two PMTs are good, the mean dEdx of both PMTs is taken, otherwise
+  * just the dEdx of the "good" PMT. If both PMTs are above saturation => dEdx=1000
+  * @param notrack Track Number
+  * @param Paddle index (0,1,...,23).
+  * @param adcflag in the paddle (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
+  * @param PadEdx dEdx from a given ToF paddle
+  * @param SatWarning 1 if the PMT ios near saturation region (adcraw ~3000)
+  */
+ void ToFLevel2::GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){
+ /*
+ Float_t  PMTsat[48] = {
+.14, 3165.48, 3153.85, 3085.73, 3089.65, 3107.64, 3097.52, 3078.37,
+.05, 3087.07, 3112.22, 3102.92, 3080.58, 3092.55, 3087.94, 3125.03,
+.09, 3143.16, 3125.51, 3181.27, 3092.09, 3124.98, 3069.3, 3095.53,
+.11, 3133.53, 3114.73, 3113.01, 3091.19, 3097.99, 3033.84, 3134.98,
+.37, 3111.04, 3066.77, 3108.17, 3133, 3111.06, 3052.52, 3140.66,
+.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] = {
+.35,3178.19,3167.38,3099.73,3117.00,3126.29,3111.44,3092.27,
+.48,3094.41,3132.13,3115.37,3099.32,3110.97,3111.80,3143.14,
+.72,3153.44,3136.00,3188.96,3104.73,3140.45,3073.18,3106.62,
+.48,3146.92,3127.24,3136.52,3109.59,3112.89,3045.15,3147.26,
+.92,3121.05,3083.25,3123.62,3150.92,3125.30,3067.60,3160.18,
+.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];
+   //
+   ToFTrkVar *trk = GetToFTrkVar(notrack);
+   if(!trk) return; //ELENA
+   //
+   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];
+ };
+ //
+ // gf Apr 07
+ /**
+  * Method to get the PMT name (like "S11_1A") if the PMT_ID is given.
+  * Indexes of corresponding  plane, paddle and  pmt are also given as output.
+  * @param ind  PMT_ID (0 - 47)
+  * @param iplane plane index (0 - 5)
+  * @param ipaddle paddle index (relative to the plane)
+  * @param ipmt pmt index (0(A), 1(B))
+  */
+ TString ToFLevel2::GetPMTName(Int_t ind, Int_t &iplane, Int_t &ipaddle,Int_t &ipmt){
+   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();
+   TString ss = pmtname(1,2);
+   iplane  = (int)(atoi(ss.Data())/10)*2-3+atoi(ss.Data())%10;
+   ss = pmtname(4);
+   ipaddle = atoi(ss.Data())-1 ;
+   if( pmtname.Contains("A") )ipmt=0;
+   if( pmtname.Contains("B") )ipmt=1;
+   return pmtname;
+ };
+ /**
+  * Method to get the PMT name (like "S11_1A") if the PMT_ID is given
+  * @param ind  PMT_ID (0 - 47)
+  */
+ TString ToFLevel2::GetPMTName(Int_t ind){
+   Int_t iplane  = -1;
+   Int_t ipaddle = -1;
+   Int_t ipmt    = -1;
+   return GetPMTName(ind,iplane,ipaddle,ipmt);
+ };
+ // wm jun 08
+ Int_t ToFLevel2::GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane){
+ return GetPaddleIdOfTrack(xtr ,ytr ,plane, 0.4);
+ }
+ // gf Apr 07
+ Int_t ToFLevel2::GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane, Float_t margin){
+   Double_t xt,yt,xl,xh,yl,yh;
+   Float_t tof11_x[8] = {-17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85};
+   Float_t tof12_y[6] = { -13.75,-8.25,-2.75,2.75,8.25,13.75};
+   Float_t tof21_y[2] = { 3.75,-3.75};
+   Float_t tof22_x[2] = { -4.5,4.5};
+   Float_t tof31_x[3] = { -6.0,0.,6.0};
+   Float_t tof32_y[3] = { -5.0,0.0,5.0};
+   //  S11 8 paddles  33.0 x 5.1 cm
+   //  S12 6 paddles  40.8 x 5.5 cm
+   //  S21 2 paddles  18.0 x 7.5 cm
+   //  S22 2 paddles  15.0 x 9.0 cm
+   //  S31 3 paddles  15.0 x 6.0 cm
+   //  S32 3 paddles  18.0 x 5.0 cm
+   Int_t paddleidoftrack=-1;
+   //
+   //--- S11 ------
+   if(plane==0){
+     xt = xtr;
+     yt = ytr;
+     paddleidoftrack=-1;
+     yl = -33.0/2. ;
+     yh =  33.0/2. ;
+     if ((yt>yl)&&(yt<yh)) {
+       for (Int_t i1=0; i1<8;i1++){
+         xl = tof11_x[i1] - (5.1-margin)/2. ;
+         xh = tof11_x[i1] + (5.1-margin)/2. ;
+         if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;
+       }
+     }
+   }
+   //      cout<<"S11  "<<paddleidoftrack[0]<<"\n";
+   //--- S12 -------
+   if(plane==1){
+     xt = xtr;
+     yt = ytr;
+     paddleidoftrack=-1;
+     xl = -40.8/2. ;
+     xh =  40.8/2. ;
+     if ((xt>xl)&&(xt<xh)) {
+       for (Int_t i1=0; i1<6;i1++){
+         yl = tof12_y[i1] - (5.5-margin)/2. ;
+         yh = tof12_y[i1] + (5.5-margin)/2. ;
+         if ((yt>yl)&&(yt<yh))  paddleidoftrack=i1;
+       }
+     }
+   }
+   //--- S21 ------
+   if(plane==2){
+     xt = xtr;
+     yt = ytr;
+     paddleidoftrack=-1;
+     xl = -18./2. ;
+     xh =  18./2. ;
+     if ((xt>xl)&&(xt<xh)) {
+       for (Int_t i1=0; i1<2;i1++){
+         yl = tof21_y[i1] - (7.5-margin)/2. ;
+         yh = tof21_y[i1] + (7.5-margin)/2. ;
+         if ((yt>yl)&&(yt<yh))  paddleidoftrack=i1;
+       }
+     }
+   }
+   //--- S22 ------
+   if(plane==3){
+     xt = xtr;
+     yt = ytr;
+     paddleidoftrack=-1;
+     yl = -15./2. ;
+     yh =  15./2. ;
+     if ((yt>yl)&&(yt<yh)) {
+       for (Int_t i1=0; i1<2;i1++){
+         xl = tof22_x[i1] - (9.0-margin)/2. ;
+         xh = tof22_x[i1] + (9.0-margin)/2. ;
+         if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;
+       }
+     }
+   }
+   //--- S31 ------
+   if(plane==4){
+     xt = xtr;
+     yt = ytr;
+     paddleidoftrack=-1;
+     yl = -15.0/2. ;
+     yh =  15.0/2. ;
+     if ((yt>yl)&&(yt<yh)) {
+       for (Int_t i1=0; i1<3;i1++){
+         xl = tof31_x[i1] - (6.0-margin)/2. ;
+         xh = tof31_x[i1] + (6.0-margin)/2. ;
+         if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;
+       }
+     }
+   }
+   //---  S32 ------
+   if(plane==5){
+     xt = xtr;
+     yt = ytr;
+     paddleidoftrack=-1;
+     xl = -18.0/2. ;
+     xh =  18.0/2. ;
+     if ((xt>xl)&&(xt<xh)) {
+       for (Int_t i1=0; i1<3;i1++){
+         yl = tof32_y[i1] - (5.0-margin)/2. ;
+         yh = tof32_y[i1] + (5.0-margin)/2. ;
+         if ((yt>yl)&&(yt<yh)) paddleidoftrack=i1;
+       }
+     }
+   }
+   return paddleidoftrack;
+ }
+ //
+ // gf Apr 07
+ void ToFLevel2::GetPMTPaddle(Int_t pmt_id, Int_t &plane, Int_t &paddle){
+   plane = GetPlaneIndex(pmt_id);
+   if(plane == 0){
+     if(pmt_id==0 || pmt_id==1)paddle=0;
+     if(pmt_id==2 || pmt_id==3)paddle=1;
+     if(pmt_id==4 || pmt_id==5)paddle=2;
+     if(pmt_id==6 || pmt_id==7)paddle=3;
+     if(pmt_id==8 || pmt_id==9)paddle=4;
+     if(pmt_id==10 || pmt_id==11)paddle=5;
+     if(pmt_id==12 || pmt_id==13)paddle=6;
+     if(pmt_id==14 || pmt_id==15)paddle=7;
+   }
+   if(plane == 1){
+     if(pmt_id==16 || pmt_id==17)paddle=0;
+     if(pmt_id==18 || pmt_id==19)paddle=1;
+     if(pmt_id==20 || pmt_id==21)paddle=2;
+     if(pmt_id==22 || pmt_id==23)paddle=3;
+     if(pmt_id==24 || pmt_id==25)paddle=4;
+     if(pmt_id==26 || pmt_id==27)paddle=5;
+   }
+   if(plane == 2){
+     if(pmt_id==28 || pmt_id==29)paddle=0;
+     if(pmt_id==30 || pmt_id==31)paddle=1;
+   }
+   if(plane == 3){
+     if(pmt_id==32 || pmt_id==33)paddle=0;
+     if(pmt_id==34 || pmt_id==35)paddle=1;
+   }
+   if(plane == 4){
+     if(pmt_id==36 || pmt_id==37)paddle=0;
+     if(pmt_id==38 || pmt_id==39)paddle=1;
+     if(pmt_id==40 || pmt_id==41)paddle=2;
+   }
+   if(plane == 5){
+     if(pmt_id==42 || pmt_id==43)paddle=0;
+     if(pmt_id==44 || pmt_id==45)paddle=1;
+     if(pmt_id==46 || pmt_id==47)paddle=2;
+   }
+   return;
+ }
+ //
+ // gf Apr 07
+ void ToFLevel2::GetPaddlePMT(Int_t paddle, Int_t &pmtleft, Int_t &pmtright){
+   if(paddle==0){
+     pmtleft=0;
+     pmtright=1;
+   }
+   if(paddle==1){
+     pmtleft=2;
+     pmtright=3;
+   }
+   if(paddle==2){
+     pmtleft=4;
+     pmtright=5;
+   }
+   if(paddle==3){
+     pmtleft=6;
+     pmtright=7;
+   }
+   if(paddle==4){
+     pmtleft=8;
+     pmtright=9;
+   }
+   if(paddle==5){
+     pmtleft=10;
+     pmtright=11;
+   }
+   if(paddle==6){
+     pmtleft=12;
+     pmtright=13;
+   }
+   if(paddle==7){
+     pmtleft=14;
+     pmtright=15;
+   }
+   if(paddle==8){
+     pmtleft=16;
+     pmtright=17;
+   }
+   if(paddle==9){
+     pmtleft=18;
+     pmtright=19;
+   }
+   if(paddle==10){
+     pmtleft=20;
+     pmtright=21;
+   }
+   if(paddle==11){
+     pmtleft=22;
+     pmtright=23;
+   }
+   if(paddle==12){
+     pmtleft=24;
+     pmtright=25;
+   }
+   if(paddle==13){
+     pmtleft=26;
+     pmtright=27;
+   }
+   if(paddle==14){
+     pmtleft=28;
+     pmtright=29;
+   }
+   if(paddle==15){
+     pmtleft=30;
+     pmtright=31;
+   }
+   if(paddle==16){
+     pmtleft=32;
+     pmtright=33;
+   }
+   if(paddle==17){
+     pmtleft=34;
+     pmtright=35;
+   }
+   if(paddle==18){
+     pmtleft=36;
+     pmtright=37;
+   }
+   if(paddle==19){
+     pmtleft=38;
+     pmtright=39;
+   }
+   if(paddle==20){
+     pmtleft=40;
+     pmtright=41;
+   }
+   if(paddle==21){
+     pmtleft=42;
+     pmtright=43;
+   }
+   if(paddle==22){
+     pmtleft=44;
+     pmtright=45;
+   }
+   if(paddle==23){
+     pmtleft=46;
+     pmtright=47;
+   }
+   return;
+ }
+ //
+ // // gf Apr 07
+ void ToFLevel2::GetPaddleGeometry(Int_t plane, Int_t paddle, Float_t &xleft, Float_t &xright, Float_t &yleft, Float_t &yright){
+   Int_t i1;
+   Float_t tof11_x[8] = {-17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85};
+   Float_t tof12_y[6] = { -13.75,-8.25,-2.75,2.75,8.25,13.75};
+   Float_t tof21_y[2] = { 3.75,-3.75};
+   Float_t tof22_x[2] = { -4.5,4.5};
+   Float_t tof31_x[3] = { -6.0,0.,6.0};
+   Float_t tof32_y[3] = { -5.0,0.0,5.0};
+   //  S11 8 paddles  33.0 x 5.1 cm
+   //  S12 6 paddles  40.8 x 5.5 cm
+   //  S21 2 paddles  18.0 x 7.5 cm
+   //  S22 2 paddles  15.0 x 9.0 cm
+   //  S31 3 paddles  15.0 x 6.0 cm
+   //  S32 3 paddles  18.0 x 5.0 cm
+   if(plane==0)
+     {
+       for (i1=0; i1<8;i1++){
+         if(i1 == paddle){
+           xleft = tof11_x[i1] - 5.1/2.;
+           xright = tof11_x[i1] + 5.1/2.;
+           yleft = -33.0/2.;
+           yright = 33.0/2.;
+         }
+       }
+     }
+   if(plane==1)
+     {
+       for (i1=0; i1<6;i1++){
+         if(i1 == paddle){
+           xleft = -40.8/2.;
+           xright = 40.8/2.;
+           yleft = tof12_y[i1] - 5.5/2.;
+           yright = tof12_y[i1] + 5.5/2.;
+         }
+       }
+     }
+   if(plane==2)
+     {
+       for (i1=0; i1<2;i1++){
+         if(i1 == paddle){
+           xleft =  -18./2.;
+           xright = 18./2.;
+           yleft = tof21_y[i1] - 7.5/2.;
+           yright = tof21_y[i1] + 7.5/2.;
+         }
+       }
+     }
+   if(plane==3)
+     {
+       for (i1=0; i1<2;i1++){
+         if(i1 == paddle){
+           xleft = tof22_x[i1] - 9.0/2.;
+           xright = tof22_x[i1] + 9.0/2.;
+           yleft = -15./2.;
+           yright = 15./2.;
+         }
+       }
+     }
+   if(plane==4)
+     {
+       for (i1=0; i1<3;i1++){
+         if(i1 == paddle){
+           xleft = tof31_x[i1] - 6.0/2.;
+           xright = tof31_x[i1] + 6.0/2.;
+           yleft = -15./2.;
+           yright = 15./2.;
+         }
+       }
+     }
+   if(plane==5)
+     {
+       for (i1=0; i1<3;i1++){
+         if(i1 == paddle){
+           xleft = -18.0/2.;
+           xright = 18.0/2.;
+           yleft = tof32_y[i1] - 5.0/2.;
+           yright = tof32_y[i1] + 5.0/2.;
+         }
+       }
+     }
+   return;
+ }
+ // gf Apr 07
+ /**
+  * Method to get the paddle index (0,...23) if the plane ID and the paddle id in the plane is given.
+  * This method is the
+  * "reverse" of method "GetPaddlePlane"
+  * @param plane    (0 - 5)
+  * @param paddle   (plane=0, paddle = 0,...5)
+  * @param padid    (0 - 23)
+  */
+ Int_t ToFLevel2::GetPaddleid(Int_t plane, Int_t paddle)
+ {
+   Int_t padid=-1;
+   Int_t pads11=8;
+   Int_t pads12=6;
+   Int_t pads21=2;
+   Int_t pads22=2;
+   Int_t pads31=3;
+   //  Int_t pads32=3;
+   if(plane == 0){
+     padid=paddle;
+   }
+   if(plane == 1){
+     padid=pads11+paddle;
+   }
+   if(plane == 2){
+     padid=pads11+pads12+paddle;
+   }
+   if(plane == 3){
+     padid=pads11+pads12+pads21+paddle;
+   }
+   if(plane == 4){
+     padid=pads11+pads12+pads21+pads22+paddle;
+   }
+   if(plane == 5){
+     padid=pads11+pads12+pads21+pads22+pads31+paddle;
+   }
+   return padid;
+ }
+ // gf Apr 07
+ /**
+  * Method to get the plane ID and the paddle id in the plane if the paddle index (0,...23) is given.
+  * This method is the
+  * "reverse" of method "GetPaddleid"
+  * @param pad      (0 - 23)
+  * @param plane    (0 - 5)
+  * @param paddle   (plane=0, paddle = 0,...5)
+  */
+ void ToFLevel2::GetPaddlePlane(Int_t pad, Int_t &plane, Int_t &paddle)
+ {
+   Int_t pads11=8;
+   Int_t pads12=6;
+   Int_t pads21=2;
+   Int_t pads22=2;
+   Int_t pads31=3;
+   // Int_t pads32=3;
+   if(pad<8){
+     plane=0;
+     paddle=pad;
+     return;
+   }
+   if(7<pad<14){
+     plane=1;
+     paddle=pad-pads11;
+     return;
+   }
+   if(13<pad<16){
+     plane=2;
+     paddle=pad-pads11-pads12;
+     return;
+   }
+   if(15<pad<18){
+     plane=3;
+     paddle=pad-pads11-pads12-pads21;
+     return;
+   }
+   if(17<pad<21){
+     plane=4;
+     paddle=pad-pads11-pads12-pads21-pads22;
+     return;
+   }
+   if(20<pad<24){
+     plane=5;
+     paddle=pad-pads11-pads12-pads21-pads22-pads31;
+     return;
+   }
+ }
+ Int_t ToFLevel2::GetNPaddle(Int_t plane){
+   Int_t npaddle=-1;
+   Int_t pads11=8;
+   Int_t pads12=6;
+   Int_t pads21=2;
+   Int_t pads22=2;
+   Int_t pads31=3;
+   Int_t pads32=3;
+   if(plane==0)npaddle=pads11;
+   if(plane==1)npaddle=pads12;
+   if(plane==2)npaddle=pads21;
+   if(plane==3)npaddle=pads22;
+   if(plane==4)npaddle=pads31;
+   if(plane==5)npaddle=pads32;
+   return npaddle;
+ }
+ /// wm feb 08
+ /**
+  * Method to calculate Beta from the 12 single measurements
+  * we check the individual weights for artificial TDC values, then calculate
+  * am mean beta for the first time. In a second step we loop again through
+  * the single measurements, checking for the residual from the mean
+  * The cut on the residual reject measurements > "x"-sigma. A chi2 value is
+  * calculated, furthermore a "quality" value by adding the weights which
+  * are finally used. If all measurements are taken, "quality" will be = 22.47.
+  * A chi2 cut around 3-4 and a quality-cut > 20 is needed for clean beta
+  * measurements like antiprotons etc.
+  * The Level2 output is derived in the fortran routines using: 10.,10.,20.
+  * @param notrack Track Number
+  * @param cut on residual: difference between single measurement and mean
+  * @param cut on "quality"
+  * @param cut on chi2
+  */
+ Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){
+ //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;
+   Float_t bxx = 100.;
+   //
+   ToFTrkVar *trk = GetToFTrkVar(notrack);
+   if(!trk) return 0; //ELENA
+   Float_t chi2,xhelp,beta_mean;
+   Float_t w_i[12],quality,sw,sxw,res,betachi,beta_mean_inv;
+   Float_t b[12],tdcfl;
+   Int_t  pmt_id,pmt_plane;
+   for (Int_t i=0; i<12; i++){
+     b[i] = trk->beta[i];
+                               }
+ //========================================================================
+ //---  Find out ToF layers with artificial TDC values & fill vector    ---
+ //========================================================================
+ Float_t  w_il[6];
+      for (Int_t jj=0; jj<6;jj++) {
+          w_il[jj] = 1000.;
+                                  }
+   for (Int_t i=0; i<trk->npmttdc; i++){
+     //
+     pmt_id = (trk->pmttdc).At(i);
+     pmt_plane = GetPlaneIndex(pmt_id);
+     tdcfl = (trk->tdcflag).At(i);
+     if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag
+                                      };
+ //========================================================================
+ //---  Set weights for the 12 measurements using information for top and bottom:
+ //---  if no measurements: weight = set to very high value=> not used
+ //---  top or bottom artificial: weight*sqrt(2)
+ //---  top and bottom artificial: weight*sqrt(2)*sqrt(2)
+ //========================================================================
+ Int_t itop[12] = {0,0,1,1,2,2,3,3,0,0,1,1};
+ Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2,3};
+      xhelp= 1E09;
+      for (Int_t jj=0; jj<12;jj++) {
+      if (jj<4)           xhelp = 0.11;    // S1-S3
+      if ((jj>3)&&(jj<8)) xhelp = 0.18;    // S2-S3
+      if (jj>7)           xhelp = 0.28;    // S1-S2
+      if ((w_il[itop[jj]] == 1000.) && (w_il[ibot[jj]] == 1000.)) xhelp = 1E09;
+      if ((w_il[itop[jj]] == 1) || (w_il[ibot[jj]] == 1.)) xhelp = xhelp*1.414 ;
+      if ((w_il[itop[jj]] == 1) && (w_il[ibot[jj]] == 1.)) xhelp = xhelp*2. ;
+      w_i[jj] = 1./xhelp;
+                                   }
+ //========================================================================
+ //--- Calculate mean beta for the first time -----------------------------
+ //--- We are using "1/beta" since its error is gaussian ------------------
+ //========================================================================
+       Int_t icount=0;
+       sw=0.;
+       sxw=0.;
+       beta_mean=100.;
+           for (Int_t jj=0; jj<12;jj++){
+         if ((fabs(1./b[jj])>0.1)&&(fabs(1./b[jj])<15.))
+          {
+             icount= icount+1;
+             sxw=sxw + (1./b[jj])*w_i[jj]*w_i[jj] ;
+             sw =sw + w_i[jj]*w_i[jj] ;
+          }
+          }
+       if (icount>0) beta_mean=1./(sxw/sw);
+       beta_mean_inv = 1./beta_mean;
+ //========================================================================
+ //--- Calculate beta for the second time, use residuals of the single
+ //--- measurements to get a chi2 value
+ //========================================================================
+       icount=0;
+       sw=0.;
+       sxw=0.;
+       betachi = 100.;
+       chi2 = 0.;
+       quality=0.;
+           for (Int_t jj=0; jj<12;jj++){
+        if ((fabs(1./b[jj])>0.1)&&(fabs(1./b[jj])<15.)&&(w_i[jj]>0.01)) {
+             res = beta_mean_inv - (1./b[jj]) ;
+             if (fabs(res*w_i[jj])<resmax)          {;
+             chi2 = chi2 + pow((res*w_i[jj]),2) ;
+             icount= icount+1;
+             sxw=sxw + (1./b[jj])*w_i[jj]*w_i[jj] ;
+             sw =sw + w_i[jj]*w_i[jj] ;
+                                                }
+                                                                         }
+                                       }
+       quality = sqrt(sw) ;
+       if (icount==0) chi2 = 1000.;
+       if (icount>0) chi2 = chi2/(icount) ;
+       if (icount>0) betachi=1./(sxw/sw);
+    bxx = 100.;
+    if ((chi2 < chi2cut)&&(quality>qualitycut)) bxx = betachi;
+   //
+   return(bxx);
+ };
+ ////////////////////////////////////////////////////
+ ////////////////////////////////////////////////////
  /**
   * Fills a struct cToFLevel2 with values from a ToFLevel2 object (to put data into a F77 common).
   */
-Line 389 
 void ToFLevel2::GetLevel2Struct(cToFLeve
+Line 1267 
 void ToFLevel2::GetLevel2Struct(cToFLeve
                l2->xtofpos[i][j]=((ToFTrkVar*)ToFTrk->At(j))->xtofpos[i];
                l2->ytofpos[i][j]=((ToFTrkVar*)ToFTrk->At(j))->ytofpos[i];
            }
+           for(Int_t i=0;i<6;i++){
+               l2->xtr_tof[i][j]=((ToFTrkVar*)ToFTrk->At(j))->xtr_tof[i];
+               l2->ytr_tof[i][j]=((ToFTrkVar*)ToFTrk->At(j))->ytr_tof[i];
+           }
        }
    } //ELENA

 Legend:



Removed from v.1.13
 


changed lines


 
Added in v.1.23
 Legend:



Removed from v.1.13
 


changed lines


 
Added in v.1.23
-Removed from v.1.13
+Added in v.1.23

	ViewVC Help
Powered by ViewVC 1.1.23