/[PAMELA software]/PadmeAmidala/src/RunGlue.cpp
ViewVC logotype

Diff of /PadmeAmidala/src/RunGlue.cpp

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

revision 1.2 by mocchiut, Wed Dec 6 07:01:11 2006 UTC revision 1.17 by mocchiut, Wed Oct 14 09:12:12 2009 UTC
# Line 11  Line 11 
11  #include <TFile.h>  #include <TFile.h>
12  #include <TTree.h>  #include <TTree.h>
13  #include <TTimeStamp.h>  #include <TTimeStamp.h>
14    #include <TTreeCloner.h>
15  #include <RunGlue.h>  #include <RunGlue.h>
16  //  //
17  ClassImp(RunGlue);  ClassImp(RunGlue);
# Line 24  RunGlue::RunGlue() { Line 25  RunGlue::RunGlue() {
25  RunGlue::RunGlue(TSQLServer *da, UInt_t ru, TString di, TString wrkdi) {  RunGlue::RunGlue(TSQLServer *da, UInt_t ru, TString di, TString wrkdi) {
26    fDBG = false;    fDBG = false;
27    li = new PamLevel2();    li = new PamLevel2();
28      //  li = 0;
29    this->Clear();    this->Clear();
30    dbc = da;    dbc = da;
31    run = ru;    run = ru;
32    dir = di;    dir = di;
33    //    //
34    TString wd = gSystem->WorkingDirectory();      wd = gSystem->WorkingDirectory();  
35    //    //
36    outdir = wrkdi;    outdir = wrkdi;
37    if ( !strcmp(outdir.Data(),"") ) outdir = wd;    if ( !strcmp(outdir.Data(),"") ) outdir = wd;
# Line 40  RunGlue::RunGlue(TSQLServer *da, UInt_t Line 42  RunGlue::RunGlue(TSQLServer *da, UInt_t
42      runmode = false;      runmode = false;
43      createlist = true;      createlist = true;
44    };    };
45      fUpgrade = true;
46    //    //
47    //  lList = new TList();    //  lList = new TList();
48    //    //
49  };  };
50    
   
