/[PAMELA software]/DarthVader/CalorimeterLevel2/src/CaloLevel0.cpp
ViewVC logotype

Diff of /DarthVader/CalorimeterLevel2/src/CaloLevel0.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.28 by mocchiut, Mon Dec 14 12:50:34 2009 UTC revision 1.30 by mocchiut, Fri Jan 29 05:49:24 2010 UTC
# Line 741  Int_t CaloLevel0::CalcCrossTalkCorrDo(TS Line 741  Int_t CaloLevel0::CalcCrossTalkCorrDo(TS
741        //        //
742        if ( !ncalibs ) return(-110);        if ( !ncalibs ) return(-110);
743        //        //
744        calo->GetEntry(calibno);        if ( calo->GetEntry(calibno) <= 0) throw -36;
745        if ( verbose ) printf(" PULSE2 using entry %u from file %s",calibno,fcalname.Data());        if ( verbose ) printf(" PULSE2 using entry %u from file %s",calibno,fcalname.Data());
746        //        //
747        // retrieve calibration table        // retrieve calibration table
# Line 822  Int_t CaloLevel0::CalcCrossTalkCorrDo(TS Line 822  Int_t CaloLevel0::CalcCrossTalkCorrDo(TS
822        //        //
823        if ( !ncalibs ) return(-110);        if ( !ncalibs ) return(-110);
824        //        //
825        calo1->GetEntry(calibno);        if ( calo1->GetEntry(calibno) <= 0 ) throw -36;
826        if ( verbose ) printf(" PULSE1 using entry %u from file %s",calibno,fcalname.Data());        if ( verbose ) printf(" PULSE1 using entry %u from file %s",calibno,fcalname.Data());
827        //        //
828        // retrieve calibration table        // retrieve calibration table
# Line 1184  Int_t CaloLevel0::Calibrate(Int_t ei){ Line 1184  Int_t CaloLevel0::Calibrate(Int_t ei){
1184    //    //
1185    // get entry ei    // get entry ei
1186    //    //
1187    l0calo->GetEntry(ei);    if ( l0calo->GetEntry(ei) <= 0 ) throw -36;
1188    //    //
1189    // if it was not a selftrigger event, could it ever been a selftrigger event? if so trigty = 3.    // if it was not a selftrigger event, could it ever been a selftrigger event? if so trigty = 3.
1190    //    //
# Line 1438  Int_t CaloLevel0::Calibrate(Int_t ei){ Line 1438  Int_t CaloLevel0::Calibrate(Int_t ei){
1438              ener = dexyc[l][m][n];              ener = dexyc[l][m][n];
1439              ener0 += ener;              ener0 += ener;
1440              clevel1->estrip[n][m][l] = 0.;              clevel1->estrip[n][m][l] = 0.;
1441              totbase += de->base[l][m][pre];              totbase += de->base[l][m][pre]/96.;
1442              totped += fabs(calped[l][m][n]);              totped += fabs(calped[l][m][n]);
1443              if ( de->base[l][m][pre] < 0 ) negbase = true;              if ( de->base[l][m][pre] < 0 ) negbase = true;
1444              if ( base0>0 && base0 < 30000. ){              if ( base0>0 && base0 < 30000. ){
# Line 1464  Int_t CaloLevel0::Calibrate(Int_t ei){ Line 1464  Int_t CaloLevel0::Calibrate(Int_t ei){
1464          //          //
1465          // check if there were problems with 5.7 or glitches in the power supply          // check if there were problems with 5.7 or glitches in the power supply
1466          //          //
1467          if ( ((ener0 == 0. && cbase0 == 0.) || negbase || totbase > 196600. || totped < 1. ) && !pproblem && clevel2->perr[se] == 0){   // check pedestal and baseline 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.]          //      if ( ((ener0 == 0. && cbase0 == 0.) || negbase || totbase > 196600. || totped < 1. ) && !pproblem && clevel2->perr[se] == 0){   // check pedestal and baseline 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.]
1468            if ( verbose ) printf(" L0 entry %i : calorimeter power problems! event marked as bad perr %f swerr %X view %i plane %i \n",ei,de->perror[se],de->stwerr[se],l,m);          if ( ((ener0 == 0. && cbase0 == 0.) || negbase || totbase > 32700. || totped < 1. ) && !pproblem && clevel2->perr[se] == 0){    // check pedestal and baseline 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.]
1469              if ( verbose ) printf(" L0 entry %i : calorimeter power problems! event marked as bad perr %f swerr %X view %i plane %i negbase %i totbase %f totped %f\n",ei,de->perror[se],de->stwerr[se],l,m, negbase, totbase, totped);
1470            pproblem = true;            pproblem = true;
1471            pe++;            pe++;
1472          };          };
# Line 2195  Int_t CaloLevel0::LoadCalib(Int_t s){ Line 2196  Int_t CaloLevel0::LoadCalib(Int_t s){
2196    //    //
2197    if ( !ncalibs ) return(-110);    if ( !ncalibs ) return(-110);
2198    //    //
2199    calo->GetEntry(calibno[s]);    if ( calo->GetEntry(calibno[s]) <= 0 ) throw -36;
2200    //    //
2201    if (ce->cstwerr[s] != 0 && ce->cperror[s] == 0 ) {    if (ce->cstwerr[s] != 0 && ce->cperror[s] == 0 ) {
2202      for ( Int_t d=0 ; d<11 ;d++  ){      for ( Int_t d=0 ; d<11 ;d++  ){
Legend:
Removed from v.1.28  
changed lines
  Added in v.1.30
  ViewVC Help
Powered by ViewVC 1.1.23