/[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.8 by mocchiut, Wed Nov 22 08:31:11 2006 UTC revision 1.11 by mocchiut, Mon Jan 22 09:16:59 2007 UTC
# Line 363  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 578  Int_t CaloProcessing::Calibrate(Int_t ei Line 593  Int_t CaloProcessing::Calibrate(Int_t ei
593            };            };
594          };          };
595          //          //
596          if ( ener0 == 0. && cbase0 == 0. && !pproblem ){          if ( ener0 == 0. && cbase0 == 0. && !pproblem && clevel2->perr[se] == 0){
597            if ( verbose ) printf(" Calorimeter power problems! event marked as bad \n");            if ( verbose ) printf(" Entry %i : calorimeter power problems! event marked as bad \n",ei);
598            pproblem = true;            pproblem = true;
599            pe++;            pe++;
600          };          };
# Line 709  void CaloProcessing::FillTrkVar(CaloLeve Line 724  void CaloProcessing::FillTrkVar(CaloLeve
724    t_ca->qlow = clevel2->qlow;    t_ca->qlow = clevel2->qlow;
725    t_ca->nlow = (Int_t)clevel2->nlow;    t_ca->nlow = (Int_t)clevel2->nlow;
726    //    //
   memcpy(t_ca->tibar,clevel2->tibar,sizeof(clevel2->tibar));  
   memcpy(t_ca->tbar,clevel2->tbar,sizeof(clevel2->tbar));  
   //  
