/[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.19 by mocchiut, Mon Mar 3 09:51:02 2008 UTC revision 1.47 by pam-fi, Tue May 19 13:46:46 2015 UTC
# Line 1  Line 1 
1  /**  /**
2   * \file ToFLevel2.cpp   * \file ToFLevel2.cpp
3   * \author Gianfranca DeRosa, Wolfgang Menn   * \author Gianfranca DeRosa, Wolfgang Menn
4     *
5     * WM dec 2008: Description of "GetdEdx" changed
6     * WM dec 2008: "GetdEdxPaddle" modified: Now includes saturation limit
7     *              PMTs higher than the saturation limit are not used for dEdx
8     * WM apr 2009: bug found by Nicola in method "GetPaddlePlane"
9   */   */
10    
 #include <TObject.h>  
11  #include <ToFLevel2.h>  #include <ToFLevel2.h>
 #include <iostream>  
12  using namespace std;  using namespace std;
13  ClassImp(ToFPMT);  ClassImp(ToFPMT);
14    ClassImp(ToFdEdx);
15    ClassImp(ToFGeom);
16  ClassImp(ToFTrkVar);  ClassImp(ToFTrkVar);
17  ClassImp(ToFLevel2);  ClassImp(ToFLevel2);
18    
# Line 16  ToFPMT::ToFPMT(){ Line 21  ToFPMT::ToFPMT(){
21    adc = 0.;    adc = 0.;
22    tdc_tw = 0.;    tdc_tw = 0.;
23    tdc = 0.;    tdc = 0.;
24      l0flag_adc = 0.;
25      l0flag_tdc = 0.;
26  }  }
27    
28  ToFPMT::ToFPMT(const ToFPMT &t){  ToFPMT::ToFPMT(const ToFPMT &t){
# Line 106  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 113  void ToFLevel2::Clear(Option_t *t){ Line 130  void ToFLevel2::Clear(Option_t *t){
130    if(PMT)PMT->Delete(); //ELENA    if(PMT)PMT->Delete(); //ELENA
131    memset(tof_j_flag, 0, 6*sizeof(Int_t));    memset(tof_j_flag, 0, 6*sizeof(Int_t));
132    unpackError = 0;    unpackError = 0;
133      unpackWarning = 0;
134    //    //
135  };  };
136    
# Line 129  void ToFLevel2::Delete(Option_t *t){ //E Line 147  void ToFLevel2::Delete(Option_t *t){ //E
147    //    //
148  }; //ELENA  }; //ELENA
149    
150    /**
151     * Retrieves the itrk-th tof track stored in the array
152     * @param itrk Array index (min 0, max ToFLevel2::ntrk())
153     *
154     */
155  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t itrk){  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t itrk){
156    //        //    
157    if(itrk >= ntrk()){    if(itrk >= ntrk()){
# Line 143  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t Line 166  ToFTrkVar *ToFLevel2::GetToFTrkVar(Int_t
166    return toftrack;    return toftrack;
167  }  }
168    
169    /**
170     * Retrieves the tof track matching the seqno-th tracker stored track.
171     * @param seqno Track sequential number
172     * (seqno = -1 for standalone tof track, seqno=0-TrkLevel2::ntrk() for tof tracks associated to a tracker track)
173     *
174     */
175    ToFTrkVar *ToFLevel2::GetToFStoredTrack(int seqno){
176    
177      if( ntrk()==0 ){
178        printf("ToFLevel2::GetToFStoredTrack(int) : requested tracker SeqNo %i but no ToFrimeter tracks are stored\n",seqno);
179        return NULL;
180      };
181      
182      ToFTrkVar *c = 0;
183      Int_t it_tof=0;
184        
185      do {
186        c = GetToFTrkVar(it_tof);
187        it_tof++;
188      } while( c && seqno != c->trkseqno && it_tof < ntrk());      
189      
190      if(!c || seqno != c->trkseqno){
191        c = 0;
192        if(seqno!=-1 ) printf("ToFLevel2::GetToFStoredTrack(int) : requested tracker SeqNo %i does not match ToFrimeter stored tracks\n",seqno);
193      };
194      return c;
195        
196    }
197    
198    
199  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit){  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit){
200    //        //    
201    if(ihit >= npmt()){    if(ihit >= npmt()){
# Line 206  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit) Line 259  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
259         false) return true;         false) return true;
260      else return false;      else return false;
261  };  };
262    
263  /**  /**
264   * Method to get the number of hit paddles on a ToF plane.   * Strict method to get the number of hit paddles on a ToF plane.
265     * The method uses "HitPaddle" which checks if there is a TDC signal
266     * from both PMTs.
267   * @param plane Plane ID (11, 12, 21, 22, 31, 32) or Plane index (0,1,2,3,4,5).   * @param plane Plane ID (11, 12, 21, 22, 31, 32) or Plane index (0,1,2,3,4,5).
268   */   */
269  Int_t ToFLevel2::GetNHitPaddles(Int_t plane){  Int_t ToFLevel2::GetNHitPaddles(Int_t plane){
# Line 216  Int_t ToFLevel2::GetNHitPaddles(Int_t pl Line 272  Int_t ToFLevel2::GetNHitPaddles(Int_t pl
272      return npad;      return npad;
273  };  };
274    
275    /**
276     * Optional method to get the number of hit paddles on a ToF plane.
277     * The method does NOT check if there is a signal from both PMTs, it only
278     * checks if there is some PMT signal in a paddle
279     * @param plane Plane ID (11, 12, 21, 22, 31, 32) or Plane index (0,1,2,3,4,5).
280     */
281    Int_t ToFLevel2::GetTrueNHitPaddles(Int_t plane){
282        Int_t npad=0;
283        TClonesArray* Pmt = this->PMT;
284        int paddle[24];
285        memset(paddle,0, 24*sizeof(int));
286        for(int i=0; i<Pmt->GetEntries(); i++) {  //loop per vedere quale TOF è colpito
287          ToFPMT* pmthit = (ToFPMT*)Pmt->At(i);
288          int pplane = -1;
289          int ppaddle = -1;
290          GetPMTPaddle(pmthit->pmt_id,pplane,ppaddle);
291          if ( pplane == plane ) paddle[ppaddle]++;
292        }
293        for(int i=0;i<24;i++) if ( paddle[i]>0 ) npad++;
294    
295        return npad;
296    };
297    
298    //new, wm Feb 15
299    //wm Nov 08
300  //gf Apr 07  //gf Apr 07
301  /**  /**
302   * Method to get the mean dEdx from a given ToF plane. This current version   * Method to get the mean dEdx from a ToF layer
303   * is just summing up all PMT signals, which will not give proper results,   * By definition there should be PMTs with dEdx values only in one paddle of a layer
304   *  and needs a revision.   * (the paddle hitted by the track), this method looks for the hitted paddle
305     * and gives the mean dEdx of that paddle as the output
306     * The method was modified for the "ToF-standalone" part in february 2015
307     * The "adcfl" option is not very useful (an artificial dEdx is per
308     * definition= 1 mip and not a real measurement), anyway left in the code
309   * @param notrack Track Number   * @param notrack Track Number
310   * @param plane Plane index (0,1,2,3,4,5)   * @param plane Plane index (0,1,2,3,4,5)
311   * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )   * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
312   */   */
313  Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){  Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane, Int_t adcfl){
314    //  ToFTrkVar *trk = GetToFTrkVar(notrack);
315        ToFTrkVar *trk = GetToFStoredTrack(notrack);//Elena 2015
316      return this->GetdEdx(trk, plane, adcfl);
317    }
318    
319    //new, wm Feb 15
320    //wm Nov 08
321    //gf Apr 07
322    /**
323     * Method to get the mean dEdx from a ToF layer
324     * By definition there should be PMTs with dEdx values only in one paddle of a layer
325     * (the paddle hitted by the track), this method looks for the hitted paddle
326     * and gives the mean dEdx of that paddle as the output
327     * The method was modified for the "ToF-standalone" part in february 2015
328     * The "adcfl" option is not very useful (an artificial dEdx is per
329     * definition= 1 mip and not a real measurement), anyway left in the code
330     * @param trk Pointer to TofTrkVar object
331     * @param plane Plane index (0,1,2,3,4,5)
332     * @param adcflag in the plane (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
333     */
334    Float_t ToFLevel2::GetdEdx(ToFTrkVar *trk, Int_t plane, Int_t adcfl){
335      
336    Float_t dedx = 0.;    Float_t dedx = 0.;
337    Float_t PadEdx =0.;    Float_t PadEdx =0.;
338    Int_t SatWarning;    Int_t SatWarning;
339    Int_t pad=-1;    Int_t pad=-1;
340    //    //
341    ToFTrkVar *trk = GetToFTrkVar(notrack);    if(!trk) cout << "ToFLevel2::GetdEdx(...) ---> NULL ToFTrkVar obj "<<endl;
342    if(!trk) return 0; //ELENA    if(!trk) return 0; //ELENA
343    //    //
344    for (Int_t ii=0; ii<GetNPaddle(plane); ii++){    // ToF standalone part
     Int_t paddleid=ii;  
     pad = GetPaddleid(plane,paddleid);  
     GetdEdxPaddle(notrack, pad, adcfl, PadEdx, SatWarning);  
     dedx += PadEdx;  
   };  
345    //    //
346      if ( trk->trkseqno == -1 ){
347        
348        //    ToFTrkVar *t_tof = trk;
349        
350        // Find the hitted paddle  (two good TDC values) using the tof_j_flag (from tofl2com.for)
351        
352        Int_t Ipaddle=-1;
353        // if tof_j_flag == 0: no paddle was hitted. Otherwise decode tof_j_flag to get the paddle
354        if (this->tof_j_flag[plane] > 0)  Ipaddle = (Int_t)log2(this->tof_j_flag[plane]) ;
355        
356        Ipaddle =  (Int_t)log2(this->tof_j_flag[plane]) ;
357        
358        // Get the dEdx of this paddle using "GetdEdxPaddle"
359        if (Ipaddle>-1) {
360          Int_t pad = GetPaddleid(plane,Ipaddle);
361          GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
362          dedx = PadEdx;
363        }
364        
365        // If there was no correct hitted paddle, but there was one (and only one) paddle with some
366        // PMT entries in the PMT-class (found with "GetTrueNHitPaddles", use the dEdx of this paddle
367        
368        if ((Ipaddle<0) && (GetTrueNHitPaddles(plane)==1)) {
369          // find the paddle by looping over the paddles in each layer
370          // since GetTrueNHitPaddles==1 this is OK
371          for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
372            Int_t paddleid=ii;
373            Int_t pad = GetPaddleid(plane,paddleid);
374            GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
375            dedx += PadEdx;
376          }
377        }
378      } else {
379        // track dependent dEdx: simple, there will be only one paddle hitted in    each layer
380        // so just loop over the paddles in each layer
381        for (Int_t ii=0; ii<GetNPaddle(plane); ii++){
382          Int_t paddleid=ii;
383          pad = GetPaddleid(plane,paddleid);
384          GetdEdxPaddle(trk, pad, adcfl, PadEdx, SatWarning);
385          dedx += PadEdx;
386        }
387      }
388      //
389    return(dedx);    return(dedx);
390  };  }
391    
392  /**  /**
393   * Method to fill the ADC_C 4x12 matrix with the dEdx values and the TDC 4x12 matrix   * Method to fill the ADC_C 4x12 matrix with the dEdx values and the TDC 4x12 matrix
# Line 319  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt Line 462  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt
462   */   */
463  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){
464    //    //
465    short tof[4][24] = {    static const short tof[4][24] = {
466      {4, 4,  4,  4,  1,  1, 2, 2,  3,  3, 3, 3,  3,  3, 1, 1,  1,  1, 2, 3,  3, 3, 3,  4},      {4, 4,  4,  4,  1,  1, 2, 2,  3,  3, 3, 3,  3,  3, 1, 1,  1,  1, 2, 3,  3, 3, 3,  4},
467      {1, 3,  5,  7, 10, 12, 2, 4,  2,  4, 6, 8, 10, 12, 1, 5,  3,  9, 7, 9, 11, 1, 5,  9},      {1, 3,  5,  7, 10, 12, 2, 4,  2,  4, 6, 8, 10, 12, 1, 5,  3,  9, 7, 9, 11, 1, 5,  9},
468      {2, 2,  2,  2,  1,  1, 1, 1,  4,  4, 4, 4,  4,  4, 2, 1,  2,  1, 2, 2,  2, 3, 3,  4},      {2, 2,  2,  2,  1,  1, 1, 1,  4,  4, 4, 4,  4,  4, 2, 1,  2,  1, 2, 2,  2, 3, 3,  4},
# Line 355  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_ Line 498  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_
498   */   */
499  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){
500    //    //
501    short tof[4][24] = {    static const short tof[4][24] = {
502      {4, 4,  4,  4,  1,  1, 2, 2,  3,  3, 3, 3,  3,  3, 1, 1,  1,  1, 2, 3,  3, 3, 3,  4},      {4, 4,  4,  4,  1,  1, 2, 2,  3,  3, 3, 3,  3,  3, 1, 1,  1,  1, 2, 3,  3, 3, 3,  4},
503      {1, 3,  5,  7, 10, 12, 2, 4,  2,  4, 6, 8, 10, 12, 1, 5,  3,  9, 7, 9, 11, 1, 5,  9},      {1, 3,  5,  7, 10, 12, 2, 4,  2,  4, 6, 8, 10, 12, 1, 5,  3,  9, 7, 9, 11, 1, 5,  9},
504      {2, 2,  2,  2,  1,  1, 1, 1,  4,  4, 4, 4,  4,  4, 2, 1,  2,  1, 2, 2,  2, 3, 3,  4},      {2, 2,  2,  2,  1,  1, 1, 1,  4,  4, 4, 4,  4,  4, 2, 1,  2,  1, 2, 2,  2, 3, 3,  4},
# Line 381  void ToFLevel2::GetPMTIndex(Int_t ind, I Line 524  void ToFLevel2::GetPMTIndex(Int_t ind, I
524    
525    
526    
527    //  wm Nov 08 revision - saturation values included
528  /// gf Apr 07  /// gf Apr 07
   
529  /**  /**
530   * Method to get the dEdx from a given ToF paddle.   * Method to get the dEdx from a given ToF paddle.
531     * If two PMTs are good, the mean dEdx of both PMTs is taken, otherwise
532     * just the dEdx of the "good" PMT. If both PMTs are above saturation => dEdx=1000
533   * @param notrack Track Number   * @param notrack Track Number
534   * @param Paddle index (0,1,...,23).   * @param Paddle index (0,1,...,23).
535   * @param adcflag in the paddle (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )   * @param adcflag in the paddle (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
# Line 394  void ToFLevel2::GetPMTIndex(Int_t ind, I Line 538  void ToFLevel2::GetPMTIndex(Int_t ind, I
538   */   */
539  void ToFLevel2::GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){  void ToFLevel2::GetdEdxPaddle(Int_t notrack, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){
540    
541    //  ToFTrkVar *trk = GetToFTrkVar(notrack);
542        ToFTrkVar *trk = GetToFStoredTrack(notrack); //Elena 2015
543      this->GetdEdxPaddle(trk, paddleid, adcfl, PadEdx, SatWarning);
544      
545    };
546    
547    //
548    //  wm Nov 08 revision - saturation values included
549    /// gf Apr 07
550    /**
551     * Method to get the dEdx from a given ToF paddle.
552     * If two PMTs are good, the mean dEdx of both PMTs is taken, otherwise
553     * just the dEdx of the "good" PMT. If both PMTs are above saturation => dEdx=1000
554     * @param notrack Track Number
555     * @param Paddle index (0,1,...,23).
556     * @param adcflag in the paddle (100<-> independent of the adcflag; !=0&&!=100 <-> at least one PMT with adcflag!=0; )
557     * @param PadEdx dEdx from a given ToF paddle
558     * @param SatWarning 1 if the PMT ios near saturation region (adcraw ~3000)
559     */
560    void ToFLevel2::GetdEdxPaddle(ToFTrkVar *trk, Int_t paddleid, Int_t adcfl, Float_t &PadEdx, Int_t &SatWarning){
561    
562      /*
563        Float_t  PMTsat[48] = {
564        3162.14, 3165.48, 3153.85, 3085.73, 3089.65, 3107.64, 3097.52, 3078.37,
565        3130.05, 3087.07, 3112.22, 3102.92, 3080.58, 3092.55, 3087.94, 3125.03,
566        3094.09, 3143.16, 3125.51, 3181.27, 3092.09, 3124.98, 3069.3, 3095.53,
567        3097.11, 3133.53, 3114.73, 3113.01, 3091.19, 3097.99, 3033.84, 3134.98,
568        3081.37, 3111.04, 3066.77, 3108.17, 3133, 3111.06, 3052.52, 3140.66,
569        3106.33, 3094.85, 3150.85, 3118.8, 3096.24, 3118.47,3111.36, 3117.11 } ;
570      */
571    
572      // new values from Napoli dec 2008
573      Float_t  PMTsat[48] = {
574        3176.35,3178.19,3167.38,3099.73,3117.00,3126.29,3111.44,3092.27,
575        3146.48,3094.41,3132.13,3115.37,3099.32,3110.97,3111.80,3143.14,
576        3106.72,3153.44,3136.00,3188.96,3104.73,3140.45,3073.18,3106.62,
577        3112.48,3146.92,3127.24,3136.52,3109.59,3112.89,3045.15,3147.26,
578        3095.92,3121.05,3083.25,3123.62,3150.92,3125.30,3067.60,3160.18,
579        3119.36,3108.92,3164.77,3133.64,3111.47,3131.98,3128.87,3135.56 };
580    
581      for (Int_t i=0; i<48;i++) PMTsat[i] = PMTsat[i] - 5.;  // safety margin
582    
583    
584    PadEdx = 0.;    PadEdx = 0.;
585    SatWarning = 1000;    //  SatWarning = 1000;
586      SatWarning = 0;   // 0=good, increase for each bad PMT
587    
588    Float_t dEdx[48] = {0};    Float_t dEdx[48] = {0};
589    Int_t pmt_id = -1;    Int_t pmt_id = -1;
590    Float_t adcraw[48];    Float_t adcraw[48];
591    //    //
592    ToFTrkVar *trk = GetToFTrkVar(notrack);    if(!trk)cout << "ToFLevel2::GetdEdxPaddle(...) ---> NULL ToFTrkVar obj "<<endl;
593    if(!trk) return; //ELENA    if(!trk) return; //ELENA
594    //    //
595    
# Line 427  void ToFLevel2::GetdEdxPaddle(Int_t notr Line 615  void ToFLevel2::GetdEdxPaddle(Int_t notr
615        adcraw[pmtright] = pmt->adc;        adcraw[pmtright] = pmt->adc;
616      }      }
617    }    }
618    
619        
620    for (Int_t i=0; i<trk->npmtadc; i++){    for (Int_t i=0; i<trk->npmtadc; i++){
621    
# Line 439  void ToFLevel2::GetdEdxPaddle(Int_t notr Line 628  void ToFLevel2::GetdEdxPaddle(Int_t notr
628      }      }
629    }    }
630    
631    if( adcraw[pmtleft] >3000 || adcraw[pmtright] >3000)SatWarning=1;  
632          //  if( adcraw[pmtleft] >3000 || adcraw[pmtright] >3000)SatWarning=1;  //old version
633    if(dEdx[pmtleft]!=0 && dEdx[pmtright]!=0){  
634      PadEdx = (dEdx[pmtleft]+dEdx[pmtright])*0.5;    // Increase SatWarning Counter for each PMT>Sat
635    }    if( adcraw[pmtleft] > PMTsat[pmtleft])SatWarning++;  
636    if(dEdx[pmtleft]==0 && dEdx[pmtright]!=0){    if( adcraw[pmtright] > PMTsat[pmtright])SatWarning++;
637      PadEdx = dEdx[pmtright];  
638    }    // if ADC  > sat set dEdx=1000
639    if(dEdx[pmtleft]!=0 && dEdx[pmtright]==0){    if( adcraw[pmtleft] > PMTsat[pmtleft]) dEdx[pmtleft] = 1000.;
640      PadEdx = dEdx[pmtleft];    if( adcraw[pmtright] > PMTsat[pmtright]) dEdx[pmtright] = 1000. ;
641    }  
642      // if two PMT are good, take mean dEdx, otherwise only the good dEdx
643      if(dEdx[pmtleft]<1000 && dEdx[pmtright]<1000) PadEdx = (dEdx[pmtleft]+dEdx[pmtright])*0.5;
644      if(dEdx[pmtleft]==1000 && dEdx[pmtright]<1000) PadEdx = dEdx[pmtright];  
645      if(dEdx[pmtleft]<1000 && dEdx[pmtright]==1000) PadEdx = dEdx[pmtleft];
646        
   return;  
