/[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.99 by mocchiut, Fri Oct 17 15:08:35 2014 UTC revision 1.100 by mocchiut, Fri Oct 17 15:58:09 2014 UTC
# Line 566  void PamLevel2::Initialize() { Line 566  void PamLevel2::Initialize() {
566    
567    run_obj = 0;//new GL_RUN();    run_obj = 0;//new GL_RUN();
568    soft_obj = 0;// Emiliano    soft_obj = 0;// Emiliano
569      proc_obj = 0;// Emiliano
570    irun = -1LL;    irun = -1LL;
571    irunt = -1LL;    irunt = -1LL;
572    totrunentry = 0LL;    totrunentry = 0LL;
# Line 592  void PamLevel2::Initialize() { Line 593  void PamLevel2::Initialize() {
593    
594    run_tree = NULL;    run_tree = NULL;
595    run_tree_clone = NULL;    run_tree_clone = NULL;
596    
597      proc_tree = NULL;
598      proc_tree_clone = NULL;
599    
600    sel_tree = NULL;    sel_tree = NULL;
601    sel_tree_clone = NULL;    sel_tree_clone = NULL;
602    
# Line 639  void PamLevel2::Initialize() { Line 644  void PamLevel2::Initialize() {
644    ND = true;    ND = true;
645    AC = true;    AC = true;
646    ORB = true;    ORB = true;
647      PROC = true;
648    GP = false;    GP = false;
649    
650    EXT = false;    EXT = false;
# Line 681  void PamLevel2::Delete() { Line 687  void PamLevel2::Delete() {
687      delete run_obj;      delete run_obj;
688    if (soft_obj)    if (soft_obj)
689      delete soft_obj; //Emiliano      delete soft_obj; //Emiliano
690      if (proc_obj)
691        delete proc_obj; //Emiliano
692    
693    //    cout << "void PamLevel2::Clear()"<<endl;    //    cout << "void PamLevel2::Clear()"<<endl;
694    if (h0_obj)    if (h0_obj)
# Line 881  void PamLevel2::Clear() { Line 889  void PamLevel2::Clear() {
889      orb2_obj->Clear();      orb2_obj->Clear();
890    if (gp_obj)    if (gp_obj)
891      gp_obj->Clear();      gp_obj->Clear();
892      if (proc_obj)
893        proc_obj->Clear();
894    
895    //    if(sorted_tracks)sorted_tracks->Clear();    //    if(sorted_tracks)sorted_tracks->Clear();
896    //    sorted_tracks.Clear();    //    sorted_tracks.Clear();
# Line 984  void PamLevel2::Reset() { Line 994  void PamLevel2::Reset() {
994    ac_obj = 0;    ac_obj = 0;
995    orb2_obj = 0;    orb2_obj = 0;
996    gp_obj = 0;    gp_obj = 0;
997      proc_obj = 0;
998    
999    trk_ext_obj     = 0;    trk_ext_obj     = 0;
1000    trk_ext_nuc_obj = 0;    trk_ext_nuc_obj = 0;
# Line 1015  void PamLevel2::Reset() { Line 1026  void PamLevel2::Reset() {
1026    //    //
1027    run_obj = 0;//new GL_RUN();    run_obj = 0;//new GL_RUN();
1028    soft_obj = 0;// Emiliano    soft_obj = 0;// Emiliano
1029      proc_obj = 0;// Emiliano
1030    irun = -1;    irun = -1;
1031    irunt = -1;    irunt = -1;
1032    totrunentry = 0LL;    totrunentry = 0LL;
# Line 1189  void *PamLevel2::GetPointerTo(const char Line 1201  void *PamLevel2::GetPointerTo(const char
1201    if (!objname.CompareTo("SoftInfo"))    if (!objname.CompareTo("SoftInfo"))
1202      return &soft_obj; // Emiliano      return &soft_obj; // Emiliano
1203    
1204      if (!objname.CompareTo("ProcessingInfo")){
1205        if (!proc_obj)
1206          proc_obj = new ProcInfo();    
1207        return &proc_obj; // Emiliano
1208      }
1209    
1210    return NULL;    return NULL;
1211  }  }
1212  ;  ;
# Line 3304  TChain *PamLevel2::GetPamTree(TList *fl, Line 3322  TChain *PamLevel2::GetPamTree(TList *fl,
3322    TChain *G = 0;    TChain *G = 0;
3323    
3324    TChain *L = 0;    TChain *L = 0;
3325      TChain *P = 0;
3326    
3327    if (TRK2 || TRK1 || TRKh)    if (TRK2 || TRK1 || TRKh)
3328      T = new TChain("Tracker");      T = new TChain("Tracker");
# Line 3323  TChain *PamLevel2::GetPamTree(TList *fl, Line 3342  TChain *PamLevel2::GetPamTree(TList *fl,
3342      B = new TChain("OrbitalInfo");      B = new TChain("OrbitalInfo");
3343    if (GP)    if (GP)
3344      G = new TChain("h20");      G = new TChain("h20");
3345      if (PROC)
3346        P = new TChain("ProcessingInfo");
3347    L = new TChain("SelectionList");    L = new TChain("SelectionList");
3348    
3349    // loop over files and create chains    // loop over files and create chains
# Line 3350  TChain *PamLevel2::GetPamTree(TList *fl, Line 3371  TChain *PamLevel2::GetPamTree(TList *fl,
3371          B->Add(name);          B->Add(name);
3372        if (GP)        if (GP)
3373          G->Add(name);          G->Add(name);
3374          if (P)
3375            P->Add(name);
3376        if (SELLI == 1)        if (SELLI == 1)
3377          L->Add(name);          L->Add(name);
3378       };       };
# Line 3475  TChain *PamLevel2::GetPamTree(TList *fl, Line 3498  TChain *PamLevel2::GetPamTree(TList *fl,
3498        L->Delete();        L->Delete();
3499    };    };
3500    
3501      //ProcessingInfo EM
3502      if ( P && P->GetEntries() ){
3503        cout << "----------------------------------------------------" << endl;
3504        cout << ">>> Found ProcessingInfo <<<" << endl;
3505        //  L->SetBranchAddress("RunEntry",&irun);
3506        P->SetBranchAddress("ProcInfo", &proc_obj);//NEWNEW
3507      }
3508    //  --------------------------------------------    //  --------------------------------------------
3509    //  return the pamela chain with all the friends    //  return the pamela chain with all the friends
3510    //  --------------------------------------------    //  --------------------------------------------
# Line 5178  void PamLevel2::CreateCloneTrees(TFile * Line 5208  void PamLevel2::CreateCloneTrees(TFile *
5208      //  cout << i<< "\t | "<<GetRunInfo()->ID<<"\t "<<GetRunInfo()->NEVENTS<< "\t "<<GetRunInfo()->RUNHEADER_TIME<<" <---> "<<GetRunInfo()->RUNTRAILER_TIME<<endl;      //  cout << i<< "\t | "<<GetRunInfo()->ID<<"\t "<<GetRunInfo()->NEVENTS<< "\t "<<GetRunInfo()->RUNHEADER_TIME<<" <---> "<<GetRunInfo()->RUNTRAILER_TIME<<endl;
5209      run_tree_clone->Fill();      run_tree_clone->Fill();
5210    }    }
5211    cout << "----------------------------------------------------" << endl;    if ( PROC ){
5212        proc_tree_clone->Branch("ProcessingInfo", "ProcInfo", GetPointerTo("ProcessingInfo"));
5213        cout << "ProcessingInfo: branch ProcessingInfo" << endl;  
5214        // ------------------
5215        // replicate processinginfo tree
5216        // ------------------
5217        //    cout << "----------------------------------------------------"<<endl;
5218        //    cout << "irun\t | RUN\t NEVENTS\t absolute time"<<endl;
5219        for (Int_t i = 0; i < proc_tree->GetEntries(); i++) {
5220          proc_tree->GetEntry(i);
5221          //        cout << i<< "\t | "<<GetRunInfo()->ID<<"\t "<<GetRunInfo()->NEVENTS<< "\t "<<GetRunInfo()->RUNHEADER_TIME<<" <---> "<<GetRunInfo()->RUNTRAILER_TIME<<endl;
5222          proc_tree_clone->Fill();
5223        }
5224        if ( SELLI != 2 || true ){      
5225          proc_obj->runID = 0;
5226          TTimeStamp *dt = new TTimeStamp();
5227          proc_obj->date = dt->AsString();
5228          delete dt;
5229          proc_obj->commandLine = Form("PamelaLevel2 was called");
5230          proc_obj->outputFilename = "";
5231          proc_obj->localDir = gSystem->WorkingDirectory();
5232          proc_obj->uname = gSystem->GetFromPipe("uname -a");
5233          proc_obj->DB = Form("mysql://%s/%s",dbc->GetHost(),dbc->GetDB());
5234          proc_tree_clone->Fill();
5235        }
5236        cout << "----------------------------------------------------" << endl;
5237      }
5238    // ------------------------------------    // ------------------------------------
5239    // add branch with dead and live times    // add branch with dead and live times
5240    // ------------------------------------    // ------------------------------------
# Line 5370  TTree* PamLevel2::GetCloneTree(TString n Line 5425  TTree* PamLevel2::GetCloneTree(TString n
5425      if (!name.CompareTo(na))      if (!name.CompareTo(na))
5426        return sel_tree_clone;        return sel_tree_clone;
5427    }    }
5428      if (proc_tree_clone && PROC) {
5429        TString na = proc_tree_clone->GetName();
5430        if (!name.CompareTo(na))
5431          return proc_tree_clone;
5432      }
5433    return NULL;    return NULL;
5434    
5435  }  }
# Line 5392  void PamLevel2::WriteCloneTrees() { Line 5452  void PamLevel2::WriteCloneTrees() {
5452        pam_tree_clone[i]->Write(pam_tree_clone[i]->GetName(),TObject::kOverwrite);        pam_tree_clone[i]->Write(pam_tree_clone[i]->GetName(),TObject::kOverwrite);
5453      };      };
5454    }    }
5455      
5456      if ( PROC ){
5457        proc_tree_clone->Write("ProcessingInfo",TObject::kOverwrite);
5458      }
5459    cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" << endl;    cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" << endl;
5460    
5461  }  }
Legend:
Removed from v.1.99  
changed lines
  Added in v.1.100
  ViewVC Help
Powered by ViewVC 1.1.23