/[PAMELA software]/PamelaLevel2/doc/examples/example2.C
ViewVC logotype

Diff of /PamelaLevel2/doc/examples/example2.C

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

revision 1.1 by pam-fi, Fri Jun 16 16:43:55 2006 UTC revision 1.2 by pam-fi, Mon Oct 9 12:37:52 2006 UTC
# Line 34  example2(TString file){ Line 34  example2(TString file){
34          //          //
35          TFile f(file);          TFile f(file);
36          //          //
37          TTree *T = pam_event->LoadPamTrees(&f);         // << load Pamela trees from file f          TTree *T = pam_event->GetPamTree(&f);         // << get Pamela trees from file f
38          Int_t nevent = T->GetEntries();          Int_t nevent = T->GetEntries();
39          // ********************          // ********************
40          // load magnetic field          // load magnetic field
# Line 57  example2(TString file){ Line 57  example2(TString file){
57          //          //
58          cout << endl<< " Start loop over events   ";          cout << endl<< " Start loop over events   ";
59          for (Int_t i=0; i<nevent;i++){          for (Int_t i=0; i<nevent;i++){
60                  //                  //
61                  T->GetEntry(i);                  T->GetEntry(i);
62                  //                  //
63                  if(pam_event->GetNTracks()==1){                 // << select events with only one track                  if(pam_event->GetNTracks()==1){                 // << select events with only one track
64                      //                      //
# Line 75  example2(TString file){ Line 75  example2(TString file){
75                                  ){                                        ){      
76                                                                    
77                              cout << endl<< "***** First trajectory"<< endl;                              cout << endl<< "***** First trajectory"<< endl;
78                              track->DoTrack2(tr1);       // << calculate the first trajectory in magnetic field                                  track->DoTrack2(tr1);   // << calculate the first trajectory in magnetic field
79  //                          tr1->Dump();                // dump the trajectory  //                          tr1->Dump();                // dump the trajectory
80                              cout << "Length: "<< tr1->GetLength()<<endl;  // << get the track length                              cout << "Length: "<< tr1->GetLength()<<endl;  // << get the track length
81                              //                              //
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.23