/[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.10 by pam-fi, Thu Jan 11 09:34:27 2007 UTC revision 1.15 by pam-de, Wed Apr 4 07:23:27 2007 UTC
# Line 1  Line 1 
1    /**
2     * \file ToFLevel2.cpp
3     * \author Gianfranca DeRosa, Wolfgang Menn
4     */
5    
6  #include <TObject.h>  #include <TObject.h>
7  #include <ToFLevel2.h>  #include <ToFLevel2.h>
8  #include <iostream>  #include <iostream>
# Line 88  ToFLevel2::ToFLevel2() {     Line 93  ToFLevel2::ToFLevel2() {    
93    //    //
94  };  };
95    
96    void ToFLevel2::Set(){//ELENA
97        if(!PMT)PMT = new TClonesArray("ToFPMT",12); //ELENA
98        if(!ToFTrk)ToFTrk = new TClonesArray("ToFTrkVar",2); //ELENA
99    }//ELENA
100    
101  void ToFLevel2::Clear(){  void ToFLevel2::Clear(){
102    //    //
103    if(ToFTrk)ToFTrk->Delete(); //ELENA    if(ToFTrk)ToFTrk->Delete(); //ELENA
# Line 143  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit) Line 153  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
153  //--------------------------------------  //--------------------------------------
154  /**  /**
155   * Method to get the plane ID (11 12 21 22 31 32) from the plane index (0 1 2 3 4 5)   * Method to get the plane ID (11 12 21 22 31 32) from the plane index (0 1 2 3 4 5)
156     * @param Plane index (0,1,2,3,4,5)
157   */   */
158    Int_t  ToFLevel2::GetToFPlaneID(Int_t ip){    Int_t  ToFLevel2::GetToFPlaneID(Int_t ip){
159        if(ip>=0 && ip<6)return 10*((int)(ip/2+1.1))+(ip%2)+1;        if(ip>=0 && ip<6)return 10*((int)(ip/2+1.1))+(ip%2)+1;
# Line 150  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit) Line 161  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
161    };    };
162  /**  /**
163   * Method to get the plane index (0 1 2 3 4 5) from the plane ID (11 12 21 22 31 32)   * Method to get the plane index (0 1 2 3 4 5) from the plane ID (11 12 21 22 31 32)
164     * @param plane Plane ID (11, 12, 21, 22, 31, 32)
165   */   */
166    Int_t  ToFLevel2::GetToFPlaneIndex(Int_t plane_id){    Int_t  ToFLevel2::GetToFPlaneIndex(Int_t plane_id){
167        if(        if(
# Line 164  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit) Line 176  ToFPMT *ToFLevel2::GetToFPMT(Int_t ihit)
176    };    };
177  /**  /**
178   * Method to know if a given ToF paddle was hit, that is there is a TDC signal   * Method to know if a given ToF paddle was hit, that is there is a TDC signal
179   * from both PMTs   * from both PMTs. The method uses the "tof_j_flag" variable.
180   * @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).
181   * @param paddle_id Paddle ID.   * @param paddle_id Paddle ID.
182   * @return 1 if the paddle was hit.   * @return 1 if the paddle was hit.
# Line 196  Int_t ToFLevel2::GetNHitPaddles(Int_t pl Line 208  Int_t ToFLevel2::GetNHitPaddles(Int_t pl
208  };  };
209    
210    
211    /**
212     * Method to get the mean dEdx from a given ToF plane. This current version
213     * is just summing up all PMT signals, which will not give proper results,
214     * and needs a revision.
215     * @param notrack Track Number
216     * @param plane Plane index (0,1,2,3,4,5).
217     */
218  Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane){  Float_t ToFLevel2::GetdEdx(Int_t notrack, Int_t plane){
219    
220    Float_t dedx = 0.;    Float_t dedx = 0.;
# Line 217  Float_t ToFLevel2::GetdEdx(Int_t notrack Line 236  Float_t ToFLevel2::GetdEdx(Int_t notrack
236      //      //
237      pl = GetPlaneIndex(pmt_id);      pl = GetPlaneIndex(pmt_id);
238      //      //
239      if ( pl == ip ){      if ( pl == ip ) dedx += (trk->dedx).At(i);
 //      // --- patch ---  
 //      bool ISNULL=false;  
 //      for(Int_t im=0; im< this->npmt(); im++){              
 //          if(this->GetToFPMT(im)->pmt_id == pmt_id && this->GetToFPMT(im)->adc == 4095){  
 //              cout << " ToFLevel2::GetdEdx(Int_t,Int_t) --> **warning** ADC(PMT="<<pmt_id<<") = 4095"<<endl;  
 //              ISNULL=true;  
 //          }  
 //      }  
 //      // --- patch ---  
 //      dedx += (trk->dedx).At(i)*(Int_t)(!ISNULL);  
         dedx += (trk->dedx).At(i);  
     }  
240      //      //
241    };    };
242    //    //
# Line 237  Float_t ToFLevel2::GetdEdx(Int_t notrack Line 244  Float_t ToFLevel2::GetdEdx(Int_t notrack
244  };  };
245    
246    
247    /**
248     * Method to fill the ADC_C 4x12 matrix with the dEdx values and the TDC 4x12 matrix
249     * with the time-walk corrected TDC values.
250     * @param notrack Track Number
251     * @param adc  ADC_C matrix with dEdx values
252     * @param tdc  TDC matrix
253     */
254  void ToFLevel2::GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]){  void ToFLevel2::GetMatrix(Int_t notrack, Float_t adc[4][12], Float_t tdc[4][12]){
255    //    //
256    for (Int_t aa=0; aa<4;aa++){    for (Int_t aa=0; aa<4;aa++){
# Line 277  void ToFLevel2::GetMatrix(Int_t notrack, Line 291  void ToFLevel2::GetMatrix(Int_t notrack,
291  };  };
292    
293    
294    /**
295     * Method to get the plane index (0 - 5) for the PMT_ID as input
296     * @param pmt_id  PMT_ID (0 - 47)
297     */
298  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt_id){  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt_id){
299    TString pmtname = GetPMTName(pmt_id);    TString pmtname = GetPMTName(pmt_id);
300    pmtname.Resize(3);    pmtname.Resize(3);
# Line 290  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt Line 307  Int_t ToFLevel2::GetPlaneIndex(Int_t pmt
307    return(-1);    return(-1);
308  };  };
309    
310    /**
311     * Method to get the PMT_ID if the index (4,12) is given. We have 4 channels on
312     * each of the 12 half-boards, this method decodes which PMT is cables to which
313     * channel.
314     * @param hh Channel
315     * @param kk HalfBoard
316     */
317  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_t kk){
318    //    //
319    short tof[4][24] = {    short tof[4][24] = {
# Line 319  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_ Line 342  Int_t ToFLevel2::GetPMTid(Int_t hh, Int_
342    return ind;    return ind;
343  };  };
344    
345    /**
346     * Method to get the PMT name (like "S11_1A") if the PMT_ID is given. The first
347     * PMT=0 is S11_1A, then S11_1B, then S11_2A, and so on, up to S32_3B which is PMT=47
348     * @param ind  PMT_ID (0 - 47)
349     */
350  TString ToFLevel2::GetPMTName(Int_t ind){  TString ToFLevel2::GetPMTName(Int_t ind){
351        
352    TString pmtname = " ";    TString pmtname = " ";
# Line 339  TString ToFLevel2::GetPMTName(Int_t ind) Line 367  TString ToFLevel2::GetPMTName(Int_t ind)
367    return pmtname;    return pmtname;
368  };  };
369    
370    /**
371     * Method to get the PMT index if the PMT ID is given. This method is the
372     * "reverse" of method "GetPMTid"
373     * @param ind  PMT_ID (0 - 47)
374     * @param hb   HalfBoard
375     * @param ch   Channel
376     */
377  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){  void ToFLevel2::GetPMTIndex(Int_t ind, Int_t &hb, Int_t &ch){
378    //    //
379    short tof[4][24] = {    short tof[4][24] = {
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.23