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

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

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

revision 1.4 by mocchiut, Wed Jun 10 13:00:23 2009 UTC revision 1.6 by mocchiut, Wed Jun 24 14:12:55 2009 UTC
# Line 459  CaloLong::CaloLong(PamLevel2 *l2p){   Line 459  CaloLong::CaloLong(PamLevel2 *l2p){  
459    //    //
460    no18x = true;    no18x = true;
461    debug = false;    debug = false;
462      maskXE = false;
463      maskXO = false;
464      maskYE = false;
465      maskYO = false;
466    //    //
467  };  };
468    
469    void CaloLong::MaskSection(TString sec){
470      sec.ToUpper();
471      if ( sec.Contains("XO") ) maskXO = true;
472      if ( sec.Contains("YO") ) maskYO = true;
473      if ( sec.Contains("XE") ) maskXE = true;
474      if ( sec.Contains("YE") ) maskYE = true;
475    }
476    
477  void CaloLong::Clear(){  void CaloLong::Clear(){
478    //    //
479    memset(eplane,0, 2*22*sizeof(Float_t));    memset(eplane,0, 2*22*sizeof(Float_t));
# Line 484  void CaloLong::Clear(){ Line 496  void CaloLong::Clear(){
496    
497  void CaloLong::Print(){  void CaloLong::Print(){
498    //    //
499    Process();    Fit();
500    //    //
501    printf("==================== Calorimeter Longitudinal Profile =======================\n");    printf("==================== Calorimeter Longitudinal Profile =======================\n");
502    printf(" OBT: %u PKT: %u ATIME: %u \n",OBT,PKT,atime);    printf(" OBT: %u PKT: %u ATIME: %u \n",OBT,PKT,atime);
# Line 597  void CaloLong::Process(){ Line 609  void CaloLong::Process(){
609    Int_t plane = 0;    Int_t plane = 0;
610    Int_t strip = 0;    Int_t strip = 0;
611    Float_t mip = 0.;    Float_t mip = 0.;
612      Bool_t gof = true;
613    for (Int_t i=0; i < L2->GetCaloLevel1()->istrip; i++){    for (Int_t i=0; i < L2->GetCaloLevel1()->istrip; i++){
614      mip = L2->GetCaloLevel1()->DecodeEstrip(i,view,plane,strip);      mip = L2->GetCaloLevel1()->DecodeEstrip(i,view,plane,strip);
615      eplane[view][plane] += mip;      gof = true;
616        if ( maskXE && (plane%2)==0 && view==1 ) gof = false;
617        if ( maskXO && (plane%2)!=0 && view==1 ) gof = false;
618        if ( maskYE && (plane%2)!=0 && view==0 ) gof = false;
619        if ( maskYO && (plane%2)==0 && view==0 ) gof = false;
620        if ( gof ) eplane[view][plane] += mip;
621    };    };
622    //    //
623    // inclination factor (stolen from Daniele's code)    // inclination factor (stolen from Daniele's code)
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.6
  ViewVC Help
Powered by ViewVC 1.1.23