/[PAMELA software]/calo/flight/CaloEnergy/src/CaloEnergy.cpp
ViewVC logotype

Diff of /calo/flight/CaloEnergy/src/CaloEnergy.cpp

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

revision 1.5 by mocchiut, Tue Aug 11 13:28:27 2009 UTC revision 1.8 by mocchiut, Wed Aug 12 14:54:52 2009 UTC
# Line 46  void CaloEnergy::UseCaloPreSampler(){ Line 46  void CaloEnergy::UseCaloPreSampler(){
46  //   cp->ForceCaloFit();  //   cp->ForceCaloFit();
47  //   cp->SetDebug(true);  //   cp->SetDebug(true);
48  //   cp->Process();  //   cp->Process();
49      if ( clong ) clong->SetCaloLevel2Pointer(cp->GetLevel2Pointer());
50  }  }
51    
52    
# Line 538  Bool_t CaloEnergy::IsInsideAcceptance(TS Line 539  Bool_t CaloEnergy::IsInsideAcceptance(TS
539      Int_t wpl = (Int_t)roundf(x0max/0.76);      Int_t wpl = (Int_t)roundf(x0max/0.76);
540      Bool_t isY = false;      Bool_t isY = false;
541      if ( ((x0max/0.76)-(Float_t)wpl) > 0. ) isY = true;      if ( ((x0max/0.76)-(Float_t)wpl) > 0. ) isY = true;
542        xomax_en = 0.;
543        yemax_en = 0.;
544        xemax_en = 0.;
545        yomax_en = 0.;
546        //
547      if ( !(wpl%2) ){      if ( !(wpl%2) ){
548        // 0, 2, 4, ...        // 0, 2, 4, ...
549        if ( isY ){        if ( isY ){
550          yomax_en = 1000.;          if ( section.Contains("YO") ) yomax_en = 1000.;
551          xemax_en = 500.;          if ( section.Contains("XE") ) xemax_en = 500.;
552          fMax_planeyo=wpl/2;          fMax_planeyo=wpl/2;
553          fMax_planexe=wpl/2;          fMax_planexe=wpl/2;
554            if ( section.Contains("XO") ) xomax_en = 10.;
555            if ( section.Contains("YE") ) yemax_en = 5.;
556        } else {        } else {
557          yomax_en = 500.;          if ( section.Contains("YO") ) yomax_en = 500.;
558          xemax_en = 1000.;          if ( section.Contains("XE") ) xemax_en = 1000.;
559          fMax_planeyo=wpl/2;          fMax_planeyo=wpl/2;
560          fMax_planexe=wpl/2;          fMax_planexe=wpl/2;
561            if ( section.Contains("XO") ) xomax_en = 5.;
562            if ( section.Contains("YE") ) yemax_en = 10.;
563        };        };
       xomax_en = 0.;  
       yemax_en = 0.;  
564      } else {      } else {
565        // 1, 3, 5, ...        // 1, 3, 5, ...
566        if ( isY ){        if ( isY ){
567          yemax_en = 1000.;          if ( section.Contains("YE") ) yemax_en = 1000.;
568          xomax_en = 500.;          if ( section.Contains("XO") ) xomax_en = 500.;
569          fMax_planeye=(wpl-1)/2;          fMax_planeye=(wpl-1)/2;
570          fMax_planexo=(wpl-1)/2;          fMax_planexo=(wpl-1)/2;
571            if ( section.Contains("XE") ) xemax_en = 10.;
572            if ( section.Contains("YO") ) yomax_en = 5.;
573        } else {        } else {
574          yemax_en = 500.;          if ( section.Contains("YE") ) yemax_en = 500.;
575          xomax_en = 1000.;          if ( section.Contains("XO") ) xomax_en = 1000.;
576          fMax_planeye=(wpl-1)/2;          fMax_planeye=(wpl-1)/2;
577          fMax_planexo=(wpl-1)/2;          fMax_planexo=(wpl-1)/2;
578            if ( section.Contains("XE") ) xemax_en = 5.;
579            if ( section.Contains("YO") ) yomax_en = 10.;
580        };        };
       xemax_en = 0.;  
       yomax_en = 0.;  
581      };      };
582      if ( debug ) printf(" x0max %f wpl %i isY %i yomax_en %f xemax_en %f yemax_en %f xomax_en %f fMaxplane %i %i %i %i\n",x0max,wpl,isY,yomax_en,xemax_en,yemax_en,xomax_en,fMax_planeyo,fMax_planexe,fMax_planeye,fMax_planexo);      if ( debug ) printf(" x0max %f x0max/0.76 %f wpl %i isY %i yomax_en %f xemax_en %f yemax_en %f xomax_en %f fMaxplane %i %i %i %i\n",x0max,(x0max/0.76),wpl,isY,yomax_en,xemax_en,yemax_en,xomax_en,fMax_planeyo,fMax_planexe,fMax_planeye,fMax_planexo);
583    };    };
584    //    //
585    Int_t nPl = fPl;    Int_t nPl = fPl;
# Line 971  void CaloEnergy::Process(TString section Line 981  void CaloEnergy::Process(TString section
981        } else {        } else {
982          clong->Fit();          clong->Fit();
983        };              };      
984        fXOen_maxplane = clong->Get_E0();        if ( clong->GetLowerLimit() != 0. || clong->GetUpperLimit() != 0. ){
985            fXOen_maxplane = clong->Get_defE0();
986          } else {
987            fXOen_maxplane = clong->Get_E0();
988          };
989        fYOen_maxplane  = 0.;        fYOen_maxplane  = 0.;
990        fYEen_maxplane = 0.;        fYEen_maxplane = 0.;
991        fXEen_maxplane = 0.;        fXEen_maxplane = 0.;
Legend:
Removed from v.1.5  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.23