/[PAMELA software]/DarthVader/OrbitalInfo/src/OrbitalInfoCore.cpp
ViewVC logotype

Diff of /DarthVader/OrbitalInfo/src/OrbitalInfoCore.cpp

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

revision 1.79 by mocchiut, Tue Oct 14 14:07:23 2014 UTC revision 1.80 by pamela, Fri Jan 16 09:14:48 2015 UTC
# Line 1482  int OrbitalInfoCore(UInt_t run, TFile *f Line 1482  int OrbitalInfoCore(UInt_t run, TFile *f
1482        Int_t errq=0;        Int_t errq=0;
1483        Int_t azim=0;        Int_t azim=0;
1484        Int_t MU=0;        Int_t MU=0;
1485        for(UInt_t mu = must;mu<RTtime2.size()-1;mu++){        for(UInt_t mu = 0;mu<RTtime2.size()-1;mu++){
1486          if(atime<=RTstart[mu+1] && atime>=RTstart[mu]){          if(atime<RTstart[mu+1] && atime>=RTstart[mu]){
1487            errq=RTerrq[mu];            errq=RTerrq[mu];
1488            azim=RTazim[mu];            azim=RTazim[mu];
1489            MU=mu;            MU=mu;
# Line 1596  int OrbitalInfoCore(UInt_t run, TFile *f Line 1596  int OrbitalInfoCore(UInt_t run, TFile *f
1596                  bank=kar*atime+bak;                  bank=kar*atime+bak;
1597                }                }
1598                if(atime>=RTstart[MU] && atime<RTpluto1[MU]){                if(atime>=RTstart[MU] && atime<RTpluto1[MU]){
1599                   Double_t s_dBdt2=(RTbpluto1[MU]-RTbank1[MU])/(RTpluto1[MU]-RTstart[MU]);                   Double_t s_dBdt2=(RTbpluto1[MU]-RTbank1[MU])/(Int_t)(RTpluto1[MU]-RTstart[MU]);
1600                   Double_t s_t2=((Double_t)RTpluto1[MU]+(Double_t)RTstart[MU])/2. - RTstart[MU];                   Double_t s_t2=((Double_t)RTpluto1[MU]+(Double_t)RTstart[MU])/2. - RTstart[MU];
1601                   Double_t s_t1=RTstart[MU]-RTstart[MU];                   Double_t s_t1=RTstart[MU]-RTstart[MU];
1602                   Double_t s_k=s_dBdt2/(s_t2-s_t1);                   Double_t s_k=s_dBdt2/(s_t2-s_t1);
# Line 1607  int OrbitalInfoCore(UInt_t run, TFile *f Line 1607  int OrbitalInfoCore(UInt_t run, TFile *f
1607                   bank=0.5*s_k*(atime-RTstart[MU])*(atime-RTstart[MU]) + s_b*(atime-RTstart[MU]) + s_c;                   bank=0.5*s_k*(atime-RTstart[MU])*(atime-RTstart[MU]) + s_b*(atime-RTstart[MU]) + s_c;
1608               }               }
1609                if(atime>RTpluto2[MU] && atime<=RTstart[MU+1]){                if(atime>RTpluto2[MU] && atime<=RTstart[MU+1]){
1610                   Double_t s_dBdt2=(RTbpluto2[MU] - RTbank2[MU])/(RTpluto2[MU]-RTstart[MU+1]);                   Double_t s_dBdt2=(RTbpluto2[MU] - RTbank2[MU])/(Int_t)(RTpluto2[MU]-RTstart[MU+1]);
1611                   Double_t s_t2=((Double_t)RTpluto2[MU]+(Double_t)RTstart[MU+1])/2. - RTstart[MU];                   Double_t s_t2=((Double_t)RTpluto2[MU]+(Double_t)RTstart[MU+1])/2. - RTstart[MU];
1612                   Double_t s_t1=RTstart[MU+1]-RTstart[MU];                   Double_t s_t1=RTstart[MU+1]-RTstart[MU];
1613                   Double_t s_k=s_dBdt2/(s_t2-s_t1);                   Double_t s_k=s_dBdt2/(s_t2-s_t1);
Legend:
Removed from v.1.79  
changed lines
  Added in v.1.80
  ViewVC Help
Powered by ViewVC 1.1.23