| 4 |
* |
* |
| 5 |
**/ |
**/ |
| 6 |
#include <TObject.h> |
#include <TObject.h> |
| 7 |
|
#include <TMath.h> |
| 8 |
#include <CaloLevel1.h> |
#include <CaloLevel1.h> |
|
|
|
| 9 |
// |
// |
| 10 |
ClassImp(CaloStrip); |
ClassImp(CaloStrip); |
| 11 |
ClassImp(CaloLevel1); |
ClassImp(CaloLevel1); |
| 115 |
**/ |
**/ |
| 116 |
void CaloStrip::Set(Float_t X, Float_t Y, Float_t Z) { |
void CaloStrip::Set(Float_t X, Float_t Y, Float_t Z) { |
| 117 |
// |
// |
| 118 |
// Float_t lShift = 0.; |
fX = X; |
| 119 |
// Float_t lPos = 0.; |
fY = Y; |
| 120 |
// extern struct shift shift; |
fZ = Z; |
| 121 |
// // |
// |
| 122 |
// fView = view++; |
Float_t zplane[22]; |
| 123 |
// fPlane = plane++; |
zplane[0] = 0.; |
| 124 |
// fStrip = strip++; |
Int_t ii = 0; |
| 125 |
// if ( fPlane%2 ){ |
for ( Int_t i = 2; i < 23; i++){ |
| 126 |
// lShift = -0.5; |
ii = i-1; |
| 127 |
// } else { |
if ( i%2 ){ |
| 128 |
// lShift = 0.5; |
zplane[ii] = zplane[ii-1] - 10.09; |
| 129 |
// }; |
} else { |
| 130 |
// // |
zplane[ii] = zplane[ii-1] - 8.09; |
| 131 |
// shift.shift = lShift; |
}; |
| 132 |
// // |
}; |
| 133 |
// Float_t zplane[22]; |
// |
| 134 |
// zplane[0] = 0.; |
Float_t dzx[22]; |
| 135 |
// Int_t ii = 0; |
Float_t dzy[22]; |
| 136 |
// for ( Int_t i = 2; i < 23; i++){ |
for ( Int_t i=0; i < 22; i++){ |
| 137 |
// ii = i-1; |
dzx[i] = fabs(fZ*10. - (zplane[i] - 5.81 + CTZ)); |
| 138 |
// if ( i%2 ){ |
dzy[i] = fabs(fZ*10. - (zplane[i] + CTZ)); |
| 139 |
// zplane[ii] = zplane[ii-1] - 10.09 + CTZ; |
}; |
| 140 |
// } else { |
// |
| 141 |
// zplane[ii] = zplane[ii-1] - 8.09 + CTZ; |
Float_t minx = TMath::MinElement(22,dzx); |
| 142 |
// }; |
Float_t miny = TMath::MinElement(22,dzy); |
| 143 |
// }; |
// |
| 144 |
// // |
// find view |
| 145 |
// millim_(&fStrip,&lPos); |
// |
| 146 |
// // |
if ( minx < miny ){ |
| 147 |
// if ( fView == 1 ){ |
fView = 1; |
| 148 |
// // |
} else { |
| 149 |
// // X view |
fView = 2; |
| 150 |
// // |
}; |
| 151 |
// fX = lPos; |
// |
| 152 |
// fY = 0.; |
// find plane |
| 153 |
// fZ = zplane[fPlane-1] - 5.81; |
// |
| 154 |
// } else { |
Float_t pos = 0.; |
| 155 |
// // |
// |
| 156 |
// // Y view |
for ( Int_t i=0; i < 22; i++){ |
| 157 |
// // |
if ( fView == 1 ){ |
| 158 |
// fX = 0.; |
if ( dzx[i] == minx ){ |
| 159 |
// fY = lPos; |
fPlane = i+1; |
| 160 |
// fZ = zplane[fPlane-1]; |
pos = fX*10. + CTX; |
| 161 |
// }; |
}; |
| 162 |
|
} else { |
| 163 |
|
if ( dzy[i] == miny ){ |
| 164 |
|
fPlane = i+1; |
| 165 |
|
pos = fY*10. + CTY; |
| 166 |
|
}; |
| 167 |
|
}; |
| 168 |
|
}; |
| 169 |
|
// |
| 170 |
|
// find strip |
| 171 |
|
// |
| 172 |
|
Float_t dxy[96]; |
| 173 |
|
Float_t stpos = 0.; |
| 174 |
|
// |
| 175 |
|
CaloStrip *ca = new CaloStrip(); |
| 176 |
|
// |
| 177 |
|
for ( Int_t i=0; i < 96; i++){ |
| 178 |
|
ca->Set(fView-1,fPlane-1,i); |
| 179 |
|
if ( fView == 1 ){ |
| 180 |
|
stpos = ca->GetX()*10. + CTX; |
| 181 |
|
} else { |
| 182 |
|
stpos = ca->GetY()*10. + CTY; |
| 183 |
|
}; |
| 184 |
|
dxy[i] = fabs(pos - stpos); |
| 185 |
|
}; |
| 186 |
|
// |
| 187 |
|
delete ca; |
| 188 |
// |
// |
| 189 |
|
Float_t mins = TMath::MinElement(96,dxy); |
| 190 |
|
// |
| 191 |
|
for ( Int_t i=0; i < 96; i++){ |
| 192 |
|
if ( dxy[i] == mins ) fStrip = i+1; |
| 193 |
|
}; |
| 194 |
}; |
}; |
|
// DISTY = PIANO(I) |
|
|
// YY(I) = DISTY * TG(2) + CY |
|
|
// BAR(M,I) = YY(I) |
|
|
// CBAR(M,I) = (-YALIG + YY(I))/10. |
|
|
// C |
|
|
// ENDIF |
|
|
// CALL LASTRISCIA(BAR(M,I),IBAR(M,I)) |
|
|
// CIBAR(M,I) = IBAR(M,I) |
|
|
// ENDDO |
|
|
// ENDIF |
|
|
// C |
|
|
// ENDDO |
|
|
|
|
|
|
|
| 195 |
|
|
| 196 |
/** |
/** |
| 197 |
* CaloLevel1 constructor |
* CaloLevel1 constructor |