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

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

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

revision 1.5 by mocchiut, Thu Sep 7 09:47:07 2006 UTC revision 1.6 by mocchiut, Tue Nov 14 14:08:50 2006 UTC
# Line 66  CaloProcessing::CaloProcessing(){ Line 66  CaloProcessing::CaloProcessing(){
66    memset(sbase, 0, 2*22*6*sizeof(Float_t));    memset(sbase, 0, 2*22*6*sizeof(Float_t));
67    calopar1 = true;    calopar1 = true;
68    calopar2 = true;    calopar2 = true;
69      calopar3 = true;
70    ftcalopar1 = 0;    ftcalopar1 = 0;
71    ttcalopar1 = 0;    ttcalopar1 = 0;
72    ftcalopar2 = 0;    ftcalopar2 = 0;
73    ttcalopar2 = 0;    ttcalopar2 = 0;
74      ftcalopar3 = 0;
75      ttcalopar3 = 0;
76  }  }
77    
78  /**  /**
# Line 170  Int_t CaloProcessing::ChkCalib(TSQLServe Line 173  Int_t CaloProcessing::ChkCalib(TSQLServe
173    
174  Int_t CaloProcessing::ChkParam(TSQLServer *dbc, UInt_t runheader){  Int_t CaloProcessing::ChkParam(TSQLServer *dbc, UInt_t runheader){
175    stringstream calfile;    stringstream calfile;
176      stringstream bmfile;
177    stringstream aligfile;    stringstream aligfile;
178    Int_t error = 0;    Int_t error = 0;
179    FILE *f = 0;    FILE *f = 0;
180      ifstream badfile;
181    GL_PARAM *glparam = new GL_PARAM();    GL_PARAM *glparam = new GL_PARAM();
182    //    //
183    if ( calopar1 || ( ttcalopar1 != 0 && ttcalopar1 < runheader ) ){    if ( calopar1 || ( ttcalopar1 != 0 && ttcalopar1 < runheader ) ){
# Line 225  Int_t CaloProcessing::ChkParam(TSQLServe Line 230  Int_t CaloProcessing::ChkParam(TSQLServe
230      ftcalopar2 = glparam->FROM_TIME;      ftcalopar2 = glparam->FROM_TIME;
231      ttcalopar2 = glparam->TO_TIME;      ttcalopar2 = glparam->TO_TIME;
232      //      //
233      if ( verbose ) printf("\n Using alignment file: \n %s \n\n",aligfile.str().c_str());      if ( verbose ) printf("\n Using alignment file: \n %s \n",aligfile.str().c_str());
234      f = fopen(aligfile.str().c_str(),"rb");      f = fopen(aligfile.str().c_str(),"rb");
235      if ( !f ){      if ( !f ){
236        if ( verbose ) printf(" CALORIMETER - ERROR: no alignement file!\n");        if ( verbose ) printf(" CALORIMETER - ERROR: no alignement file!\n");
# Line 244  Int_t CaloProcessing::ChkParam(TSQLServe Line 249  Int_t CaloProcessing::ChkParam(TSQLServe
249      fclose(f);      fclose(f);
250    };    };
251    //    //
252      // Load offline bad strip mask
253      //
254      if ( calopar3 || ( ttcalopar3 != 0 && ttcalopar3 < runheader ) ){
255        calopar3 = false;
256        //
257        // determine where I can find calorimeter alignment file  
258        //
259        //
260        error = 0;
261        error = glparam->Query_GL_PARAM(runheader,103,dbc);
262        if ( error < 0 ) return(error);
263        //
264        bmfile.str("");
265        bmfile << glparam->PATH.Data() << "/";
266        bmfile << glparam->NAME.Data();
267        ftcalopar3 = glparam->FROM_TIME;
268        ttcalopar3 = glparam->TO_TIME;
269        //
270        if ( verbose ) printf("\n Using bad strip offline mask file: \n %s \n\n",bmfile.str().c_str());
271        badfile.open(bmfile.str().c_str());
272        if ( !badfile ){
273          if ( verbose ) printf(" CALORIMETER - ERROR: no bad strip offline mask file!\n");
274          return(-115);
275        };
276        //
277        Bool_t isdone = false;
278        Int_t bad = 0;
279        Int_t view = 1;
280        Int_t strip = 0;
281        Int_t plane = 21;
282        while ( !isdone ) {
283          badfile >> bad;
284          obadmask[view][plane][strip] = bad;
285          if ( debug && bad ) printf(" SETTING view %i plane %i strip %i  BAD = %i \n",view,plane,strip,bad);
286          strip++;
287          if ( strip > 95 ){
288            strip = 0;
289            plane--;
290            if ( plane < 0 ){
291              plane = 21;
292              view--;
293            };
294            if ( view < 0 ) isdone = true;
295          };
296        };
297        //
298        badfile.close();
299      };
300      //
301    delete glparam;    delete glparam;
302    //    //
303    return(0);    return(0);
# Line 256  void CaloProcessing::FindBaseRaw(Int_t l Line 310  void CaloProcessing::FindBaseRaw(Int_t l
310      Float_t rms = 0.;      Float_t rms = 0.;
311      base[l][m][pre] = 0.;      base[l][m][pre] = 0.;
312      for (Int_t e = pre*16; e < (pre+1)*16 ; e++){      for (Int_t e = pre*16; e < (pre+1)*16 ; e++){
313          if ( calgood[l][m][e] == 0. &&  dexy[l][m][e]-calped[l][m][e] < minstrip &&  dexy[l][m][e] > 0.) {          if ( calgood[l][m][e] == 0. && obadmask[l][m][e] == 0 &&  dexy[l][m][e]-calped[l][m][e] < minstrip &&  dexy[l][m][e] > 0.) {
314              minstrip = dexy[l][m][e]-calped[l][m][e];              minstrip = dexy[l][m][e]-calped[l][m][e];
315              rms = calthr[l][m][pre];              rms = calthr[l][m][pre];
316          };          };
# Line 339  Int_t CaloProcessing::Calibrate(Int_t ei Line 393  Int_t CaloProcessing::Calibrate(Int_t ei
393    Int_t chdone[4] = {0,0,0,0};    Int_t chdone[4] = {0,0,0,0};
394    Int_t pe = 0;    Int_t pe = 0;
395    //    //
396      Float_t ener0 = 0.;
397      Float_t cbase0 = 0.;
398      Bool_t pproblem = false;
399      //
400      Float_t tim = 0.;
401      Int_t plo = 0;
402      Int_t fbi = 0;
403      Int_t cle = 0;
404      //
405    // run over views and planes    // run over views and planes
406    //    //
407    for (Int_t l = 0; l < 2; l++){    for (Int_t l = 0; l < 2; l++){
# Line 395  Int_t CaloProcessing::Calibrate(Int_t ei Line 458  Int_t CaloProcessing::Calibrate(Int_t ei
458          // run over preamplifiers          // run over preamplifiers
459          //          //
460          pre = -1;          pre = -1;
461            cbase0 = 0.;
462          for (Int_t i = 0; i < 3; i++){          for (Int_t i = 0; i < 3; i++){
463            for (Int_t j = 0; j < 2; j++){            for (Int_t j = 0; j < 2; j++){
464              pre = j + i*2;              pre = j + i*2;
# Line 403  Int_t CaloProcessing::Calibrate(Int_t ei Line 467  Int_t CaloProcessing::Calibrate(Int_t ei
467              //              //
468              if ( !isRAW ) {              if ( !isRAW ) {
469                base[l][m][pre] = de->base[l][m][pre] ;                  base[l][m][pre] = de->base[l][m][pre] ;  
470                  cbase0 += base[l][m][pre];
471              } else {              } else {
472                //                //
473                // if it is a raw event and we haven't checked                // if it is a raw event and we haven't checked
474                // yet, calculate the baseline.                // yet, calculate the baseline.
475                //                //
476                FindBaseRaw(l,m,pre);                FindBaseRaw(l,m,pre);
477                  cbase0 += base[l][m][pre];
478              };              };
479            };            };
480          };          };
# Line 416  Int_t CaloProcessing::Calibrate(Int_t ei Line 482  Int_t CaloProcessing::Calibrate(Int_t ei
482          // run over strips          // run over strips
483          //          //
484          pre = -1;          pre = -1;
485            ener0 = 0.;
486          for (Int_t i = 0 ; i < 3 ; i++){          for (Int_t i = 0 ; i < 3 ; i++){
487            ip[i] = 0;            ip[i] = 0;
488            for (Int_t n = i*32 ; n < (i+1)*32 ; n++){                            for (Int_t n = i*32 ; n < (i+1)*32 ; n++){                
# Line 476  Int_t CaloProcessing::Calibrate(Int_t ei Line 543  Int_t CaloProcessing::Calibrate(Int_t ei
543                doneb = 1;                doneb = 1;
544              };              };
545              ener = dexyc[l][m][n];              ener = dexyc[l][m][n];
546                ener0 += ener;
547              clevel1->estrip[n][m][l] = 0.;              clevel1->estrip[n][m][l] = 0.;
548              if ( base0>0 && base0 < 30000. ){              if ( base0>0 && base0 < 30000. ){
549                if ( !donec && (base0 - base1 + base2) != 0. ){                if ( !donec && (base0 - base1 + base2) != 0. ){
# Line 507  Int_t CaloProcessing::Calibrate(Int_t ei Line 575  Int_t CaloProcessing::Calibrate(Int_t ei
575                clevel1->estrip[j][m][l] +=  qpre[ipre] * 0.00478;                clevel1->estrip[j][m][l] +=  qpre[ipre] * 0.00478;
576              };              };
577            };            };
578          };          };
579            //
580            if ( ener0 == 0. && cbase0 == 0. && !pproblem ){
581              if ( verbose ) printf(" Calorimeter power problems! event marked as bad \n");
582              pproblem = true;
583              pe++;
584            };
585          //          //
586          Int_t j4 = -4;          Int_t j4 = -4;
587          Int_t jjj = -3;          Int_t jjj = -3;
# Line 527  Int_t CaloProcessing::Calibrate(Int_t ei Line 601  Int_t CaloProcessing::Calibrate(Int_t ei
601            };            };
602            if ( j4 >= 0 && j4 < 96 ){            if ( j4 >= 0 && j4 < 96 ){
603              //              //
604              // NOTICE: THE FOLLOWING LINE EXCLUDE ALL STRIPS FOR WHICH THE RMS*4 IS GREATER THAN 26 !!! <===V============ IMPORTANT! =================>              // NOTICE: THE FOLLOWING LINE EXCLUDE ALL STRIPS FOR WHICH THE RMS*4 IS GREATER THAN 26 !!! <=============== IMPORTANT! =================>
605              //              //
606              if ( clevel1->estrip[j4][m][l]!=0. && ( clevel1->estrip[j4][m][l] < clevel1->emin || calrms[l][m][j4] > 26 )){              if ( obadmask[l][m][j4] == 1 || clevel1->estrip[j4][m][l] <= clevel1->emin || calrms[l][m][j4] > 26 ){
607                clevel1->estrip[j4][m][l] = 0.;                clevel1->estrip[j4][m][l] = 0.;
608              };              };
609                //
610                // code and save the energy for each strip in svstrip
611                //
612              if ( clevel1->estrip[j4][m][l] > clevel1->emin ){              if ( clevel1->estrip[j4][m][l] > clevel1->emin ){
613                  //
614                  tim = 100000.;
615                  plo = m;
616                  fbi = 0;
617                  if ( clevel1->estrip[j4][m][l] > 1. ){
618                    tim = 10000.;
619                    plo = m;
620                    fbi = 1;
621                  };
622                  if ( clevel1->estrip[j4][m][l] > 10. ){
623                    tim = 1000.;
624                    plo = 22 + m;
625                    fbi = 1;
626                  };
627                  if ( clevel1->estrip[j4][m][l] > 100. ){
628                    tim = 100.;
629                    plo = 22 + m;
630                    fbi = 0;
631                  };
632                  if ( clevel1->estrip[j4][m][l] > 1000. ){
633                    tim = 10.;
634                    plo = 44 + m;
635                    fbi = 0;
636                  };
637                  if ( clevel1->estrip[j4][m][l] > 10000. ){
638                    tim = 1.;
639                    plo = 66 + m;
640                    fbi = 0;
641                  };
642                  //
643                  cle = (Int_t)lroundf(tim*clevel1->estrip[j4][m][l]);
644                  //
645                if ( l == 0 ){                if ( l == 0 ){
646                  //                  //
647                  // +-PPSSmmmm.mmmm                  // +-PPSSmmmm.mmmm
648                  //                  //
649                  svstrip[istrip] = ((Float_t)m)*1000000. + ((Float_t)j4)*10000. + clevel1->estrip[j4][m][l];                  svstrip[istrip] = fbi*1000000000 + plo*10000000 + j4*100000 + cle;
650                } else {                } else {
651                  svstrip[istrip] = -(((Float_t)m)*1000000. + ((Float_t)j4)*10000. + clevel1->estrip[j4][m][l]);                  svstrip[istrip] = -(fbi*1000000000 + plo*10000000 + j4*100000 + cle);
652                };                };
653                  //              if ( ei == 770 ) printf(" j %i l %i m %i estrip %f \n",j4,l,m,clevel1->estrip[j4][m][l]);
654                  //              if ( ei == 770 ) printf(" num lim %i fbi %i tim %f plo %i cle %i \n",numeric_limits<Int_t>::max(),fbi,tim,plo,cle);
655                  //              if ( ei == 770 ) printf(" svstrip %i \n",svstrip[istrip]);
656                  //
657                istrip++;                istrip++;
658              };              };
659            };            };
# Line 634  void CaloProcessing::FillCommonVar(CaloL Line 747  void CaloProcessing::FillCommonVar(CaloL
747    memcpy(ca->perr,clevel2->perr,sizeof(clevel2->perr));    memcpy(ca->perr,clevel2->perr,sizeof(clevel2->perr));
748    memcpy(ca->swerr,clevel2->swerr,sizeof(clevel2->swerr));    memcpy(ca->swerr,clevel2->swerr,sizeof(clevel2->swerr));
749    memcpy(ca->crc,clevel2->crc,sizeof(clevel2->crc));    memcpy(ca->crc,clevel2->crc,sizeof(clevel2->crc));
   ca->estrip = TArrayF(0,svstrip);  
750    ca->nstrip = (Int_t)clevel2->nstrip;    ca->nstrip = (Int_t)clevel2->nstrip;
751    ca->qtot = clevel2->qtot;    ca->qtot = clevel2->qtot;
752    ca->impx = clevel2->impx;    ca->impx = clevel2->impx;
# Line 646  void CaloProcessing::FillCommonVar(CaloL Line 758  void CaloProcessing::FillCommonVar(CaloL
758    ca->qmax = clevel2->qmax;    ca->qmax = clevel2->qmax;
759    ca->elen = clevel2->elen;    ca->elen = clevel2->elen;
760    ca->selen = clevel2->selen;    ca->selen = clevel2->selen;
761    ca->estrip = TArrayF(ca->nstrip,svstrip);    ca->estrip = TArrayI(ca->nstrip,svstrip);
762    memcpy(ca->qq,clevel2->qq,sizeof(clevel2->qq));    memcpy(ca->qq,clevel2->qq,sizeof(clevel2->qq));
763    memcpy(ca->planemax,clevel2->planemax,sizeof(clevel2->planemax));    memcpy(ca->planemax,clevel2->planemax,sizeof(clevel2->planemax));
764    memcpy(ca->varcfit,clevel2->varcfit,sizeof(clevel2->varcfit));    memcpy(ca->varcfit,clevel2->varcfit,sizeof(clevel2->varcfit));
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.23