/[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.9 by mocchiut, Fri May 4 10:05:38 2007 UTC revision 1.15 by mocchiut, Fri Jan 30 09:37:18 2009 UTC
# Line 25  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 46  RunGlue::RunGlue(TSQLServer *da, UInt_t Line 47  RunGlue::RunGlue(TSQLServer *da, UInt_t
47    //    //
48  };  };
49    
   
50  void RunGlue::Clear() {  void RunGlue::Clear() {
51    fEnd = false;    fEnd = false;
52    run = 0;    run = 0;
# Line 75  void RunGlue::SetDList(TString list) { Line 75  void RunGlue::SetDList(TString list) {
75    if (list.Contains("+ALL", TString::kIgnoreCase)) RUN = true;    if (list.Contains("+ALL", TString::kIgnoreCase)) RUN = true;
76    //    //
77    fDList = list;    fDList = list;
78    li->SetWhichTrees(fDList);    //  li->SetWhichTrees(fDList);
79    //    //
80    if ( DebugMode() ) printf(" Detector list is %s \n",fDList.Data());    if ( DebugMode() ) printf(" Detector list is %s \n",fDList.Data());
81    //    //
82  };  };
83    
84  TList *RunGlue::GetRunList(){  TList *RunGlue::GetRunList(){
85    //    //  
86    lList = new TList();    lList = new TList();
87    lList->Clear();    lList->Clear();
88    TString thisrun;    TString thisrun;
# Line 173  TList *RunGlue::GetRunList(){ Line 173  TList *RunGlue::GetRunList(){
173        if ( li->CheckLevel2File(thisrun) ){        if ( li->CheckLevel2File(thisrun) ){
174          lList->Add(su);          lList->Add(su);
175        } else {        } else {
176          if ( DebugMode() ) printf(" RUN %s DISCARDED \n",Row->GetField(0));                  if ( DebugMode() ) printf(" RUN %s DISCARDED \n",Row->GetField(0));
177        };        };
178        //        //
179      };      };
# Line 394  Bool_t RunGlue::OpenFile(){ Line 394  Bool_t RunGlue::OpenFile(){
394      return(false);      return(false);
395    } else {    } else {
396      //      //
397        //    Target = TFile::Open((this->GetFilename()).Data(), "RECREATE" );
398        Target = TFile::Open((this->GetFilename()).Data(), "NEW" );    
399        //    
400      printf(" Output file does not exist, creating it\n");      printf(" Output file does not exist, creating it\n");
401      //      //
402      Long64_t maxsize = 99900000000LL;        Long64_t maxsize = 99900000000LL;  
403      //      //
404      Target = TFile::Open((this->GetFilename()).Data(), "RECREATE" );      if ( !Target || Target->IsZombie()) {
405          printf("Error opening target file (does %s exist? are the permissions ok?)\n",(this->GetFilename()).Data());
406          //      exit(1);
407          return(false);
408        }
409        //
410        Target->Write();
411        //
412      //fastMethod = kTRUE;      //fastMethod = kTRUE;
413      fastMethod = kFALSE;      fastMethod = kFALSE;
414      //          //    
     //      
     if ( !Target || Target->IsZombie()) {  
       printf("Error opening target file (does %s exist?)\n",(this->GetFilename()).Data());  
       exit(1);  
     }  
415      //      //
416      TTree::SetMaxTreeSize(maxsize);      TTree::SetMaxTreeSize(maxsize);
417      Target->SetCompressionLevel(2);      //    Target->SetCompressionLevel(2);
418        Target->SetCompressionLevel(3);
419        //Target->SetCompressionLevel(5);
420      //      //
421      fOpen = true;      fOpen = true;
422      return(true);      return(true);
# Line 528  void RunGlue::UpdateDB(TList *dlist){ Line 535  void RunGlue::UpdateDB(TList *dlist){
535  //  //
536  void RunGlue::MergeRootfile(TList *sourcelist) {  void RunGlue::MergeRootfile(TList *sourcelist) {
537    //    //
538    merged = false;    merged = true;
   //  
   TDirectory *target = Target;  
   //  
   if ( DebugMode() ) printf("\nTarget path is: %s \n",target->GetPath());  
   //  
   TString path( (char*)strstr( target->GetPath(), ":" ) );  
   path.Remove(0,2);  
   //  
   TDirectory *first_source = (TDirectory*)sourcelist->First();  
   THashList allNames;  
   //  
   while ( first_source ){  
     //  
     TDirectory *current_sourcedir = first_source->GetDirectory(path);  
     //  
     if ( !current_sourcedir ){  
       first_source = (TDirectory*)sourcelist->After(first_source);  
       continue;  
     }  
     //  
     // loop over all keys in this directory  
     //  
     TChain *globChain = 0;      
     TIter nextkey( current_sourcedir->GetListOfKeys() );  
     TKey *key = 0;  
     TKey *oldkey = 0;  
     TH1::AddDirectory(kFALSE); // gain time, do not add the objects in the list in memory  
     //  
     while ( (key = (TKey*)nextkey()) ) {  
       //  
       //      printf(" target ls \n");  
       //      Target->ls();  
       //  
       if ( current_sourcedir == target ) break;  
       //  
       if ( oldkey && !strcmp(oldkey->GetName(),key->GetName()) ) continue; //keep only the highest cycle number for each key  
       //  
       if ( allNames.FindObject(key->GetName()) ) continue;  
       //  
       if ( DebugMode() ) printf(" Key name is -%s- \n",key->GetName());  
       //  
       //  
       //  
       if ( !strcmp(key->GetName(),"Tracker") && !li->IsTRK2() && !li->IsTRK1() && !li->IsTRKh() ) continue;  
       //  
       if ( !strcmp(key->GetName(),"Calorimeter") && !li->IsCAL2() && !li->IsCAL1() ) continue;  
       //  
       if ( !strcmp(key->GetName(),"ToF") && !li->IsTOF() ) continue;  
       //  
       if ( !strcmp(key->GetName(),"Trigger") && !li->IsTRG() ) continue;  
       //  
       if ( !strcmp(key->GetName(),"Anticounter") && !li->IsAC() ) continue;  
       //  
       if ( !strcmp(key->GetName(),"S4") && !li->IsS4() ) continue;  
       //  
       if ( !strcmp(key->GetName(),"NeutronD") && !li->IsND() ) continue;  
       //  
       if ( !strcmp(key->GetName(),"OrbitalInfo") && !li->IsORB() ) continue;  
       //  
       if ( !strcmp(key->GetName(),"Run") && !li->IsRUN() ) continue;  
       //  
       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") ){  
         if ( DebugMode() ) printf(" ERROR UNKNOWN KEY %s !\n",key->GetName());  
         continue;  
       };  
       //  
       allNames.Add(new TObjString(key->GetName()));  
       //  
       // read object from first source file  
       //  
       current_sourcedir->cd();  
       TObject *obj = key->ReadObj();  
       //  
       if ( obj->IsA()->InheritsFrom("TTree") ){  
         //  
         // loop over all source files create a chain of Trees "globChain"  
         //  
         TString obj_name;  
         //  
         if ( path.Length() ) {  
           obj_name = path + "/" + obj->GetName();  
         } else {  
           obj_name = obj->GetName();  
         };  
         //  
         globChain = new TChain(obj_name);  
         //  
         globChain->SetCacheSize(0);  
         //  
         globChain->Add(first_source->GetName());  
         //  
         TFile *nextsource = (TFile*)sourcelist->After( first_source );  
         //  
         while ( nextsource ) {              
           //  
           //do not add to the list a file that does not contain this Tree  
           //  
           TFile *curf = TFile::Open(nextsource->GetName());  
           //  
           if ( curf ) {  
             //  
             Bool_t mustAdd = kFALSE;  
             //  
             if (curf->FindKey(obj_name)) {  
               //  
               mustAdd = kTRUE;  
               //  
             } else {  
               //  
               //we could be more clever here. No need to import the object  
               //we are missing a function in TDirectory  
               //  
               TObject *aobj = curf->Get(obj_name);  
               //  
               if ( aobj ){  
                 mustAdd = kTRUE;  
                 delete aobj;  
               };  
             };  
             if (mustAdd) {  
               globChain->Add(nextsource->GetName());  
             };  
           };  
           delete curf;  
           nextsource = (TFile*)sourcelist->After(nextsource);  
         };  
         //  
         delete nextsource;  
         //  
       } else {  
         //  
         // object is of no type that we know or can handle  
         //  
         if ( DebugMode() ) cout << "Unknown object type, name: "  
              << obj->GetName() << " title: " << obj->GetTitle() << endl;  
       };  
       //  
       // now write the merged histogram (which is "in" obj) to the target file  
       // note that this will just store obj in the current directory level,  
       // which is not persistent until the complete directory itself is stored  
       // by "target->Write()" below  
       //  
       if ( obj ){  
         //  
         target->cd();  
         //  
         if( obj->IsA()->InheritsFrom("TTree") ) {  
           //  
           Long64_t nfiles = 0;  
           if ( fastMethod ){  
             //      globChain->Merge(target->GetFile(),0,"C keep fast");  
             nfiles=this->Mergy((TChain*)globChain,target->GetFile(),0,"C keep fast");  
           } else {  
             //globChain->Merge(target->GetFile(),0,"C keep");  
             nfiles=this->Mergy((TChain*)globChain,target->GetFile(),0,"C keep");  
           };  
           //  
           merged = true;  
           //  
           if ( DebugMode() ){  
             printf(" Merged %i files\n",(int)nfiles);  
             globChain->ls();  
           };  
           //  
           delete globChain;  
           //      return;  
           //  
         } else {  
           //  
           //      obj->Write( key->GetName() ); // <==================================================================================  
           //  
         };  
       };  
       delete obj;  
       oldkey = key;  
     };  
     //  
     first_source = (TDirectory*)sourcelist->After(first_source);  
     //  
   };  
   // save modifications to target file  
   //  
   target->SaveSelf(kTRUE);  
539    //    //
540      if ( li ) li->Delete();
541      li = new PamLevel2(wd,sourcelist,fDList);
542      li->SetSELLI(2);
543      //
544      Target->cd();
545      //
546      li->CreateCloneTrees(Target);
547      ULong64_t nevents = li->GetEntries();
548      printf(" NEVENTS %llu \n",nevents);
549      for(ULong64_t iev=0; iev<nevents; iev++){
550        li->Clear();
551        if( li->GetEntry(iev) ){
552          li->FillCloneTrees();
553        };
554      };
555      Target->cd();
556      li->WriteCloneTrees();
557      printf("Written file %s \n",Target->GetName());
558      //  Target->Write();
559      //  TTree *slist = (TTree*)Target->Get("SelectionList");
560      //  slist->Delete("all");
561      Target->Close();
562    
563    //   //
564    //   TDirectory *target = Target;
565    //   //
566    //   if ( DebugMode() ) printf("\nTarget path is: %s \n",target->GetPath());
567    //   //
568    //   TString path( (char*)strstr( target->GetPath(), ":" ) );
569    //   path.Remove(0,2);
570    //   //
571    //   TDirectory *first_source = (TDirectory*)sourcelist->First();
572    //   THashList allNames;
573    //   //
574    //   while ( first_source ){
575    //     //
576    //     TDirectory *current_sourcedir = first_source->GetDirectory(path);
577    //     //
578    //     if ( !current_sourcedir ){
579    //       first_source = (TDirectory*)sourcelist->After(first_source);
580    //       continue;
581    //     }
582    //     //
583    //     // loop over all keys in this directory
584    //     //
585    //     TChain *globChain = 0;    
586    //     TIter nextkey( current_sourcedir->GetListOfKeys() );
587    //     TKey *key = 0;
588    //     TKey *oldkey = 0;
589    //     TH1::AddDirectory(kFALSE); // gain time, do not add the objects in the list in memory
590    //     //
591    //     while ( (key = (TKey*)nextkey()) ) {
592    //       //
593    //       //      printf(" target ls \n");
594    //       //      Target->ls();
595    //       //
596    //       if ( current_sourcedir == target ) break;
597    //       //
598    //       if ( oldkey && !strcmp(oldkey->GetName(),key->GetName()) ) continue; //keep only the highest cycle number for each key
599    //       //
600    //       if ( allNames.FindObject(key->GetName()) ) continue;
601    //       //
602    //       if ( DebugMode() ) printf(" Key name is -%s- \n",key->GetName());
603    //       //
604    //       //
605    //       //
606    //       if ( !strcmp(key->GetName(),"Tracker") && !li->IsTRK2() && !li->IsTRK1() && !li->IsTRKh() ) continue;
607    //       //
608    //       if ( !strcmp(key->GetName(),"Calorimeter") && !li->IsCAL2() && !li->IsCAL1() ) continue;
609    //       //
610    //       if ( !strcmp(key->GetName(),"ToF") && !li->IsTOF() ) continue;
611    //       //
612    //       if ( !strcmp(key->GetName(),"Trigger") && !li->IsTRG() ) continue;
613    //       //
614    //       if ( !strcmp(key->GetName(),"Anticounter") && !li->IsAC() ) continue;
615    //       //
616    //       if ( !strcmp(key->GetName(),"S4") && !li->IsS4() ) continue;
617    //       //
618    //       if ( !strcmp(key->GetName(),"NeutronD") && !li->IsND() ) continue;
619    //       //
620    //       if ( !strcmp(key->GetName(),"OrbitalInfo") && !li->IsORB() ) continue;
621    //       //
622    //       if ( !strcmp(key->GetName(),"Run") && !li->IsRUN() ) continue;
623    //       //
624    //       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") ){
625    //      if ( DebugMode() ) printf(" ERROR UNKNOWN KEY %s !\n",key->GetName());
626    //      continue;
627    //       };
628    //       //
629    //       printf(" 1 \n");
630    //       allNames.Add(new TObjString(key->GetName()));
631    //       printf(" 2 \n");
632    //       //
633    //       // read object from first source file
634    //       //
635    //       current_sourcedir->cd();
636    //       TObject *obj = key->ReadObj();
637    //       printf(" 3 \n");
638    //       //
639    //       if ( obj->IsA()->InheritsFrom("TTree") ){
640    //      //
641    //              // loop over all source files create a chain of Trees "globChain"
642    //      //
643    //      TString obj_name;
644    //      printf(" 4 \n");
645    //      //
646    //      if ( path.Length() ) {
647    //        obj_name = path + "/" + obj->GetName();
648    //      } else {
649    //        obj_name = obj->GetName();
650    //      };
651    //      //
652    //      globChain = new TChain(obj_name);
653    //      //
654    //      globChain->SetCacheSize(0);
655    //      //
656    //      printf(" 5 \n");
657    //      globChain->Add(first_source->GetName());
658    //      printf(" 6 \n");
659    //      //
660    //      TFile *nextsource = (TFile*)sourcelist->After( first_source );
661    //      //
662    //      while ( nextsource ) {            
663    //        //
664    //        //do not add to the list a file that does not contain this Tree
665    //        //
666    //        printf(" 7 \n");
667    //        TFile *curf = TFile::Open(nextsource->GetName());
668    //        //
669    //        if ( curf ) {
670    //          //
671    //          Bool_t mustAdd = kFALSE;
672    //          //
673    //          if (curf->FindKey(obj_name)) {
674    //            //
675    //            mustAdd = kTRUE;
676    //            //
677    //          } else {
678    //            //
679    //            //we could be more clever here. No need to import the object
680    //            //we are missing a function in TDirectory
681    //            //
682    //            TObject *aobj = curf->Get(obj_name);
683    //            //
684    //            if ( aobj ){
685    //              mustAdd = kTRUE;
686    //              delete aobj;
687    //            };
688    //          };
689    //          if (mustAdd) {
690    //            printf(" 8 \n");
691    //            globChain->Add(nextsource->GetName());
692    //            printf(" 9 \n");
693    //          };
694    //        };
695    //        delete curf;
696    //        nextsource = (TFile*)sourcelist->After(nextsource);
697    //      };
698    //      //
699    //      printf(" 10 \n");      
700    //      delete nextsource;
701    //      printf(" 11 \n");
702    //      //
703    //       } else {
704    //      //
705    //      // object is of no type that we know or can handle
706    //      //
707    //      if ( DebugMode() ) cout << "Unknown object type, name: "
708    //           << obj->GetName() << " title: " << obj->GetTitle() << endl;
709    //       };
710    //       //
711    //       // now write the merged histogram (which is "in" obj) to the target file
712    //       // note that this will just store obj in the current directory level,
713    //       // which is not persistent until the complete directory itself is stored
714    //       // by "target->Write()" below
715    //       //
716    //       printf(" 12 \n");
717    //       if ( obj ){
718    //      //
719    //      target->cd();
720    //      //
721    //      if( obj->IsA()->InheritsFrom("TTree") ) {
722    //        //
723    //        Long64_t nfiles = 0;
724    //        if ( fastMethod ){
725    //          //      globChain->Merge(target->GetFile(),0,"C keep fast");
726    //          nfiles=this->Mergy((TChain*)globChain,target->GetFile(),0,"C keep fast");
727    //        } else {
728    //          //globChain->Merge(target->GetFile(),0,"C keep");
729    //          printf(" 13 %llu \n",globChain->GetEntries());          
730    //          nfiles=this->Mergy((TChain*)globChain,target->GetFile(),0,"C keep");
731    //          printf(" 14 \n");
732    //        };
733    //        //
734    //        merged = true;
735    //        //
736    //        if ( DebugMode() ){
737    //          printf(" Merged %i files\n",(int)nfiles);
738    //          globChain->ls();
739    //        };
740    //        //
741    //        delete globChain;
742    //        //      return;
743    //        //
744    //      } else {
745    //        //
746    //        //      obj->Write( key->GetName() ); // <==================================================================================
747    //        //
748    //      };
749    //       };
750    //       delete obj;
751    //       oldkey = key;
752    //     };
753    //     //
754    //     first_source = (TDirectory*)sourcelist->After(first_source);
755    //     //
756    //   };
757    //   // save modifications to target file
758    //   //
759    //   target->SaveSelf(kTRUE);
760    //   //
761  };  };
762    
763    
764  Long64_t RunGlue::Mergy(TChain *mychain, TFile* file, Int_t basketsize, Option_t* option){  // Long64_t RunGlue::Mergy(TChain *mychain, TFile* file, Int_t basketsize, Option_t* option){
765     // We must have been passed a file, we will use it  //    // We must have been passed a file, we will use it
766     // later to reset the compression level of the branches.  //    // later to reset the compression level of the branches.
767     if (!file) {  //    if (!file) {
768        // FIXME: We need an error message here.  //       // FIXME: We need an error message here.
769       return 0;  //      printf(" 19 \n");
770     }  //      return 0;
771    //    }
772     // Options  //    printf(" 20 \n");
773     Bool_t fastClone = kFALSE;  
774     TString opt = option;  //    // Options
775     opt.ToLower();  //    Bool_t fastClone = kFALSE;
776     if (opt.Contains("fast")) {  //    TString opt = option;
777        fastClone = kTRUE;  //    opt.ToLower();
778     }  //    if (opt.Contains("fast")) {
779    //       fastClone = kTRUE;
780     // The chain tree must have a list of branches  //       printf(" 21 \n");
781     // because we may try to change their basket  //    }
782     // size later.  //    // The chain tree must have a list of branches
783     TObjArray* lbranches = mychain->GetListOfBranches();  //    // because we may try to change their basket
784     if (!lbranches) {  //    // size later.
785        // FIXME: We need an error message here.  //    TObjArray* lbranches = mychain->GetListOfBranches();
786        return 0;  //    if (!lbranches) {
787     }  //       // FIXME: We need an error message here.
788    //      printf(" 22 \n");
789     // The chain must have a current tree because  //       return 0;
790     // that is the one we will clone.  //    }
791     //   if (!fTree) {  
792        // -- LoadTree() has not yet been called, no current tree.  //    //   file->cd();
793        // FIXME: We need an error message here.  //    // The chain must have a current tree because
794     //      return 0;  //    // that is the one we will clone.
795     //   }  //    //   if (!fTree) {
796    //       // -- LoadTree() has not yet been called, no current tree.
797     // Copy the chain's current tree without  //       // FIXME: We need an error message here.
798     // copying any entries, we will do that later.  //    //      return 0;
799     TTree* newTree = mychain->CloneTree(0);  //    //   }
800     if (!newTree) {  
801        // FIXME: We need an error message here.  //    // Copy the chain's current tree without
802        return 0;  //    // copying any entries, we will do that later.
803     }  //    printf(" 23 \n");
804    //    TTree* newTree = mychain->CloneTree(0);
805     // Strip out the (potential) directory name.  //    if (!newTree) {
806     // FIXME: The merged chain may or may not have the  //       // FIXME: We need an error message here.
807     //        same name as the original chain.  This is  //      printf(" 24 \n");
808     //        bad because the chain name determines the  //       return 0;
809     //        names of the trees in the chain by default.  //    }
810     newTree->SetName(gSystem->BaseName(mychain->GetName()));  //    printf(" 25 \n");
811    
812     // FIXME: Why do we do this?  
813     //   newTree->SetAutoSave(-1);  //    // Strip out the (potential) directory name.
814     newTree->SetAutoSave(900000000000000LL);  //    // FIXME: The merged chain may or may not have the
815    //    //        same name as the original chain.  This is
816     // Circularity is incompatible with merging, it may  //    //        bad because the chain name determines the
817     // force us to throw away entries, which is not what  //    //        names of the trees in the chain by default.
818     // we are supposed to do.  //    printf(" 26 \n");
819     newTree->SetCircular(0);  //    newTree->SetName(gSystem->BaseName(mychain->GetName()));
820    //    printf(" 27 \n");
821     // Reset the compression level of the branches.  
822     if (opt.Contains("c")) {  //    // FIXME: Why do we do this?
823        TBranch* branch = 0;  //    //   newTree->SetAutoSave(-1);
824        TIter nextb(newTree->GetListOfBranches());  //    newTree->SetAutoSave(900000000000000LL);
825        while ((branch = (TBranch*) nextb())) {  
826           branch->SetCompressionLevel(file->GetCompressionLevel());  //    printf(" 28 \n");
827        }  //    // Circularity is incompatible with merging, it may
828     }  //    // force us to throw away entries, which is not what
829    //    // we are supposed to do.
830     // Reset the basket size of the branches.  //    newTree->SetCircular(0);
831     if (basketsize > 1000) {  
832        TBranch* branch = 0;  //    // Reset the compression level of the branches.
833        TIter nextb(newTree->GetListOfBranches());  //    if (opt.Contains("c")) {
834        while ((branch = (TBranch*) nextb())) {  //       TBranch* branch = 0;
835           branch->SetBasketSize(basketsize);  //       TIter nextb(newTree->GetListOfBranches());
836        }  //       while ((branch = (TBranch*) nextb())) {
837     }  //      printf(" 29 \n");
838    //          branch->SetCompressionLevel(file->GetCompressionLevel());
839     Long64_t nentries = mychain->GetEntriesFast();  //       }
840    //    }
841     // Copy the entries.  
842     if (fastClone) {  //    printf(" 30 \n");
843        // For each tree in the chain.  //    // Reset the basket size of the branches.
844        for (Long64_t i = 0; i < nentries; i += mychain->GetTree()->GetEntries()) {  //    if (basketsize > 1000) {
845           if (mychain->LoadTree(i) < 0) {  //       TBranch* branch = 0;
846              break;  //       TIter nextb(newTree->GetListOfBranches());
847           }  //       while ((branch = (TBranch*) nextb())) {
848           TTreeCloner cloner(mychain->GetTree(), newTree, option);  //      printf(" 31 \n");
849           if (cloner.IsValid()) {  //          branch->SetBasketSize(basketsize);
850              newTree->SetEntries(newTree->GetEntries() + mychain->GetTree()->GetEntries());  //       }
851              cloner.Exec();  //    }
852           } else {  
853              if (mychain->GetFile()) {  //    printf(" 32 \n");
854                printf("Merge Skipped file %s\n", mychain->GetFile()->GetName());  //    Long64_t nentries = mychain->GetEntriesFast();
855                 //            } else {  //    //Long64_t nentries = mychain->GetEntries();
856                 //               Warning("Merge", "Skipped file number %d\n", fTreeNumber);  
857              }  //    // Copy the entries.
858           }  //    if (fastClone) {
859        }  //       // For each tree in the chain.
860     } else {  //       for (Long64_t i = 0; i < nentries; i += mychain->GetTree()->GetEntries()) {
861        for (Long64_t i = 0; i < nentries; i++) {  //          if (mychain->LoadTree(i) < 0) {
862           if (mychain->GetEntry(i) <= 0) {  //             break;
863              break;  //          }
864           }  //          TTreeCloner cloner(mychain->GetTree(), newTree, option);
865           newTree->Fill();  //          if (cloner.IsValid()) {
866        }  //             newTree->SetEntries(newTree->GetEntries() + mychain->GetTree()->GetEntries());
867     }  //             cloner.Exec();
868    //          } else {
869     // Write the new tree header.  //             if (mychain->GetFile()) {
870     newTree->Write();  //            printf("Merge Skipped file %s\n", mychain->GetFile()->GetName());
871    //             //            } else {
872     // Get our return value.  //             //               Warning("Merge", "Skipped file number %d\n", fTreeNumber);
873     Int_t nfiles = newTree->GetFileNumber() + 1;  //             }
874    //          }
875     // Close and delete the current file of the new tree.  //       }
876     if (!opt.Contains("keep")) {  //    } else {
877        // FIXME: What happens to fDirectory in newTree here?  //      printf(" 33 %llu \n",nentries);
878        delete newTree->GetCurrentFile();  //      mychain->Print();
879     }  //      for (Long64_t i = 0; i < nentries; i++) {
880     return nfiles;  //        printf(" i %llu \n",i);
881  }  //        // for (Long64_t i = 0; i < 1; i++) {
882    //          if (mychain->GetEntry(i) <= 0) {
883    //             break;
884    //          }
885    //          newTree->Fill();
886    //       }
887    //       printf(" 34 \n");
888    //    }
889    
890    //    // Write the new tree header.
891    //    printf(" 35 \n");
892    //    newTree->Write();
893    
894    //    printf(" 36 \n");
895    //    // Get our return value.
896    //    Int_t nfiles = newTree->GetFileNumber() + 1;
897    
898    //    // Close and delete the current file of the new tree.
899    //    if (!opt.Contains("keep")) {
900    //       // FIXME: What happens to fDirectory in newTree here?
901    //       delete newTree->GetCurrentFile();
902    //    }
903    //    return nfiles;
904    // }
Legend:
Removed from v.1.9  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.23