| 799 |
|
const std::vector<DTChamber*> & chs = dtGeom->chambers(); |
| 800 |
|
|
| 801 |
|
|
| 802 |
< |
int qual[5][4][14]; |
| 802 |
> |
int qual[5][5][15]; |
| 803 |
|
for (int iw=0; iw<5; iw++) // Wheel 1-5 |
| 804 |
|
for (int sec=1; sec<15; ++sec) // section 1 to 14 |
| 805 |
< |
for (int ist=0; ist<5; ++ist) { |
| 805 |
> |
for (int ist=1; ist<5; ++ist) { |
| 806 |
|
qual[iw][ist-1][sec-1] = -1; |
| 807 |
|
} |
| 808 |
|
|
| 881 |
|
|
| 882 |
|
const DTChamber* ch = dtGeom->chamber(seg->chamberId()); |
| 883 |
|
GlobalPoint glbPoint = ch->toGlobal((*seg).localPosition()); |
| 884 |
< |
GlobalVector glbDir = ch->toGlobal((*seg).localDirection()); |
| 884 |
> |
//GlobalVector glbDir = ch->toGlobal((*seg).localDirection()); |
| 885 |
|
|
| 886 |
|
float radtodeg = 57.296; |
| 887 |
|
|
| 987 |
|
} |
| 988 |
|
// end trigger efficiency analysis ------------ |
| 989 |
|
// prepare vectors for re-fitting |
| 990 |
< |
float xfit[12], yfit [12]; |
| 991 |
< |
float sfit[12]; // flag (-1 or +1) for left/right semicell, input to fitline_t0 function |
| 990 |
> |
//float xfit[12], yfit [12]; |
| 991 |
> |
//float sfit[12]; // flag (-1 or +1) for left/right semicell, input to fitline_t0 function |
| 992 |
|
DTSuperLayerId slid1(phiSeg->chamberId(),1); |
| 993 |
|
/// Mean timer analysis |
| 994 |
|
DTMeanTimer meanTimer1(dtGeom->superLayer(slid1), phiHits, eventSetup, theSync); |
| 1013 |
|
float xhit = (specificHits[sphit]).localPosition().x() + laylocal.x() ; // localPosition is in the layer (not in the cell) |
| 1014 |
|
DTWireId wireId = (specificHits[sphit]).wireId(); |
| 1015 |
|
float ttrig = theSync->offset(wireId); |
| 1016 |
< |
int ill = il; |
| 1017 |
< |
if ( isl > 1) ill = il + 4; |
| 1016 |
> |
//int ill = il; |
| 1017 |
> |
//if ( isl > 1) ill = il + 4; |
| 1018 |
|
if (il == 1 || il == 3 ) tmax123[isl] = tmax123[isl] + 0.5*( (specificHits[sphit]).digiTime() - ttrig ); |
| 1019 |
|
if (il == 2) tmax123[isl] = tmax123[isl] + (specificHits[sphit]).digiTime() - ttrig; |
| 1020 |
< |
xfit[nhitfit] = laylocal.z(); // input to fitline function |
| 1021 |
< |
yfit[nhitfit] = xhit; // input to fitline function |
| 1022 |
< |
sfit[nhitfit] = 1.; // mark left/right semi-cell (1=Right) input to fitline_t0 function ) |
| 1023 |
< |
if( (specificHits[sphit]).lrSide() == 2 ) sfit[nhitfit] = -1.; // -1 = Left |
| 1020 |
> |
//xfit[nhitfit] = laylocal.z(); // input to fitline function |
| 1021 |
> |
//yfit[nhitfit] = xhit; // input to fitline function |
| 1022 |
> |
//sfit[nhitfit] = 1.; // mark left/right semi-cell (1=Right) input to fitline_t0 function ) |
| 1023 |
> |
//if( (specificHits[sphit]).lrSide() == 2 ) sfit[nhitfit] = -1.; // -1 = Left |
| 1024 |
|
nhitfit++; |
| 1025 |
|
// ============== histograms of residuals ====================== |
| 1026 |
|
stringstream hTitle; |
| 1349 |
|
|
| 1350 |
|
int SCsect=0; int SCst=0; |
| 1351 |
|
|
| 1352 |
< |
int bx[5][4][5]; |
| 1352 |
> |
// int bx[5][4][5]; |
| 1353 |
|
float bxbest[5][4][14]; |
| 1354 |
|
// float qual[]={-1,-1,-1,-1}; |
| 1355 |
|
int qual[5][4][14]; |
| 1372 |
|
|
| 1373 |
|
if(doWheel[iw]) |
| 1374 |
|
{ |
| 1375 |
< |
bx[iw][SCst-1][0]=0; bx[iw][SCst-1][1]=0; |
| 1375 |
> |
//bx[iw][SCst-1][0]=0; bx[iw][SCst-1][1]=0; |
| 1376 |
|
|
| 1377 |
|
int countsec=-1; |
| 1378 |
|
for (int ns=0; ns<12; ns++) { |
| 1399 |
|
hasTr[iw][SCsect-1][SCst-1]=true; |
| 1400 |
|
} |
| 1401 |
|
if (trigtrack->quality()<7) { |
| 1402 |
< |
bx[iw][SCst-1][ntrCh]=trigtrack->bx(); // dice a che bx c'e' stato trigger |
| 1402 |
> |
//bx[iw][SCst-1][ntrCh]=trigtrack->bx(); // dice a che bx c'e' stato trigger |
| 1403 |
|
float ibx=trigtrack->bx(); |
| 1404 |
|
ntrCh++; |
| 1405 |
|
if (ntrCh>4) break; |