51  void RunGlue::Clear() {  void RunGlue::Clear() {
52    fEnd = false;    fEnd = false;
53    run = 0;    run = 0;
# Line 74  void RunGlue::SetDList(TString list) { Line 76  void RunGlue::SetDList(TString list) {
76    if (list.Contains("+ALL", TString::kIgnoreCase)) RUN = true;    if (list.Contains("+ALL", TString::kIgnoreCase)) RUN = true;
77    //    //
78    fDList = list;    fDList = list;
79    li->SetWhichTrees(fDList);    //  li->SetWhichTrees(fDList);
80    //    //
81    if ( DebugMode() ) printf(" Detector list is %s \n",fDList.Data());    if ( DebugMode() ) printf(" Detector list is %s \n",fDList.Data());
82    //    //
83  };  };
84    
85  TList *RunGlue::GetRunList(){  TList *RunGlue::GetRunList(){
86    //    //  
87    lList = new TList();    lList = new TList();
88    lList->Clear();    lList->Clear();
89    TString thisrun;    TString thisrun;
# Line 115  TList *RunGlue::GetRunList(){ Line 117  TList *RunGlue::GetRunList(){
117      UInt_t lowerlimit = llim->GetSec();      UInt_t lowerlimit = llim->GetSec();
118      UInt_t upperlimit = lowerlimit + 86401;      UInt_t upperlimit = lowerlimit + 86401;
119      //      //
     stringstream myquery;  
     myquery.str("");  
     myquery << "SELECT ID FROM GL_RUN WHERE RUNHEADER_TIME>=" << (UInt_t)lowerlimit  
             << " AND RUNHEADER_TIME<" << (UInt_t)upperlimit << " ORDER BY RUNHEADER_TIME ASC;";  
     //  
     pResult = dbc->Query(myquery.str().c_str());  
     for( UInt_t r=0; r < 1000; r++){  
       Row = pResult->Next();        
       if( Row == NULL ) break;  
       if ( DebugMode() ) printf(" %u RUN %s \n",r,Row->GetField(0));  
       //  
       thisrun=dir+(TString)Row->GetField(0)+".Level2.root";        
       //  
       if ( DebugMode() ) printf(" Filename is %s \n",thisrun.Data());  
       //  
       su = TFile::Open(thisrun);  
       if ( li->CheckLevel2File(thisrun) ){  
         lList->Add(su);  
       } else {  
         if ( DebugMode() ) printf(" RUN %s DISCARDED \n",Row->GetField(0));  
       };  
       //  
     };  
     //  
120      if ( DebugMode() ) printf(" YY %u MM %u DD %u ll %u ul %u \n",YY,MM,DD,lowerlimit,upperlimit);      if ( DebugMode() ) printf(" YY %u MM %u DD %u ll %u ul %u \n",YY,MM,DD,lowerlimit,upperlimit);
121      //      //
122      YY -= 2000;      YY -= 2000;
# Line 172  TList *RunGlue::GetRunList(){ Line 150  TList *RunGlue::GetRunList(){
150      //      //
151      printf(" Output filename is %s \n",fFilename.Data());      printf(" Output filename is %s \n",fFilename.Data());
152      //      //
153        if ( !this->OpenFile() ){
154             fEnd = true;
155             return(0);
156        };
157        //
158        stringstream myquery;
159        myquery.str("");
160        myquery << "SELECT ID FROM GL_RUN WHERE RUNHEADER_TIME>=" << (UInt_t)lowerlimit
161                << " AND RUNHEADER_TIME<" << (UInt_t)upperlimit << " ORDER BY RUNHEADER_TIME ASC;";
162        //
163        pResult = dbc->Query(myquery.str().c_str());
164        for( UInt_t r=0; r < 1000; r++){
165          Row = pResult->Next();      
166          if( Row == NULL ) break;
167          if ( DebugMode() ) printf(" %u RUN %s \n",r,Row->GetField(0));
168          //
169          thisrun=dir+(TString)Row->GetField(0)+".Level2.root";      
170          //
171          if ( DebugMode() ) printf(" Filename is %s \n",thisrun.Data());
172          //
173          su = TFile::Open(thisrun);
174          if ( li->CheckLevel2File(thisrun) ){
175            lList->Add(su);
176          } else {
177                    if ( DebugMode() ) printf(" RUN %s DISCARDED \n",Row->GetField(0));
178          };
179          //
180        };
181        //
182        //
183      delete glrun;      delete glrun;
184      delete dbtime;      delete dbtime;
185      //      //
# Line 276  TList *RunGlue::GetRunList(){ Line 284  TList *RunGlue::GetRunList(){
284      UInt_t lowerlimit = llim->GetSec();      UInt_t lowerlimit = llim->GetSec();
285      UInt_t upperlimit = lowerlimit + 86401;      UInt_t upperlimit = lowerlimit + 86401;
286      //      //
     stringstream myquery;  
     myquery.str("");  
     myquery << "SELECT ID FROM GL_RUN WHERE RUNHEADER_TIME>=" << (UInt_t)lowerlimit  
             << " AND RUNHEADER_TIME<" << (UInt_t)upperlimit << " ORDER BY RUNHEADER_TIME ASC;";  
     //  
     pResult = dbc->Query(myquery.str().c_str());  
     for( UInt_t r=0; r < 1000; r++){  
       Row = pResult->Next();        
       if( Row == NULL ) break;  
       if ( DebugMode() ) printf(" %u RUN %s \n",r,Row->GetField(0));  
       //  
       thisrun=dir+(TString)Row->GetField(0)+".Level2.root";        
       //  
       if ( DebugMode() ) printf(" Filename is %s \n",thisrun.Data());  
       //  
       su = TFile::Open(thisrun);  
       if ( su ){  
         TFile *su0 = TFile::Open(thisrun);  
         fDoneList->Add(su0);  
         fNlistdone++;  
       };  
       if ( li->CheckLevel2File(thisrun) ){  
         lList->Add(su);  
         if ( DebugMode() ) printf(" RUN %s ADDED \n",Row->GetField(0));  
       } else {  
         if ( DebugMode() ) printf(" RUN %s DISCARDED \n",Row->GetField(0));  
       };  
       //  
     };  
     //  
287      if ( DebugMode() ) printf(" YY %u MM %u DD %u ll %u ul %u \n",YY,MM,DD,lowerlimit,upperlimit);      if ( DebugMode() ) printf(" YY %u MM %u DD %u ll %u ul %u \n",YY,MM,DD,lowerlimit,upperlimit);
288      //      //
289      YY -= 2000;      YY -= 2000;
# Line 339  TList *RunGlue::GetRunList(){ Line 317  TList *RunGlue::GetRunList(){
317      //      //
318      printf(" Output filename is %s \n",fFilename.Data());      printf(" Output filename is %s \n",fFilename.Data());
319      //      //
320        if ( !this->OpenFile() ) return(0);
321        //
322        stringstream myquery;
323        myquery.str("");
324        myquery << "SELECT ID FROM GL_RUN WHERE RUNHEADER_TIME>=" << (UInt_t)lowerlimit
325                << " AND RUNHEADER_TIME<" << (UInt_t)upperlimit << " ORDER BY RUNHEADER_TIME ASC;";
326        //
327        pResult = dbc->Query(myquery.str().c_str());
328        for( UInt_t r=0; r < 1000; r++){
329          Row = pResult->Next();      
330          if( Row == NULL ) break;
331          if ( DebugMode() ) printf(" %u RUN %s \n",r,Row->GetField(0));
332          //
333          thisrun=dir+(TString)Row->GetField(0)+".Level2.root";      
334          //
335          if ( DebugMode() ) printf(" Filename is %s \n",thisrun.Data());
336          //
337          su = TFile::Open(thisrun);
338          if ( su ){
339            TFile *su0 = TFile::Open(thisrun);
340            fDoneList->Add(su0);
341            fNlistdone++;
342          };
343          if ( li->CheckLevel2File(thisrun) ){
344            lList->Add(su);
345            if ( DebugMode() ) printf(" RUN %s ADDED \n",Row->GetField(0));
346          } else {
347            if ( DebugMode() ) printf(" RUN %s DISCARDED \n",Row->GetField(0));
348          };
349          //
350        };
351        //
352      if ( DebugMode() ){      if ( DebugMode() ){
353        UInt_t ll = 0;        UInt_t ll = 0;
354        while ( (TFile*)lList->At(ll) ){        while ( (TFile*)lList->At(ll) ){
# Line 371  TList *RunGlue::GetRunList(){ Line 381  TList *RunGlue::GetRunList(){
381    
382  Bool_t RunGlue::OpenFile(){  Bool_t RunGlue::OpenFile(){
383    //    //
384      fOpen = false;
385    printf(" Check if output file already exists \n");    printf(" Check if output file already exists \n");
386    Target = TFile::Open((this->GetFilename()).Data(), "READ" );    ifstream myfile;
387    //    myfile.open((this->GetFilename()).Data());
388    if ( Target ){    if ( myfile ){
389      Target->Close();      //  Target = TFile::Open((this->GetFilename()).Data(), "READ" );
390        //
391        //  if ( Target ){
392        //    Target->Close();
393        myfile.close();
394      printf("Error opening target file, %s  already exist!\n",(this->GetFilename()).Data());      printf("Error opening target file, %s  already exist!\n",(this->GetFilename()).Data());
395      return(false);      return(false);
396    } else {    } else {
397      //      //
398        // try to lock tables and work on a db level to avoid multiple file simultaneous access
399        //
400        stringstream oss;
401        //  
402        if ( fUpgrade ){
403          oss.str("");
404          oss << "lock table GL_ROOT write;";
405          TSQLResult *result = 0;
406          TSQLRow *row = 0;
407          result = dbc->Query(oss.str().c_str());
408          oss.str("");
409          oss << "select ID from GL_ROOT where NAME='" << ((TString)gSystem->BaseName(this->GetFilename())).Data() << "';";
410          result = dbc->Query(oss.str().c_str());
411          row = result->Next();
412          if ( row ){
413            printf("Error opening target file (does %s exist? already in the DB? are the permissions ok?)\n",(this->GetFilename()).Data());
414            return(false);
415          } else {
416            oss.str("");
417            oss << "insert into GL_ROOT (PATH,NAME) values ('$PAM_L2','" << ((TString)gSystem->BaseName(this->GetFilename())).Data() << "');";      
418            result = dbc->Query(oss.str().c_str());
419          };
420          oss.str("");
421          oss << "unlock tables;";
422          result = dbc->Query(oss.str().c_str());
423        };
424        //    Target = TFile::Open((this->GetFilename()).Data(), "RECREATE" );
425        Target = TFile::Open((this->GetFilename()).Data(), "NEW" );    
426        //    
427      printf(" Output file does not exist, creating it\n");      printf(" Output file does not exist, creating it\n");
428      //      //
429      Long64_t maxsize = 10000000000LL;        Long64_t maxsize = 99900000000LL;  
430      //      //
     Target = TFile::Open((this->GetFilename()).Data(), "RECREATE" );  
     fastMethod = kTRUE;  
     //fastMethod = kFALSE;  
     //      
     //      
431      if ( !Target || Target->IsZombie()) {      if ( !Target || Target->IsZombie()) {
432        printf("Error opening target file (does %s exist?)\n",(this->GetFilename()).Data());        printf("Error opening target file (does %s exist? are the permissions ok?)\n",(this->GetFilename()).Data());
433        exit(1);        //      exit(1);
434          return(false);
435      }      }
436      //      //
437        Target->Write();
438        //
439        //fastMethod = kTRUE;
440        fastMethod = kFALSE;
441        //    
442        //
443      TTree::SetMaxTreeSize(maxsize);      TTree::SetMaxTreeSize(maxsize);
444        //    Target->SetCompressionLevel(2);
445        Target->SetCompressionLevel(3);
446        //Target->SetCompressionLevel(5);
447      //      //
448        fOpen = true;
449      return(true);      return(true);
450      //      //
451    };    };
# Line 455  void RunGlue::UpdateDB(TList *dlist){ Line 505  void RunGlue::UpdateDB(TList *dlist){
505    };    };
506    //    //
507    myquery.str("");    myquery.str("");
   myquery << "insert into GL_ROOT (ID_RAW,PATH,NAME) values (" << idr << ",'" << outdir.Data() << "','" << ((TString)gSystem->BaseName(this->GetFilename())).Data() << "');";  
   if ( DebugMode() ) printf(" query is %s \n",myquery.str().c_str());  
   //  
   pResult = dbc->Query(myquery.str().c_str());  
   if ( !pResult ){  
     printf(" ERROR WRITING ON DB!\n");  
     return;  
   };  
   myquery.str("");  
508    myquery << "select ID from GL_ROOT where NAME='" << ((TString)gSystem->BaseName(this->GetFilename())).Data() << "';";    myquery << "select ID from GL_ROOT where NAME='" << ((TString)gSystem->BaseName(this->GetFilename())).Data() << "';";
509    if ( DebugMode() ) printf(" query is %s \n",myquery.str().c_str());    if ( DebugMode() ) printf(" query is %s \n",myquery.str().c_str());
510    //    //
# Line 479  void RunGlue::UpdateDB(TList *dlist){ Line 520  void RunGlue::UpdateDB(TList *dlist){
520    } else {    } else {
521      idl2 = (UInt_t)atoll(Row->GetField(0));      idl2 = (UInt_t)atoll(Row->GetField(0));
522    };    };
523      myquery.str("");
524      //  myquery << "insert into GL_ROOT (ID_RAW,PATH,NAME) values (" << idr << ",'" << outdir.Data() << "','" << ((TString)gSystem->BaseName(this->GetFilename())).Data() << "');";
525      //  myquery << "insert into GL_ROOT (ID_RAW,PATH,NAME) values (" << idr << ",'$PAM_L2','" << ((TString)gSystem->BaseName(this->GetFilename())).Data() << "');";
526      myquery << "update GL_ROOT set ID_RAW=" << idr << " where ID=" << idl2 << ";";
527      if ( DebugMode() ) printf(" query is %s \n",myquery.str().c_str());
528      //
529      pResult = dbc->Query(myquery.str().c_str());
530      if ( !pResult ){
531        printf(" ERROR WRITING ON DB!\n");
532        return;
533      };
534    //    //
535    //    //
536    //    //
# Line 511  void RunGlue::UpdateDB(TList *dlist){ Line 563  void RunGlue::UpdateDB(TList *dlist){
563  //  //
564  void RunGlue::MergeRootfile(TList *sourcelist) {  void RunGlue::MergeRootfile(TList *sourcelist) {
565    //    //
566    merged = false;    merged = true;
   //  
   TDirectory *target = Target;  
567    //    //
568    if ( DebugMode() ) printf("\nTarget path is: %s \n",target->GetPath());    //  if ( li ) li->Delete();
569    //    PamLevel2 *nli = new PamLevel2(wd,sourcelist,fDList);
570    TString path( (char*)strstr( target->GetPath(), ":" ) );    nli->SetSELLI(2);
571    path.Remove(0,2);    //
572    //    Target->cd();
573    TDirectory *first_source = (TDirectory*)sourcelist->First();    //
574    THashList allNames;    nli->CreateCloneTrees(Target);
575    //    ULong64_t nevents = nli->GetEntries();        
576    while ( first_source ){    printf(" NEVENTS %llu \n",nevents);
577      //    for(ULong64_t iev=0; iev<nevents; iev++){
578      TDirectory *current_sourcedir = first_source->GetDirectory(path);      nli->Clear();
579      //      if( nli->GetEntry(iev) ){
580      if ( !current_sourcedir ){        nli->FillCloneTrees();
581        first_source = (TDirectory*)sourcelist->After(first_source);      };
582        continue;    };
583      }    Target->cd();
584      //    nli->WriteCloneTrees();
585      // loop over all keys in this directory    printf("Written file %s \n",Target->GetName());
586      //    //  Target->Write();
587      TChain *globChain = 0;    //  TTree *slist = (TTree*)Target->Get("SelectionList");
588      TIter nextkey( current_sourcedir->GetListOfKeys() );    //  slist->Delete("all");
589      TKey *key, *oldkey=0;    Target->Close();
590      TH1::AddDirectory(kFALSE); // gain time, do not add the objects in the list in memory  
591      //  //   //
592      while ( (key = (TKey*)nextkey()) ) {  //   TDirectory *target = Target;
593        //  //   //
594        if ( current_sourcedir == target ) break;  //   if ( DebugMode() ) printf("\nTarget path is: %s \n",target->GetPath());
595        //  //   //
596        if (oldkey && !strcmp(oldkey->GetName(),key->GetName())) continue; //keep only the highest cycle number for each key  //   TString path( (char*)strstr( target->GetPath(), ":" ) );
597        //  //   path.Remove(0,2);
598        if (allNames.FindObject(key->GetName())) continue;  //   //
599        //  //   TDirectory *first_source = (TDirectory*)sourcelist->First();
600        if ( DebugMode() ) printf(" Key name is -%s- \n",key->GetName());  //   THashList allNames;
601        //  //   //
602        //  //   while ( first_source ){
603        //  //     //
604        if ( !strcmp(key->GetName(),"Tracker") && !li->TRK2 && !li->TRK1 && !li->TRKh ) continue;  //     TDirectory *current_sourcedir = first_source->GetDirectory(path);
605        //  //     //
606        if ( !strcmp(key->GetName(),"Calorimeter") && !li->CAL2 && !li->CAL1 ) continue;  //     if ( !current_sourcedir ){
607        //  //       first_source = (TDirectory*)sourcelist->After(first_source);
608        if ( !strcmp(key->GetName(),"ToF") && !li->TOF ) continue;  //       continue;
609        //  //     }
610        if ( !strcmp(key->GetName(),"Trigger") && !li->TRG ) continue;  //     //
611        //  //     // loop over all keys in this directory
612        if ( !strcmp(key->GetName(),"Anticounter") && !li->AC ) continue;  //     //
613        //  //     TChain *globChain = 0;    
614        if ( !strcmp(key->GetName(),"S4") && !li->S4 ) continue;  //     TIter nextkey( current_sourcedir->GetListOfKeys() );
615        //  //     TKey *key = 0;
616        if ( !strcmp(key->GetName(),"NeutronD") && !li->ND ) continue;  //     TKey *oldkey = 0;
617        //  //     TH1::AddDirectory(kFALSE); // gain time, do not add the objects in the list in memory
618        if ( !strcmp(key->GetName(),"OrbitalInfo") && !li->ORB ) continue;  //     //
619        //  //     while ( (key = (TKey*)nextkey()) ) {
620        if ( !strcmp(key->GetName(),"Run") && !RUN ) continue;  //       //
621        //  //       //      printf(" target ls \n");
622        if ( strcmp(key->GetName(),"Calorimeter") && strcmp(key->GetName(),"Tracker") && strcmp(key->GetName(),"ToF") && strcmp(key->GetName(),"Trigger") && strcmp(key->GetName(),"Anticounter") && strcmp(key->GetName(),"S4") && strcmp(key->GetName(),"NeutronD") && strcmp(key->GetName(),"OrbitalInfo") && strcmp(key->GetName(),"Run") && strcmp(key->GetName(),"ProcessID0") ){  //       //      Target->ls();
623          if ( DebugMode() ) printf(" ERROR UNKNOWN KEY %s !\n",key->GetName());  //       //
624          continue;  //       if ( current_sourcedir == target ) break;
625        };  //       //
626        //  //       if ( oldkey && !strcmp(oldkey->GetName(),key->GetName()) ) continue; //keep only the highest cycle number for each key
627        allNames.Add(new TObjString(key->GetName()));  //       //
628        //  //       if ( allNames.FindObject(key->GetName()) ) continue;
629        // read object from first source file  //       //
630        //  //       if ( DebugMode() ) printf(" Key name is -%s- \n",key->GetName());
631        current_sourcedir->cd();  //       //
632        TObject *obj = key->ReadObj();  //       //
633        //  //       //
634        if ( obj->IsA()->InheritsFrom("TTree") ){  //       if ( !strcmp(key->GetName(),"Tracker") && !li->IsTRK2() && !li->IsTRK1() && !li->IsTRKh() ) continue;
635          //  //       //
636          // loop over all source files create a chain of Trees "globChain"  //       if ( !strcmp(key->GetName(),"Calorimeter") && !li->IsCAL2() && !li->IsCAL1() ) continue;
637          //  //       //
638          TString obj_name;  //       if ( !strcmp(key->GetName(),"ToF") && !li->IsTOF() ) continue;
639          //  //       //
640          if ( path.Length() ) {  //       if ( !strcmp(key->GetName(),"Trigger") && !li->IsTRG() ) continue;
641            obj_name = path + "/" + obj->GetName();  //       //
642          } else {  //       if ( !strcmp(key->GetName(),"Anticounter") && !li->IsAC() ) continue;
643            obj_name = obj->GetName();  //       //
644          };  //       if ( !strcmp(key->GetName(),"S4") && !li->IsS4() ) continue;
645          //  //       //
646          globChain = new TChain(obj_name);  //       if ( !strcmp(key->GetName(),"NeutronD") && !li->IsND() ) continue;
647          //  //       //
648          globChain->Add(first_source->GetName());  //       if ( !strcmp(key->GetName(),"OrbitalInfo") && !li->IsORB() ) continue;
649          //  //       //
650          TFile *nextsource = (TFile*)sourcelist->After( first_source );  //       if ( !strcmp(key->GetName(),"Run") && !li->IsRUN() ) continue;
651          //  //       //
652          while ( nextsource ) {              //       if ( strcmp(key->GetName(),"Calorimeter") && strcmp(key->GetName(),"Tracker") && strcmp(key->GetName(),"ToF") && strcmp(key->GetName(),"Trigger") && strcmp(key->GetName(),"Anticounter") && strcmp(key->GetName(),"S4") && strcmp(key->GetName(),"NeutronD") && strcmp(key->GetName(),"OrbitalInfo") && strcmp(key->GetName(),"Run") && strcmp(key->GetName(),"ProcessID0") ){
653            //  //      if ( DebugMode() ) printf(" ERROR UNKNOWN KEY %s !\n",key->GetName());
654            //do not add to the list a file that does not contain this Tree  //      continue;
655            //  //       };
656            TFile *curf = TFile::Open(nextsource->GetName());  //       //
657            //  //       printf(" 1 \n");
658            if ( curf ) {  //       allNames.Add(new TObjString(key->GetName()));
659              //  //       printf(" 2 \n");
660              Bool_t mustAdd = kFALSE;  //       //
661              //  //       // read object from first source file
662              if (curf->FindKey(obj_name)) {  //       //
663                //  //       current_sourcedir->cd();
664                mustAdd = kTRUE;  //       TObject *obj = key->ReadObj();
665                //  //       printf(" 3 \n");
666              } else {  //       //
667                //  //       if ( obj->IsA()->InheritsFrom("TTree") ){
668                //we could be more clever here. No need to import the object  //      //
669                //we are missing a function in TDirectory  //              // loop over all source files create a chain of Trees "globChain"
670                //  //      //
671                TObject *aobj = curf->Get(obj_name);  //      TString obj_name;
672                //  //      printf(" 4 \n");
673                if ( aobj ){  //      //
674                  mustAdd = kTRUE;  //      if ( path.Length() ) {
675                  delete aobj;  //        obj_name = path + "/" + obj->GetName();
676                };  //      } else {
677              };  //        obj_name = obj->GetName();
678              if (mustAdd) {  //      };
679                globChain->Add(nextsource->GetName());  //      //
680              };  //      globChain = new TChain(obj_name);
681            };  //      //
682            delete curf;  //      globChain->SetCacheSize(0);
683            nextsource = (TFile*)sourcelist->After(nextsource);  //      //
684          };  //      printf(" 5 \n");
685          //  //      globChain->Add(first_source->GetName());
686        } else {  //      printf(" 6 \n");
687          //  //      //
688          // object is of no type that we know or can handle  //      TFile *nextsource = (TFile*)sourcelist->After( first_source );
689          //  //      //
690          if ( DebugMode() ) cout << "Unknown object type, name: "  //      while ( nextsource ) {            
691               << obj->GetName() << " title: " << obj->GetTitle() << endl;  //        //
692        };  //        //do not add to the list a file that does not contain this Tree
693        //  //        //
694        // now write the merged histogram (which is "in" obj) to the target file  //        printf(" 7 \n");
695        // note that this will just store obj in the current directory level,  //        TFile *curf = TFile::Open(nextsource->GetName());
696        // which is not persistent until the complete directory itself is stored  //        //
697        // by "target->Write()" below  //        if ( curf ) {
698        //  //          //
699        if ( obj ){  //          Bool_t mustAdd = kFALSE;
700          //  //          //
701          target->cd();  //          if (curf->FindKey(obj_name)) {
702          //  //            //
703          if( obj->IsA()->InheritsFrom("TTree") ) {  //            mustAdd = kTRUE;
704            //  //            //
705            //  //          } else {
706            if ( fastMethod ){  //            //
707              globChain->Merge(target->GetFile(),0,"keep fast");  //            //we could be more clever here. No need to import the object
708            } else {  //            //we are missing a function in TDirectory
709              globChain->Merge(target->GetFile(),0,"keep");  //            //
710            };  //            TObject *aobj = curf->Get(obj_name);
711            //  //            //
712            merged = true;  //            if ( aobj ){
713            //  //              mustAdd = kTRUE;
714            if ( DebugMode() ) globChain->ls();  //              delete aobj;
715            //  //            };
716            delete globChain;  //          };
717            //  //          if (mustAdd) {
718          } else {  //            printf(" 8 \n");
719            //  //            globChain->Add(nextsource->GetName());
720            //      obj->Write( key->GetName() ); // <==================================================================================  //            printf(" 9 \n");
721            //  //          };
722          };  //        };
723        };  //        delete curf;
724        oldkey = key;  //        nextsource = (TFile*)sourcelist->After(nextsource);
725      };  //      };
726      //  //      //
727      first_source = (TDirectory*)sourcelist->After(first_source);  //      printf(" 10 \n");      
728      //  //      delete nextsource;
729    };  //      printf(" 11 \n");
730    // save modifications to target file  //      //
731    target->SaveSelf(kTRUE);  //       } else {
732    //  //      //
733  }  //      // object is of no type that we know or can handle
734    //      //
735    //      if ( DebugMode() ) cout << "Unknown object type, name: "
736    //           << obj->GetName() << " title: " << obj->GetTitle() << endl;
737    //       };
738    //       //
739    //       // now write the merged histogram (which is "in" obj) to the target file
740    //       // note that this will just store obj in the current directory level,
741    //       // which is not persistent until the complete directory itself is stored
742    //       // by "target->Write()" below
743    //       //
744    //       printf(" 12 \n");
745    //       if ( obj ){
746    //      //
747    //      target->cd();
748    //      //
749    //      if( obj->IsA()->InheritsFrom("TTree") ) {
750    //        //
751    //        Long64_t nfiles = 0;
752    //        if ( fastMethod ){
753    //          //      globChain->Merge(target->GetFile(),0,"C keep fast");
754    //          nfiles=this->Mergy((TChain*)globChain,target->GetFile(),0,"C keep fast");
755    //        } else {
756    //          //globChain->Merge(target->GetFile(),0,"C keep");
757    //          printf(" 13 %llu \n",globChain->GetEntries());          
758    //          nfiles=this->Mergy((TChain*)globChain,target->GetFile(),0,"C keep");
759    //          printf(" 14 \n");
760    //        };
761    //        //
762    //        merged = true;
763    //        //
764    //        if ( DebugMode() ){
765    //          printf(" Merged %i files\n",(int)nfiles);
766    //          globChain->ls();
767    //        };
768    //        //
769    //        delete globChain;
770    //        //      return;
771    //        //
772    //      } else {
773    //        //
774    //        //      obj->Write( key->GetName() ); // <==================================================================================
775    //        //
776    //      };
777    //       };
778    //       delete obj;
779    //       oldkey = key;
780    //     };
781    //     //
782    //     first_source = (TDirectory*)sourcelist->After(first_source);
783    //     //
784    //   };
785    //   // save modifications to target file
786    //   //
787    //   target->SaveSelf(kTRUE);
788    //   //
789    };
790    
791    
792    // Long64_t RunGlue::Mergy(TChain *mychain, TFile* file, Int_t basketsize, Option_t* option){
793    //    // We must have been passed a file, we will use it
794    //    // later to reset the compression level of the branches.
795    //    if (!file) {
796    //       // FIXME: We need an error message here.
797    //      printf(" 19 \n");
798    //      return 0;
799    //    }
800    //    printf(" 20 \n");
801    
802    //    // Options
803    //    Bool_t fastClone = kFALSE;
804    //    TString opt = option;
805    //    opt.ToLower();
806    //    if (opt.Contains("fast")) {
807    //       fastClone = kTRUE;
808    //       printf(" 21 \n");
809    //    }
810    //    // The chain tree must have a list of branches
811    //    // because we may try to change their basket
812    //    // size later.
813    //    TObjArray* lbranches = mychain->GetListOfBranches();
814    //    if (!lbranches) {
815    //       // FIXME: We need an error message here.
816    //      printf(" 22 \n");
817    //       return 0;
818    //    }
819    
820    //    //   file->cd();
821    //    // The chain must have a current tree because
822    //    // that is the one we will clone.
823    //    //   if (!fTree) {
824    //       // -- LoadTree() has not yet been called, no current tree.
825    //       // FIXME: We need an error message here.
826    //    //      return 0;
827    //    //   }
828    
829    //    // Copy the chain's current tree without
830    //    // copying any entries, we will do that later.
831    //    printf(" 23 \n");
832    //    TTree* newTree = mychain->CloneTree(0);
833    //    if (!newTree) {
834    //       // FIXME: We need an error message here.
835    //      printf(" 24 \n");
836    //       return 0;
837    //    }
838    //    printf(" 25 \n");
839    
840    
841    //    // Strip out the (potential) directory name.
842    //    // FIXME: The merged chain may or may not have the
843    //    //        same name as the original chain.  This is
844    //    //        bad because the chain name determines the
845    //    //        names of the trees in the chain by default.
846    //    printf(" 26 \n");
847    //    newTree->SetName(gSystem->BaseName(mychain->GetName()));
848    //    printf(" 27 \n");
849    
850    //    // FIXME: Why do we do this?
851    //    //   newTree->SetAutoSave(-1);
852    //    newTree->SetAutoSave(900000000000000LL);
853    
854    //    printf(" 28 \n");
855    //    // Circularity is incompatible with merging, it may
856    //    // force us to throw away entries, which is not what
857    //    // we are supposed to do.
858    //    newTree->SetCircular(0);
859    
860    //    // Reset the compression level of the branches.
861    //    if (opt.Contains("c")) {
862    //       TBranch* branch = 0;
863    //       TIter nextb(newTree->GetListOfBranches());
864    //       while ((branch = (TBranch*) nextb())) {
865    //      printf(" 29 \n");
866    //          branch->SetCompressionLevel(file->GetCompressionLevel());
867    //       }
868    //    }
869    
870    //    printf(" 30 \n");
871    //    // Reset the basket size of the branches.
872    //    if (basketsize > 1000) {
873    //       TBranch* branch = 0;
874    //       TIter nextb(newTree->GetListOfBranches());
875    //       while ((branch = (TBranch*) nextb())) {
876    //      printf(" 31 \n");
877    //          branch->SetBasketSize(basketsize);
878    //       }
879    //    }
880    
881    //    printf(" 32 \n");
882    //    Long64_t nentries = mychain->GetEntriesFast();
883    //    //Long64_t nentries = mychain->GetEntries();
884    
885    //    // Copy the entries.
886    //    if (fastClone) {
887    //       // For each tree in the chain.
888    //       for (Long64_t i = 0; i < nentries; i += mychain->GetTree()->GetEntries()) {
889    //          if (mychain->LoadTree(i) < 0) {
890    //             break;
891    //          }
892    //          TTreeCloner cloner(mychain->GetTree(), newTree, option);
893    //          if (cloner.IsValid()) {
894    //             newTree->SetEntries(newTree->GetEntries() + mychain->GetTree()->GetEntries());
895    //             cloner.Exec();
896    //          } else {
897    //             if (mychain->GetFile()) {
898    //            printf("Merge Skipped file %s\n", mychain->GetFile()->GetName());
899    //             //            } else {
900    //             //               Warning("Merge", "Skipped file number %d\n", fTreeNumber);
901    //             }
902    //          }
903    //       }
904    //    } else {
905    //      printf(" 33 %llu \n",nentries);
906    //      mychain->Print();
907    //      for (Long64_t i = 0; i < nentries; i++) {
908    //        printf(" i %llu \n",i);
909    //        // for (Long64_t i = 0; i < 1; i++) {
910    //          if (mychain->GetEntry(i) <= 0) {
911    //             break;
912    //          }
913    //          newTree->Fill();
914    //       }
915    //       printf(" 34 \n");
916    //    }
917    
918    //    // Write the new tree header.
919    //    printf(" 35 \n");
920    //    newTree->Write();
921    
922    //    printf(" 36 \n");
923    //    // Get our return value.
924    //    Int_t nfiles = newTree->GetFileNumber() + 1;
925    
926    //    // Close and delete the current file of the new tree.
927    //    if (!opt.Contains("keep")) {
928    //       // FIXME: What happens to fDirectory in newTree here?
929    //       delete newTree->GetCurrentFile();
930    //    }
931    //    return nfiles;
932    // }
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.17
  ViewVC Help
Powered by ViewVC 1.1.23