| 26 |
|
|
| 27 |
void CaloNuclei::Clear(){ |
void CaloNuclei::Clear(){ |
| 28 |
// |
// |
| 29 |
|
tr = 0; |
| 30 |
interplane = 0; |
interplane = 0; |
| 31 |
preq = 0.; |
preq = 0.; |
| 32 |
postq = 0.; |
postq = 0.; |
| 44 |
// |
// |
| 45 |
Process(); |
Process(); |
| 46 |
// |
// |
| 47 |
printf("===================================================================\n"); |
printf("========================================================================\n"); |
| 48 |
printf(" OBT: %u PKT: %u ATIME: %u \n",OBT,PKT,atime); |
printf(" OBT: %u PKT: %u ATIME: %u Track %i \n",OBT,PKT,atime,tr); |
| 49 |
printf(" interplane [number of available dE/dx before interaction]: %i\n",interplane); |
printf(" interplane [number of available dE/dx before interaction]: %i\n",interplane); |
| 50 |
printf(" ethr [threshold used to determine interplane]:............ %f \n",ethr); |
printf(" ethr [threshold used to determine interplane]:............ %f \n",ethr); |
| 51 |
printf(" dedx1 [dE/dx from the first calorimeter plane]:........... %f \n",dedx1); |
printf(" dedx1 [dE/dx from the first calorimeter plane]:........... %f \n",dedx1); |
| 58 |
printf(" N [no of used plane]:..................................... %i \n",N); |
printf(" N [no of used plane]:..................................... %i \n",N); |
| 59 |
printf(" R [no strip used per plane ]:............................. %i \n",R); |
printf(" R [no strip used per plane ]:............................. %i \n",R); |
| 60 |
printf(" qpremeanN [truncated mean using N planes and R strips]:... %f \n",qpremeanN); |
printf(" qpremeanN [truncated mean using N planes and R strips]:... %f \n",qpremeanN); |
| 61 |
printf("===================================================================\n"); |
printf("========================================================================\n"); |
| 62 |
// |
// |
| 63 |
}; |
}; |
| 64 |
|
|
| 69 |
|
|
| 70 |
|
|
| 71 |
void CaloNuclei::Process(){ |
void CaloNuclei::Process(){ |
| 72 |
// |
Process(0); |
| 73 |
|
}; |
| 74 |
|
|
| 75 |
|
void CaloNuclei::Process(Int_t ntr){ |
| 76 |
|
// |
| 77 |
if ( !L2 ){ |
if ( !L2 ){ |
| 78 |
printf(" ERROR: cannot find PamLevel2 object, use the correct constructor or check your program!\n"); |
printf(" ERROR: cannot find PamLevel2 object, use the correct constructor or check your program!\n"); |
| 79 |
printf(" ERROR: CaloNuclei variables not filled \n"); |
printf(" ERROR: CaloNuclei variables not filled \n"); |
| 95 |
// |
// |
| 96 |
if ( !newentry ) return; |
if ( !newentry ) return; |
| 97 |
// |
// |
| 98 |
|
tr = ntr; |
| 99 |
|
// |
| 100 |
// printf(" Processing event at OBT %u PKT %u time %u \n",OBT,PKT,atime); |
// printf(" Processing event at OBT %u PKT %u time %u \n",OBT,PKT,atime); |
| 101 |
// |
// |
| 102 |
Clear(); |
Clear(); |
| 103 |
// |
// |
| 104 |
PamTrack *track = 0; |
PamTrack *track = 0; |
| 105 |
track = L2->GetTrack(0); |
track = L2->GetTrack(ntr); |
| 106 |
// |
// |
| 107 |
Int_t view = 0; |
Int_t view = 0; |
| 108 |
Int_t plane = 0; |
Int_t plane = 0; |