/[PAMELA software]/PamelaLevel2/src/PamLevel2.cpp
ViewVC logotype

Diff of /PamelaLevel2/src/PamLevel2.cpp

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

revision 1.16 by pam-fi, Tue Jan 16 13:19:58 2007 UTC revision 1.32 by pam-fi, Tue Mar 20 09:57:40 2007 UTC
# Line 100  PamLevel2::PamLevel2(){ Line 100  PamLevel2::PamLevel2(){
100  //     nd_obj    = NDLevel2::GetNDLevel2();  //     nd_obj    = NDLevel2::GetNDLevel2();
101  //     ac_obj    = AcLevel2::GetAcLevel2();  //     ac_obj    = AcLevel2::GetAcLevel2();
102  //     orb_obj   = OrbitalInfo::GetOrbitalInfo();  //     orb_obj   = OrbitalInfo::GetOrbitalInfo();
103        
104        h0_obj    = 0;
105        trk0_obj  = 0;
106        calo0_obj  = 0;
107    
108      trk2_obj  = 0;      trk2_obj  = 0;
109      trk1_obj  = 0;      trk1_obj  = 0;
110      trkh_obj  = 0;      trkh_obj  = 0;
# Line 114  PamLevel2::PamLevel2(){ Line 119  PamLevel2::PamLevel2(){
119    
120      run_obj   = 0;//new GL_RUN();      run_obj   = 0;//new GL_RUN();
121      soft_obj   = 0;// Emiliano      soft_obj   = 0;// Emiliano
122            irun = -1;
123        runfirstentry = 0ULL;
124        runlastentry = 0ULL;
125    
126        l0_file = NULL;
127        l0_tree = NULL;
128        iroot   = -1;
129        dbc     = 0;
130    
131        prevshift = 0;
132        
133        irun = -1;
134        run_tree = NULL;
135        run_tree_clone = NULL;
136        sel_tree = NULL;
137        sel_tree_clone = NULL;
138        
139        irunentry = -1;
140        pam_tree = NULL;
141        for(Int_t i=0; i<NCLONES; i++ )pam_tree_clone[i]=NULL;
142    
143        host = "mysql://localhost/pamelaprod";
144        user = "anonymous";
145        psw = "";
146        const char *pamdbhost=gSystem->Getenv("PAM_DBHOST");
147        const char *pamdbuser=gSystem->Getenv("PAM_DBUSER");
148        const char *pamdbpsw=gSystem->Getenv("PAM_DBPSW");
149        if ( !pamdbhost ) pamdbhost = "";
150        if ( !pamdbuser ) pamdbuser = "";
151        if ( !pamdbpsw ) pamdbpsw = "";
152        if ( strcmp(pamdbhost,"") ) host = pamdbhost;
153        if ( strcmp(pamdbuser,"") ) user = pamdbuser;
154        if ( strcmp(pamdbpsw,"") ) psw = pamdbpsw;
155    
156      for(Int_t i=0; i<8; i++ )tree_clone[i]=NULL;      
157    
158  //    sorted_tracks = 0;//new TRefArray();  //    sorted_tracks = 0;//new TRefArray();
159            
160        CAL0 = false;
161      CAL1 = true;      CAL1 = true;
162      CAL2 = true;      CAL2 = true;
163      TRK2 = true;      TRK2 = true;
164      TRK1 = false;      TRK1 = false;
165      TRKh = false;      TRKh = false;
166      TRG = true;      TRKh = false;
167      TOF = true;      TRG  = true;
168      S4  = true;      TOF  = true;
169      ND  = true;      TOF0 = false;
170      AC  = true;      S4   = true;
171      ORB = true;      ND   = true;
172        AC   = true;
173        ORB  = true;
174            
175      RUN = true;      RUN  = true;
176    
177        SELLI = -1;
178    
179      tsorted=0;      tsorted=0;
180      timage=0;      timage=0;
# Line 147  void PamLevel2::Delete(){ Line 189  void PamLevel2::Delete(){
189      if(soft_obj)delete soft_obj; //Emiliano      if(soft_obj)delete soft_obj; //Emiliano
190    
191  //    cout << "void PamLevel2::Clear()"<<endl;  //    cout << "void PamLevel2::Clear()"<<endl;
192        if(h0_obj)   delete h0_obj;
193        if(trk0_obj) delete trk0_obj;      
194        if(calo0_obj) delete calo0_obj;    
195      if(trk1_obj) delete trk1_obj;            if(trk1_obj) delete trk1_obj;      
196      if(trk2_obj) delete trk2_obj;      if(trk2_obj) delete trk2_obj;
197      if(trkh_obj) delete trkh_obj;      if(trkh_obj) delete trkh_obj;
# Line 167  void PamLevel2::Delete(){ Line 212  void PamLevel2::Delete(){
212          timage->Delete();          timage->Delete();
213          delete timage;          delete timage;
214      }      }
215    
216        if(dbc){
217            dbc->Close();
218            delete dbc;
219        }
220        
221        if(l0_file)l0_file->Close();
222        if(pam_tree)pam_tree->Delete();;
223        if(run_tree)run_tree->Delete();;
224        if(sel_tree)sel_tree->Delete();;
225        for(Int_t i=0; i<NCLONES; i++ )if(pam_tree_clone[i])pam_tree_clone[i]->Delete();;
226        if(run_tree_clone)run_tree_clone->Delete();;
227        if(sel_tree_clone)sel_tree_clone->Delete();;
228        
229    
230  };  };
231    
232  /**  /**
# Line 177  void PamLevel2::Clear(){ Line 237  void PamLevel2::Clear(){
237  //    cout << "void PamLevel2::Clear()"<<endl;  //    cout << "void PamLevel2::Clear()"<<endl;
238    
239  //  //
240    // This method is called once for every entry but RunInfo and SoftInfo do not change until the next run so we cannot clear them here unless we don't
241    // want to load them for each event even if they are the same...
242    //
243  //    if(run_obj)delete run_obj;  //    if(run_obj)delete run_obj;
244      if(run_obj) run_obj->Clear();  // Emiliano: Do not deallocate run_obj here, it will give segmentation fault! call clear instead  //    if(run_obj) run_obj->Clear();  // Emiliano: Do not deallocate run_obj here, it will give segmentation fault! call clear instead
245      if(soft_obj) soft_obj->Clear();  //    if(soft_obj) soft_obj->Clear();
246    
247        if(h0_obj)   h0_obj->Clear();      
248        if(trk0_obj) trk0_obj->Clear();    
249        if(calo0_obj) calo0_obj->Clear();  
250      if(trk1_obj) trk1_obj->Clear();          if(trk1_obj) trk1_obj->Clear();    
251      if(trk2_obj) trk2_obj->Clear();      if(trk2_obj) trk2_obj->Clear();
252      if(trkh_obj) trkh_obj->Clear();      if(trkh_obj) trkh_obj->Clear();
# Line 493  void PamLevel2::SortTracks(TString how){ Line 559  void PamLevel2::SortTracks(TString how){
559                  use_CAL            &&                  use_CAL            &&
560                  calo2_obj->npcfit[1] > 15     &&   //no. of fit planes on Y view                  calo2_obj->npcfit[1] > 15     &&   //no. of fit planes on Y view
561                  calo2_obj->varcfit[1] < 1000. &&  //fit variance on Y view                  calo2_obj->varcfit[1] < 1000. &&  //fit variance on Y view
562                    cp && ci &&
563                  true){                  true){
564    
565                                    
# Line 514  void PamLevel2::SortTracks(TString how){ Line 581  void PamLevel2::SortTracks(TString how){
581                  cout << "void PamLevel2::SortTracks(TString how): how= "<<how<<" but ToFLevel2 not loaded !!!";                  cout << "void PamLevel2::SortTracks(TString how): how= "<<how<<" but ToFLevel2 not loaded !!!";
582                  return;                  return;
583              };              };
584              if( use_TOF ){              if( use_TOF && op && oi ){
585                                    
586                  Int_t nphit_p =0;                  Int_t nphit_p =0;
587                  Int_t nphit_i =0;                  Int_t nphit_i =0;
# Line 770  TTree *PamLevel2::GetPamTree(TFile *f, T Line 837  TTree *PamLevel2::GetPamTree(TFile *f, T
837    
838  //     if( !detlist.IsNull() ) SetWhichTrees(detlist);  //     if( !detlist.IsNull() ) SetWhichTrees(detlist);
839  //     else                    GetWhichTrees(f);  //     else                    GetWhichTrees(f);
840        
841        cout << "TTree *PamLevel2::GetPamTree(TFile *f, TString detlist ) -- obsolte "<<endl;
842    
843      SetWhichTrees(detlist);      SetWhichTrees(detlist);
844    
845      TTree *Trout =0;      TTree *Trout =0;
# Line 777  TTree *PamLevel2::GetPamTree(TFile *f, T Line 847  TTree *PamLevel2::GetPamTree(TFile *f, T
847      TString fname = f->GetName();      TString fname = f->GetName();
848      if(!CheckLevel2File(fname))return NULL;      if(!CheckLevel2File(fname))return NULL;
849    
850      UInt_t *found=0;      //    UInt_t *found=0;    
851    
852      cout<< "GetPamTree(TFile*,TString): detector list --> ";      cout<< "GetPamTree(TFile*,TString): detector list --> ";
853      if(TRK1)cout<<"TRK1 ";      if(TRK1)cout<<"TRK1 ";
# Line 792  TTree *PamLevel2::GetPamTree(TFile *f, T Line 862  TTree *PamLevel2::GetPamTree(TFile *f, T
862      if(S4)cout<<"S4 ";      if(S4)cout<<"S4 ";
863      if(ORB)cout<<"ORB ";      if(ORB)cout<<"ORB ";
864      cout << endl;      cout << endl;
865        if(SELLI)cout<<">>> Found selection-list <<<"<<endl;
866    
867      f->cd();      f->cd();
868    
# Line 799  TTree *PamLevel2::GetPamTree(TFile *f, T Line 870  TTree *PamLevel2::GetPamTree(TFile *f, T
870      TTree *T = (TTree*)f->Get("Tracker");      TTree *T = (TTree*)f->Get("Tracker");
871      if(T && (TRK2||TRK1||TRKh)) {      if(T && (TRK2||TRK1||TRKh)) {
872          if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));          if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
873          else    T->SetBranchStatus("TrkLevel2",0,found);  //      else    T->SetBranchStatus("TrkLevel2",0,found);
874          if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;          if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;
875          if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));          if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));
876          else    T->SetBranchStatus("TrkLevel1",0,found);  //      else    T->SetBranchStatus("TrkLevel1",0,found);
877          if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;          if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;
878          if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));          if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));
879          else    T->SetBranchStatus("TrkHough",0,found);  //      else    T->SetBranchStatus("TrkHough",0,found);
880          if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;          if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;
881          if(!Trout)Trout=T;          if(!Trout)Trout=T;
882          else Trout->AddFriend(T);          else Trout->AddFriend(T);
# Line 816  TTree *PamLevel2::GetPamTree(TFile *f, T Line 887  TTree *PamLevel2::GetPamTree(TFile *f, T
887      TTree *C = (TTree*)f->Get("Calorimeter");      TTree *C = (TTree*)f->Get("Calorimeter");
888      if(C && (CAL2||CAL1)) {      if(C && (CAL2||CAL1)) {
889          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));
890          else    C->SetBranchStatus("CaloLevel2",0,found);  //      else    C->SetBranchStatus("CaloLevel2",0,found);
891          if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;          if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;
892          if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));          if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));
893          else    C->SetBranchStatus("CaloLevel1",0,found);  //      else    C->SetBranchStatus("CaloLevel1",0,found);
894          if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;          if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;
895          if(!Trout)Trout=C;          if(!Trout)Trout=C;
896          else Trout->AddFriend(C);          else Trout->AddFriend(C);
# Line 887  TTree *PamLevel2::GetPamTree(TFile *f, T Line 958  TTree *PamLevel2::GetPamTree(TFile *f, T
958      }else{      }else{
959          cout << "OrbitalInfo  : missing tree"<<endl;          cout << "OrbitalInfo  : missing tree"<<endl;
960      };      };
961    
962        TTree *L = (TTree*)f->Get("SelectionList");
963        if(L && SELLI==1) {
964    //      L->SetBranchAddress("RunEntry",&irun);
965    //      cout << "SelectionList: set branch address RunEntry"<<endl;
966    //      L->SetBranchAddress("EventEntry",&irunentry);
967    //      cout << "SelectionList: set branch address EventEntry"<<endl;
968    //      if(!Trout)Trout=O;
969    //      else Trout->AddFriend("SelectionList");
970            cout << " TTree *PamLevel2::GetPamTree(TFile, TString) >>> SelectionList not implemente!!!!"<<endl;
971            sel_tree = 0;
972        }else{
973            cout << "SelectionList  : missing tree"<<endl;
974        };
975            
976      cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;      cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;
977    
978        pam_tree = (TChain*)Trout;
979            
980      return Trout;      return Trout;
981            
# Line 916  TList*  PamLevel2::GetListOfLevel2Files( Line 1003  TList*  PamLevel2::GetListOfLevel2Files(
1003      contents->SetOwner();      contents->SetOwner();
1004            
1005  //    char *fullpath;  //    char *fullpath;
1006      const char *fullpath;  //    const char *fullpath;
1007            
1008      // if no input file list is given:        // if no input file list is given:  
1009      if ( flisttxt != "" ){      if ( flisttxt != "" ){
1010                    
1011  //      if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) ){          //      if( !gSystem->IsFileInIncludePath(flisttxt,&fullpath) ){        
1012          if( !(fullpath = gSystem->FindFile("./",flisttxt)) ){  //      if( !(fullpath = gSystem->FindFile("./",flisttxt)) ){
1013              cout <<"File "<<flisttxt<<" not found"<<endl;  //          cout <<"File "<<flisttxt<<" not found"<<endl;
1014              return 0;  //          return 0;
1015          }        //      }      
1016  //              flisttxt = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt));  //      flisttxt = fullpath;
1017          flisttxt = fullpath;        if ( !flisttxt.EndsWith(".root") ){
1018    
1019            flisttxt = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt));
1020    
1021          if( !gSystem->ChangeDirectory(ddir) ){          if( !gSystem->ChangeDirectory(ddir) ){
1022              cout << "Cannot change directory : "<<ddir<<endl;              cout << "Cannot change directory : "<<ddir<<endl;
# Line 937  TList*  PamLevel2::GetListOfLevel2Files( Line 1026  TList*  PamLevel2::GetListOfLevel2Files(
1026          cout <<"Input file list : " << flisttxt <<endl;          cout <<"Input file list : " << flisttxt <<endl;
1027          ifstream in;          ifstream in;
1028          in.open(flisttxt, ios::in); //open input file list          in.open(flisttxt, ios::in); //open input file list
1029            if(!in.good()){
1030                cout <<" ERROR opening the file "<<endl;
1031                gSystem->ChangeDirectory(wdir); // back to the working directory
1032                return 0;
1033            }      
1034          int line=0;          int line=0;
1035          while (1) {          while (1) {
1036              TString file;              TString file;
# Line 951  TList*  PamLevel2::GetListOfLevel2Files( Line 1045  TList*  PamLevel2::GetListOfLevel2Files(
1045              if(file.Contains("#"))file = file(0,file.First("#"));              if(file.Contains("#"))file = file(0,file.First("#"));
1046  //          cout <<"(2) " << file << endl;  //          cout <<"(2) " << file << endl;
1047  //          if( gSystem->IsFileInIncludePath(file,&fullpath) ){  //          if( gSystem->IsFileInIncludePath(file,&fullpath) ){
1048              if( (fullpath = gSystem->FindFile(ddir,file)) ){  //          if( (fullpath = gSystem->FindFile(ddir,file)) ){
1049                if( file.EndsWith(".root") ){
1050                    char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(file));
1051                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));// add file to the list                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));// add file to the list
1052                    delete fullpath;
1053              }else{              }
1054    //          }else{
1055  //              if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl;  //              if(file.Data()!="")cout << "File: "<<file<<" ---> missing "<< endl;
1056              };  //          };
1057          };          };
1058          in.close();          in.close();
1059                  } else {
1060              char *fullpath = gSystem->ConcatFileName(gSystem->DirName(flisttxt),gSystem->BaseName(flisttxt));
1061              contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));// add file to the list
1062              delete fullpath;
1063          };        
1064      }else{      }else{
1065                    
1066          cout << "No input file list given."<<endl;          cout << "No input file list given."<<endl;
# Line 979  TList*  PamLevel2::GetListOfLevel2Files( Line 1079  TList*  PamLevel2::GetListOfLevel2Files(
1079          while ( (questo = (TSystemFile*) next()) ) {          while ( (questo = (TSystemFile*) next()) ) {
1080              TString name =  questo-> GetName();              TString name =  questo-> GetName();
1081              if( name.EndsWith(".root") ){              if( name.EndsWith(".root") ){
1082                  const char *fullpath = gSystem->FindFile(ddir,name);  //              const char *fullpath = gSystem->FindFile(ddir,name);
1083  //              char *fullpath;  //              char *fullpath;
1084  //              gSystem->IsFileInIncludePath(name,&fullpath);  //              gSystem->IsFileInIncludePath(name,&fullpath);
1085                    char *fullpath = gSystem->ConcatFileName(gSystem->DirName(ddir),gSystem->BaseName(name));
1086                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));                  contents->Add(new TSystemFile(fullpath,gSystem->DirName(fullpath)));
1087                    delete fullpath;
1088              };              };
1089          }          }
1090          delete temp;          delete temp;
# Line 992  TList*  PamLevel2::GetListOfLevel2Files( Line 1094  TList*  PamLevel2::GetListOfLevel2Files(
1094      gSystem->ChangeDirectory(wdir); // back to the working directory      gSystem->ChangeDirectory(wdir); // back to the working directory
1095  //      cout << endl << "Selected files:" << endl;  //      cout << endl << "Selected files:" << endl;
1096  //      contents->Print();  //      contents->Print();
1097      cout << contents->GetEntries()<<" files selected\n";      cout << contents->GetEntries()<<" files \n";
1098  //      cout << endl;  //      cout << endl;
1099  //      cout << "Working directory: "<< gSystem->WorkingDirectory()<< endl;  //      cout << "Working directory: "<< gSystem->WorkingDirectory()<< endl;
1100      return contents;      return contents;
# Line 1036  TChain *PamLevel2::GetPamTree(TList *fl, Line 1138  TChain *PamLevel2::GetPamTree(TList *fl,
1138      TChain *N = 0;      TChain *N = 0;
1139      TChain *A = 0;      TChain *A = 0;
1140      TChain *B = 0;      TChain *B = 0;
1141    
1142        TChain *L = 0;
1143            
1144      if(TRK2||TRK1||TRKh) T = new TChain("Tracker");          if(TRK2||TRK1||TRKh) T = new TChain("Tracker");    
1145      if(CAL2||CAL1)       C = new TChain("Calorimeter");      if(CAL2||CAL1)       C = new TChain("Calorimeter");
# Line 1045  TChain *PamLevel2::GetPamTree(TList *fl, Line 1149  TChain *PamLevel2::GetPamTree(TList *fl,
1149      if(ND)               N = new TChain("NeutronD");      if(ND)               N = new TChain("NeutronD");
1150      if(AC)               A = new TChain("Anticounter");      if(AC)               A = new TChain("Anticounter");
1151      if(ORB)              B = new TChain("OrbitalInfo");      if(ORB)              B = new TChain("OrbitalInfo");
1152    
1153         L = new TChain("SelectionList");
1154            
1155      // loop over files and create chains              // loop over files and create chains        
1156      TIter next(fl);      TIter next(fl);
# Line 1061  TChain *PamLevel2::GetPamTree(TList *fl, Line 1167  TChain *PamLevel2::GetPamTree(TList *fl,
1167              if(ND)               N->Add(name);              if(ND)               N->Add(name);
1168              if(AC)               A->Add(name);              if(AC)               A->Add(name);
1169              if(ORB)              B->Add(name);              if(ORB)              B->Add(name);
1170                if(SELLI==1)         L->Add(name);
1171          };          };
1172      }      }
1173            
1174      cout << "done chain \n";      cout << "done chain \n";
1175    
1176      UInt_t *found=0;  //    UInt_t *found=0;
1177  // Tracker  // Tracker
1178      if(T && (TRK2||TRK1||TRKh)) {      if(T && (TRK2||TRK1||TRKh)) {
1179          if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));          if(TRK2)T->SetBranchAddress("TrkLevel2", GetPointerTo("TrkLevel2"));
1180          else    T->SetBranchStatus("TrkLevel2",0,found);  //      else    T->SetBranchStatus("TrkLevel2",0,found);
1181          if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;          if(TRK2)cout << "Tracker      : set branch address TrkLevel2"<<endl;
1182          if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));          if(TRK1)T->SetBranchAddress("TrkLevel1", GetPointerTo("TrkLevel1"));
1183          else    T->SetBranchStatus("TrkLevel1",0,found);  //      else    T->SetBranchStatus("TrkLevel1",0,found);
1184          if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;          if(TRK1)cout << "Tracker      : set branch address TrkLevel1"<<endl;
1185          if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));          if(TRKh)T->SetBranchAddress("TrkHough", GetPointerTo("TrkHough"));
1186          else    T->SetBranchStatus("TrkHough",0,found);  //      else    T->SetBranchStatus("TrkHough",0,found);
1187          if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;          if(TRKh)cout << "Tracker      : set branch address TrkHough"<<endl;
1188          if(!Trout)Trout=T;          if(!Trout)Trout=T;
1189          else Trout->AddFriend("Tracker");          else Trout->AddFriend("Tracker");
# Line 1086  TChain *PamLevel2::GetPamTree(TList *fl, Line 1193  TChain *PamLevel2::GetPamTree(TList *fl,
1193      // Calorimeter      // Calorimeter
1194      if(C && (CAL2||CAL1)) {      if(C && (CAL2||CAL1)) {
1195          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));          if(CAL2)C->SetBranchAddress("CaloLevel2", GetPointerTo("CaloLevel2"));
1196          else    C->SetBranchStatus("CaloLevel2",0,found);  //      else    C->SetBranchStatus("CaloLevel2",0,found);
1197          if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;          if(CAL2)cout << "Calorimeter  : set branch address CaloLevel2"<<endl;
1198          if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));          if(CAL1)C->SetBranchAddress("CaloLevel1", GetPointerTo("CaloLevel1"));
1199          else    C->SetBranchStatus("CaloLevel1",0,found);  //      else    C->SetBranchStatus("CaloLevel1",0,found);
1200          if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;          if(CAL1)cout << "Calorimeter  : set branch address CaloLevel1"<<endl;
1201          if(!Trout)Trout=C;          if(!Trout)Trout=C;
1202          else Trout->AddFriend("Calorimeter");          else Trout->AddFriend("Calorimeter");
# Line 1150  TChain *PamLevel2::GetPamTree(TList *fl, Line 1257  TChain *PamLevel2::GetPamTree(TList *fl,
1257      }else{      }else{
1258          cout << "OrbitalInfo  : missing tree"<<endl;          cout << "OrbitalInfo  : missing tree"<<endl;
1259      };      };
1260    
1261        // Selection List
1262        if(L && SELLI==1) {
1263            cout<<">>> Found selection-list <<<"<<endl;
1264            L->SetBranchAddress("RunEntry",&irun);
1265            cout << "SelectionList: set branch address RunEntry"<<endl;
1266            L->SetBranchAddress("EventEntry",&irunentry);
1267            cout << "SelectionList: set branch address EventEntry"<<endl;
1268            sel_tree = L;
1269    //      if(!Trout)Trout=O;
1270    //      else Trout->AddFriend("SelectionList");
1271        }else{
1272    //      cout << "SelectionList  : missing tree"<<endl;
1273            if(L)L->Delete();
1274        };
1275            
1276  //    cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;  //    cout<<endl<<" Number of entries: "<<Trout->GetEntries()<<endl<<endl;
1277    
1278            pam_tree = Trout;
1279  //     // Tracker  
 //     if(TRK2||TRK1||TRKh){  
 //      if(!Tout)Tout=T;  
 //      else Tout->AddFriend("Tracker");  
 //     }  
 //     // Calorimeter  
 //     if(CAL2||CAL1) {          
 //      if(!Tout)Tout=C;  
 //      else Tout->AddFriend("Calorimeter");  
 //     };  
 //     // ToF      
 //     if(TOF) {  
 //      if(!Tout)Tout=O;  
 //      else Tout->AddFriend("ToF");  
 //     };  
 //     // Trigger  
 //     if(TRG) {  
 //      if(!Tout)Tout=R;  
 //      else Tout->AddFriend("Trigger");  
 //     };  
 //     // S4  
 //     if(S4) {  
 //      if(!Tout)Tout=S;  
 //      else Tout->AddFriend("S4");  
 //     };  
 //     // Neutron Detector  
 //     if(ND) {  
 //      if(!Tout)Tout=N;  
 //      else Tout->AddFriend("NeutronD");  
 //     };  
 //     // Anticounters  
 //     if(AC) {  
 //      if(!Tout)Tout=A;  
 //      else Tout->AddFriend("Anticounter");  
 //     };  
 //     // OrbitalInfo  
 //     if(ORB) {  
 //      if(!Tout)Tout=B;  
 //      else Tout->AddFriend("OrbitalInfo");  
 //     };  
       
 //    cout<<endl<<" Number of entries: "<<Tout->GetEntries()<<endl<<endl;  
       
 //    if( Tout->GetEntries() )PamLevel2::SetBranchAddress( Tout );  
       
1280      return Trout;      return Trout;
1281  }  }
1282    
# Line 1280  void PamLevel2::SetBranchAddress(TTree * Line 1358  void PamLevel2::SetBranchAddress(TTree *
1358          t->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));          t->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));
1359          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;
1360      };      };
1361        // SelectionList
1362        if(SELLI==1) {
1363            t->SetBranchAddress("RunEntry", &irun);
1364            cout << "SelectionList: set branch address RunEntry"<<endl;
1365            t->SetBranchAddress("EventEntry", &irunentry);
1366            cout << "SelectionList: set branch address EventEntry"<<endl;
1367        };
1368            
1369  }  }
1370  /**  /**
# Line 1353  void PamLevel2::SetBranchAddress(TChain Line 1438  void PamLevel2::SetBranchAddress(TChain
1438          t->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));          t->SetBranchAddress("OrbitalInfo", GetPointerTo("OrbitalInfo"));
1439          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;          cout << "OrbitalInfo  : set branch address OrbitalInfo"<<endl;
1440      };      };
1441        // SelectionList
1442        if(SELLI==1) {
1443            t->SetBranchAddress("RunEntry", &irun);
1444            cout << "SelectionList: set branch address RunEntry"<<endl;
1445            t->SetBranchAddress("EventEntry", &irunentry);
1446            cout << "SelectionList: set branch address EventEntry"<<endl;
1447        };
1448            
1449  }  }
1450    
# Line 1380  TChain *PamLevel2::GetRunTree(TList *fl) Line 1472  TChain *PamLevel2::GetRunTree(TList *fl)
1472              R->Add(name);              R->Add(name);
1473          };          };
1474      }      }
1475    
1476        if(R->GetNtrees()){
1477            
1478      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));          R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));
1479      cout << "Run         : set branch address RunInfo"<<endl;          cout << "Run         : set branch address RunInfo"<<endl;
1480      R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano          R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano
1481      cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano          cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano
1482            }else{
1483            delete R;
1484            R=0;
1485        }
1486    
1487        run_tree = R;
1488    
1489      return R;      return R;
1490            
1491  }  }
# Line 1400  TChain *PamLevel2::GetRunTree(TList *fl) Line 1500  TChain *PamLevel2::GetRunTree(TList *fl)
1500   */   */
1501  TTree *PamLevel2::GetRunTree(TFile *f){  TTree *PamLevel2::GetRunTree(TFile *f){
1502    
1503        cout << "TTree *PamLevel2::GetRunTree(TFile *f) -- obsolte "<<endl;
1504    
1505      TTree *R = (TTree*)f->Get("Run");      TTree *R = (TTree*)f->Get("Run");
1506            
1507      R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));      if(R){
1508      cout << "Run         : set branch address RunInfo"<<endl;          R->SetBranchAddress("RunInfo", GetPointerTo("RunInfo"));
1509      R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano          cout << "Run         : set branch address RunInfo"<<endl;
1510      cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano          R->SetBranchAddress("SoftInfo", GetPointerTo("SoftInfo")); // Emiliano
1511                cout << "Software    : set branch address SoftInfo"<<endl; // Emiliano
1512        }
1513    
1514        run_tree = (TChain*)R;
1515    
1516      return R;      return R;
1517            
1518  }  }
1519    /**
1520     * Update the runinfo informations (to be used to have Run infos event by event basis)
1521     * @param run Pointer to the chain/tree which contains run infos
1522     * @return true if a new run has been read, false if it is still the same run
1523     */
1524    Bool_t PamLevel2::UpdateRunInfo(TChain *run, ULong64_t iev){
1525      //
1526      // check if we have already called once GetEntry, if not call it
1527      //
1528      if ( run->GetEntries() <= 0 ) return(false);
1529      //
1530      
1531      Int_t oldrun = irun;
1532      // --------------------------------------
1533      // if it is a full file (not preselected)
1534      // --------------------------------------
1535      if(SELLI==0){
1536          if ( irun < 0 ){
1537              irun = 0;
1538              run->GetEntry(irun);
1539              runfirstentry = 0ULL;
1540              runlastentry += (ULong64_t)(this->GetRunInfo()->NEVENTS);
1541              if ( (ULong64_t)(this->GetRunInfo()->NEVENTS) > 0ULL ) runlastentry -= 1ULL;
1542              return(true);
1543          };      
1544    
1545          if( !GetOrbitalInfo() ){
1546              cout << "Bool_t PamLevel2::UpdateRunInfo(TChain *run, ULong64_t iev) -- ERROR -- missing OrbitalInfo "<<endl;
1547              return(false);
1548          }
1549          // modificato il controllo sull'aggiornamento del run, per evitare problemi
1550          // dovuti agli eventi annidati (NB! NEVENTS conta anche questi!!)
1551          //
1552          bool fromfirst = true;
1553          //
1554          while ( !(GetOrbitalInfo()->absTime >= GetRunInfo()->RUNHEADER_TIME && GetOrbitalInfo()->absTime <= GetRunInfo()->RUNTRAILER_TIME) && irun < run->GetEntries() ){
1555    //      while ( iev > (runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS-1)) && irun < run->GetEntries() ){
1556              irun++;
1557              run->GetEntry(irun);
1558              runfirstentry = runlastentry;
1559              if ( (ULong64_t)(this->GetRunInfo()->NEVENTS) > 0ULL ) runfirstentry += 1ULL;
1560              runlastentry += (ULong64_t)(this->GetRunInfo()->NEVENTS);
1561              //      printf(" iev %llu %u %llu \n",iev,this->GetRunInfo()->NEVENTS,(ULong64_t)(runfirstentry+(ULong64_t)(this->GetRunInfo()->NEVENTS)));
1562              //      printf(" abstime %u trailertime %u \n",GetOrbitalInfo()->absTime,GetRunInfo()->RUNTRAILER_TIME);
1563              //      printf(" IDRUN %u \n",GetRunInfo()->ID);
1564              //
1565              prevshift = 0;
1566              //
1567              if ( irun == run->GetEntries()-1 && fromfirst && !(GetOrbitalInfo()->absTime >= GetRunInfo()->RUNHEADER_TIME && GetOrbitalInfo()->absTime <= GetRunInfo()->RUNTRAILER_TIME)){
1568                  printf(" resetting irun \n");
1569                  fromfirst = false;
1570                  irun = 0;
1571                  run->GetEntry(irun);
1572                  runfirstentry = 0ULL;
1573                  runlastentry += (ULong64_t)(this->GetRunInfo()->NEVENTS);
1574                  if ( (ULong64_t)(this->GetRunInfo()->NEVENTS) > 0ULL ) runlastentry -= 1ULL;
1575              };
1576              //
1577          };
1578          //
1579          if ( !fromfirst && !(GetOrbitalInfo()->absTime >= GetRunInfo()->RUNHEADER_TIME && GetOrbitalInfo()->absTime <= GetRunInfo()->RUNTRAILER_TIME) ) printf(" Something very wrong here: cannot find RUN containing absolute time %u \n",GetOrbitalInfo()->absTime);
1580          //
1581          if ( irun == oldrun || irun >= run->GetEntries() ) return(false);
1582          //
1583          //  printf(" iev %llu irun %i nevents %u 1st %llu last %llu \n",iev,irun,this->GetRunInfo()->NEVENTS,(ULong64_t)runfirstentry,(ULong64_t)runlastentry);
1584          //
1585          return(true);    
1586      };
1587      // ----------------------------------------------------
1588      // if it is a preselected file (there is SelectionList)
1589      // NBNB - the event tree MUST be read first
1590      // ----------------------------------------------------
1591      if(SELLI==1){      
1592          sel_tree->GetEntry(iev);
1593          if(irun != oldrun){
1594              run->GetEntry(irun);
1595              prevshift = 0;
1596              return true;
1597          }
1598          return false;
1599      }
1600    
1601      return false;
1602      //
1603    };
1604    /**
1605     * Update the runinfo informations (to be used to have Run infos event by event basis)
1606     * @param run Pointer to the chain/tree which contains run infos
1607     * @return true if a new run has been read, false if it is still the same run
1608     */
1609    Bool_t PamLevel2::UpdateRunInfo(TTree *run, ULong64_t iev){
1610      return(UpdateRunInfo((TChain*)run,iev));
1611    };
1612    
1613  //--------------------------------------  //--------------------------------------
1614  //  //
1615  //  //
# Line 1422  TTree *PamLevel2::GetRunTree(TFile *f){ Line 1621  TTree *PamLevel2::GetRunTree(TFile *f){
1621  void PamLevel2::SetWhichTrees(TString detlist){  void PamLevel2::SetWhichTrees(TString detlist){
1622                    
1623      if(detlist.IsNull() || detlist.Contains("+ALL", TString::kIgnoreCase)){      if(detlist.IsNull() || detlist.Contains("+ALL", TString::kIgnoreCase)){
1624            CAL0 = false;
1625          CAL1 = true;          CAL1 = true;
1626          CAL2 = true;          CAL2 = true;
1627          TRK2 = true;          TRK2 = true;
1628          TRK1 = false;          TRK1 = false;
1629          TRKh = false;          TRKh = false;
1630            TRK0 = false;
1631          TRG = true;          TRG = true;
1632          TOF = true;          TOF = true;
1633            TOF0 = false;
1634          S4  = true;          S4  = true;
1635          ND  = true;          ND  = true;
1636          AC  = true;          AC  = true;
1637          ORB = true;          ORB = true;
1638      }else if( detlist.Contains("-ALL", TString::kIgnoreCase) ){      }else if( detlist.Contains("-ALL", TString::kIgnoreCase) ){
1639            CAL0 = false;
1640          CAL1 = false;          CAL1 = false;
1641          CAL2 = false;          CAL2 = false;
1642          TRK2 = false;          TRK2 = false;
1643          TRK1 = false;          TRK1 = false;
1644          TRKh = false;          TRKh = false;
1645            TRK0 = false;
1646          TRG = false;          TRG = false;
1647          TOF = false;          TOF = false;
1648            TOF0 = false;
1649          S4  = false;          S4  = false;
1650          ND  = false;          ND  = false;
1651          AC  = false;          AC  = false;
# Line 1452  void PamLevel2::SetWhichTrees(TString de Line 1657  void PamLevel2::SetWhichTrees(TString de
1657          if ( detlist.Contains("-CAL1", TString::kIgnoreCase) )CAL1=false;          if ( detlist.Contains("-CAL1", TString::kIgnoreCase) )CAL1=false;
1658          if ( detlist.Contains("+CAL1", TString::kIgnoreCase) )CAL1=true;          if ( detlist.Contains("+CAL1", TString::kIgnoreCase) )CAL1=true;
1659      };      };
1660    
1661        if( detlist.Contains("CAL0", TString::kIgnoreCase) ){
1662            if ( detlist.Contains("-CAL0", TString::kIgnoreCase) )CAL0=false;
1663            if ( detlist.Contains("+CAL0", TString::kIgnoreCase) )CAL0=true;
1664        };
1665                    
1666      if( detlist.Contains("CAL2", TString::kIgnoreCase)){      if( detlist.Contains("CAL2", TString::kIgnoreCase)){
1667          if ( detlist.Contains("-CAL2", TString::kIgnoreCase) )CAL2=false;          if ( detlist.Contains("-CAL2", TString::kIgnoreCase) )CAL2=false;
# Line 1464  void PamLevel2::SetWhichTrees(TString de Line 1674  void PamLevel2::SetWhichTrees(TString de
1674          CAL1=false;          CAL1=false;
1675      }      }
1676  //  -------------------------------------------------------------------------  //  -------------------------------------------------------------------------
1677        if( detlist.Contains("TRK0", TString::kIgnoreCase) ){
1678            if ( detlist.Contains("-TRK0", TString::kIgnoreCase) )TRK0=false;
1679            if ( detlist.Contains("+TRK0", TString::kIgnoreCase) )TRK0=true;
1680        };
1681    
1682      if( detlist.Contains("TRK1", TString::kIgnoreCase) ){      if( detlist.Contains("TRK1", TString::kIgnoreCase) ){
1683          if ( detlist.Contains("-TRK1", TString::kIgnoreCase) )TRK1=false;          if ( detlist.Contains("-TRK1", TString::kIgnoreCase) )TRK1=false;
1684          if ( detlist.Contains("+TRK1", TString::kIgnoreCase) )TRK1=true;          if ( detlist.Contains("+TRK1", TString::kIgnoreCase) )TRK1=true;
# Line 1492  void PamLevel2::SetWhichTrees(TString de Line 1707  void PamLevel2::SetWhichTrees(TString de
1707            
1708      if( detlist.Contains("-TOF", TString::kIgnoreCase) )TOF = false;      if( detlist.Contains("-TOF", TString::kIgnoreCase) )TOF = false;
1709      else if( detlist.Contains("+TOF", TString::kIgnoreCase) )TOF = true;      else if( detlist.Contains("+TOF", TString::kIgnoreCase) )TOF = true;
1710    
1711        if( detlist.Contains("-TOF0", TString::kIgnoreCase) )TOF0 = false;
1712        else if( detlist.Contains("+TOF0", TString::kIgnoreCase) )TOF0 = true;
1713            
1714      if( detlist.Contains("-S4",  TString::kIgnoreCase) )S4  = false;      if( detlist.Contains("-S4",  TString::kIgnoreCase) )S4  = false;
1715      else if( detlist.Contains("+S4",  TString::kIgnoreCase) )S4  = true;      else if( detlist.Contains("+S4",  TString::kIgnoreCase) )S4  = true;
# Line 1526  void PamLevel2::SetWhichTrees(TString de Line 1744  void PamLevel2::SetWhichTrees(TString de
1744   * Set tree/branch detector flags from the content of a tree   * Set tree/branch detector flags from the content of a tree
1745   */   */
1746  void  PamLevel2::GetWhichTrees(TFile* f){  void  PamLevel2::GetWhichTrees(TFile* f){
1747                
1748    
1749    
1750        cout << "void  PamLevel2::GetWhichTrees(TFile* f) --- WARNING!! --- ...potrebbe non funzionare "<<endl;
1751      // -----------      // -----------
1752      // reset flags      // reset flags
1753      // -----------      // -----------
# Line 1698  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1919  Bool_t  PamLevel2::CheckLevel2File(TStri
1919            
1920      Bool_t RUN__ok    = false;      Bool_t RUN__ok    = false;
1921            
1922            Bool_t SELLI__ok  = false;
1923    
1924      cout << "Checking file: "<<name<<endl;      cout << "Checking file: "<<name<<endl;
1925      TFile *f = new TFile(name.Data());      TFile *f = new TFile(name.Data());
# Line 1715  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1936  Bool_t  PamLevel2::CheckLevel2File(TStri
1936    
1937      while( (key = (TKey*)next()) ){      while( (key = (TKey*)next()) ){
1938                    
1939    //      cout << key->GetName() << endl;
1940  //      cout << key->GetName() << ""<<key->GetClassName()<<endl;  //      cout << key->GetName() << ""<<key->GetClassName()<<endl;
1941  //              cout << " Get tree: " << f->Get(key->GetName())<<endl;  //              cout << " Get tree: " << f->Get(key->GetName())<<endl;
1942  //      nev_previous = nev;  //      nev_previous = nev;
# Line 1725  Bool_t  PamLevel2::CheckLevel2File(TStri Line 1947  Bool_t  PamLevel2::CheckLevel2File(TStri
1947  //          return false;  //          return false;
1948  //      };  //      };
1949    
1950    
1951          if( !strcmp(key->GetName(),"Run"        ) )RUN__ok = true;                if( !strcmp(key->GetName(),"Run"        ) )RUN__ok = true;      
1952    
1953          //=========================================================              //=========================================================    
1954            if( !strcmp(key->GetName(),"SelectionList"    ) ){
1955                SELLI__ok = true;
1956                if(SELLI==1){
1957                    Int_t nevt = ((TTree*)f->Get("SelectionList"))->GetEntries();
1958                    if( nev && nevt!=nev){
1959                        cout << "File: "<< f->GetName() <<" discarded ---- SelectionList tree has "<<nevt<<" events instead of "<<nev<< endl;
1960                        return false;
1961                    }
1962                    nev=nevt;
1963                }
1964            }
1965    
1966            //=========================================================    
1967          if( !strcmp(key->GetName(),"Trigger"    ) ){          if( !strcmp(key->GetName(),"Trigger"    ) ){
1968              TRG__ok = true;              TRG__ok = true;
1969              if(TRG){              if(TRG){
# Line 1839  Bool_t  PamLevel2::CheckLevel2File(TStri Line 2075  Bool_t  PamLevel2::CheckLevel2File(TStri
2075          };                };      
2076    
2077      };      };
2078    
2079        if( SELLI==-1 )SELLI = (Int_t)SELLI__ok;
2080        if( SELLI==0 && SELLI__ok ){
2081            cout << "File: "<< f->GetName() <<" discarded ---- found SelectionList (it is not a full-event file)" << endl;
2082            return false;  
2083        }
2084        if( SELLI==1 && !SELLI__ok ){
2085            cout << "File: "<< f->GetName() <<" discarded ---- SelectionList missing" << endl;
2086            return false;
2087        }
2088            
2089    //    cout << "SELLI "<<SELLI<<endl;
2090    
2091  //     cout<< "CheckLevel2File(TString): detector list --> ";  //     cout<< "CheckLevel2File(TString): detector list --> ";
2092  //     if(TRK1__ok)cout<<"TRK1 ";  //     if(TRK1__ok)cout<<"TRK1 ";
2093  //     if(TRK2__ok)cout<<"TRK2 ";  //     if(TRK2__ok)cout<<"TRK2 ";
# Line 1948  void PamLevel2::CreateCloneTrees0( TChai Line 2196  void PamLevel2::CreateCloneTrees0( TChai
2196      cout << "Create clones of PAMELA trees "<<endl;      cout << "Create clones of PAMELA trees "<<endl;
2197            
2198      Int_t i=0;      Int_t i=0;
2199      tree_clone[i] = fChain->GetTree()->CloneTree(0);      pam_tree_clone[i] = fChain->GetTree()->CloneTree(0);
2200      TString name =  tree_clone[i]->GetName();      TString name =  pam_tree_clone[i]->GetName();
2201      name.Append("_clone");      name.Append("_clone");
2202  //    tree_clone[i]->SetName(name.Data());  //    pam_tree_clone[i]->SetName(name.Data());
2203      cout << tree_clone[i]->GetName() <<endl;      cout << pam_tree_clone[i]->GetName() <<endl;
2204      i++;      i++;
2205    
2206      TList *li = fChain->GetListOfFriends();      TList *li = fChain->GetListOfFriends();
# Line 1962  void PamLevel2::CreateCloneTrees0( TChai Line 2210  void PamLevel2::CreateCloneTrees0( TChai
2210      while( (T_friend = (TFriendElement*)next()) ){      while( (T_friend = (TFriendElement*)next()) ){
2211  //      cout<<T_friend->IsA()->GetName()<<" "<<T_friend->GetName()<<hex << T_friend->GetTree() << dec<<endl;  //      cout<<T_friend->IsA()->GetName()<<" "<<T_friend->GetName()<<hex << T_friend->GetTree() << dec<<endl;
2212  //      cout<<T_friend->GetTree()->GetName()<< endl;  //      cout<<T_friend->GetTree()->GetName()<< endl;
2213          tree_clone[i] = T_friend->GetTree()->CloneTree(0);          pam_tree_clone[i] = T_friend->GetTree()->CloneTree(0);
2214          tree_clone[i]->SetAutoSave(1000000);          pam_tree_clone[i]->SetAutoSave(1000000);
2215          name =  tree_clone[i]->GetName();          name =  pam_tree_clone[i]->GetName();
2216          name.Append("_clone");          name.Append("_clone");
2217  //      tree_clone[i]->SetName(name.Data());  //      pam_tree_clone[i]->SetName(name.Data());
2218          cout << tree_clone[i]->GetName() << endl;          cout << pam_tree_clone[i]->GetName() << endl;
2219          i++;          i++;
2220      }      }
2221            
# Line 1987  void PamLevel2::CreateCloneTrees(TFile * Line 2235  void PamLevel2::CreateCloneTrees(TFile *
2235      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2236      cout << "Create new PAMELA trees "<<endl;      cout << "Create new PAMELA trees "<<endl;
2237                    
     Int_t i=0;  
2238    
2239        run_tree_clone = new TTree("Run","PAMELA Level2 data from the GL_RUN table ");
2240        run_tree_clone->Branch("RunInfo","GL_RUN",GetPointerTo("RunInfo"));
2241        cout << "Run          : branch RunInfo"<<endl;
2242        run_tree_clone->Branch("SoftInfo","SoftInfo",GetPointerTo("SoftInfo"));
2243        cout << "Run          : branch SoftInfo"<<endl;
2244                
2245    
2246    
2247        sel_tree_clone = new TTree("SelectionList","List of selected events ");
2248        sel_tree_clone->Branch("RunEntry",&irun,"runentry/I");
2249        sel_tree_clone->Branch("EventEntry",&irunentry,"eventry/I");
2250        
2251    
2252        Int_t i=0;
2253      if(TRK1||TRK2||TRKh){      if(TRK1||TRK2||TRKh){
2254          tree_clone[i] = new TTree("Tracker","PAMELA tracker level2 data ");              pam_tree_clone[i] = new TTree("Tracker","PAMELA tracker level2 data ");
2255          if(TRK1) {          if(TRK1) {
2256              tree_clone[i]->Branch("TrkLevel1","TrkLevel1", GetPointerTo("TrkLevel1"));              pam_tree_clone[i]->Branch("TrkLevel1","TrkLevel1", GetPointerTo("TrkLevel1"));
2257              tree_clone[i]->BranchRef();              pam_tree_clone[i]->BranchRef();
2258              cout << "Tracker      : branch TrkLevel1"<<endl;              cout << "Tracker      : branch TrkLevel1"<<endl;
2259                cout << "CreateCloneTrees " << GetTrkLevel1()<<endl;
2260          };          };
2261          if(TRK2) {          if(TRK2) {
2262              tree_clone[i]->Branch("TrkLevel2", "TrkLevel2",GetPointerTo("TrkLevel2"));              pam_tree_clone[i]->Branch("TrkLevel2", "TrkLevel2",GetPointerTo("TrkLevel2"));
2263              cout << "Tracker      : branch TrkLevel2"<<endl;                      cout << "Tracker      : branch TrkLevel2"<<endl;        
2264          };          };
2265          if(TRKh) {          if(TRKh) {
2266              tree_clone[i]->Branch("TrkHough","TrkHough", GetPointerTo("TrkHough"));              pam_tree_clone[i]->Branch("TrkHough","TrkHough", GetPointerTo("TrkHough"));
2267              cout << "Tracker      : branch TrkHough"<<endl;              cout << "Tracker      : branch TrkHough"<<endl;
2268          };          };
2269          i++;          i++;
# Line 2009  void PamLevel2::CreateCloneTrees(TFile * Line 2271  void PamLevel2::CreateCloneTrees(TFile *
2271    
2272      // Calorimeter      // Calorimeter
2273      if(CAL1||CAL2){      if(CAL1||CAL2){
2274          tree_clone[i] = new TTree("Calorimeter","PAMELA calorimeter level2 data ");              pam_tree_clone[i] = new TTree("Calorimeter","PAMELA calorimeter level2 data ");
2275          if(CAL1) {          if(CAL1) {
2276              tree_clone[i]->Branch("CaloLevel1", "CaloLevel1", GetPointerTo("CaloLevel1"));              pam_tree_clone[i]->Branch("CaloLevel1", "CaloLevel1", GetPointerTo("CaloLevel1"));
2277              cout << "Calorimeter  : branch CaloLevel1"<<endl;              cout << "Calorimeter  : branch CaloLevel1"<<endl;
2278          };          };
2279          if(CAL2) {          if(CAL2) {
2280              tree_clone[i]->Branch("CaloLevel2","CaloLevel2", GetPointerTo("CaloLevel2"));              pam_tree_clone[i]->Branch("CaloLevel2","CaloLevel2", GetPointerTo("CaloLevel2"));
2281              cout << "Calorimeter  : branch CaloLevel2"<<endl;              cout << "Calorimeter  : branch CaloLevel2"<<endl;
2282          };          };
2283          i++;          i++;
# Line 2023  void PamLevel2::CreateCloneTrees(TFile * Line 2285  void PamLevel2::CreateCloneTrees(TFile *
2285    
2286      // ToF          // ToF    
2287      if(TOF) {      if(TOF) {
2288          tree_clone[i] = new TTree("ToF","PAMELA ToF level2 data ");              pam_tree_clone[i] = new TTree("ToF","PAMELA ToF level2 data ");
2289          tree_clone[i]->Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2"));          pam_tree_clone[i]->Branch("ToFLevel2","ToFLevel2", GetPointerTo("ToFLevel2"));
2290          cout << "ToF          : branch ToFLevel2"<<endl;          cout << "ToF          : branch ToFLevel2"<<endl;
2291          i++;          i++;
2292      };      };
2293      // Trigger      // Trigger
2294      if(TRG) {      if(TRG) {
2295          tree_clone[i] = new TTree("Trigger","PAMELA trigger level2 data ");              pam_tree_clone[i] = new TTree("Trigger","PAMELA trigger level2 data ");
2296          tree_clone[i]->Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2"));          pam_tree_clone[i]->Branch("TrigLevel2","TrigLevel2", GetPointerTo("TrigLevel2"));
2297          cout << "Trigger      : branch TrigLevel2"<<endl;          cout << "Trigger      : branch TrigLevel2"<<endl;
2298          i++;          i++;
2299      };      };
2300      // S4      // S4
2301      if(S4) {      if(S4) {
2302          tree_clone[i] = new TTree("S4","PAMELA S4 level2 data ");                pam_tree_clone[i] = new TTree("S4","PAMELA S4 level2 data ");  
2303          tree_clone[i]->Branch("S4Level2","S4Level2", GetPointerTo("S4Level2"));          pam_tree_clone[i]->Branch("S4Level2","S4Level2", GetPointerTo("S4Level2"));
2304          cout << "S4           : branch S4Level2"<<endl;          cout << "S4           : branch S4Level2"<<endl;
2305          i++;          i++;
2306      };      };
2307      // Neutron Detector      // Neutron Detector
2308      if(ND) {      if(ND) {
2309          tree_clone[i] = new TTree("NeutronD","PAMELA neutron detector level2 data ");            pam_tree_clone[i] = new TTree("NeutronD","PAMELA neutron detector level2 data ");      
2310          tree_clone[i]->Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2"));          pam_tree_clone[i]->Branch("NDLevel2","NDLevel2", GetPointerTo("NDLevel2"));
2311          cout << "NeutronD     : branch NDLevel2"<<endl;          cout << "NeutronD     : branch NDLevel2"<<endl;
2312          i++;          i++;
2313      };      };
2314      // Anticounters      // Anticounters
2315      if(AC) {      if(AC) {
2316          tree_clone[i] = new TTree("Anticounter","PAMELA anticounter detector level2 data ");              pam_tree_clone[i] = new TTree("Anticounter","PAMELA anticounter detector level2 data ");        
2317          tree_clone[i]->Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2"));          pam_tree_clone[i]->Branch("AcLevel2","AcLevel2", GetPointerTo("AcLevel2"));
2318          cout << "Anticounter  : branch AcLevel2"<<endl;          cout << "Anticounter  : branch AcLevel2"<<endl;
2319          i++;          i++;
2320      };      };
2321      // OrbitalInfo      // OrbitalInfo
2322      if(ORB) {      if(ORB) {
2323          tree_clone[i] = new TTree("OrbitalInfo","PAMELA oribital info  ");                pam_tree_clone[i] = new TTree("OrbitalInfo","PAMELA oribital info  ");  
2324          tree_clone[i]->Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo"));          pam_tree_clone[i]->Branch("OrbitalInfo","OrbitalInfo", GetPointerTo("OrbitalInfo"));
2325          cout << "OrbitalInfo  : branch OrbitalInfo"<<endl;          cout << "OrbitalInfo  : branch OrbitalInfo"<<endl;
2326          i++;          i++;
2327      };      };
# Line 2074  void PamLevel2::CreateCloneTrees(TFile * Line 2336  void PamLevel2::CreateCloneTrees(TFile *
2336   */   */
2337  //void PamLevel2::FillNewPamTree(TTree *T){  //void PamLevel2::FillNewPamTree(TTree *T){
2338  void PamLevel2::FillCloneTrees(){  void PamLevel2::FillCloneTrees(){
2339                
2340      for(Int_t i=0; i<8; i++){  //    cout << "PamLevel2::FillCloneTrees()" << irunentry << endl;
2341          if(tree_clone[i])tree_clone[i]->Fill();  
2342        for(Int_t i=0; i<NCLONES; i++){
2343            if(pam_tree_clone[i])pam_tree_clone[i]->Fill();
2344      }      }
2345        if(sel_tree_clone)sel_tree_clone->Fill();    
2346    
2347  }  }
2348    
2349    
2350  TTree* PamLevel2::GetCloneTree(TString name){  TTree* PamLevel2::GetCloneTree(TString name){
2351    
2352      for(Int_t i=0; i<8; i++){        for(Int_t i=0; i<NCLONES; i++){    
2353          if(tree_clone[i]){          if(pam_tree_clone[i]){
2354              TString na = tree_clone[i]->GetName();              TString na = pam_tree_clone[i]->GetName();
2355              if(!name.CompareTo(na))return tree_clone[i];              if(!name.CompareTo(na))return pam_tree_clone[i];
2356          };          };
2357      }      }
2358        if(run_tree_clone){
2359            TString na = run_tree_clone->GetName();
2360            if(!name.CompareTo(na))return run_tree_clone;
2361        }
2362        if(sel_tree_clone){
2363            TString na = sel_tree_clone->GetName();
2364            if(!name.CompareTo(na))return sel_tree_clone;
2365        }
2366      return NULL;      return NULL;
2367    
2368  }  }
2369  void PamLevel2::WriteCloneTrees(){  void PamLevel2::WriteCloneTrees(){
2370      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
2371      cout << "Write clones of PAMELA trees "<<endl;      cout << "Write clones of PAMELA trees "<<endl;
2372      for(Int_t i=0; i<8; i++){        cout << run_tree_clone->GetName()<<endl;    
2373          if(tree_clone[i]){      run_tree_clone->Write();
2374              cout << tree_clone[i]->GetName()<<endl;      cout << sel_tree_clone->GetName()<<endl;    
2375              tree_clone[i]->Write();      sel_tree_clone->Write();
2376        for(Int_t i=0; i<NCLONES; i++){    
2377            if(pam_tree_clone[i]){
2378                cout << pam_tree_clone[i]->GetName()<<endl;
2379                pam_tree_clone[i]->Write();
2380          };          };
2381      }      }
2382      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;
# Line 2106  void PamLevel2::WriteCloneTrees(){ Line 2384  void PamLevel2::WriteCloneTrees(){
2384  }  }
2385    
2386  /**  /**
2387   * Create a new (empty) Pamela trees   * Method to get level2-trees entry, the corresponding run entry and (if required) the level0 entry.
2388   */   */
2389  // TTree* PamLevel2::GetNewPamTree(){  Int_t PamLevel2::GetEntry(Int_t iee){
2390        
2391        if(!pam_tree){
2392            cout << " Int_t PamLevel2::GetEntry(Int_t) -- ERROR -- level2 trees not loaded"<<endl;
2393            return 0;
2394        }
2395    
 //     if(tree_clone)return tree_clone;  
2396    
2397  //     TTree *Tout = 0;      //
2398        // This is a sort of bug: if you don't have the run tree you don't want to exit here you want to have loaded the event anyway...
2399        //
2400        //    if(!run_tree ){
2401        //  cout << " Int_t PamLevel2::GetEntry(Int_t) -- ERROR -- run tree not loeaded"<<endl;
2402        //  return 0;
2403        //    }
2404    
2405        Int_t ii=0;
2406        //-------------------------------
2407        ii = iee;
2408        if( !pam_tree->GetEntry(ii) ) return 0;
2409    
2410        //
2411        // ... that's way I put it here. Notice that nothing change in the code (is backward compatible) since in any case you return with 0.
2412        // in theory one would like to return 1 if run is not loaded but now I don't have the will to add that 2 lines of code and it is not
2413        // a problem if you don't check the return code of getentry.
2414        //
2415        if(!run_tree ){
2416            if ( TRK0 || CAL0 || TOF0 || RUN ) { //forse cosi` va bene per tornare 1?
2417                    cout << " Int_t PamLevel2::GetEntry(Int_t) -- ERROR -- run tree not loaded"<<endl;
2418                    return 0;
2419            }  else {
2420                    return 1; //cosi` se non c'e` run esce qua...
2421            }
2422        }
2423    
2424        //-------------------------------
2425        ii = iee;
2426        Bool_t UPDATED = UpdateRunInfo(run_tree,ii);
2427        if(SELLI==0)irunentry = iee-runfirstentry;
2428        if(UPDATED && run_tree_clone)run_tree_clone->Fill();
2429    
2430    //    cout << "PamLevel2::GetEntry("<<iee<<") "<<irun<<" "<<runfirstentry<<" "<<irunentry<<endl;
2431    
2432        if( TRK0 || CAL0 || TOF0 )GetYodaEntry( );
2433    
2434        return 1;
2435    
2436    }
2437    
2438    /**
2439     * Method to retrieve the level0 tree (YODA tree) that contains the current event.
2440     * Given the run ID (...), if needed it query the DB and load the proper file.
2441     * @return Pointer to the tree
2442     */
2443    
2444    
2445    TTree* PamLevel2::GetYodaTree( ){
2446        
2447    //    cout << "TTree* PamLevel2::GetYodaTree( )"<<endl;
2448        //===================================
2449        // check if iroot has changed
2450        //===================================
2451        if( irun<0 ){
2452            cout << "PamLevel2::GetYodaTree() -- ERROR "<<endl;
2453            cout << "In order to use this method you have to load the RunInfo tree "<<endl;
2454            cout << "with the method TChain* PamLevel2::GetRunTree(TString, TString)"<<endl;
2455    //      cout << " - read the event with PamLevel2::GetEntry(Int_t)"<<endl;
2456            return NULL;
2457        }    
2458        Int_t irootnew = run_obj->ID_ROOT_L0;
2459    //     cout << "iroot    "<<iroot<<endl;
2460    //     cout << "irootnew "<<irootnew<<endl;
2461    
2462        //===================================
2463        // load the level0 file
2464        // (if not already loaded)
2465        //===================================
2466        if( iroot != irootnew || !l0_tree){
2467            iroot = irootnew;
2468            //===================================
2469            // open the DB connection
2470            // (if not already opened)
2471            //===================================
2472            if(!dbc || (dbc && !dbc->IsConnected())){
2473                cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<<endl;
2474                cout<<"Connecting to DB"<<endl;
2475                cout<<"HOST "<<host<<endl;
2476                cout<<"USER "<<user<<endl;
2477                cout<<"PSW  "<<psw<<endl;
2478                dbc = TSQLServer::Connect(host.Data(),user.Data(),psw.Data());
2479                if( !dbc )return NULL;
2480                if( !dbc->IsConnected() )return NULL;    
2481    //          cout<<"...done"<<endl;
2482                cout<<"~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~"<<endl;
2483            }else{
2484    //          cout<<"DB already connected"<<endl;
2485            }
2486            GL_ROOT glroot = GL_ROOT();
2487            if( glroot.Query_GL_ROOT(iroot,dbc) ){
2488                cout << "TTree* PamLevel2::GetYodaTree( ) -- ERROR -- level0 file iroot = "<<iroot<< " does not exists"<<endl;
2489                return NULL;
2490            };
2491            TString filename = glroot.PATH + glroot.NAME;
2492            if(l0_file){
2493                l0_file->Close();
2494                l0_file->Delete();
2495            }
2496            cout << "Opening LEVEL0 file: "<< filename << endl;
2497            FileStat_t t;
2498            if( gSystem->GetPathInfo(filename.Data(),t) ){
2499                cout << " PamLevel2::GetYodaTree() -- ERROR opening file "<<endl;
2500                return NULL;
2501            }
2502            l0_file = new TFile(filename);
2503            if( !l0_file )return NULL;
2504            l0_tree = (TTree*)l0_file->Get("Physics");
2505            if(!h0_obj)h0_obj = new EventHeader();
2506            l0_tree->SetBranchAddress("Header" ,&h0_obj);
2507            prevshift = 0;
2508            //---------------------------------------------------
2509            // TRACKER:
2510            if(TRK0){
2511                if(!trk0_obj){
2512                    trk0_obj = new TrkLevel0();
2513                    trk0_obj->Set();
2514                };
2515                l0_tree->SetBranchAddress("Tracker" ,trk0_obj->GetPointerToTrackerEvent());
2516            }
2517            //---------------------------------------------------
2518            // CALORIMETER:
2519            if(CAL0){
2520                if(!calo0_obj){
2521                    calo0_obj = new CaloLevel0();
2522                    calo0_obj->Set();
2523                };
2524                l0_tree->SetBranchAddress("Calorimeter" ,calo0_obj->GetPointerToCalorimeterEvent());
2525                //      cout << "PamLevel2::GetYodaTree() --- level0 calorimeter not implemented "<<endl;
2526            }
2527            //---------------------------------------------------
2528            // TOF:
2529            if(TOF0){
2530                cout << "PamLevel2::GetYodaTree() --- level0 TOF not implemented "<<endl;
2531            }
2532    
2533        };
2534    
2535        if(!dbc || (dbc && !dbc->IsConnected())){
2536            cout << " TTree* PamLevel2::GetYodaTree( ) -- no DB connected... hai fatto qualche cazzata "<<endl;
2537        }
2538    
2539        if ( TRK0 ) TrkParams::SetCalib(run_obj,dbc);
2540    
2541    
2542    //    cout << l0_tree << endl;
2543            
2544  //     cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;      return l0_tree;
2545  //     cout << "Create new PAMELA trees "<<endl;  
2546    }
2547    
2548    /**
2549     * Method to retrieve the level0 tree (YODA tree) that contains the current event.
2550     */
2551    Int_t PamLevel2::GetYodaEntry(){
2552    
2553    //    cout << "Int_t PamLevel2::GetYodaEntry()"<<endl;
2554        if(!GetYodaTree())return 0;
2555            
2556        // patch
2557        if( (Int_t)irunentry < 0){
2558            cout << "Int_t PamLevel2::GetYodaEntry() -- ATTENZIONE -- irunentry negativo?!?! "<<(Int_t)irunentry<<endl;
2559            irunentry=0;
2560        }
2561        //  ---------------------------------
2562        //  if file is NOT a preselected file
2563        //  ---------------------------------
2564        UInt_t quellagiusta = irunentry + run_obj->EV_FROM;
2565            
 //     if(TRK1||TRK2||TRKh){  
 //      TTree *T = new TTree("Tracker_clone","PAMELA tracker level2 data ");      
 //      if(TRK1) {  
 //          trk1_clone = new TrkLevel1();  
 //          T->Branch("TrkLevel1","TrkLevel1", &trk1_clone);  
 //          T->BranchRef();  
 //          cout << "Tracker      : branch TrkLevel1"<<endl;  
 //      };  
 //      if(TRK2) {  
 //          trk2_clone = new TrkLevel2();  
 //          T->Branch("TrkLevel2", "TrkLevel2",&trk2_clone);  
 //          cout << "Tracker      : branch TrkLevel2"<<endl;          
 //      };  
 //      if(TRKh) {  
 //          trkh_clone = new TrkHough();  
 //          T->Branch("TrkHough","TrkHough", &trkh_clone);  
 //          cout << "Tracker      : branch TrkHough"<<endl;  
 //      };  
 //      if(!Tout)Tout=T;  
 //      else Tout->AddFriend("Tracker_clone");  
 //     }  
   
 //     // Calorimeter  
 //     if(CAL1||CAL2){  
 //      TTree *C = new TTree("Calorimeter_clone","PAMELA calorimeter level2 data ");      
 //      if(CAL1) {  
 //          calo1_clone = new CaloLevel1();  
 //          C->Branch("CaloLevel1", "CaloLevel1", &calo1_clone);  
 //          cout << "Calorimeter  : branch CaloLevel1"<<endl;  
 //      };  
 //      if(CAL2) {  
 //          calo2_clone = new CaloLevel2();  
 //          C->Branch("CaloLevel2","CaloLevel2", &calo2_clone);  
 //          cout << "Calorimeter  : branch CaloLevel2"<<endl;  
 //      };  
 //      if(!Tout)Tout=C;  
 //      else Tout->AddFriend("Calorimeter_clone");  
 //     }      
   
 //     // ToF      
 //     if(TOF) {  
 //      TTree *O = new TTree("ToF_clone","PAMELA ToF level2 data ");      
 //      tof_clone = new ToFLevel2();  
 //      O->Branch("ToFLevel2","ToFLevel2", &tof_clone);  
 //      cout << "ToF          : branch ToFLevel2"<<endl;  
 //      if(!Tout)Tout=O;  
 //      else Tout->AddFriend("ToF_clone");  
 //     };  
 //     // Trigger  
 //     if(TRG) {  
 //      TTree *R = new TTree("Trigger_clone","PAMELA trigger level2 data ");      
 //      trig_clone = new TrigLevel2();  
 //      R->Branch("TrigLevel2","TrigLevel2", &trig_clone);  
 //      cout << "Trigger      : branch TrigLevel2"<<endl;  
 //      if(!Tout)Tout=R;  
 //      else Tout->AddFriend("Trigger_clone");  
 //     };  
 //     // S4  
 //     if(S4) {  
 //      TTree *S = new TTree("S4_clone","PAMELA S4 level2 data ");        
 //      s4_clone = new S4Level2();  
 //      S->Branch("S4Level2","S4Level2", &s4_clone);  
 //      cout << "S4           : branch S4Level2"<<endl;  
 //      if(!Tout)Tout=S;  
 //      else Tout->AddFriend("S4_clone");  
 //     };  
 //     // Neutron Detector  
 //     if(ND) {  
 //      TTree *N = new TTree("NeutronD_clone","PAMELA neutron detector level2 data ");    
 //      nd_clone = new NDLevel2();  
 //      N->Branch("NDLevel2","NDLevel2", &nd_clone);  
 //      cout << "NeutronD     : branch NDLevel2"<<endl;  
 //      if(!Tout)Tout=N;  
 //      else Tout->AddFriend("NeutronD_clone");  
 //     };  
 //     // Anticounters  
 //     if(AC) {  
 //      TTree *A = new TTree("Anticounter_clone","PAMELA anticounter detector level2 data ");    
 //      ac_clone = new AcLevel2();  
 //      A->Branch("AcLevel2","AcLevel2", &ac_clone);  
 //      cout << "Anticounter  : branch AcLevel2"<<endl;  
 //      if(!Tout)Tout=A;  
 //      else Tout->AddFriend("Anticounter_clone");  
 //     };  
 //     // OrbitalInfo  
 //     if(ORB) {  
 //      TTree *B = new TTree("OrbitalInfo_clone","PAMELA oribital info  ");      
 //      orb_clone = new OrbitalInfo();  
 //      B->Branch("OrbitalInfo","OrbitalInfo", &orb_clone);  
 //      cout << "OrbitalInfo  : branch OrbitalInfo"<<endl;  
 //      if(!Tout)Tout=B;  
 //      else Tout->AddFriend("OrbitalInfo_clone");  
 //     };  
 //     cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" <<endl;  
2566    
2567  //     tree_clone = Tout;  //     cout << " irun "<< irun << " irunentry "<< irunentry<<" run_obj->EV_FROM "<<run_obj->EV_FROM <<" quella giusta "<<quellagiusta << endl;
2568  //     tree_clone->SetDirectory(0);  //     cout << " iroot "<<iroot<<" run_obj->ID_ROOT_L0 "<<run_obj->ID_ROOT_L0<<endl;
2569    //     cout << " time "<< GetOrbitalInfo()->absTime << endl;
2570    //     cout << " trk_calib_used "<<run_obj->TRK_CALIB_USED<< endl;
2571        
2572        if( !GetOrbitalInfo() ){
2573            cout << "Int_t PamLevel2::GetYodaEntry() -- ERROR -- missing OrbitalInfo "<<endl;
2574            return 0;
2575        }
2576    
2577  //     return Tout;      // ---------------------------------------------------------------------
2578  // }      // ATTENTION!!!
2579  // /**      // If data are not pre-processed with cleaner, the level0 tree may contain
2580  //  * Fill a tree (created with GetNewPamTree)      // spurious nested physics packets.
2581  //  *      // The GL_RUN variables EV_FROM, EV_TO and NEVENTS counts also these entries
2582  //  */      // while level2 tree DOES NOT!!
2583  // //void PamLevel2::FillNewPamTree(TTree *T){      // This means that "quellagiusta" in these cases is not correct.
2584  // void PamLevel2::FillNewPamTree(){      // In order to retrieve the correct level0 event, I implemented a check
2585                // of the OBT and pkt-number. In case of mismatch, the level0 entry number
2586        // is shift forward until when the packets match.
2587        // ---------------------------------------------------------------------
2588        Int_t answer = 0;
2589        Int_t shift =0;
2590        //    printf(" siamo qui %i %i \n",shift,prevshift);
2591        do{
2592            if(shift>0){    
2593              cout << " level0 <--> level2 mismatch ( irun "<<irun<<" irunentry "<<irunentry<<" shift "<<shift<<" )"<<endl;
2594            }
2595            answer = l0_tree->GetEntry(quellagiusta+shift+prevshift);
2596            shift++;
2597            if( !GetEventHeader() ){
2598                cout << "Int_t PamLevel2::GetYodaEntry() -- ERROR -- missing EventHeader "<<endl;
2599                return 0;
2600            }
2601    //      if( (quellagiusta+shift) == l0_tree->GetEntries()+1 )cout << ">>> end of level0 tree <<<"<<endl;
2602    //      cout << " GetOrbitalInfo()->OBT "<< GetOrbitalInfo()->OBT << endl;
2603    //      cout << " GetEventHeader()->GetPscuHeader()->GetOrbitalTime() "<< GetEventHeader()->GetPscuHeader()->GetOrbitalTime() << endl;
2604    //      cout << " GetOrbitalInfo()->pkt_num "<< GetOrbitalInfo()->pkt_num << endl;
2605    //      cout << " GetEventHeader()->GetPscuHeader()->GetCounter() "<< GetEventHeader()->GetPscuHeader()->GetCounter() << endl;
2606    //      printf(" IDRUN %u \n",GetRunInfo()->ID);
2607    //
2608            if ( prevshift != 0 && (quellagiusta+shift) == GetYodaTree()->GetEntries() ){
2609              prevshift = 0;
2610              shift = -1;
2611            };
2612    
2613        }while( ( GetOrbitalInfo()->OBT != (UInt_t)(GetEventHeader()->GetPscuHeader()->GetOrbitalTime()) || GetOrbitalInfo()->pkt_num != (UInt_t)(GetEventHeader()->GetPscuHeader()->GetCounter())) && (quellagiusta+shift) < GetYodaTree()->GetEntries() );
2614    
2615        if ( (quellagiusta+shift) == GetYodaTree()->GetEntries() ) cout << " Big trouble here, no such event in Level0 data! " <<endl;
2616    
2617    //    cout << "LA ENTRY GIUSTA E`: "<<quellagiusta<<" (spero...)"<<endl;
2618    //    return GetYodaTree()->GetEntry(quellagiusta);
2619        if ( shift > 1 ) prevshift += (shift-1);
2620        
2621        return answer;
2622        
2623    }
2624    
 //     if(trk1_clone)  *trk1_clone = *trk1_obj;  
 //     if(trk2_clone){  
 //      trk2_clone->Clear();  
 //      trk2_obj->Copy(*trk2_clone);  
 // //   *trk2_clone = *trk2_obj;  
 //     }  
 //     if(trkh_clone)  *trkh_clone = *trkh_obj;  
 //     if(calo1_clone){  
 // //   *calo1_clone = *calo1_obj;  
 //      calo1_clone->Clear();  
 //      calo1_obj->Copy(*calo1_clone);  
 //     }  
 //     if(calo2_clone){  
 // //   *calo2_clone = *calo2_obj;  
 //      calo2_clone->Clear();  
 //      calo2_obj->Copy(*calo2_clone);  
 //     }  
 //     if(tof_clone)   *tof_clone  = *tof_obj;  
 //     if(trig_clone)  *trig_clone = *trig_obj;  
 //     if(s4_clone)    *s4_clone   = *s4_obj;  
 //     if(nd_clone)    *nd_clone   = *nd_obj;  
 //     if(ac_clone)    *ac_clone   = *ac_obj;  
 //     if(orb_clone)   *orb_clone  = *orb_obj;  
   
 //     TTree *T = tree_clone;  
   
 //     T->Fill(); //fill main tree  
 // //    cout<<T->IsA()->GetName()<<" "<<T->GetName()<<endl;  
 //     TList *li = T->GetListOfFriends();  
 //     TIter next(li);  
 //     TFriendElement* T_friend=0;  
 //     while( (T_friend = (TFriendElement*)next()) ){  
 // //   cout<<T_friend->IsA()->GetName()<<" "<<T_friend->GetName()<<hex << T_friend->GetTree() << dec<<endl;  
 //      T_friend->GetTree()->Fill();//fill friends  
 //     }  
2625    
 // }  
Legend:
Removed from v.1.16  
changed lines
  Added in v.1.32
  ViewVC Help
Powered by ViewVC 1.1.23