/[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.13 by mocchiut, Wed Aug 19 07:25:34 2009 UTC revision 1.14 by mocchiut, Thu Sep 10 12:53:57 2009 UTC
# Line 676  CaloLong::CaloLong(PamLevel2 *l2p){   Line 676  CaloLong::CaloLong(PamLevel2 *l2p){  
676    suf = "";    suf = "";
677    xyaverage = true;    xyaverage = true;
678    //    //
679      letmax = -1;
680      heavytail = false;
681      lmipth = 100.;
682      //
683  };  };
684    
685  void CaloLong::MaskSection(TString sec){  void CaloLong::MaskSection(TString sec){
# Line 974  void CaloLong::SetEnergies(Float_t myene Line 978  void CaloLong::SetEnergies(Float_t myene
978    //    //
979    if ( vmax == 0 ) pmax++;    if ( vmax == 0 ) pmax++;
980    etmax = pmax * X0pl;    etmax = pmax * X0pl;
981      if ( letmax < 0. ) letmax = etmax;
982    //    //
983    if ( debug ) this->Print();    if ( debug ) this->Print();
984    if ( debug ) printf(" exit \n");    if ( debug ) printf(" exit \n");
# Line 1155  void CaloLong::Fit(Bool_t draw){ Line 1160  void CaloLong::Fit(Bool_t draw){
1160          exx[numpo] = 0.1;          exx[numpo] = 0.1;
1161          yyy[numpo] = enemip;          yyy[numpo] = enemip;
1162          eyy[numpo] = sqrt(enemip*3.)+sqrt(5.);          eyy[numpo] = sqrt(enemip*3.)+sqrt(5.);
1163            if ( xpos > letmax && enemip > lmipth && heavytail) eyy[numpo] = (sqrt(enemip*3.)+sqrt(5.))/numpo;
1164          numpo++;          numpo++;
1165          //      th->Fill(xpos,enemip);          //      th->Fill(xpos,enemip);
1166          if ( debug ) printf(" AVE Filling: st %i xpos %f energy %f \n",st,xpos,enemip);          if ( debug ) printf(" AVE Filling: st %i xpos %f energy %f \n",st,xpos,enemip);
# Line 1187  void CaloLong::Fit(Bool_t draw){ Line 1193  void CaloLong::Fit(Bool_t draw){
1193            exx[numpo] = 0.1;            exx[numpo] = 0.1;
1194            yyy[numpo] = enemip;            yyy[numpo] = enemip;
1195            eyy[numpo] = sqrt(enemip*3.)+sqrt(5.);            eyy[numpo] = sqrt(enemip*3.)+sqrt(5.);
1196            if ( xpos > letmax && enemip > lmipth && heavytail ) eyy[numpo] = (sqrt(enemip*3.)+sqrt(5.))/numpo;
1197            //      eyy[numpo] = sqrt(enemip)/(st*0.95);            //      eyy[numpo] = sqrt(enemip)/(st*0.95);
1198            numpo++;            numpo++;
1199            //      th->Fill(xpos,enemip);            //      th->Fill(xpos,enemip);
# Line 1237  void CaloLong::Fit(Bool_t draw){ Line 1244  void CaloLong::Fit(Bool_t draw){
1244    //  lfit->SetParLimits(0,0.,1000.);    //  lfit->SetParLimits(0,0.,1000.);
1245    //  lfit->SetParLimits(1,-1.,80.);    //  lfit->SetParLimits(1,-1.,80.);
1246    //  lfit->SetParLimits(2,-1.,10.);    //  lfit->SetParLimits(2,-1.,10.);
1247    //  TString optio = "ROW"; // "RO"    // TString optio = "ROW"; // "RO"
1248    TString optio = "RO"; // "RO"    TString optio = "RO"; // "RO"
1249    if ( debug ){    if ( debug ){
1250      optio += "NV";      optio += "NV";
# Line 1247  void CaloLong::Fit(Bool_t draw){ Line 1254  void CaloLong::Fit(Bool_t draw){
1254      if ( draw ) optio += "Q";      if ( draw ) optio += "Q";
1255    };    };
1256    //    //
1257    if ( debug ) printf(" OK, start the fitting procedure...\n");    if ( debug ) printf(" OK, start the fitting procedure... options: %s \n",optio.Data());
1258    //    //
1259    //  fitresult = th->Fit("lfit",optio);    //  fitresult = th->Fit("lfit",optio);
1260    fitresult = gh->Fit("lfit",optio);    fitresult = gh->Fit("lfit",optio);
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.23