647  };  };
 //  
   
648    
649  // gf Apr 07  // gf Apr 07
650    
# Line 470  TString ToFLevel2::GetPMTName(Int_t ind, Line 660  TString ToFLevel2::GetPMTName(Int_t ind,
660        
661    TString pmtname = " ";    TString pmtname = " ";
662        
663    TString photoS[48] = {    static const TString photoS[48] = {
664      "S11_1A", "S11_1B", "S11_2A", "S11_2B", "S11_3A", "S11_3B", "S11_4A",      "S11_1A", "S11_1B", "S11_2A", "S11_2B", "S11_3A", "S11_3B", "S11_4A",
665      "S11_4B",      "S11_4B",
666      "S11_5A", "S11_5B", "S11_6A", "S11_6B", "S11_7A", "S11_7B", "S11_8A",      "S11_5A", "S11_5B", "S11_6A", "S11_6B", "S11_7A", "S11_7B", "S11_8A",
# Line 508  TString ToFLevel2::GetPMTName(Int_t ind) Line 698  TString ToFLevel2::GetPMTName(Int_t ind)
698        
699  };  };
700    
701    // wm jun 08
 // gf Apr 07  
702  Int_t ToFLevel2::GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane){  Int_t ToFLevel2::GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane){
703    return GetPaddleIdOfTrack(xtr ,ytr ,plane, 0.4);
704    }
705    
706    // gf Apr 07
707    Int_t ToFLevel2::GetPaddleIdOfTrack(Float_t xtr, Float_t ytr, Int_t plane, Float_t margin){
708      
709    Double_t xt,yt,xl,xh,yl,yh;    Double_t xt,yt,xl,xh,yl,yh;
710        
711    Float_t tof11_x[8] = {-17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85};    Float_t tof11_x[8] = {-17.85,-12.75,-7.65,-2.55,2.55,7.65,12.75,17.85};
# Line 541  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa Line 735  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa
735      yh =  33.0/2. ;      yh =  33.0/2. ;
736      if ((yt>yl)&&(yt<yh)) {      if ((yt>yl)&&(yt<yh)) {
737        for (Int_t i1=0; i1<8;i1++){        for (Int_t i1=0; i1<8;i1++){
738          xl = tof11_x[i1] - (5.1-0.4)/2. ;          xl = tof11_x[i1] - (5.1-margin)/2. ;
739          xh = tof11_x[i1] + (5.1-0.4)/2. ;          xh = tof11_x[i1] + (5.1-margin)/2. ;
740          if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;          if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;
741        }        }
742      }      }
# Line 559  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa Line 753  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa
753            
754      if ((xt>xl)&&(xt<xh)) {      if ((xt>xl)&&(xt<xh)) {
755        for (Int_t i1=0; i1<6;i1++){        for (Int_t i1=0; i1<6;i1++){
756          yl = tof12_y[i1] - (5.5-0.4)/2. ;          yl = tof12_y[i1] - (5.5-margin)/2. ;
757          yh = tof12_y[i1] + (5.5-0.4)/2. ;          yh = tof12_y[i1] + (5.5-margin)/2. ;
758          if ((yt>yl)&&(yt<yh))  paddleidoftrack=i1;          if ((yt>yl)&&(yt<yh))  paddleidoftrack=i1;
759        }        }
760      }      }
# Line 577  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa Line 771  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa
771            
772      if ((xt>xl)&&(xt<xh)) {      if ((xt>xl)&&(xt<xh)) {
773        for (Int_t i1=0; i1<2;i1++){        for (Int_t i1=0; i1<2;i1++){
774          yl = tof21_y[i1] - (7.5-0.4)/2. ;          yl = tof21_y[i1] - (7.5-margin)/2. ;
775          yh = tof21_y[i1] + (7.5-0.4)/2. ;          yh = tof21_y[i1] + (7.5-margin)/2. ;
776          if ((yt>yl)&&(yt<yh))  paddleidoftrack=i1;          if ((yt>yl)&&(yt<yh))  paddleidoftrack=i1;
777        }        }
778      }      }
# Line 594  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa Line 788  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa
788            
789      if ((yt>yl)&&(yt<yh)) {      if ((yt>yl)&&(yt<yh)) {
790        for (Int_t i1=0; i1<2;i1++){        for (Int_t i1=0; i1<2;i1++){
791          xl = tof22_x[i1] - (9.0-0.4)/2. ;          xl = tof22_x[i1] - (9.0-margin)/2. ;
792          xh = tof22_x[i1] + (9.0-0.4)/2. ;          xh = tof22_x[i1] + (9.0-margin)/2. ;
793          if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;          if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;
794        }        }
795      }      }
# Line 611  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa Line 805  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa
805            
806      if ((yt>yl)&&(yt<yh)) {      if ((yt>yl)&&(yt<yh)) {
807        for (Int_t i1=0; i1<3;i1++){        for (Int_t i1=0; i1<3;i1++){
808          xl = tof31_x[i1] - (6.0-0.4)/2. ;          xl = tof31_x[i1] - (6.0-margin)/2. ;
809          xh = tof31_x[i1] + (6.0-0.4)/2. ;          xh = tof31_x[i1] + (6.0-margin)/2. ;
810          if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;          if ((xt>xl)&&(xt<xh))  paddleidoftrack=i1;
811        }        }
812      }      }
# Line 628  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa Line 822  Int_t ToFLevel2::GetPaddleIdOfTrack(Floa
822            
823      if ((xt>xl)&&(xt<xh)) {      if ((xt>xl)&&(xt<xh)) {
824        for (Int_t i1=0; i1<3;i1++){        for (Int_t i1=0; i1<3;i1++){
825          yl = tof32_y[i1] - (5.0-0.4)/2. ;          yl = tof32_y[i1] - (5.0-margin)/2. ;
826          yh = tof32_y[i1] + (5.0-0.4)/2. ;          yh = tof32_y[i1] + (5.0-margin)/2. ;
827          if ((yt>yl)&&(yt<yh)) paddleidoftrack=i1;          if ((yt>yl)&&(yt<yh)) paddleidoftrack=i1;
828        }        }
829      }      }
# Line 696  void ToFLevel2::GetPMTPaddle(Int_t pmt_i Line 890  void ToFLevel2::GetPMTPaddle(Int_t pmt_i
890  // gf Apr 07  // gf Apr 07
891    
892  void ToFLevel2::GetPaddlePMT(Int_t paddle, Int_t &pmtleft, Int_t &pmtright){  void ToFLevel2::GetPaddlePMT(Int_t paddle, Int_t &pmtleft, Int_t &pmtright){
893      pmtleft=paddle*2;
894    if(paddle==0){    pmtright= pmtleft+1;  
     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;  
   }  
     
895    return;    return;
896  }  }
897    
# Line 930  void ToFLevel2::GetPaddleGeometry(Int_t Line 1005  void ToFLevel2::GetPaddleGeometry(Int_t
1005   */   */
1006  Int_t ToFLevel2::GetPaddleid(Int_t plane, Int_t paddle)  Int_t ToFLevel2::GetPaddleid(Int_t plane, Int_t paddle)
1007  {  {
   
1008    Int_t padid=-1;    Int_t padid=-1;
1009    Int_t pads11=8;    Int_t pads[6]={8,6,2,2,3,3};
   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;  
   }  
1010    
1011    if(plane == 5){    int somma=0;
1012      padid=pads11+pads12+pads21+pads22+pads31+paddle;    int np=plane;
1013      for(Int_t j=0; j<np; j++){
1014        somma+=pads[j];
1015    }    }
1016      padid=paddle+somma;
1017    return padid;    return padid;
1018    
1019  }  }
# Line 994  void ToFLevel2::GetPaddlePlane(Int_t pad Line 1044  void ToFLevel2::GetPaddlePlane(Int_t pad
1044      return;      return;
1045    }    }
1046    
1047    if(7<pad<14){    if((7<pad)&&(pad<14)){
1048      plane=1;      plane=1;
1049      paddle=pad-pads11;      paddle=pad-pads11;
1050      return;      return;
1051    }    }
1052        
1053    if(13<pad<16){    if((13<pad)&&(pad<16)){
1054      plane=2;      plane=2;
1055      paddle=pad-pads11-pads12;      paddle=pad-pads11-pads12;
1056      return;      return;
1057    }    }
1058    
1059    if(15<pad<18){    if((15<pad)&&(pad<18)){
1060      plane=3;      plane=3;
1061      paddle=pad-pads11-pads12-pads21;      paddle=pad-pads11-pads12-pads21;
1062      return;      return;
1063    }    }
1064    
1065    if(17<pad<21){    if((17<pad)&&(pad<21)){
1066      plane=4;      plane=4;
1067      paddle=pad-pads11-pads12-pads21-pads22;      paddle=pad-pads11-pads12-pads21-pads22;
1068      return;      return;
1069    }    }
1070    
1071    if(20<pad<24){    if((20<pad)&&(pad<24)){
1072      plane=5;      plane=5;
1073      paddle=pad-pads11-pads12-pads21-pads22-pads31;      paddle=pad-pads11-pads12-pads21-pads22-pads31;
1074      return;      return;
# Line 1070  Int_t ToFLevel2::GetNPaddle(Int_t plane) Line 1120  Int_t ToFLevel2::GetNPaddle(Int_t plane)
1120   * @param cut on chi2   * @param cut on chi2
1121   */   */
1122    
 Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){  
1123    
1124  //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;  Float_t ToFTrkVar::CalcBeta( Float_t resmax, Float_t qualitycut, Float_t chi2cut){
1125    
1126    
1127    Float_t bxx = 100.;    Float_t bxx = 100.;
1128    //    //
1129    ToFTrkVar *trk = GetToFTrkVar(notrack);    ToFTrkVar *trk = this;
   if(!trk) return 0; //ELENA  
1130    
1131    
1132    Float_t chi2,xhelp,beta_mean;    Float_t chi2,xhelp,beta_mean;
# Line 1104  Float_t  w_il[6]; Line 1153  Float_t  w_il[6];
1153    for (Int_t i=0; i<trk->npmttdc; i++){    for (Int_t i=0; i<trk->npmttdc; i++){
1154      //      //
1155      pmt_id = (trk->pmttdc).At(i);      pmt_id = (trk->pmttdc).At(i);
1156      pmt_plane = GetPlaneIndex(pmt_id);      pmt_plane = ToFLevel2::GetPlaneIndex(pmt_id);
1157      tdcfl = (trk->tdcflag).At(i);      tdcfl = (trk->tdcflag).At(i);
1158      if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag      if (w_il[pmt_plane] != 1.) w_il[pmt_plane] = tdcfl; //tdcflag
1159                                       };                                       };
# Line 1191  Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2, Line 1240  Int_t ibot[12] = {4,5,4,5,4,5,4,5,2,3,2,
1240    //    //
1241    return(bxx);    return(bxx);
1242  };  };
1243    ////////////////////////////////////////////////////
1244    ////////////////////////////////////////////////////
1245    /**
1246     * See ToFTrkVar::CalcBeta(Float_t,Float_t, Float_t).
1247     */
1248    Float_t ToFLevel2::CalcBeta(Int_t notrack, Float_t resmax, Float_t qualitycut, Float_t chi2cut){
1249    
1250    //  cout<<" in CalcBeta "<<resmax<<" "<<chi2cut<<" "<<qualitycut<<endl;
1251    
1252    //  ToFTrkVar *trk = GetToFTrkVar(notrack); //wrong!
1253        ToFTrkVar *trk = GetToFStoredTrack(notrack);//Elena Apr 2015
1254      if(!trk) return 0; //ELENA
1255    
1256      return trk->CalcBeta(resmax,qualitycut,chi2cut);
1257    
1258    };
1259    
1260    
1261  ////////////////////////////////////////////////////  ////////////////////////////////////////////////////
# Line 1243  void ToFLevel2::GetLevel2Struct(cToFLeve Line 1308  void ToFLevel2::GetLevel2Struct(cToFLeve
1308        }        }
1309    } //ELENA    } //ELENA
1310  }  }
1311    
1312    
1313    //
1314    // Reprocessing tool // Emiliano 08/04/07
1315    //
1316    Int_t ToFLevel2::Process(TrkLevel2 *trk, TrigLevel2 *trg, GL_RUN *run, OrbitalInfo *orb, Bool_t force){
1317      //
1318      // Copiare qui qualcosa di simile a calonuclei per evitare di riprocessare sempre tutto
1319      //
1320      printf("\n\n\n ERROR: NOT IMPLEMENTED ANYMORE, write Emiliano if you need this method (Emiliano.Mocchiutti@ts.infn.it) \n\n\n");
1321      return(-1);
1322      //   //
1323      //   // structures to communicate with F77
1324      //   //
1325      //   extern struct ToFInput  tofinput_;
1326    //   extern struct ToFOutput tofoutput_;
1327    //   //
1328    //   // DB connection
1329    //   //
1330    //   TString host;
1331    //   TString user;
1332    //   TString psw;
1333    //   const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");
1334    //   const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");
1335    //   const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");
1336    //   if ( !pamdbhost ) pamdbhost = "";
1337    //   if ( !pamdbuser ) pamdbuser = "";
1338    //   if ( !pamdbpsw ) pamdbpsw = "";
1339    //   if ( strcmp(pamdbhost,"") ) host = pamdbhost;
1340    //   if ( strcmp(pamdbuser,"") ) user = pamdbuser;
1341    //   if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;
1342    //   //
1343    //   //
1344    //   TSQLServer *dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());
1345    //   if ( !dbc->IsConnected() ) return 1;
1346    //   stringstream myquery;
1347    //   myquery.str("");
1348    //   myquery << "SET time_zone='+0:00';";
1349    //   dbc->Query(myquery.str().c_str());
1350    //   delete dbc->Query("SET sql_mode = 'NO_UNSIGNED_SUBTRACTION';");
1351    //   GL_PARAM *glparam = new GL_PARAM();
1352    //   glparam->Query_GL_PARAM(1,1,dbc); // parameters stored in DB in GL_PRAM table
1353    //   trk->LoadField(glparam->PATH+glparam->NAME);
1354    //   //
1355    //   Bool_t defcal = true;
1356    //   Int_t error=glparam->Query_GL_PARAM(run->RUNHEADER_TIME,201,dbc); // parameters stored in DB in GL_PRAM table
1357    //   if ( error<0 ) {
1358    //     return(1);
1359    //   };
1360    //   printf(" Reading ToF parameter file: %s \n",(glparam->PATH+glparam->NAME).Data());
1361    //   if ( (UInt_t)glparam->TO_TIME != (UInt_t)4294967295UL ) defcal = false;
1362    //   //
1363    //   Int_t nlen = (Int_t)(glparam->PATH+glparam->NAME).Length();
1364    //   rdtofcal((char *)(glparam->PATH+glparam->NAME).Data(),&nlen);
1365    //   //
1366    //   Int_t adc[4][12];
1367    //   Int_t tdc[4][12];
1368    //   Float_t tdcc[4][12];
1369    //   //
1370    //   // process tof data
1371    //   //
1372    //   for (Int_t hh=0; hh<12;hh++){
1373    //     for (Int_t kk=0; kk<4;kk++){
1374    //            adc[kk][hh] = 4095;
1375    //            tdc[kk][hh] = 4095;
1376    //            tdcc[kk][hh] = 4095.;
1377    //            tofinput_.adc[hh][kk] = 4095;
1378    //            tofinput_.tdc[hh][kk] = 4095;
1379    //     };
1380    //   };
1381    //   Int_t ntrkentry = 0;
1382    //   Int_t npmtentry = 0;
1383    //   Int_t gg = 0;
1384    //   Int_t hh = 0;
1385    //   Int_t adcf[48];
1386    //   memset(adcf, 0, 48*sizeof(Int_t));
1387    //   Int_t tdcf[48];
1388    //   memset(tdcf, 0, 48*sizeof(Int_t));
1389    //   for (Int_t pm=0; pm < this->ntrk() ; pm++){
1390    //      ToFTrkVar *ttf = this->GetToFTrkVar(pm);
1391    //      for ( Int_t nc=0; nc < ttf->npmttdc; nc++){
1392    //             if ( (ttf->tdcflag).At(nc) != 0 ) tdcf[(ttf->pmttdc).At(nc)] = 1;
1393    //      };
1394    //      for ( Int_t nc=0; nc < ttf->npmtadc; nc++){
1395    //             if ( (ttf->adcflag).At(nc) != 0 ) adcf[(ttf->pmtadc).At(nc)] = 1;
1396    //      };
1397    //   };
1398    //   //
1399    //   for (Int_t pm=0; pm < this->npmt() ; pm++){
1400    //      ToFPMT *pmt = this->GetToFPMT(pm);
1401    //      this->GetPMTIndex(pmt->pmt_id, gg, hh);
1402    //      if ( adcf[pmt->pmt_id] == 0 ){
1403    //              tofinput_.adc[gg][hh] = (int)pmt->adc;
1404    //              adc[hh][gg] = (int)pmt->adc;
1405    //      };
1406    //      if ( tdcf[pmt->pmt_id] == 0 ){
1407    //              tofinput_.tdc[gg][hh] = (int)pmt->tdc;
1408    //              tdc[hh][gg] = (int)pmt->tdc;
1409    //      };
1410    //      tdcc[hh][gg] = (float)pmt->tdc_tw;
1411    //      // Int_t pppid = this->GetPMTid(hh,gg);
1412    //      //      printf(" pm %i pmt_id %i pppid %i hh %i gg %i tdcc %f tdc %f adc %f \n",pm,pmt->pmt_id,pppid,hh,gg,pmt->tdc_tw,pmt->tdc,pmt->adc);
1413    //   };
1414    //   //
1415    //   Int_t unpackError = this->unpackError;
1416    //   //
1417    //   for (Int_t hh=0; hh<5;hh++){
1418    //      tofinput_.patterntrig[hh]=trg->patterntrig[hh];
1419    //   };
1420    //   //
1421    //   this->Clear();
1422    //   //
1423    //       Int_t pmt_id = 0;
1424    //       ToFPMT *t_pmt = new ToFPMT();
1425    //       if(!(this->PMT)) this->PMT = new TClonesArray("ToFPMT",12); //ELENA
1426    //       TClonesArray &tpmt = *this->PMT;
1427    //       ToFTrkVar *t_tof = new ToFTrkVar();
1428    //       if(!(this->ToFTrk)) this->ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA
1429    //       TClonesArray &t = *this->ToFTrk;
1430    //       //
1431    //       //
1432    //       // Here we have calibrated data, ready to be passed to the FORTRAN routine which will extract common and track-related  variables.
1433    //       //
1434    //       npmtentry = 0;
1435    //       //
1436    //       ntrkentry = 0;
1437    //       //
1438    //       // Calculate tracks informations from ToF alone
1439    //       //
1440    //       tofl2com();
1441    //       //
1442    //       memcpy(this->tof_j_flag,tofoutput_.tof_j_flag,6*sizeof(Int_t));
1443    //       //
1444    //       t_tof->trkseqno = -1;
1445    //       //
1446    //       // and now we must copy from the output structure to the level2 class:
1447    //       //
1448    //       t_tof->npmttdc = 0;
1449    //       //
1450    //       for (Int_t hh=0; hh<12;hh++){
1451    //         for (Int_t kk=0; kk<4;kk++){
1452    //           if ( tofoutput_.tofmask[hh][kk] != 0 ){
1453    //             pmt_id = this->GetPMTid(kk,hh);
1454    //             t_tof->pmttdc.AddAt(pmt_id,t_tof->npmttdc);
1455    //             t_tof->tdcflag.AddAt(tofoutput_.tdcflagtof[hh][kk],t_tof->npmttdc); // gf: Jan 09/07
1456    //             t_tof->npmttdc++;
1457    //           };
1458    //         };
1459    //       };
1460    //       for (Int_t kk=0; kk<13;kk++){
1461    //         t_tof->beta[kk] = tofoutput_.betatof_a[kk];
1462    //       }
1463    //       //
1464    //       t_tof->npmtadc = 0;
1465    //       for (Int_t hh=0; hh<12;hh++){
1466    //         for (Int_t kk=0; kk<4;kk++){
1467    //           if ( tofoutput_.adctof_c[hh][kk] < 1000 ){
1468    //             t_tof->dedx.AddAt(tofoutput_.adctof_c[hh][kk],t_tof->npmtadc);
1469    //             pmt_id = this->GetPMTid(kk,hh);
1470    //             t_tof->pmtadc.AddAt(pmt_id,t_tof->npmtadc);
1471    //             t_tof->adcflag.AddAt(tofoutput_.adcflagtof[hh][kk],t_tof->npmtadc); // gf: Jan 09/07
1472    //             t_tof->npmtadc++;
1473    //           };
1474    //         };
1475    //       };
1476    //       //
1477    //       memcpy(t_tof->xtofpos,tofoutput_.xtofpos,sizeof(t_tof->xtofpos));
1478    //       memcpy(t_tof->ytofpos,tofoutput_.ytofpos,sizeof(t_tof->ytofpos));
1479    //       memcpy(t_tof->xtr_tof,tofoutput_.xtr_tof,sizeof(t_tof->xtr_tof));
1480    //       memcpy(t_tof->ytr_tof,tofoutput_.ytr_tof,sizeof(t_tof->ytr_tof));
1481    //       //
1482    //       new(t[ntrkentry]) ToFTrkVar(*t_tof);
1483    //       ntrkentry++;
1484    //       t_tof->Clear();
1485    //       //
1486    //       //
1487    //       //
1488    //       t_pmt->Clear();
1489    //       //
1490    //       for (Int_t hh=0; hh<12;hh++){
1491    //         for (Int_t kk=0; kk<4;kk++){
1492    //          // new WM
1493    //           if ( tofoutput_.tdc_c[hh][kk] < 4095 || adc[kk][hh] < 4095  || tdc[kk][hh] < 4095 ){
1494    // //          if ( tdcc[kk][hh] < 4095. || adc[kk][hh] < 4095  || tdc[kk][hh] < 4095 ){
1495    //             //
1496    //             t_pmt->pmt_id = this->GetPMTid(kk,hh);
1497    //             t_pmt->tdc_tw = tofoutput_.tdc_c[hh][kk];
1498    //             t_pmt->adc = (Float_t)adc[kk][hh];
1499    //             t_pmt->tdc = (Float_t)tdc[kk][hh];
1500    //             //
1501    //             new(tpmt[npmtentry]) ToFPMT(*t_pmt);
1502    //             npmtentry++;
1503    //             t_pmt->Clear();
1504    //           };
1505    //         };
1506    //       };
1507    //       //
1508    //       // Calculate track-related variables
1509    //       //
1510    //       if ( trk->ntrk() > 0 ){
1511    //         //
1512    //         // We have at least one track
1513    //         //
1514    //         //
1515    //         // Run over tracks
1516    //         //
1517    //         for(Int_t nt=0; nt < trk->ntrk(); nt++){
1518    //           //
1519    //           TrkTrack *ptt = trk->GetStoredTrack(nt);
1520    //           //
1521    //           // Copy the alpha vector in the input structure
1522    //           //
1523    //           for (Int_t e = 0; e < 5 ; e++){
1524    //             tofinput_.al_pp[e] = ptt->al[e];
1525    //           };
1526    //           //
1527    //           // Get tracker related variables for this track
1528    //           //
1529    //           toftrk();
1530    //           //
1531    //           // Copy values in the class from the structure (we need to use a temporary class to store variables).
1532    //           //
1533    //           t_tof->npmttdc = 0;
1534    //           for (Int_t hh=0; hh<12;hh++){
1535    //             for (Int_t kk=0; kk<4;kk++){
1536    //               if ( tofoutput_.tofmask[hh][kk] != 0 ){
1537    //                 pmt_id = this->GetPMTid(kk,hh);
1538    //                 t_tof->pmttdc.AddAt(pmt_id,t_tof->npmttdc);
1539    //                 t_tof->tdcflag.AddAt(tofoutput_.tdcflag[hh][kk],t_tof->npmttdc); // gf: Jan 09/07
1540    //                 t_tof->npmttdc++;
1541    //               };
1542    //             };
1543    //           };
1544    //           for (Int_t kk=0; kk<13;kk++){
1545    //             t_tof->beta[kk] = tofoutput_.beta_a[kk];
1546    //           };
1547    //           //
1548    //           t_tof->npmtadc = 0;
1549    //           for (Int_t hh=0; hh<12;hh++){
1550    //             for (Int_t kk=0; kk<4;kk++){
1551    //               if ( tofoutput_.adc_c[hh][kk] < 1000 ){
1552    //                 t_tof->dedx.AddAt(tofoutput_.adc_c[hh][kk],t_tof->npmtadc);
1553    //                 pmt_id = this->GetPMTid(kk,hh);
1554    //                 t_tof->pmtadc.AddAt(pmt_id,t_tof->npmtadc);
1555    //                 t_tof->adcflag.AddAt(tofoutput_.adcflag[hh][kk],t_tof->npmtadc); // gf: Jan 09/07
1556    //                 t_tof->npmtadc++;
1557    //               };
1558    //             };
1559    //           };
1560    //           //
1561    //           memcpy(t_tof->xtofpos,tofoutput_.xtofpos,sizeof(t_tof->xtofpos));
1562    //           memcpy(t_tof->ytofpos,tofoutput_.ytofpos,sizeof(t_tof->ytofpos));
1563    //           memcpy(t_tof->xtr_tof,tofoutput_.xtr_tof,sizeof(t_tof->xtr_tof));
1564    //           memcpy(t_tof->ytr_tof,tofoutput_.ytr_tof,sizeof(t_tof->ytr_tof));
1565    //           //
1566    //           // Store the tracker track number in order to be sure to have shyncronized data during analysis
1567    //           //
1568    //           t_tof->trkseqno = nt;
1569    //           //
1570    //           // create a new object for this event with track-related variables
1571    //           //
1572    //           new(t[ntrkentry]) ToFTrkVar(*t_tof);
1573    //           ntrkentry++;
1574    //           t_tof->Clear();
1575    //           //
1576    //         }; // loop on all the tracks
1577    //       //
1578    //       this->unpackError = unpackError;
1579    //       if ( defcal ){
1580    //         this->default_calib = 1;
1581    //       } else {
1582    //         this->default_calib = 0;
1583    //       };
1584    //};
1585    //  return(0);
1586    }
1587    
1588    bool ToFLevel2::bit(int decimal, char pos){
1589      return( (decimal>>pos)%2 );
1590    }
1591    
1592    bool ToFLevel2::checkPMT(TString givenpmt){
1593      TClonesArray* Pmt = this->PMT;
1594      //  printf(" ou %s entries %i \n",givenpmt.Data(),Pmt->GetEntries());
1595      for(int i=0; i<Pmt->GetEntries(); i++) {  
1596        ToFPMT* pmthit = (ToFPMT*)Pmt->At(i);
1597        TString pmtname = this->GetPMTName(pmthit->pmt_id);
1598        //    printf(" name %s \n",pmtname.Data());
1599        if ( !strcmp(pmtname.Data(),givenpmt.Data()) )
1600          return true;
1601      }
1602      //  printf(" PMT %s missing \n",givenpmt.Data());
1603      return false;
1604    }
1605    
1606    bool ToFLevel2::checkPMTpatternPMThit(TrigLevel2 *trg, int &pmtpattern, int &pmtnosignal){
1607      UInt_t *patterntrig = trg->patterntrig;
1608      pmtpattern = 0;
1609      pmtnosignal = 0;
1610      bool good = true;
1611      //S3
1612      if ( this->bit(patterntrig[2],0) ){ pmtpattern++;  if ( !this->checkPMT("S31_1A")){ pmtnosignal++; good = false;}}
1613      if ( this->bit(patterntrig[2],1) ){ pmtpattern++;  if ( !this->checkPMT("S31_2A")){ pmtnosignal++; good = false;}}
1614      if ( this->bit(patterntrig[2],2) ){ pmtpattern++;  if ( !this->checkPMT("S31_3A")){ pmtnosignal++; good = false;}}
1615      if ( this->bit(patterntrig[2],3) ){ pmtpattern++;  if ( !this->checkPMT("S31_1B")){ pmtnosignal++; good = false;}}
1616      if ( this->bit(patterntrig[2],4) ){ pmtpattern++;  if ( !this->checkPMT("S31_2B")){ pmtnosignal++; good = false;}}
1617      if ( this->bit(patterntrig[2],5) ){ pmtpattern++;  if ( !this->checkPMT("S31_3B")){ pmtnosignal++; good = false;}}      
1618      if ( this->bit(patterntrig[2],6) ){ pmtpattern++;  if ( !this->checkPMT("S32_1A")){ pmtnosignal++; good = false;}}
1619      if ( this->bit(patterntrig[2],7) ){ pmtpattern++;  if ( !this->checkPMT("S32_2A")){ pmtnosignal++; good = false;}}
1620      if ( this->bit(patterntrig[2],8) ){ pmtpattern++;  if ( !this->checkPMT("S32_3A")){ pmtnosignal++; good = false;}}
1621      if ( this->bit(patterntrig[2],9) ){ pmtpattern++;  if ( !this->checkPMT("S32_1B")){ pmtnosignal++; good = false;}}
1622      if ( this->bit(patterntrig[2],10) ){ pmtpattern++;  if ( !this->checkPMT("S32_2B")){ pmtnosignal++; good = false;}}
1623      if ( this->bit(patterntrig[2],11) ){ pmtpattern++;  if ( !this->checkPMT("S32_3B")){ pmtnosignal++; good = false;}}      
1624      //S2
1625      if ( this->bit(patterntrig[3],0) ){ pmtpattern++;  if ( !this->checkPMT("S21_1A")){ pmtnosignal++; good = false;}}
1626      if ( this->bit(patterntrig[3],1) ){ pmtpattern++;  if ( !this->checkPMT("S21_2A")){ pmtnosignal++; good = false;}}
1627      if ( this->bit(patterntrig[3],2) ){ pmtpattern++;  if ( !this->checkPMT("S21_1B")){ pmtnosignal++; good = false;}}
1628      if ( this->bit(patterntrig[3],3) ){ pmtpattern++;  if ( !this->checkPMT("S21_2B")){ pmtnosignal++; good = false;}}      
1629      if ( this->bit(patterntrig[3],4) ){ pmtpattern++;  if ( !this->checkPMT("S22_1A")){ pmtnosignal++; good = false;}}
1630      if ( this->bit(patterntrig[3],5) ){ pmtpattern++;  if ( !this->checkPMT("S22_2A")){ pmtnosignal++; good = false;}}
1631      if ( this->bit(patterntrig[3],6) ){ pmtpattern++;  if ( !this->checkPMT("S22_1B")){ pmtnosignal++; good = false;}}
1632      if ( this->bit(patterntrig[3],7) ){ pmtpattern++;  if ( !this->checkPMT("S22_2B")){ pmtnosignal++; good = false;}}      
1633      //S12
1634      if ( this->bit(patterntrig[4],0) ){ pmtpattern++;  if ( !this->checkPMT("S12_1A")){ pmtnosignal++; good = false;}}
1635      if ( this->bit(patterntrig[4],1) ){ pmtpattern++;  if ( !this->checkPMT("S12_2A")){ pmtnosignal++; good = false;}}
1636      if ( this->bit(patterntrig[4],2) ){ pmtpattern++;  if ( !this->checkPMT("S12_3A")){ pmtnosignal++; good = false;}}
1637      if ( this->bit(patterntrig[4],3) ){ pmtpattern++;  if ( !this->checkPMT("S12_4A")){ pmtnosignal++; good = false;}}
1638      if ( this->bit(patterntrig[4],4) ){ pmtpattern++;  if ( !this->checkPMT("S12_5A")){ pmtnosignal++; good = false;}}
1639      if ( this->bit(patterntrig[4],5) ){ pmtpattern++;  if ( !this->checkPMT("S12_6A")){ pmtnosignal++; good = false;}}      
1640      if ( this->bit(patterntrig[4],6) ){ pmtpattern++;  if ( !this->checkPMT("S12_1A")){ pmtnosignal++; good = false;}}
1641      if ( this->bit(patterntrig[4],7) ){ pmtpattern++;  if ( !this->checkPMT("S12_2A")){ pmtnosignal++; good = false;}}
1642      if ( this->bit(patterntrig[4],8) ){ pmtpattern++;  if ( !this->checkPMT("S12_3A")){ pmtnosignal++; good = false;}}
1643      if ( this->bit(patterntrig[4],9) ){ pmtpattern++;  if ( !this->checkPMT("S12_4B")){ pmtnosignal++; good = false;}}
1644      if ( this->bit(patterntrig[4],10) ){ pmtpattern++; if ( !this->checkPMT("S12_5B")){ pmtnosignal++; good = false;}}
1645      if ( this->bit(patterntrig[4],11) ){ pmtpattern++; if ( !this->checkPMT("S12_6B")){ pmtnosignal++; good = false;}}      
1646      //S11
1647      if ( this->bit(patterntrig[5],0) ){ pmtpattern++;  if ( !this->checkPMT("S11_1A")){ pmtnosignal++; good = false;}}
1648      if ( this->bit(patterntrig[5],1) ){ pmtpattern++;  if ( !this->checkPMT("S11_2A")){ pmtnosignal++; good = false;}}
1649      if ( this->bit(patterntrig[5],2) ){ pmtpattern++;  if ( !this->checkPMT("S11_3A")){ pmtnosignal++; good = false;}}
1650      if ( this->bit(patterntrig[5],3) ){ pmtpattern++;  if ( !this->checkPMT("S11_4A")){ pmtnosignal++; good = false;}}
1651      if ( this->bit(patterntrig[5],4) ){ pmtpattern++;  if ( !this->checkPMT("S11_5A")){ pmtnosignal++; good = false;}}
1652      if ( this->bit(patterntrig[5],5) ){ pmtpattern++;  if ( !this->checkPMT("S11_6A")){ pmtnosignal++; good = false;}}
1653      if ( this->bit(patterntrig[5],6) ){ pmtpattern++;  if ( !this->checkPMT("S11_7A")){ pmtnosignal++; good = false;}}
1654      if ( this->bit(patterntrig[5],7) ){ pmtpattern++;  if ( !this->checkPMT("S11_8A")){ pmtnosignal++; good = false;}}      
1655      if ( this->bit(patterntrig[5],8) ){ pmtpattern++;  if ( !this->checkPMT("S11_1B")){ pmtnosignal++; good = false;}}
1656      if ( this->bit(patterntrig[5],9) ){ pmtpattern++;  if ( !this->checkPMT("S11_2B")){ pmtnosignal++; good = false;}}
1657      if ( this->bit(patterntrig[5],10) ){ pmtpattern++; if ( !this->checkPMT("S11_3B")){ pmtnosignal++; good = false;}}
1658      if ( this->bit(patterntrig[5],11) ){ pmtpattern++; if ( !this->checkPMT("S11_4B")){ pmtnosignal++; good = false;}}
1659      if ( this->bit(patterntrig[5],12) ){ pmtpattern++; if ( !this->checkPMT("S11_5B")){ pmtnosignal++; good = false;}}
1660      if ( this->bit(patterntrig[5],13) ){ pmtpattern++; if ( !this->checkPMT("S11_6B")){ pmtnosignal++; good = false;}}
1661      if ( this->bit(patterntrig[5],14) ){ pmtpattern++; if ( !this->checkPMT("S11_7B")){ pmtnosignal++; good = false;}}
1662      if ( this->bit(patterntrig[5],15) ){ pmtpattern++; if ( !this->checkPMT("S11_8B")){ pmtnosignal++; good = false;}}
1663    
1664      return good;
1665    }
1666    
1667    bool ToFLevel2::checkPMTpmttrig(TrigLevel2 *trg){
1668      //  UInt_t *patterntrig = trg->patterntrig;
1669      int rS11 = 0;
1670      int rS12 = 0;
1671      int rS21 = 0;
1672      int rS22 = 0;
1673      int rS31 = 0;
1674      int rS32 = 0;
1675    
1676      // trigger configuration for the event from saved pmts
1677      TClonesArray* Pmt = this->PMT;
1678      for(int i=0; i<Pmt->GetEntries(); i++) {  
1679        ToFPMT* pmthit = (ToFPMT*)Pmt->At(i);
1680        TString pmtname = this->GetPMTName(pmthit->pmt_id);
1681        if ( pmtname.Contains("S11") ) rS11++;
1682        if ( pmtname.Contains("S12") ) rS12++;
1683        if ( pmtname.Contains("S21") ) rS21++;
1684        if ( pmtname.Contains("S22") ) rS22++;
1685        if ( pmtname.Contains("S31") ) rS31++;
1686        if ( pmtname.Contains("S32") ) rS32++;
1687      }
1688      int rTOF1 = (rS11 + rS12) * (rS21 + rS22) * (rS31 + rS32);
1689      int rTOF2 = (rS11 * rS12) * (rS21 * rS22) * (rS31 * rS32);
1690    
1691      int rTOF3 = (rS21 + rS22) * (rS31 + rS32);
1692      int rTOF4 = (rS21 * rS22) * (rS31 * rS32);
1693    
1694      int rTOF5 = rS12 * (rS21 * rS22);
1695    
1696      int rTOF6 = (rS11 + rS12) * (rS31 + rS32);
1697      int rTOF7 = (rS11 * rS12) * (rS31 * rS32);
1698    
1699    
1700      // trigger configuration of the run
1701      bool TCTOF1 = false;
1702      bool TCTOF2 = false;
1703      bool TCTOF3 = false;
1704      bool TCTOF4 = false;
1705      bool TCTOF5 = false;
1706      bool TCTOF6 = false;
1707      bool TCTOF7 = false;
1708      if ( trg->trigconf & (1<<0) ) TCTOF1 = true;
1709      if ( trg->trigconf & (1<<1) ) TCTOF2 = true;
1710      if ( trg->trigconf & (1<<2) ) TCTOF3 = true;
1711      if ( trg->trigconf & (1<<3) ) TCTOF4 = true;
1712      if ( trg->trigconf & (1<<4) ) TCTOF5 = true;
1713      if ( trg->trigconf & (1<<5) ) TCTOF6 = true;
1714      if ( trg->trigconf & (1<<6) ) TCTOF7 = true;
1715    
1716      // do patterntrig pmts match the trigger configuration?
1717      bool pmtsconf_trigconf_match = true;
1718      if ( rTOF1 == 0 && TCTOF1 ) pmtsconf_trigconf_match = false;
1719      if ( rTOF2 == 0 && TCTOF2 ) pmtsconf_trigconf_match = false;
1720      if ( rTOF3 == 0 && TCTOF3 ) pmtsconf_trigconf_match = false;
1721      if ( rTOF4 == 0 && TCTOF4 ) pmtsconf_trigconf_match = false;
1722      if ( rTOF5 == 0 && TCTOF5 ) pmtsconf_trigconf_match = false;
1723      if ( rTOF6 == 0 && TCTOF6 ) pmtsconf_trigconf_match = false;
1724      if ( rTOF7 == 0 && TCTOF7 ) pmtsconf_trigconf_match = false;
1725    
1726      return pmtsconf_trigconf_match;
1727    }
1728    
1729    void ToFLevel2::printPMT(){
1730      TClonesArray* Pmt = this->PMT;
1731      for(int i=0; i<Pmt->GetEntries(); i++) {  
1732        ToFPMT* pmthit = (ToFPMT*)Pmt->At(i);
1733        TString pmtname = this->GetPMTName(pmthit->pmt_id);
1734        printf(" PMT hit: %s \n",pmtname.Data());
1735      }
1736    }
1737    
1738    
1739    ToFdEdx::ToFdEdx()
1740    {
1741      memset(conn,0,12*sizeof(Bool_t));
1742      memset(ts,0,12*sizeof(UInt_t));
1743      memset(te,0,12*sizeof(UInt_t));
1744      eDEDXpmt = new TArrayF(48);
1745      Define_PMTsat();
1746      Clear();
1747    }
1748    
1749    ToFdEdx::~ToFdEdx(){
1750      Clear();
1751      Delete();
1752    }
1753    
1754    void ToFdEdx::Delete(Option_t *option){
1755      if ( eDEDXpmt ){
1756        eDEDXpmt->Set(0);
1757        if ( eDEDXpmt) delete eDEDXpmt;
1758      }
1759    }
1760    
1761    //------------------------------------------------------------------------
1762    void ToFdEdx::CheckConnectors(UInt_t atime, GL_PARAM *glparam, TSQLServer *dbc)
1763    {
1764      for(int i=0; i<12; i++){
1765        if(atime<=ts[i] || atime>te[i]){
1766          Int_t error=glparam->Query_GL_PARAM(atime,210+i,dbc); // parameters stored in DB in GL_PRAM table
1767          if ( error<0 ) {
1768            conn[i]=false;
1769            ts[i]=0;
1770            te[i]=numeric_limits<UInt_t>::max();
1771          };
1772          if ( !error ){
1773            conn[i]=true;
1774            ts[i]=glparam->FROM_TIME;
1775            te[i]=glparam->TO_TIME;
1776          }
1777          if ( error>0 ){
1778            conn[i]=false;
1779            ts[i]=glparam->TO_TIME;
1780            TSQLResult *pResult;
1781            TSQLRow *row;
1782            TString query= Form("SELECT FROM_TIME FROM GL_PARAM WHERE TYPE=%i AND FROM_TIME>=%i ORDER BY FROM_TIME ASC LIMIT 1;",210+i,atime);
1783            pResult=dbc->Query(query.Data());
1784            if(!pResult->GetRowCount()){
1785              te[i]=numeric_limits<UInt_t>::max();
1786            }else{
1787              row=pResult->Next();
1788              te[i]=(UInt_t)atoll(row->GetField(0));
1789            }
1790          }
1791          //
1792          
1793        }
1794      }
1795    
1796    }
1797    //------------------------------------------------------------------------
1798    void ToFdEdx::Clear(Option_t *option)
1799    {
1800      //
1801      // Set arrays and initialize structure
1802      //  eDEDXpmt.Set(48);    eDEDXpmt.Reset(-1);   // Set array size  and reset structure
1803      eDEDXpmt->Set(48);    eDEDXpmt->Reset(-1);   // Set array size  and reset structure
1804      //
1805    };
1806    
1807    //------------------------------------------------------------------------
1808    void ToFdEdx::Print(Option_t *option)
1809    {
1810      //
1811      printf("========================================================================\n");
1812    
1813    };
1814    
1815    //------------------------------------------------------------------------
1816    void ToFdEdx::Init(pamela::tof::TofEvent *tofl0)
1817    {
1818      //
1819      ToFLevel2 tf;
1820      for (Int_t gg=0; gg<4;gg++){
1821        for (Int_t hh=0; hh<12;hh++){
1822          //          tofinput_.tdc[hh][gg]=tofEvent->tdc[gg][hh];          
1823          int mm = tf.GetPMTid(gg,hh);        
1824          adc[mm]= (0xFFF & tofl0->adc[gg][hh]); // EM, exclude warning bits
1825        };      
1826      };
1827      
1828    };
1829    
1830    //------------------------------------------------------------------------
1831    void ToFdEdx::Init(Int_t gg, Int_t hh, Float_t adce)
1832    {
1833      //
1834      ToFLevel2 tf;
1835      //  for (Int_t gg=0; gg<4;gg++){
1836      //    for (Int_t hh=0; hh<12;hh++){
1837      int mm = tf.GetPMTid(gg,hh);    
1838      adc[mm]=adce;
1839      
1840    };
1841    //------------------------------------------------------------------------
1842    void ToFdEdx::Process(UInt_t atime, Float_t betamean, Float_t *xtr_tof, Float_t *ytr_tof, Int_t exitat)
1843    {
1844      bool debug = false;
1845      if ( debug ) printf(" INSIDE TOFDEDX PROCESS \n");
1846      // the parameters should be already initialised by InitPar()
1847      //  printf(" in process \n");
1848      Clear();
1849    
1850     // define angle:  
1851      double dx   = xtr_tof[1] - xtr_tof[5];
1852      double dy   = ytr_tof[0] - ytr_tof[4];
1853      double dr   = sqrt(dx*dx+dy*dy);
1854      double theta=atan(dr/76.81);
1855      //
1856      if ( xtr_tof[1] > 99. ||  xtr_tof[5] > 99. || ytr_tof[0] > 99. ||  ytr_tof[4] > 99. ) theta = 0.;
1857      for (Int_t ii=0; ii<6; ii++){
1858        if ( xtr_tof[ii] > 99. ) xtr_tof[ii] = 0.;
1859        if ( ytr_tof[ii] > 99. ) ytr_tof[ii] = 0.;
1860      };
1861      //
1862      if ( debug ) printf(" theta %f \n",theta);
1863      if ( debug ) printf(" xtr_tof %.1f %.1f %.1f %.1f %.1f %.1f \n",xtr_tof[0],xtr_tof[1],xtr_tof[2],xtr_tof[3],xtr_tof[4],xtr_tof[5]);
1864      if ( debug ) printf(" ytr_tof %.1f %.1f %.1f %.1f %.1f %.1f \n",ytr_tof[0],ytr_tof[1],ytr_tof[2],ytr_tof[3],ytr_tof[4],ytr_tof[5]);
1865      //--------------------- TABLE OF PERIODS WITH HV PROBLEMS ----------------------------
1866      
1867      int Aconn=conn[0];    // PMT 0,20,22,24
1868      int Bconn=conn[1];    // PMT 6,12,26,34
1869      int Cconn=conn[2];    // PMT 4,14,28,32
1870      int Dconn=conn[3];    // PMT 2,8,10,30
1871      int Econn=conn[4];    // PMT 42,43,44,47
1872      int Fconn=conn[5];    // PMT 7,19,23,27
1873      int Gconn=conn[6];    // PMT 3,11,25,33
1874      int Hconn=conn[7];    // PMT 1,9,13,21
1875      int Iconn=conn[8];    // PMT 5,29,31,35
1876      int Lconn=conn[9];    // PMT 37,40,45,46
1877      int Mconn=conn[10];    // PMT 15,16,17,18
1878      int Nconn=conn[11];    // PMT 36,38,39,41
1879      if( false ) cout << Gconn << Iconn << Lconn <<endl; // to avoid compilation warnings
1880        
1881      //  printf(" size %i \n",eDEDXpmt.GetSize());
1882      for( int ii=0; ii<48; ii++ ) {
1883        //
1884        //    eDEDXpmt.SetAt(-1.,ii);
1885        //    printf(" ii %i beta %f atime %u xtr 1 %f ytr 1 %f adc %f \n",ii,betamean,atime,xtr_tof[0],ytr_tof[0],adc[ii]);
1886        if ( debug ) printf("II %i adc %f \n",ii,adc[ii]);
1887    
1888        if( adc[ii] >= 4095. ){
1889          //      eDEDXpmt[ii] = 0.;
1890          eDEDXpmt->AddAt(0.,ii);
1891          if ( debug ) printf(" %i adc>4095 \n",ii);
1892          continue; // EMILIANO
1893        };
1894    
1895        if( adc[ii] >= (PMTsat[ii]-5.) && adc[ii] < 4095. ){
1896          eDEDXpmt->AddAt(1000.,ii);
1897          if ( debug ) printf(" %i adc> pmtsat && adc<4095 \n",ii);
1898          continue; // EMILIANO
1899        };
1900    
1901        if( adc[ii] <= 0. ) {
1902          eDEDXpmt->AddAt(1500.,ii);
1903          if ( debug ) printf(" %i adc<=0 \n",ii);
1904          continue;
1905        };
1906        //
1907        double adcpC   = f_adcPC( adc[ii] );    // - adc conversion in pC
1908        if ( exitat == 0 ){
1909          eDEDXpmt->AddAt((Float_t)adcpC,ii);
1910          continue;
1911        }
1912        //    printf(" e qua? \n");
1913    
1914        double adccorr = adcpC*fabs(cos(theta));    
1915        if ( debug ) printf(" adccorr %f \n",adccorr);
1916        if(adccorr<=0.){
1917          if ( debug ) printf(" %i adccorr<=0 \n",ii);
1918          //      eDEDXpmt->AddAt((Float_t)adcpC,ii);//?
1919          continue;
1920        }
1921        if ( exitat == 1 ){
1922          eDEDXpmt->AddAt((Float_t)adccorr,ii);
1923          continue;
1924        }
1925        //    printf(" e quo? \n");
1926    
1927        //    int standard=0;
1928        int S115B_ok=0;
1929        int S115B_break=0;
1930    
1931        if(atime<1158720000)S115B_ok=1;
1932        else S115B_break=1;
1933    
1934    
1935        //------------------------------------------------------------------------
1936        //    printf(" e qui? \n");
1937        //---------------------------------------------------- Z reconstruction
1938    
1939        double adcHe, adcnorm, adclin, dEdx;//, Zeta; // EM GCC4.7
1940    
1941        adcHe=-2;
1942        adcnorm=-2;
1943        adclin=-2;
1944        dEdx=-2;
1945        //    Zeta=-2;//EM GCC4.7
1946        Double_t correction = 1.;
1947    
1948        if(Aconn==1 && (ii==0 || ii==20 || ii==22 || ii==24)){
1949          correction = 1.675;
1950        }
1951        else if(Bconn==1 && (ii==6 || ii==12 || ii==26 || ii==34)){
1952          correction = 2.482;
1953        }
1954        else if(Cconn==1 && (ii==4 || ii==14 || ii==28 || ii==32)){
1955          correction = 1.464;
1956        }
1957        else if(Dconn==1 && (ii==2 || ii==8 || ii==10 || ii==30)){
1958          correction = 1.995;
1959        }
1960        else if(Econn==1 && (ii==42 || ii==43 || ii==44 || ii==47)){
1961          correction = 1.273;
1962        }
1963        else if(Fconn==1 && (ii==7 || ii==19 || ii==23 || ii==27)){
1964          correction = 1.565;
1965        }
1966        else if(Mconn==1 && (ii==15 || ii==16 || ii==17 || ii==18)){
1967          correction = 1.565;
1968        }
1969        else if(Nconn==1 && (ii==36 || ii==38 || ii==39 || ii==41)){
1970          correction = 1.018;
1971        }
1972        else if(Hconn==1 && (ii==1 || ii==13 || ii==21 || (ii==9&&S115B_ok==1))){
1973          correction = 1.84;
1974        }
1975        else if(S115B_break==1 && ii==9 && Hconn==1){
1976          correction = 1.64;
1977        }
1978        else correction = 1.;
1979        
1980        if( ii==9 && S115B_break==1 ){
1981          adcHe   = f_att5B( ytr_tof[0] )/correction;
1982        } else {
1983          adcHe   = Get_adc_he(ii, xtr_tof, ytr_tof)/correction;
1984        };
1985        if(adcHe<=0){
1986          if ( debug ) printf(" %i adcHe<=0 \n",ii);
1987          //      eDEDXpmt->AddAt((Float_t)adccorr,ii); //?
1988          continue;
1989        }
1990        if ( exitat == 2 ){
1991          if(ii==9 && S115B_break==1)  eDEDXpmt->AddAt(36.*(Float_t)adccorr/adcHe,ii);
1992          else  adclin  = 4.*(Float_t)adccorr/adcHe;
1993          continue;
1994        }
1995    
1996        if(ii==9 && S115B_break==1)  adcnorm = f_pos5B(adccorr);
1997        else adcnorm = f_pos( (parPos[ii]), adccorr);
1998        if(adcnorm<=0){
1999          if ( debug ) printf(" %i adcnorm<=0 \n",ii);
2000          //      eDEDXpmt->AddAt((Float_t)adccorr,ii);//?
2001          continue;
2002        }
2003        if ( debug ) printf(" adcnorm %f \n",adcnorm);
2004    
2005        if(ii==9 && S115B_break==1)  adclin  = 36.*adcnorm/adcHe;
2006        else  adclin  = 4.*adcnorm/adcHe;
2007        if ( debug ) printf(" adclin %f \n",adclin);
2008        if(adclin<=0){
2009          if ( debug ) printf(" %i adclin<=0 \n",ii);
2010          //      eDEDXpmt->AddAt((Float_t)adccorr,ii);//?
2011          continue;
2012        }
2013        if ( exitat == 3 ){
2014          if(ii==9 && S115B_break==1)  eDEDXpmt->AddAt((Float_t)adclin,ii);
2015          else  eDEDXpmt->AddAt((Float_t)adclin,ii);
2016          continue;
2017        }
2018        //
2019        if ( betamean > 99. ){
2020          //      eDEDXpmt.AddAt((Float_t)adclin,ii);
2021          eDEDXpmt->AddAt((Float_t)adclin,ii);
2022          //      printf(" AAPMT IS %i dedx is %f vector is %f \n",ii,adclin,eDEDXpmt[ii]);
2023          if ( debug ) printf(" %i betamean > 99 \n",ii);
2024          continue;
2025        };
2026        //
2027        double dEdxHe=-2;
2028        if(ii==9 && S115B_break==1){
2029          if( betamean <1. ) dEdxHe = f_BB5B( betamean );
2030          else                       dEdxHe = 33;
2031        } else {
2032          if( betamean <1. ) dEdxHe = f_BB( (parBBneg[ii]), betamean );
2033          else                       dEdxHe = parBBpos[ii];
2034        }
2035        
2036        if ( debug ) printf(" dEdxHe %f \n",dEdxHe);
2037        
2038        if(dEdxHe<=0){
2039          eDEDXpmt->AddAt((Float_t)adclin,ii);
2040          if ( debug ) printf(" %i dEdxHe<=0 \n",ii);
2041          continue;
2042        };
2043    
2044        if(ii==9 && S115B_break==1)  dEdx = f_desatBB5B( adclin );
2045        else  dEdx = f_desatBB((parDesatBB[ii]), adclin );
2046    
2047        if(dEdx<=0){
2048          eDEDXpmt->AddAt((Float_t)adclin,ii);
2049          if ( debug ) printf(" %i dEdx<=0 \n",ii);
2050          continue;
2051        };
2052    
2053        if ( debug ) printf(" dEdx %f \n",dEdx);
2054        eDEDXpmt->AddAt((Float_t)dEdx,ii);
2055        //    eDEDXpmt.AddAt((Float_t)dEdx,ii);
2056    
2057        //    printf(" PMT IS %i dedx is %f vector is %f \n",ii,dEdx,eDEDXpmt[ii]);
2058    
2059      }  //end loop on 48 PMT
2060    
2061    };
2062    
2063    
2064    //------------------------------------------------------------------------
2065    void ToFdEdx::Define_PMTsat()
2066    {
2067      Float_t  sat[48] = {
2068        3176.35,3178.19,3167.38,3099.73,3117.00,3126.29,3111.44,3092.27,
2069        3146.48,3094.41,3132.13,3115.37,3099.32,3110.97,3111.80,3143.14,
2070        3106.72,3153.44,3136.00,3188.96,3104.73,3140.45,3073.18,3106.62,
2071        3112.48,3146.92,3127.24,3136.52,3109.59,3112.89,3045.15,3147.26,
2072        3095.92,3121.05,3083.25,3123.62,3150.92,3125.30,3067.60,3160.18,
2073        3119.36,3108.92,3164.77,3133.64,3111.47,3131.98,3128.87,3135.56 };
2074      PMTsat.Set(48,sat);
2075    }
2076    
2077    //------------------------------------------------------------------------
2078    void ToFdEdx::ReadParBBpos( const char *fname )
2079    {
2080      //  printf("read %s\n",fname);
2081      parBBpos.Set(48);
2082      FILE *fattin = fopen( fname , "r" );
2083      for (int i=0; i<48; i++) {
2084        int   tid=0;
2085        float  tp;
2086        if(fscanf(fattin,"%d %f",
2087                  &tid, &tp )!=2) break;
2088        parBBpos[i]=tp;
2089      }
2090      fclose(fattin);
2091    }
2092    
2093    //------------------------------------------------------------------------
2094    void ToFdEdx::ReadParDesatBB( const char *fname )
2095    {
2096      //  printf("read %s\n",fname);
2097      FILE *fattin = fopen( fname , "r" );
2098      for (int i=0; i<48; i++) {
2099        int   tid=0;
2100        float  tp[3];
2101        if(fscanf(fattin,"%d %f %f %f",
2102                  &tid, &tp[0], &tp[1], &tp[2] )!=4) break;
2103        parDesatBB[i].Set(3,tp);
2104      }
2105      fclose(fattin);
2106    }
2107    
2108    
2109    //------------------------------------------------------------------------
2110    void ToFdEdx::ReadParBBneg( const char *fname )
2111    
2112    {
2113      //  printf("read %s\n",fname);
2114      FILE *fattin = fopen( fname , "r" );
2115      for (int i=0; i<48; i++) {
2116        int   tid=0;
2117        float  tp[3];
2118        if(fscanf(fattin,"%d %f %f %f",
2119                  &tid, &tp[0], &tp[1], &tp[2] )!=4) break;
2120        parBBneg[i].Set(3,tp);
2121      }
2122      fclose(fattin);
2123    }
2124    
2125    //------------------------------------------------------------------------
2126    void ToFdEdx::ReadParPos( const char *fname )
2127    {
2128      //  printf("read %s\n",fname);
2129      FILE *fattin = fopen( fname , "r" );
2130      for (int i=0; i<48; i++) {
2131        int   tid=0;
2132        float  tp[4];
2133        if(fscanf(fattin,"%d %f %f %f %f",
2134                  &tid, &tp[0], &tp[1], &tp[2], &tp[3])!=5) break;
2135        parPos[i].Set(4,tp);
2136      }
2137      fclose(fattin);
2138    }
2139    
2140    //------------------------------------------------------------------------
2141    void ToFdEdx::ReadParAtt( const char *fname )
2142    {
2143      //  printf("read %s\n",fname);
2144      FILE *fattin = fopen( fname , "r" );
2145      for (int i=0; i<48; i++) {
2146        int   tid=0;
2147        float  tp[6];
2148        if(fscanf(fattin,"%d %f %f %f %f %f %f",
2149                  &tid, &tp[0], &tp[1], &tp[2], &tp[3], &tp[4], &tp[5] )!=7) break;
2150        parAtt[i].Set(6,tp);
2151      }
2152      fclose(fattin);
2153    }
2154    
2155    
2156    
2157    
2158    
2159    
2160    double ToFdEdx::f_att( TArrayF &p, float x )
2161    {
2162      return
2163        p[0] +
2164        p[1]*x +
2165        p[2]*x*x +
2166        p[3]*x*x*x +
2167        p[4]*x*x*x*x +
2168        p[5]*x*x*x*x*x;
2169    }
2170    //------------------------------------------------------------------------
2171    double ToFdEdx::f_att5B( float x )
2172    {
2173      return
2174        101.9409 +
2175        6.643781*x +
2176        0.2765518*x*x +
2177        0.004617647*x*x*x +
2178        0.0006195132*x*x*x*x +
2179        0.00002813734*x*x*x*x*x;
2180    }
2181    
2182    
2183    double ToFdEdx::f_pos( TArrayF &p, float x )
2184    {
2185      return
2186        p[0] +
2187        p[1]*x +
2188        p[2]*x*x +
2189        p[3]*x*x*x;
2190    }
2191    
2192    double ToFdEdx::f_pos5B( float x )
2193    {
2194      return
2195        15.45132 +
2196        0.8369721*x +
2197        0.0005*x*x;
2198    }
2199    
2200    
2201    
2202    double ToFdEdx::f_adcPC( float x )
2203    {
2204      return 28.12+0.6312*x-5.647e-05*x*x+3.064e-08*x*x*x;
2205    }
2206    
2207    
2208    float ToFdEdx::Get_adc_he( int id, float pl_x[6], float pl_y[6])
2209    {
2210    
2211      //
2212      // input: id - pmt [0:47}
2213      //             pl_x - coord x of the tof plane
2214      //             pl_y - coord y
2215    
2216      adc_he = 0;
2217      if( eGeom.GetXY(id)==1 )  adc_he = f_att( (parAtt[id]), pl_x[eGeom.GetPlane(id)] );
2218      if( eGeom.GetXY(id)==2 )  adc_he = f_att( (parAtt[id]), pl_y[eGeom.GetPlane(id)] );
2219      return adc_he;
2220    }
2221    
2222    //------------------------------------------------------------------------
2223    double ToFdEdx::f_BB( TArrayF &p, float x )
2224    {
2225      return  p[0]/(x*x)*(log(x*x/(1-x*x)) - p[1]*x*x - p[2]);
2226    }
2227    
2228    //------------------------------------------------------------------------
2229    double ToFdEdx::f_BB5B( float x )
2230    {
2231      return  0.165797/(x*x)*(log(x*x/(1-x*x)) + 140.481*x*x + 52.9258);
2232    }
2233    //------------------------------------------------------------------------
2234    double ToFdEdx::f_desatBB( TArrayF &p, float x )
2235    {
2236      return
2237        p[0] +
2238        p[1]*x +
2239        p[2]*x*x;
2240    }
2241    
2242    //------------------------------------------------------------------------
2243    double ToFdEdx::f_desatBB5B( float x )
2244    {
2245      return
2246        -2.4 +
2247        0.75*x +
2248        0.009*x*x;
2249    }
2250    
Legend:
Removed from v.1.19  
changed lines
  Added in v.1.47
  ViewVC Help
Powered by ViewVC 1.1.23