45 |
delete clong; |
delete clong; |
46 |
clong = new CaloLong(L2); |
clong = new CaloLong(L2); |
47 |
clong->SplitInto(0,22); |
clong->SplitInto(0,22); |
48 |
clong->HeavyTail(true); |
// clong->HeavyTail(true); |
49 |
|
clong->HeavyTail(false); |
50 |
}; |
}; |
51 |
if ( cp ) delete cp; |
if ( cp ) delete cp; |
52 |
cp = NULL; |
cp = NULL; |
74 |
clong = new CaloLong(L2); |
clong = new CaloLong(L2); |
75 |
if ( cp ) clong->SetCaloLevel2Pointer(cp->GetCaloLevel2Pointer()); |
if ( cp ) clong->SetCaloLevel2Pointer(cp->GetCaloLevel2Pointer()); |
76 |
clong->SplitInto(0,22); |
clong->SplitInto(0,22); |
77 |
clong->HeavyTail(true); |
// clong->HeavyTail(true); |
78 |
|
clong->HeavyTail(false); |
79 |
}; |
}; |
80 |
// |
// |
81 |
} |
} |
95 |
asntr = "start"; |
asntr = "start"; |
96 |
// |
// |
97 |
debug = false; |
debug = false; |
98 |
|
usepl18x = false; |
99 |
// |
// |
100 |
indep = false; |
indep = false; |
101 |
// |
// |
531 |
// |
// |
532 |
for (Int_t ch=0; ch< L2->GetCaloLevel1()->istrip; ch++){ |
for (Int_t ch=0; ch< L2->GetCaloLevel1()->istrip; ch++){ |
533 |
en = L2->GetCaloLevel1()->DecodeEstrip(ch,view,plane,strip); |
en = L2->GetCaloLevel1()->DecodeEstrip(ch,view,plane,strip); |
534 |
|
// |
535 |
|
if ( !usepl18x && view==0 && plane==18 ) en = 0.; |
536 |
|
// |
537 |
enstrip[view][plane][strip]=en; |
enstrip[view][plane][strip]=en; |
538 |
}; |
}; |
539 |
// |
// |
1140 |
// |
// |
1141 |
// check if the cylinder of integration can go out of the sensor given the frame which has been set (if we use all the calorimeter fRad is < 0 and the printout is suppressed) |
// check if the cylinder of integration can go out of the sensor given the frame which has been set (if we use all the calorimeter fRad is < 0 and the printout is suppressed) |
1142 |
// |
// |
1143 |
if ( (fM1+0.122-0.244*(Float_t)fRad) < 0. ) printf("Error: (fM1+0.122-0.244*(Float_t)fRad) < 0. fM1 %f fRad %i %f \n",fM1,fRad,(fM1+0.122-0.244*(Float_t)fRad)); |
if ( (fM1+0.122-0.244*(Float_t)fRad) < 0.001 ) printf("Error: (fM1+0.122-0.244*(Float_t)fRad) < 0. fM1 %f fRad %i %f \n",fM1,fRad,(fM1+0.122-0.244*(Float_t)fRad)); |
1144 |
// |
// |
1145 |
if ( fLong ){ |
if ( fLong ){ |
1146 |
if ( debug ) printf(" ==================================================================> LONGITUDINAL FIT! \n"); |
if ( debug ) printf(" ==================================================================> LONGITUDINAL FIT! \n"); |