/[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.59 by mocchiut, Sat May 19 06:25:27 2012 UTC revision 1.60 by pam-mep, Thu Jun 7 14:04:00 2012 UTC
# Line 9  Line 9 
9  // ROOT headers  // ROOT headers
10  //  //
11  //#include <TCanvas.h>  //#include <TCanvas.h>
12  //#include <TH2F.h> //for test only. Vitaly.  #include <TH2F.h> //for test only. Vitaly.
13    #include <TVector3.h>
14  //#include <TF1.h>  //#include <TF1.h>
15    
16  #include <TTree.h>  #include <TTree.h>
# Line 290  int OrbitalInfoCore(UInt_t run, TFile *f Line 291  int OrbitalInfoCore(UInt_t run, TFile *f
291    }    }
292    if ( verbose ) cout<<"We have read recovered data"<<endl;    if ( verbose ) cout<<"We have read recovered data"<<endl;
293    
   
294    // IGRF stuff moved inside run loop!      // IGRF stuff moved inside run loop!  
295    
296    for (Int_t ip=0;ip<nz;ip++){    for (Int_t ip=0;ip<nz;ip++){
# Line 1064  int OrbitalInfoCore(UInt_t run, TFile *f Line 1064  int OrbitalInfoCore(UInt_t run, TFile *f
1064                    
1065          if ( debug ) printf(" fuffi \n");          if ( debug ) printf(" fuffi \n");
1066    
         //sineparam(q0sine,qtime,q0,qRoll,qPitch,0.60);  
         //sineparam(q1sine,qtime,q1,qRoll,qPitch,0.82);  
         //sineparam(q2sine,qtime,q2,qRoll,qPitch,0.82);  
         //sineparam(q3sine,qtime,q3,qRoll,qPitch,0.60);  
         //sineparam(Yawsine,qtime,qYaw,qRoll,qPitch,4);  
   
1067          if ( debug ) printf(" puffi \n");          if ( debug ) printf(" puffi \n");
1068          Double_t tmin = 9999999999.;          Double_t tmin = 9999999999.;
1069          Double_t tmax = 0.;          Double_t tmax = 0.;
# Line 1079  int OrbitalInfoCore(UInt_t run, TFile *f Line 1073  int OrbitalInfoCore(UInt_t run, TFile *f
1073          }          }
1074          if ( debug ) printf(" gnfuffi \n");          if ( debug ) printf(" gnfuffi \n");
1075    
         //q0testing->SetName("q0testing");  
         //q1testing->SetName("q1testing");  
         //q2testing->SetName("q2testing");  
         //q3testing->SetName("q3testing");  
           
 //      Int_t ss=10.*(tmax-tmin);  
         //q0testing->SetBins(ss,tmin,tmax,1000,-1.,1.);  
         //Pitchtesting->SetBins(ss,tmin,tmax,1000,-40.,40.);  
   
 //      for(Int_t tre = 0;tre<qtime.size();tre++){  
           //cout<<"q0["<<tre<<" = "<<q0[tre]<<endl;  
           //q0testing->Fill(qtime[tre],q0[tre]);  
           //q1testing->Fill(qtime[tre],q1[tre]);  
           //Pitchtesting->Fill(qtime[tre],qPitch[tre],100);  
           //if(qmode[tre] == -10)Pitchtesting->Fill(qtime[tre],10,100);  
           //q2testing->Fill(qtime[tre],q2[tre],100);  
           //q3testing->Fill(qtime[tre],q3[tre],100);  
 //      }  
           
         //for(Int_t tre=0;tre<q0sine.size();tre++)cout<<q1sine[tre].A<<"*sin("<<q1sine[tre].b<<"x+"<<q1sine[tre].c<<")\t time start: "<<q1sine[tre].startPoint<<"\ttime end: "<<q1sine[tre].finishPoint<<endl;  
         //for(Int_t tre=0;tre<q0sine.size();tre++)cout<<q1sine[tre].A<<"*sin("<<q1sine[tre].b<<"x+"<<q1sine[tre].c<<")\t time start: "<<q0sine[tre].startPoint<<"\ttime end: "<<q0sine[tre].finishPoint<<endl;  
1076        } // if we processed first event        } // if we processed first event
1077    
1078                
1079        //Filling Inclination information        //Filling Inclination information
1080        Double_t incli = 0;        Double_t incli = 0;
# Line 1228  int OrbitalInfoCore(UInt_t run, TFile *f Line 1202  int OrbitalInfoCore(UInt_t run, TFile *f
1202          orbitalinfo->Beast = beast;          orbitalinfo->Beast = beast;
1203          orbitalinfo->Bdown = bdown;          orbitalinfo->Bdown = bdown;
1204          orbitalinfo->Babs = babs;          orbitalinfo->Babs = babs;
1205            orbitalinfo->M = dimo;
1206          orbitalinfo->BB0 = babs/bequ;          orbitalinfo->BB0 = babs/bequ;
1207          orbitalinfo->L = xl;                orbitalinfo->L = xl;      
1208          // Set Stormer vertical cutoff using L shell.          // Set Stormer vertical cutoff using L shell.
Legend:
Removed from v.1.59  
changed lines
  Added in v.1.60
  ViewVC Help
Powered by ViewVC 1.1.23