| 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); |
| 49 |
|
clong->HeavyTail(false); |
| 50 |
}; |
}; |
| 51 |
if ( cp ) delete cp; |
if ( cp ) delete cp; |
| 52 |
cp = NULL; |
cp = NULL; |
| 63 |
// cp->ForceCaloFit(); |
// cp->ForceCaloFit(); |
| 64 |
// cp->SetDebug(true); |
// cp->SetDebug(true); |
| 65 |
// cp->Process(); |
// cp->Process(); |
| 66 |
if ( clong ) clong->SetCaloLevel2Pointer(cp->GetLevel2Pointer()); |
if ( clong ) clong->SetCaloLevel2Pointer(cp->GetCaloLevel2Pointer()); |
| 67 |
} |
} |
| 68 |
|
|
| 69 |
|
|
| 72 |
fLong = true; |
fLong = true; |
| 73 |
if ( !clong ){ |
if ( !clong ){ |
| 74 |
clong = new CaloLong(L2); |
clong = new CaloLong(L2); |
| 75 |
if ( cp ) clong->SetCaloLevel2Pointer(cp->GetLevel2Pointer()); |
if ( cp ) clong->SetCaloLevel2Pointer(cp->GetCaloLevel2Pointer()); |
| 76 |
clong->SplitInto(0,22); |
clong->SplitInto(0,22); |
| 77 |
|
// 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.000000001 ) 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"); |