727    if ( trkseqno == -1 ){    if ( trkseqno == -1 ){
728      ca->impx = clevel2->impx;      //    ca->impx = clevel2->impx;
729      ca->impy = clevel2->impy;      //    ca->impy = clevel2->impy;
730      ca->tanx = clevel2->tanx;      ca->tanx[1] = clevel2->tanx;
731      ca->tany = clevel2->tany;      ca->tany[1] = clevel2->tany;
732      ca->elen = clevel2->elen;      ca->elen = clevel2->elen;
733      ca->selen = clevel2->selen;      ca->selen = clevel2->selen;
734      memcpy(ca->cibar,clevel2->cibar,sizeof(clevel2->cibar));      //    memcpy(ca->cibar,clevel2->cibar,sizeof(clevel2->cibar));
735      memcpy(ca->cbar,clevel2->cbar,sizeof(clevel2->cbar));      //    memcpy(ca->cbar,clevel2->cbar,sizeof(clevel2->cbar));
736        memcpy(t_ca->tibar,clevel2->cibar,sizeof(clevel2->cibar));
737        memcpy(t_ca->tbar,clevel2->cbar,sizeof(clevel2->cbar));
738      memcpy(ca->planemax,clevel2->planemax,sizeof(clevel2->planemax));      memcpy(ca->planemax,clevel2->planemax,sizeof(clevel2->planemax));
739      memcpy(ca->varcfit,clevel2->varcfit,sizeof(clevel2->varcfit));      ca->varcfit[2] = clevel2->varcfit[0];
740      memcpy(ca->npcfit,clevel2->npcfit,sizeof(clevel2->npcfit));      ca->varcfit[3] = clevel2->varcfit[1];
741        ca->npcfit[2] = clevel2->npcfit[0];
742        ca->npcfit[3] = clevel2->npcfit[1];    
743        //    memcpy(ca->varcfit,clevel2->varcfit,sizeof(clevel2->varcfit));
744        //    memcpy(ca->npcfit,clevel2->npcfit,sizeof(clevel2->npcfit));
745      } else {
746        memcpy(t_ca->tibar,clevel2->tibar,sizeof(clevel2->tibar));
747        memcpy(t_ca->tbar,clevel2->tbar,sizeof(clevel2->tbar));
748    };    };
749      //
750    //    //
751      if ( !(ca->CaloTrk) ) ca->CaloTrk = new TClonesArray("CaloTrkVar",1); //ELENA
752    TClonesArray &t = *ca->CaloTrk;    TClonesArray &t = *ca->CaloTrk;
753    new(t[nutrk]) CaloTrkVar(*t_ca);    new(t[nutrk]) CaloTrkVar(*t_ca);
754    //    //
# Line 746  void CaloProcessing::FillCommonVar(CaloL Line 769  void CaloProcessing::FillCommonVar(CaloL
769    } else {    } else {
770      ca->selftrigger = 0;      ca->selftrigger = 0;
771    };    };
772      //
773      ca->selftrigger += (Int_t)clevel2->wartrig;
774      //
775    memcpy(ca->perr,clevel2->perr,sizeof(clevel2->perr));    memcpy(ca->perr,clevel2->perr,sizeof(clevel2->perr));
776    memcpy(ca->swerr,clevel2->swerr,sizeof(clevel2->swerr));    memcpy(ca->swerr,clevel2->swerr,sizeof(clevel2->swerr));
777    memcpy(ca->crc,clevel2->crc,sizeof(clevel2->crc));    memcpy(ca->crc,clevel2->crc,sizeof(clevel2->crc));
778    ca->nstrip = (Int_t)clevel2->nstrip;    ca->nstrip = (Int_t)clevel2->nstrip;
779    ca->qtot = clevel2->qtot;    ca->qtot = clevel2->qtot;
780    ca->impx = clevel2->impx;    //  ca->impx = clevel2->impx;
781    ca->impy = clevel2->impy;    //  ca->impy = clevel2->impy;
782    ca->tanx = clevel2->tanx;    ca->tanx[0] = clevel2->tanx;
783    ca->tany = clevel2->tany;    ca->tany[0] = clevel2->tany;
784    ca->nx22 = (Int_t)clevel2->nx22;    ca->nx22 = (Int_t)clevel2->nx22;
785    ca->qx22 = clevel2->qx22;    ca->qx22 = clevel2->qx22;
786    ca->qmax = clevel2->qmax;    ca->qmax = clevel2->qmax;
# Line 762  void CaloProcessing::FillCommonVar(CaloL Line 788  void CaloProcessing::FillCommonVar(CaloL
788    ca->selen = clevel2->selen;    ca->selen = clevel2->selen;
789    memcpy(ca->qq,clevel2->qq,sizeof(clevel2->qq));    memcpy(ca->qq,clevel2->qq,sizeof(clevel2->qq));
790    memcpy(ca->planemax,clevel2->planemax,sizeof(clevel2->planemax));    memcpy(ca->planemax,clevel2->planemax,sizeof(clevel2->planemax));
791    memcpy(ca->varcfit,clevel2->varcfit,sizeof(clevel2->varcfit));    ca->varcfit[0] = clevel2->varcfit[0];
792    memcpy(ca->npcfit,clevel2->npcfit,sizeof(clevel2->npcfit));    ca->varcfit[1] = clevel2->varcfit[1];
793      ca->npcfit[0] = clevel2->npcfit[0];
794      ca->npcfit[1] = clevel2->npcfit[1];
795      ca->fitmode[0] = clevel2->fmode[0];
796      ca->fitmode[1] = clevel2->fmode[1];
797      //  memcpy(ca->varcfit,clevel2->varcfit,sizeof(clevel2->varcfit));
798      //  memcpy(ca->npcfit,clevel2->npcfit,sizeof(clevel2->npcfit));
799    memcpy(ca->cibar,clevel2->cibar,sizeof(clevel2->cibar));    memcpy(ca->cibar,clevel2->cibar,sizeof(clevel2->cibar));
800    memcpy(ca->cbar,clevel2->cbar,sizeof(clevel2->cbar));    memcpy(ca->cbar,clevel2->cbar,sizeof(clevel2->cbar));
801    //    //
# Line 828  void CaloProcessing::ClearTrkVar(){ Line 860  void CaloProcessing::ClearTrkVar(){
860  void CaloProcessing::ClearCommonVar(){  void CaloProcessing::ClearCommonVar(){
861    istrip = 0;    istrip = 0;
862    clevel2->trigty = -1.;    clevel2->trigty = -1.;
863      clevel2->wartrig = 0.;
864    clevel2->good = 0;    clevel2->good = 0;
865    clevel2->nstrip = 0.;    clevel2->nstrip = 0.;
866    clevel2->qtot = 0.;    clevel2->qtot = 0.;
867    clevel2->impx = 0.;    //  clevel2->impx = 0.;
868    clevel2->impy = 0.;    //  clevel2->impy = 0.;
869    clevel2->tanx = 0.;    clevel2->tanx = 0.; // this is correct since it refers to the fortran structure
870    clevel2->tany = 0.;    clevel2->tany = 0.; // this is correct since it refers to the fortran structure
871    clevel2->qmax = 0.;    clevel2->qmax = 0.;
872    clevel2->nx22 = 0.;    clevel2->nx22 = 0.;
873    clevel2->qx22 = 0.;    clevel2->qx22 = 0.;
# Line 842  void CaloProcessing::ClearCommonVar(){ Line 875  void CaloProcessing::ClearCommonVar(){
875    memset(clevel2->swerr, 0, 4*sizeof(Int_t));    memset(clevel2->swerr, 0, 4*sizeof(Int_t));
876    memset(clevel2->crc, 0, 4*sizeof(Int_t));    memset(clevel2->crc, 0, 4*sizeof(Int_t));
877    memset(clevel2->qq, 0, 4*sizeof(Int_t));    memset(clevel2->qq, 0, 4*sizeof(Int_t));
878    memset(clevel2->varcfit, 0, 2*sizeof(Float_t));    memset(clevel2->varcfit, 0, 4*sizeof(Float_t));
879    memset(clevel2->npcfit, 0, 2*sizeof(Int_t));    memset(clevel2->npcfit, 0, 4*sizeof(Int_t));
880    memset(clevel2->planemax, 0, 2*sizeof(Int_t));    memset(clevel2->planemax, 0, 2*sizeof(Int_t));
881      memset(clevel2->fmode, 0, 2*sizeof(Int_t));
882    memset(clevel2->cibar, 0, 2*22*sizeof(Int_t));    memset(clevel2->cibar, 0, 2*22*sizeof(Int_t));
883    memset(clevel2->cbar, 0, 2*22*sizeof(Float_t));    memset(clevel2->cbar, 0, 2*22*sizeof(Float_t));
884  }  }
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.23