/[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.43 by mocchiut, Thu Feb 26 11:27:35 2015 UTC revision 1.45 by pam-fi, Tue Mar 24 11:09:54 2015 UTC
# Line 319  Float_t ToFLevel2::GetdEdx(Int_t notrack Line 319  Float_t ToFLevel2::GetdEdx(Int_t notrack
319   * The method was modified for the "ToF-standalone" part in february 2015   * The method was modified for the "ToF-standalone" part in february 2015
320   * The "adcfl" option is not very useful (an artificial dEdx is per   * The "adcfl" option is not very useful (an artificial dEdx is per
321   * 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
322   * @param trk Pointer to TofTrkVar class   * @param trk Pointer to TofTrkVar object
323   * @param plane Plane index (0,1,2,3,4,5)   * @param plane Plane index (0,1,2,3,4,5)
324   * @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; )
325   */   */
# Line 336  Float_t ToFLevel2::GetdEdx(ToFTrkVar *tr Line 336  Float_t ToFLevel2::GetdEdx(ToFTrkVar *tr
336    //    //
337    if ( trk->trkseqno == -1 ){    if ( trk->trkseqno == -1 ){
338            
339      ToFTrkVar *t_tof = trk;      //    ToFTrkVar *t_tof = trk;
340            
341      // Find the hitted paddle  (two good TDC values) using the tof_j_flag (from tofl2com.for)      // Find the hitted paddle  (two good TDC values) using the tof_j_flag (from tofl2com.for)
342            
# Line 349  Float_t ToFLevel2::GetdEdx(ToFTrkVar *tr Line 349  Float_t ToFLevel2::GetdEdx(ToFTrkVar *tr
349      // Get the dEdx of this paddle using "GetdEdxPaddle"      // Get the dEdx of this paddle using "GetdEdxPaddle"
350      if (Ipaddle>-1) {      if (Ipaddle>-1) {
351        Int_t pad = GetPaddleid(plane,Ipaddle);        Int_t pad = GetPaddleid(plane,Ipaddle);
352        GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);        GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
353        dedx = PadEdx;        dedx = PadEdx;
354      }      }
355            
# Line 362  Float_t ToFLevel2::GetdEdx(ToFTrkVar *tr Line 362  Float_t ToFLevel2::GetdEdx(ToFTrkVar *tr
362        for (Int_t ii=0; ii<GetNPaddle(plane); ii++){        for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
363          Int_t paddleid=ii;          Int_t paddleid=ii;
364          Int_t pad = GetPaddleid(plane,paddleid);          Int_t pad = GetPaddleid(plane,paddleid);
365          GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);          GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
366          dedx += PadEdx;          dedx += PadEdx;
367        }        }
368      }      }
# Line 372  Float_t ToFLevel2::GetdEdx(ToFTrkVar *tr Line 372  Float_t ToFLevel2::GetdEdx(ToFTrkVar *tr
372      for (Int_t ii=0; ii<GetNPaddle(plane); ii++){      for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
373        Int_t paddleid=ii;        Int_t paddleid=ii;
374        pad = GetPaddleid(plane,paddleid);        pad = GetPaddleid(plane,paddleid);
375        GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);        GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
376        dedx += PadEdx;        dedx += PadEdx;
377      }      }
378    }    }
# Line 453  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt Line 453  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt
453   */   */
454  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){
455    //    //
456    short tof[4][24] = {    static const short tof[4][24] = {
457      {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},
458      {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},
459      {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 489  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_ Line 489  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_
489   */   */
490  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){
491    //    //
492    short tof[4][24] = {    static const short tof[4][24] = {
493      {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},
494      {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},
495      {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 529  void ToFLevel2::GetPMTIndex(Int_t ind, I Line 529  void ToFLevel2::GetPMTIndex(Int_t ind, I
529   */   */
530  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){
531    
   /*  
     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];  
   //  
532    ToFTrkVar *trk = GetToFTrkVar(notrack);    ToFTrkVar *trk = GetToFTrkVar(notrack);
533    if(!trk) return; //ELENA    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];  
534        
535  };  };
536    
# Line 731  TString ToFLevel2::GetPMTName(Int_t ind, Line 649  TString ToFLevel2::GetPMTName(Int_t ind,
649        
650    TString pmtname = " ";    TString pmtname = " ";
651        
652    TString photoS[48] = {    static const TString photoS[48] = {
653      "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",
654      "S11_4B",      "S11_4B",
655      "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 1191  Int_t ToFLevel2::GetNPaddle(Int_t plane) Line 1109  Int_t ToFLevel2::GetNPaddle(Int_t plane)
1109   * @param cut on chi2   * @param cut on chi2
1110   */   */
1111    
 Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){  
1112    
1113  //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;  Float_t ToFTrkVar::CalcBeta( Float_t resmax, Float_t qualitycut, Float_t chi2cut){
1114    
1115    
1116    Float_t bxx = 100.;    Float_t bxx = 100.;
1117    //    //
1118    ToFTrkVar *trk = GetToFTrkVar(notrack);    ToFTrkVar *trk = this;
   if(!trk) return 0; //ELENA  
1119    
1120    
1121    Float_t chi2,xhelp,beta_mean;    Float_t chi2,xhelp,beta_mean;
# Line 1225  Float_t  w_il[6]; Line 1142  Float_t  w_il[6];
1142    for (Int_t i=0; i<trk->npmttdc; i++){    for (Int_t i=0; i<trk->npmttdc; i++){
1143      //      //
1144      pmt_id = (trk->pmttdc).At(i);      pmt_id = (trk->pmttdc).At(i);
1145      pmt_plane = GetPlaneIndex(pmt_id);      pmt_plane = ToFLevel2::GetPlaneIndex(pmt_id);
1146      tdcfl = (trk->tdcflag).At(i);      tdcfl = (trk->tdcflag).At(i);
1147      if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag      if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag
1148                                       };                                       };
# Line 1312  Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2, Line 1229  Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2,
1229    //    //
1230    return(bxx);    return(bxx);
1231  };  };
1232    ////////////////////////////////////////////////////
1233    ////////////////////////////////////////////////////
1234    /**
1235     * See ToFTrkVar::CalcBeta(Float_t,Float_t, Float_t).
1236     */
1237    Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){
1238    
1239    //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;
1240    
1241      ToFTrkVar *trk = GetToFTrkVar(notrack);
1242      if(!trk) return 0; //ELENA
1243    
1244      return trk->CalcBeta(resmax,qualitycut,chi2cut);
1245    
1246    };
1247    
1248    
1249  ////////////////////////////////////////////////////  ////////////////////////////////////////////////////
Legend:
Removed from v.1.43  
changed lines
  Added in v.1.45
  ViewVC Help
Powered by ViewVC 1.1.23