/[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.4 by mocchiut, Wed Sep 6 11:03:31 2006 UTC revision 1.10 by mocchiut, Wed Jan 17 09:58:17 2007 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 91  void CaloProcessing::ProcessingInit(TSQL Line 94  void CaloProcessing::ProcessingInit(TSQL
94    GL_CALO_CALIB *glcalo = new GL_CALO_CALIB();    GL_CALO_CALIB *glcalo = new GL_CALO_CALIB();
95    //    //
96    sgnl = 0;    sgnl = 0;
97      UInt_t uptime = 0;
98    //    //
99    for (Int_t s = 0; s < 4; s++){    for (Int_t s = 0; s < 4; s++){
100      idcalib[s] = 0;      idcalib[s] = 0;
# Line 99  void CaloProcessing::ProcessingInit(TSQL Line 103  void CaloProcessing::ProcessingInit(TSQL
103      calibno[s] = 0;      calibno[s] = 0;
104      ClearCalibVals(s);      ClearCalibVals(s);
105      //      //
106      sgnl = glcalo->Query_GL_CALO_CALIB(hs,s,dbc);      sgnl = glcalo->Query_GL_CALO_CALIB(hs,uptime,s,dbc);
107      if ( sgnl < 0 ){      if ( sgnl < 0 ){
108          if ( verbose ) printf(" CALORIMETER - ERROR: error from GLTables\n");          if ( verbose ) printf(" CALORIMETER - ERROR: error from GLTables\n");
109          return;          return;
# Line 107  void CaloProcessing::ProcessingInit(TSQL Line 111  void CaloProcessing::ProcessingInit(TSQL
111      //        //  
112      idcalib[s] = glcalo->ID_ROOT_L0;      idcalib[s] = glcalo->ID_ROOT_L0;
113      fromtime[s] = glcalo->FROM_TIME;      fromtime[s] = glcalo->FROM_TIME;
114      totime[s] = glcalo->TO_TIME;      if ( glcalo->TO_TIME < hs ){ // calibration is corrupted and we are using the one that preceed the good one
115          totime[s] = uptime;
116        } else {
117          totime[s] = glcalo->TO_TIME;
118        };
119      calibno[s] = glcalo->EV_ROOT;      calibno[s] = glcalo->EV_ROOT;
120      //      //
121      if ( totime[s] == 0 ){      if ( totime[s] == 0 ){
# Line 156  Int_t CaloProcessing::ChkCalib(TSQLServe Line 164  Int_t CaloProcessing::ChkCalib(TSQLServe
164    Int_t sgnl = 0;    Int_t sgnl = 0;
165    for ( Int_t s = 0; s < 4; s++){    for ( Int_t s = 0; s < 4; s++){
166      if ( atime > totime[s] ){      if ( atime > totime[s] ){
167          if ( !dbc->IsConnected() ) throw -116;
168        sgnl = Update(dbc,atime,s);        sgnl = Update(dbc,atime,s);
169        if ( sgnl < 0 ) return(sgnl);        if ( sgnl < 0 ) return(sgnl);
170      };      };
# Line 165  Int_t CaloProcessing::ChkCalib(TSQLServe Line 174  Int_t CaloProcessing::ChkCalib(TSQLServe
174    
175  Int_t CaloProcessing::ChkParam(TSQLServer *dbc, UInt_t runheader){  Int_t CaloProcessing::ChkParam(TSQLServer *dbc, UInt_t runheader){
176    stringstream calfile;    stringstream calfile;
177      stringstream bmfile;
178    stringstream aligfile;    stringstream aligfile;
179    Int_t error = 0;    Int_t error = 0;
180    FILE *f = 0;    FILE *f = 0;
181      ifstream badfile;
182    GL_PARAM *glparam = new GL_PARAM();    GL_PARAM *glparam = new GL_PARAM();
183    //    //
184    if ( calopar1 || ( ttcalopar1 != 0 && ttcalopar1 < runheader ) ){    if ( calopar1 || ( ttcalopar1 != 0 && ttcalopar1 < runheader ) ){
# Line 220  Int_t CaloProcessing::ChkParam(TSQLServe Line 231  Int_t CaloProcessing::ChkParam(TSQLServe
231      ftcalopar2 = glparam->FROM_TIME;      ftcalopar2 = glparam->FROM_TIME;
232      ttcalopar2 = glparam->TO_TIME;      ttcalopar2 = glparam->TO_TIME;
233      //      //
234      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());
235      f = fopen(aligfile.str().c_str(),"rb");      f = fopen(aligfile.str().c_str(),"rb");
236      if ( !f ){      if ( !f ){
237        if ( verbose ) printf(" CALORIMETER - ERROR: no alignement file!\n");        if ( verbose ) printf(" CALORIMETER - ERROR: no alignement file!\n");
# Line 239  Int_t CaloProcessing::ChkParam(TSQLServe Line 250  Int_t CaloProcessing::ChkParam(TSQLServe
250      fclose(f);      fclose(f);
251    };    };
252    //    //
253      // Load offline bad strip mask
254      //
255      if ( calopar3 || ( ttcalopar3 != 0 && ttcalopar3 < runheader ) ){
256        calopar3 = false;
257        //
258        // determine where I can find calorimeter alignment file  
259        //
260        //
261        error = 0;
262        error = glparam->Query_GL_PARAM(runheader,103,dbc);
263        if ( error < 0 ) return(error);
264        //
265        bmfile.str("");
266        bmfile << glparam->PATH.Data() << "/";
267        bmfile << glparam->NAME.Data();
268        ftcalopar3 = glparam->FROM_TIME;
269        ttcalopar3 = glparam->TO_TIME;
270        //
271        if ( verbose ) printf("\n Using bad strip offline mask file: \n %s \n\n",bmfile.str().c_str());
272        badfile.open(bmfile.str().c_str());
273        if ( !badfile ){
274          if ( verbose ) printf(" CALORIMETER - ERROR: no bad strip offline mask file!\n");
275          return(-115);
276        };
277        //
278        Bool_t isdone = false;
279        Int_t bad = 0;
280        Int_t view = 1;
281        Int_t strip = 0;
282        Int_t plane = 21;
283        while ( !isdone ) {
284          badfile >> bad;
285          obadmask[view][plane][strip] = bad;
286          if ( debug && bad ) printf(" SETTING view %i plane %i strip %i  BAD = %i \n",view,plane,strip,bad);
287          strip++;
288          if ( strip > 95 ){
289            strip = 0;
290            plane--;
291            if ( plane < 0 ){
292              plane = 21;
293              view--;
294            };
295            if ( view < 0 ) isdone = true;
296          };
297        };
298        //
299        badfile.close();
300      };
301      //
302    delete glparam;    delete glparam;
303    //    //
304    return(0);    return(0);
# Line 251  void CaloProcessing::FindBaseRaw(Int_t l Line 311  void CaloProcessing::FindBaseRaw(Int_t l
311      Float_t rms = 0.;      Float_t rms = 0.;
312      base[l][m][pre] = 0.;      base[l][m][pre] = 0.;
313      for (Int_t e = pre*16; e < (pre+1)*16 ; e++){      for (Int_t e = pre*16; e < (pre+1)*16 ; e++){
314          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.) {
315              minstrip = dexy[l][m][e]-calped[l][m][e];              minstrip = dexy[l][m][e]-calped[l][m][e];
316              rms = calthr[l][m][pre];              rms = calthr[l][m][pre];
317          };          };
# Line 303  Int_t CaloProcessing::Calibrate(Int_t ei Line 363  Int_t CaloProcessing::Calibrate(Int_t ei
363    Int_t val = 0;    Int_t val = 0;
364    Int_t del = 1100;    Int_t del = 1100;
365    if ( clevel2->trigty != 2. ){    if ( clevel2->trigty != 2. ){
366        Bool_t ck = false;
367      for (Int_t sec = 0; sec < 4; sec++){      for (Int_t sec = 0; sec < 4; sec++){
368        val = (Int_t)de->calselftrig[sec][6];        val = (Int_t)de->calselftrig[sec][6];
369        del = delay(val);        del = delay(val);
370        if ( del < 1100 ){        if ( del < 1100 ){
371            clevel2->wartrig = 0.;      
372          clevel2->trigty = 3.;          clevel2->trigty = 3.;
373            ck = true;
374          break;          break;
375        };        };
376      };      };
377        if ( !ck ) clevel2->wartrig = 100.;      
378      } else {
379        Bool_t ck = false;
380        for (Int_t sec = 0; sec < 4; sec++){
381          val = (Int_t)de->calselftrig[sec][6];
382          del = delay(val);
383          if ( del < 1100 ){
384            clevel2->wartrig = 0.;      
385            ck = true;
386          };
387        };
388        if ( !ck ) clevel2->wartrig = 100.;      
389    };    };
390    //    //
391    Int_t se = 5;    Int_t se = 5;
# Line 334  Int_t CaloProcessing::Calibrate(Int_t ei Line 409  Int_t CaloProcessing::Calibrate(Int_t ei
409    Int_t chdone[4] = {0,0,0,0};    Int_t chdone[4] = {0,0,0,0};
410    Int_t pe = 0;    Int_t pe = 0;
411    //    //
412      Float_t ener0 = 0.;
413      Float_t cbase0 = 0.;
414      Bool_t pproblem = false;
415      //
416      Float_t tim = 0.;
417      Int_t plo = 0;
418      Int_t fbi = 0;
419      Int_t cle = 0;
420      //
421    // run over views and planes    // run over views and planes
422    //    //
423    for (Int_t l = 0; l < 2; l++){    for (Int_t l = 0; l < 2; l++){
# Line 390  Int_t CaloProcessing::Calibrate(Int_t ei Line 474  Int_t CaloProcessing::Calibrate(Int_t ei
474          // run over preamplifiers          // run over preamplifiers
475          //          //
476          pre = -1;          pre = -1;
477            cbase0 = 0.;
478          for (Int_t i = 0; i < 3; i++){          for (Int_t i = 0; i < 3; i++){
479            for (Int_t j = 0; j < 2; j++){            for (Int_t j = 0; j < 2; j++){
480              pre = j + i*2;              pre = j + i*2;
# Line 398  Int_t CaloProcessing::Calibrate(Int_t ei Line 483  Int_t CaloProcessing::Calibrate(Int_t ei
483              //              //
484              if ( !isRAW ) {              if ( !isRAW ) {
485                base[l][m][pre] = de->base[l][m][pre] ;                  base[l][m][pre] = de->base[l][m][pre] ;  
486                  cbase0 += base[l][m][pre];
487              } else {              } else {
488                //                //
489                // if it is a raw event and we haven't checked                // if it is a raw event and we haven't checked
490                // yet, calculate the baseline.                // yet, calculate the baseline.
491                //                //
492                FindBaseRaw(l,m,pre);                FindBaseRaw(l,m,pre);
493                  cbase0 += base[l][m][pre];
494              };              };
495            };            };
496          };          };
# Line 411  Int_t CaloProcessing::Calibrate(Int_t ei Line 498  Int_t CaloProcessing::Calibrate(Int_t ei
498          // run over strips          // run over strips
499          //          //
500          pre = -1;          pre = -1;
501            ener0 = 0.;
502          for (Int_t i = 0 ; i < 3 ; i++){          for (Int_t i = 0 ; i < 3 ; i++){
503            ip[i] = 0;            ip[i] = 0;
504            for (Int_t n = i*32 ; n < (i+1)*32 ; n++){                            for (Int_t n = i*32 ; n < (i+1)*32 ; n++){                
# Line 471  Int_t CaloProcessing::Calibrate(Int_t ei Line 559  Int_t CaloProcessing::Calibrate(Int_t ei
559                doneb = 1;                doneb = 1;
560              };              };
561              ener = dexyc[l][m][n];              ener = dexyc[l][m][n];
562                ener0 += ener;
563              clevel1->estrip[n][m][l] = 0.;              clevel1->estrip[n][m][l] = 0.;
564              if ( base0>0 && base0 < 30000. ){              if ( base0>0 && base0 < 30000. ){
565                if ( !donec && (base0 - base1 + base2) != 0. ){                if ( !donec && (base0 - base1 + base2) != 0. ){
# Line 502  Int_t CaloProcessing::Calibrate(Int_t ei Line 591  Int_t CaloProcessing::Calibrate(Int_t ei
591                clevel1->estrip[j][m][l] +=  qpre[ipre] * 0.00478;                clevel1->estrip[j][m][l] +=  qpre[ipre] * 0.00478;
592              };              };
593            };            };
594          };          };
595            //
596            if ( ener0 == 0. && cbase0 == 0. && !pproblem ){
597              if ( verbose ) printf(" Calorimeter power problems! event marked as bad \n");
598              pproblem = true;
599              pe++;
600            };
601          //          //
602          Int_t j4 = -4;          Int_t j4 = -4;
603          Int_t jjj = -3;          Int_t jjj = -3;
# Line 522  Int_t CaloProcessing::Calibrate(Int_t ei Line 617  Int_t CaloProcessing::Calibrate(Int_t ei
617            };            };
618            if ( j4 >= 0 && j4 < 96 ){            if ( j4 >= 0 && j4 < 96 ){
619              //              //
620              // 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! =================>
621              //              //
622              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 ){
623                clevel1->estrip[j4][m][l] = 0.;                clevel1->estrip[j4][m][l] = 0.;
624              };              };
625                //
626                // code and save the energy for each strip in svstrip
627                //
628              if ( clevel1->estrip[j4][m][l] > clevel1->emin ){              if ( clevel1->estrip[j4][m][l] > clevel1->emin ){
629                  //
630                  tim = 100000.;
631                  plo = m;
632                  fbi = 0;
633                  if ( clevel1->estrip[j4][m][l] > 0.99995 ){
634                    tim = 10000.;
635                    plo = m;
636                    fbi = 1;
637                  };
638                  if ( clevel1->estrip[j4][m][l] > 9.9995 ){
639                    tim = 1000.;
640                    plo = 22 + m;
641                    fbi = 1;
642                  };
643                  if ( clevel1->estrip[j4][m][l] > 99.995 ){
644                    tim = 100.;
645                    plo = 22 + m;
646                    fbi = 0;
647                  };
648                  if ( clevel1->estrip[j4][m][l] > 999.95 ){
649                    tim = 10.;
650                    plo = 44 + m;
651                    fbi = 0;
652                  };
653                  if ( clevel1->estrip[j4][m][l] > 9999.5 ){
654                    tim = 1.;
655                    plo = 66 + m;
656                    fbi = 0;
657                  };
658                  //
659                  cle = (Int_t)lroundf(tim*clevel1->estrip[j4][m][l]);
660                  //
661                if ( l == 0 ){                if ( l == 0 ){
662                  //                  //
663                  // +-PPSSmmmm.mmmm                  // +-PPSSmmmm.mmmm
664                  //                  //
665                  svstrip[istrip] = ((Float_t)m)*1000000. + ((Float_t)j4)*10000. + clevel1->estrip[j4][m][l];                  svstrip[istrip] = fbi*1000000000 + plo*10000000 + j4*100000 + cle;
666                } else {                } else {
667                  svstrip[istrip] = -(((Float_t)m)*1000000. + ((Float_t)j4)*10000. + clevel1->estrip[j4][m][l]);                  svstrip[istrip] = -(fbi*1000000000 + plo*10000000 + j4*100000 + cle);
668                };                };
669                  //
670                  //              if ( ei >= -770 ) printf(" j %i l %i m %i estrip %f \n",j4,l,m,clevel1->estrip[j4][m][l]);
671                  //              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);
672                  //              if ( ei >= -770 ) printf(" svstrip %i \n",svstrip[istrip]);
673                  //
674                istrip++;                istrip++;
675              };              };
676            };            };
# Line 606  void CaloProcessing::FillTrkVar(CaloLeve Line 741  void CaloProcessing::FillTrkVar(CaloLeve
741      memcpy(ca->npcfit,clevel2->npcfit,sizeof(clevel2->npcfit));      memcpy(ca->npcfit,clevel2->npcfit,sizeof(clevel2->npcfit));
742    };    };
743    //    //
744      if(!(ca->CaloTrk))ca->CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA
745    TClonesArray &t = *ca->CaloTrk;    TClonesArray &t = *ca->CaloTrk;
746    new(t[nutrk]) CaloTrkVar(*t_ca);    new(t[nutrk]) CaloTrkVar(*t_ca);
747    //    //
# Line 618  void CaloProcessing::GetCommonVar(){ Line 754  void CaloProcessing::GetCommonVar(){
754    calol2cm();    calol2cm();
755  }  }
756    
757  void CaloProcessing::FillCommonVar(CaloLevel2 *ca){  void CaloProcessing::FillCommonVar(CaloLevel1 *c1, CaloLevel2 *ca){
758    //    //
759    ca->good = clevel2->good;    ca->good = clevel2->good;
760    if ( clevel2->trigty == 2. ){    if ( clevel2->trigty == 2. ){
# Line 626  void CaloProcessing::FillCommonVar(CaloL Line 762  void CaloProcessing::FillCommonVar(CaloL
762    } else {    } else {
763      ca->selftrigger = 0;      ca->selftrigger = 0;
764    };    };
765      //
766      ca->selftrigger += (Int_t)clevel2->wartrig;
767      //
768    memcpy(ca->perr,clevel2->perr,sizeof(clevel2->perr));    memcpy(ca->perr,clevel2->perr,sizeof(clevel2->perr));
769    memcpy(ca->swerr,clevel2->swerr,sizeof(clevel2->swerr));    memcpy(ca->swerr,clevel2->swerr,sizeof(clevel2->swerr));
770    memcpy(ca->crc,clevel2->crc,sizeof(clevel2->crc));    memcpy(ca->crc,clevel2->crc,sizeof(clevel2->crc));
   ca->estrip = TArrayF(0,svstrip);  
771    ca->nstrip = (Int_t)clevel2->nstrip;    ca->nstrip = (Int_t)clevel2->nstrip;
772    ca->qtot = clevel2->qtot;    ca->qtot = clevel2->qtot;
773    ca->impx = clevel2->impx;    ca->impx = clevel2->impx;
# Line 641  void CaloProcessing::FillCommonVar(CaloL Line 779  void CaloProcessing::FillCommonVar(CaloL
779    ca->qmax = clevel2->qmax;    ca->qmax = clevel2->qmax;
780    ca->elen = clevel2->elen;    ca->elen = clevel2->elen;
781    ca->selen = clevel2->selen;    ca->selen = clevel2->selen;
   ca->estrip = TArrayF(ca->nstrip,svstrip);  
782    memcpy(ca->qq,clevel2->qq,sizeof(clevel2->qq));    memcpy(ca->qq,clevel2->qq,sizeof(clevel2->qq));
783    memcpy(ca->planemax,clevel2->planemax,sizeof(clevel2->planemax));    memcpy(ca->planemax,clevel2->planemax,sizeof(clevel2->planemax));
784    memcpy(ca->varcfit,clevel2->varcfit,sizeof(clevel2->varcfit));    memcpy(ca->varcfit,clevel2->varcfit,sizeof(clevel2->varcfit));
# Line 649  void CaloProcessing::FillCommonVar(CaloL Line 786  void CaloProcessing::FillCommonVar(CaloL
786    memcpy(ca->cibar,clevel2->cibar,sizeof(clevel2->cibar));    memcpy(ca->cibar,clevel2->cibar,sizeof(clevel2->cibar));
787    memcpy(ca->cbar,clevel2->cbar,sizeof(clevel2->cbar));    memcpy(ca->cbar,clevel2->cbar,sizeof(clevel2->cbar));
788    //    //
789      if ( c1 ){
790        c1->istrip = istrip;
791        c1->estrip = TArrayI(istrip,svstrip);
792      };
793      //
794  }  }
795    
796  void CaloProcessing::ClearStructs(){  void CaloProcessing::ClearStructs(){
# Line 705  void CaloProcessing::ClearTrkVar(){ Line 847  void CaloProcessing::ClearTrkVar(){
847  void CaloProcessing::ClearCommonVar(){  void CaloProcessing::ClearCommonVar(){
848    istrip = 0;    istrip = 0;
849    clevel2->trigty = -1.;    clevel2->trigty = -1.;
850      clevel2->wartrig = 0.;
851    clevel2->good = 0;    clevel2->good = 0;
852    clevel2->nstrip = 0.;    clevel2->nstrip = 0.;
853    clevel2->qtot = 0.;    clevel2->qtot = 0.;
# Line 791  Int_t CaloProcessing::Update(TSQLServer Line 934  Int_t CaloProcessing::Update(TSQLServer
934    calibno[s] = 0;    calibno[s] = 0;
935    ClearCalibVals(s);    ClearCalibVals(s);
936    //    //
937    sgnl = glcalo->Query_GL_CALO_CALIB(atime,s,dbc);    UInt_t uptime = 0;
938      //
939      sgnl = glcalo->Query_GL_CALO_CALIB(atime,uptime,s,dbc);
940    if ( sgnl < 0 ){    if ( sgnl < 0 ){
941      if ( verbose ) printf(" CALORIMETER - ERROR: error from GLTables\n");      if ( verbose ) printf(" CALORIMETER - ERROR: error from GLTables\n");
942      return(sgnl);      return(sgnl);
# Line 799  Int_t CaloProcessing::Update(TSQLServer Line 944  Int_t CaloProcessing::Update(TSQLServer
944    //      //  
945    idcalib[s] = glcalo->ID_ROOT_L0;    idcalib[s] = glcalo->ID_ROOT_L0;
946    fromtime[s] = glcalo->FROM_TIME;    fromtime[s] = glcalo->FROM_TIME;
947    totime[s] = glcalo->TO_TIME;    if ( glcalo->TO_TIME < atime ){ // calibration is corrupted and we are using the one that preceed the good one
948        totime[s] = uptime;
949      } else {
950        totime[s] = glcalo->TO_TIME;
951      };
952      //  totime[s] = glcalo->TO_TIME;
953    calibno[s] = glcalo->EV_ROOT;    calibno[s] = glcalo->EV_ROOT;
954    //    //
955    if ( totime[s] == 0 ){    if ( totime[s] == 0 ){
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.10
  ViewVC Help
Powered by ViewVC 1.1.23