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> |
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; |
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( |
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. |
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.; |
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++){ |
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); |
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] = { |
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 = " "; |
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] = { |