4236 |
UInt_t totime = 0; |
UInt_t totime = 0; |
4237 |
UInt_t obt = 0; |
UInt_t obt = 0; |
4238 |
UInt_t pkt = 0; |
UInt_t pkt = 0; |
4239 |
|
Float_t totped = 0.; |
4240 |
// |
// |
4241 |
tr = (TTree*)file->Get("CalibCalPed"); |
tr = (TTree*)file->Get("CalibCalPed"); |
4242 |
if ( !tr || tr->IsZombie() ) throw -21; |
if ( !tr || tr->IsZombie() ) throw -21; |
4254 |
if ( calibCalPed->cstwerr[section] ){ |
if ( calibCalPed->cstwerr[section] ){ |
4255 |
valid = 1; |
valid = 1; |
4256 |
if ( calibCalPed->cperror[section] ) valid = 0; |
if ( calibCalPed->cperror[section] ) valid = 0; |
4257 |
|
// |
4258 |
|
// check pedestal values for one plane, if all zeros calibration is not valid (calorimeter power problems) [8th data reduction bug, fixed on 25/11/2009 by E.M.] |
4259 |
|
// |
4260 |
|
totped = 0.; |
4261 |
|
Int_t ns = 0; |
4262 |
|
if ( section == 2 ) ns = 3; |
4263 |
|
if ( section == 3 ) ns = 1; |
4264 |
|
if ( section == 1 ) ns = 2; |
4265 |
|
for (UInt_t ss=0; ss<96; ss++){ |
4266 |
|
totped += fabs(calibCalPed->calped[ns][0][ss]); |
4267 |
|
} |
4268 |
|
if ( totped < 1. ){ |
4269 |
|
if ( IsDebug() ) printf(" Section %i totped %f - No calibration data! Calorimeter power problems? \n",section,totped); |
4270 |
|
valid = 0; |
4271 |
|
}; |
4272 |
|
// |
4273 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
4274 |
obt = ph->GetOrbitalTime(); |
obt = ph->GetOrbitalTime(); |
4275 |
pkt = ph->GetCounter(); |
pkt = ph->GetCounter(); |
4437 |
tr->SetBranchAddress("Header", &eh); |
tr->SetBranchAddress("Header", &eh); |
4438 |
nevents = tr->GetEntries(); |
nevents = tr->GetEntries(); |
4439 |
// |
// |
4440 |
|
Float_t totpul = 0.; |
4441 |
|
// |
4442 |
if ( nevents > 0 ){ |
if ( nevents > 0 ){ |
4443 |
// |
// |
4444 |
for (UInt_t i=0; i < nevents; i++){ |
for (UInt_t i=0; i < nevents; i++){ |
4448 |
if ( cp1->pstwerr[section] && cp1->unpackError == 0 ){ |
if ( cp1->pstwerr[section] && cp1->unpackError == 0 ){ |
4449 |
valid = 1; |
valid = 1; |
4450 |
if ( cp1->pperror[section] ) valid = 0; |
if ( cp1->pperror[section] ) valid = 0; |
4451 |
|
// |
4452 |
|
// check pulse values for one plane, if all zeros calibration is not valid (calorimeter power problems) [8th data reduction bug, fixed on 25/11/2009 by E.M.] |
4453 |
|
// |
4454 |
|
totpul = 0.; |
4455 |
|
Int_t ns = 0; |
4456 |
|
if ( section == 2 ) ns = 3; |
4457 |
|
if ( section == 3 ) ns = 1; |
4458 |
|
if ( section == 1 ) ns = 2; |
4459 |
|
for (UInt_t ss=0; ss<96; ss++){ |
4460 |
|
totpul += cp1->calpuls[ns][0][ss]; |
4461 |
|
} |
4462 |
|
if ( totpul >= 3145632. ){ |
4463 |
|
if ( IsDebug() ) printf(" PULSE1 Section %i totpul %f - No calibration data! Calorimeter power problems? \n",section,totpul); |
4464 |
|
valid = 0; |
4465 |
|
}; |
4466 |
|
// |
4467 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
4468 |
obt = ph->GetOrbitalTime(); |
obt = ph->GetOrbitalTime(); |
4469 |
pkt = ph->GetCounter(); |
pkt = ph->GetCounter(); |
4621 |
if ( cp2->pstwerr[section] && cp2->unpackError == 0 ){ |
if ( cp2->pstwerr[section] && cp2->unpackError == 0 ){ |
4622 |
valid = 1; |
valid = 1; |
4623 |
if ( cp2->pperror[section] ) valid = 0; |
if ( cp2->pperror[section] ) valid = 0; |
4624 |
|
// |
4625 |
|
// check pulse values for one plane, if all zeros calibration is not valid (calorimeter power problems) [8th data reduction bug, fixed on 25/11/2009 by E.M.] |
4626 |
|
// |
4627 |
|
totpul = 0.; |
4628 |
|
Int_t ns = 0; |
4629 |
|
if ( section == 2 ) ns = 3; |
4630 |
|
if ( section == 3 ) ns = 1; |
4631 |
|
if ( section == 1 ) ns = 2; |
4632 |
|
for (UInt_t ss=0; ss<96; ss++){ |
4633 |
|
totpul += cp2->calpuls[ns][0][ss]; |
4634 |
|
} |
4635 |
|
if ( totpul >= 3145632. ){ |
4636 |
|
if ( IsDebug() ) printf(" PULSE2 Section %i totpul %f - No calibration data! Calorimeter power problems? \n",section,totpul); |
4637 |
|
valid = 0; |
4638 |
|
}; |
4639 |
|
// |
4640 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
4641 |
obt = ph->GetOrbitalTime(); |
obt = ph->GetOrbitalTime(); |
4642 |
pkt = ph->GetCounter(); |
pkt = ph->GetCounter(); |