/[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.15 by mocchiut, Mon Sep 14 13:56:50 2009 UTC revision 1.16 by mocchiut, Fri Nov 27 10:30:29 2009 UTC
# Line 678  CaloLong::CaloLong(PamLevel2 *l2p){   Line 678  CaloLong::CaloLong(PamLevel2 *l2p){  
678    //    //
679    letmax = -1;    letmax = -1;
680    heavytail = false;    heavytail = false;
681    lmipth = 100.;    //  lmipth = 100.;
682      lmipth = 0.;
683    //    //
684  };  };
685    
# Line 1165  void CaloLong::Fit(Bool_t draw){ Line 1166  void CaloLong::Fit(Bool_t draw){
1166          exx[numpo] = 0.1;          exx[numpo] = 0.1;
1167          yyy[numpo] = enemip;          yyy[numpo] = enemip;
1168          eyy[numpo] = sqrt(enemip*3.)+sqrt(5.);          eyy[numpo] = sqrt(enemip*3.)+sqrt(5.);
1169          if ( xpos > letmax && enemip > lmipth && heavytail) eyy[numpo] = (sqrt(enemip*3.)+sqrt(5.))/numpo;          //      if ( xpos > letmax && enemip > lmipth && heavytail) eyy[numpo] = (sqrt(enemip*3.)+sqrt(5.))/numpo;
1170            if ( xpos > letmax && enemip > lmipth && heavytail) eyy[numpo] = sqrt(enemip)/5.;
1171            if ( xpos > letmax-1 && xpos < letmax+1 && heavytail ) eyy[numpo] /= 5.;
1172            //if ( xpos > letmax-2 && xpos < letmax+1 && heavytail ) eyy[numpo] /= 7.;
1173            if ( xpos < 3. && heavytail ) eyy[numpo] /= 5.;
1174          numpo++;          numpo++;
1175          //      th->Fill(xpos,enemip);          //      th->Fill(xpos,enemip);
1176          if ( debug ) printf(" AVE Filling: st %i xpos %f energy %f error %f letmax %f lmpith %f \n",st,xpos,enemip,eyy[numpo-1],letmax,lmipth);          if ( debug ) printf(" AVE Filling: st %i xpos %f energy %f error %f letmax %f lmpith %f \n",st,xpos,enemip,eyy[numpo-1],letmax,lmipth);
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.23