/[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.96 by pam-ts, Wed Oct 15 10:02:51 2014 UTC revision 1.107 by mocchiut, Fri Oct 17 16:25:38 2014 UTC
# Line 374  PamTrack::PamTrack(ExtTrack* t, CaloTrkV Line 374  PamTrack::PamTrack(ExtTrack* t, CaloTrkV
374  PamTrack::PamTrack(const PamTrack& track) {  PamTrack::PamTrack(const PamTrack& track) {
375    
376    TrkTrack *t = track.trk_track;    TrkTrack *t = track.trk_track;
377      ExtTrack *e = track.trk_ext_track;
378    CaloTrkVar *c = track.calo_track;    CaloTrkVar *c = track.calo_track;
379    ToFTrkVar *o = track.tof_track;    ToFTrkVar *o = track.tof_track;
380    OrbitalInfoTrkVar *r = track.orb_track;    OrbitalInfoTrkVar *r = track.orb_track;
381    
382    trk_track = 0;    trk_ext_track = 0;
383    calo_track = 0;    trk_track     = 0;
384    tof_track = 0;    calo_track    = 0;
385    orb_track = 0;    tof_track     = 0;
386      orb_track     = 0;
387      if(e)
388          trk_ext_track = new ExtTrack(*e);
389    if (t)    if (t)
390      trk_track = new TrkTrack(*t);      trk_track = new TrkTrack(*t);
391    if (c)    if (c)
# Line 390  PamTrack::PamTrack(const PamTrack& track Line 394  PamTrack::PamTrack(const PamTrack& track
394      tof_track = new ToFTrkVar(*o);      tof_track = new ToFTrkVar(*o);
395    if (r)    if (r)
396      orb_track = new OrbitalInfoTrkVar(*r);      orb_track = new OrbitalInfoTrkVar(*r);
397    
398    candeleteobj = 1;    candeleteobj = 1;
399    pscore = 0;    pscore = 0;
400    iscore = 0;    iscore = 0;
401    
402  }  }
403  void PamTrack::Clear() {  void PamTrack::Clear(Option_t *option) {
404    
405    //    cout << "PamTrack::Clear() "<<candeleteobj<<endl;  //    cout << "PamTrack::Clear( "<<option<<" ) "<<candeleteobj<<endl;
406    if (candeleteobj) {    if (candeleteobj) {
407                
408      if (trk_ext_track)      if (trk_ext_track)
409        trk_ext_track->ExtTrack::Clear();        trk_ext_track->ExtTrack::Clear(option);
410      if (trk_track)      if (trk_track)
411        trk_track->TrkTrack::Clear();        trk_track->TrkTrack::Clear();
412      if (calo_track)      if (calo_track)
# Line 426  void PamTrack::Delete() { Line 431  void PamTrack::Delete() {
431    //    cout << "PamTrack::Delete() "<<candeleteobj<<endl;    //    cout << "PamTrack::Delete() "<<candeleteobj<<endl;
432    if (candeleteobj) {    if (candeleteobj) {
433      if (trk_ext_track) {      if (trk_ext_track) {
434        trk_ext_track->ExtTrack::Clear();        trk_ext_track->ExtTrack::Clear("C");
435        delete trk_ext_track;        delete trk_ext_track;
436      }      }
437      if (trk_track) {      if (trk_track) {
# Line 561  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 571  void PamLevel2::Initialize() { Line 577  void PamLevel2::Initialize() {
577    gltsync = 0; // Emiliano    gltsync = 0; // Emiliano
578    fUpdateRunInfo = true; // Emiliano    fUpdateRunInfo = true; // Emiliano
579    fUseDBinRunInfo = true; // Emiliano    fUseDBinRunInfo = true; // Emiliano
580      fDiscarded = false; //EM
581    isSync = false; // by default assume that the level2 file(s) is(are) not sinchronized between L0/DB and L2, that is we miss some packets in L2 due to nested/DV-skipped events    isSync = false; // by default assume that the level2 file(s) is(are) not sinchronized between L0/DB and L2, that is we miss some packets in L2 due to nested/DV-skipped events
582    il0entry = 0LL;    il0entry = 0LL;
583    //  hasL0EE = true;    //  hasL0EE = true;
# Line 586  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 633  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 675  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 875  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 955  void PamLevel2::Reset() { Line 971  void PamLevel2::Reset() {
971    if (sel_tree)    if (sel_tree)
972      sel_tree->Delete();;      sel_tree->Delete();;
973    sel_tree = NULL;    sel_tree = NULL;
974    
975      if (proc_tree)
976        proc_tree->Delete();
977      proc_tree = NULL;
978    //    //
979    // Close file    // Close file
980    //    //
# Line 978  void PamLevel2::Reset() { Line 998  void PamLevel2::Reset() {
998    ac_obj = 0;    ac_obj = 0;
999    orb2_obj = 0;    orb2_obj = 0;
1000    gp_obj = 0;    gp_obj = 0;
1001      proc_obj = 0;
1002    
1003    trk_ext_obj     = 0;    trk_ext_obj     = 0;
1004    trk_ext_nuc_obj = 0;    trk_ext_nuc_obj = 0;
# Line 1009  void PamLevel2::Reset() { Line 1030  void PamLevel2::Reset() {
1030    //    //
1031    run_obj = 0;//new GL_RUN();    run_obj = 0;//new GL_RUN();
1032    soft_obj = 0;// Emiliano    soft_obj = 0;// Emiliano
1033      proc_obj = 0;// Emiliano
1034    irun = -1;    irun = -1;
1035    irunt = -1;    irunt = -1;
1036    totrunentry = 0LL;    totrunentry = 0LL;
# Line 1183  void *PamLevel2::GetPointerTo(const char Line 1205  void *PamLevel2::GetPointerTo(const char
1205    if (!objname.CompareTo("SoftInfo"))    if (!objname.CompareTo("SoftInfo"))
1206      return &soft_obj; // Emiliano      return &soft_obj; // Emiliano
1207    
1208      if (!objname.CompareTo("ProcInfo")){
1209        if (!proc_obj)
1210          proc_obj = new ProcInfo();    
1211        return &proc_obj; // Emiliano
1212      }
1213    
1214    return NULL;    return NULL;
1215  }  }
1216  ;  ;
# Line 1448  void PamLevel2::SortTracks() { Line 1476  void PamLevel2::SortTracks() {
1476    // create TCloneArrays to store tracks and its images    // create TCloneArrays to store tracks and its images
1477    if (!tsorted)    if (!tsorted)
1478      tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());      tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
1479    tsorted->Delete();    tsorted->Clear("C+C");//Delete();
1480    TClonesArray &ttsorted = *tsorted;    TClonesArray &ttsorted = *tsorted;
1481    
1482    if (!timage)    if (!timage)
1483      timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());      timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
1484    timage->Delete();    timage->Clear("C+C");//Delete();
1485    TClonesArray &ttimage = *timage;    TClonesArray &ttimage = *timage;
1486    
1487    
# Line 1989  void PamLevel2::SortTracksNew() { Line 2017  void PamLevel2::SortTracksNew() {
2017    //-----------------------------------------------------------    //-----------------------------------------------------------
2018    // create/reset TCloneArrays to store tracks and their images    // create/reset TCloneArrays to store tracks and their images
2019    //-----------------------------------------------------------    //-----------------------------------------------------------
2020    
2021    //  cout << " PamLevel2::SortTracksNew() --- Clear TClonesArray objects"<<endl;
2022        
2023    // main tracks from standard alg    // main tracks from standard alg
2024    if (!tsorted)  //   if (!tsorted)
2025        tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());  //       tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
2026    tsorted->Delete();  //   tsorted->Clear("C+C");//Delete();
2027    // track images from standard alg  //   // track images from standard alg
2028    if (!timage)  //   if (!timage)
2029        timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());  //       timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
2030    timage->Delete();  //   timage->Clear("C+C");//Delete();
2031    // tracks from extended algorythm  //   // tracks from extended algorythm
2032    if(EXT && !text)  //   if(EXT && !text)
2033        text = new TClonesArray("PamTrack",trk_ext_obj->GetEntries());  //       text = new TClonesArray("PamTrack",trk_ext_obj->GetEntries());
2034    if(text)text->Delete();  //   if(text)text->Clear("C+C");//Delete();
2035    
2036      if(tsorted)delete tsorted;
2037      if(timage) delete timage;
2038      if(text)   delete text;
2039      tsorted = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
2040      timage = new TClonesArray("PamTrack", trk2_obj->GetNTracks());
2041      text = new TClonesArray("PamTrack",trk_ext_obj->GetEntries());
2042    
2043    //-----------------------------------------------------------    //-----------------------------------------------------------
2044    // create/reset TCloneArrays to store tracks and their images    // create/reset TCloneArrays to store tracks and their images
2045    //-----------------------------------------------------------    //-----------------------------------------------------------
2046    if(NUC){    if(NUC){
2047    
2048          
2049          if(tsorted_nuc)delete tsorted_nuc;
2050          if(timage_nuc) delete timage_nuc;
2051          if(text_nuc)   delete text_nuc;
2052          tsorted_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
2053          timage_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
2054          text_nuc = new TClonesArray("PamTrack",trk_ext_nuc_obj->GetEntries());
2055    
2056    // main tracks from standard alg    // main tracks from standard alg
2057        if (!tsorted_nuc)  //       if (!tsorted_nuc)
2058            tsorted_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());  //        tsorted_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
2059        tsorted_nuc->Delete();  //       tsorted_nuc->Clear("C+C");//Delete();
2060        // track images from standard alg  //       // track images from standard alg
2061        if (!timage_nuc)  //       if (!timage_nuc)
2062            timage_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());  //        timage_nuc = new TClonesArray("PamTrack", trk2_nuc_obj->GetNTracks());
2063        timage_nuc->Delete();  //       timage_nuc->Clear("C+C");//Delete();
2064        // tracks from extended algorythm  //       // tracks from extended algorythm
2065        if(EXT && !text_nuc)  //       if(EXT && !text_nuc)
2066            text_nuc = new TClonesArray("PamTrack",trk_ext_nuc_obj->GetEntries());  //        text_nuc = new TClonesArray("PamTrack",trk_ext_nuc_obj->GetEntries());
2067        if(text_nuc)text_nuc->Delete();  //       if(text_nuc)text_nuc->Clear("C+C");//Delete();
2068    
2069    }    }
2070    //--------------------------------------------------    //--------------------------------------------------
2071    // retrieve sorting method    // retrieve sorting method
# Line 3279  TChain *PamLevel2::GetPamTree(TList *fl, Line 3326  TChain *PamLevel2::GetPamTree(TList *fl,
3326    TChain *G = 0;    TChain *G = 0;
3327    
3328    TChain *L = 0;    TChain *L = 0;
3329      TChain *P = 0;
3330    
3331    if (TRK2 || TRK1 || TRKh)    if (TRK2 || TRK1 || TRKh)
3332      T = new TChain("Tracker");      T = new TChain("Tracker");
# Line 3298  TChain *PamLevel2::GetPamTree(TList *fl, Line 3346  TChain *PamLevel2::GetPamTree(TList *fl,
3346      B = new TChain("OrbitalInfo");      B = new TChain("OrbitalInfo");
3347    if (GP)    if (GP)
3348      G = new TChain("h20");      G = new TChain("h20");
3349      if (PROC)
3350        P = new TChain("ProcessingInfo");
3351    L = new TChain("SelectionList");    L = new TChain("SelectionList");
3352    
3353    // loop over files and create chains    // loop over files and create chains
# Line 3325  TChain *PamLevel2::GetPamTree(TList *fl, Line 3375  TChain *PamLevel2::GetPamTree(TList *fl,
3375          B->Add(name);          B->Add(name);
3376        if (GP)        if (GP)
3377          G->Add(name);          G->Add(name);
3378          if (P)
3379            P->Add(name);
3380        if (SELLI == 1)        if (SELLI == 1)
3381          L->Add(name);          L->Add(name);
3382       };       };
# Line 3450  TChain *PamLevel2::GetPamTree(TList *fl, Line 3502  TChain *PamLevel2::GetPamTree(TList *fl,
3502        L->Delete();        L->Delete();
3503    };    };
3504    
3505      //ProcessingInfo EM
3506      if ( P && P->GetEntries() ){
3507        cout << "----------------------------------------------------" << endl;
3508        cout << ">>> Found ProcessingInfo <<<" << endl;
3509        //  L->SetBranchAddress("RunEntry",&irun);
3510        P->SetBranchAddress("ProcInfo", &proc_obj);//NEWNEW
3511        proc_tree = P;
3512      }
3513    //  --------------------------------------------    //  --------------------------------------------
3514    //  return the pamela chain with all the friends    //  return the pamela chain with all the friends
3515    //  --------------------------------------------    //  --------------------------------------------
3516    
3517    pam_tree = Trout;    pam_tree = Trout;
   
3518    return Trout;    return Trout;
3519  }  }
3520    
# Line 4547  void PamLevel2::GetWhichTrees(TFile* f) Line 4606  void PamLevel2::GetWhichTrees(TFile* f)
4606    //    cout << "Checking file: "<<f->GetName()<<endl;    //    cout << "Checking file: "<<f->GetName()<<endl;
4607    if (!f || f->IsZombie()) {    if (!f || f->IsZombie()) {
4608      cout << "File: " << f->GetName() << " Non valid root file" << endl;      cout << "File: " << f->GetName() << " Non valid root file" << endl;
4609        fDiscarded = true;
4610      return;      return;
4611    }    }
4612    
# Line 4753  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4813  Bool_t PamLevel2::CheckLevel2File(TStrin
4813    TFile *f = new TFile(name.Data());    TFile *f = new TFile(name.Data());
4814    if (!f || f->IsZombie()) {    if (!f || f->IsZombie()) {
4815      cout << "File: " << f->GetName() << " discarded ---- Non valid root file" << endl;      cout << "File: " << f->GetName() << " discarded ---- Non valid root file" << endl;
4816        fDiscarded = true;
4817      return false;      return false;
4818    }    }
4819    //    cout << "Get list of keys: "<<f<<endl;    //    cout << "Get list of keys: "<<f<<endl;
# Line 4791  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4852  Bool_t PamLevel2::CheckLevel2File(TStrin
4852          if (nev && nevt != nev) {          if (nev && nevt != nev) {
4853            cout << "File: " << f->GetName() << " discarded ---- SelectionList tree has " << nevt            cout << "File: " << f->GetName() << " discarded ---- SelectionList tree has " << nevt
4854                << " events instead of " << nev << endl;                << " events instead of " << nev << endl;
4855              fDiscarded = true;
4856            return false;            return false;
4857          }          }
4858          nev = nevt;          nev = nevt;
# Line 4805  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4867  Bool_t PamLevel2::CheckLevel2File(TStrin
4867          if (nev && nevt != nev) {          if (nev && nevt != nev) {
4868            cout << "File: " << f->GetName() << " discarded ---- Trigger tree has " << nevt << " events instead of "            cout << "File: " << f->GetName() << " discarded ---- Trigger tree has " << nevt << " events instead of "
4869                << nev << endl;                << nev << endl;
4870              fDiscarded = true;
4871            return false;            return false;
4872          }          }
4873          nev = nevt;          nev = nevt;
# Line 4818  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4881  Bool_t PamLevel2::CheckLevel2File(TStrin
4881          if (nev && nevt != nev) {          if (nev && nevt != nev) {
4882            cout << "File: " << f->GetName() << " discarded ---- ToF tree has " << nevt << " events instead of " << nev            cout << "File: " << f->GetName() << " discarded ---- ToF tree has " << nevt << " events instead of " << nev
4883                << endl;                << endl;
4884              fDiscarded = true;
4885            return false;            return false;
4886          }          }
4887          nev = nevt;          nev = nevt;
# Line 4831  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4895  Bool_t PamLevel2::CheckLevel2File(TStrin
4895          if (nev && nevt != nev) {          if (nev && nevt != nev) {
4896            cout << "File: " << f->GetName() << " discarded ---- S4 tree has " << nevt << " events instead of " << nev            cout << "File: " << f->GetName() << " discarded ---- S4 tree has " << nevt << " events instead of " << nev
4897                << endl;                << endl;
4898              fDiscarded = true;
4899            return false;            return false;
4900          }          }
4901          nev = nevt;          nev = nevt;
# Line 4845  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4910  Bool_t PamLevel2::CheckLevel2File(TStrin
4910          if (nev && nevt != nev) {          if (nev && nevt != nev) {
4911            cout << "File: " << f->GetName() << " discarded ---- NeutronD tree has " << nevt << " events instead of "            cout << "File: " << f->GetName() << " discarded ---- NeutronD tree has " << nevt << " events instead of "
4912                << nev << endl;                << nev << endl;
4913              fDiscarded = true;
4914            return false;            return false;
4915          }          }
4916          nev = nevt;          nev = nevt;
# Line 4858  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4924  Bool_t PamLevel2::CheckLevel2File(TStrin
4924          if (nev && nevt != nev) {          if (nev && nevt != nev) {
4925            cout << "File: " << f->GetName() << " discarded ---- Anticounter tree has " << nevt << " events instead of "            cout << "File: " << f->GetName() << " discarded ---- Anticounter tree has " << nevt << " events instead of "
4926                << nev << endl;                << nev << endl;
4927              fDiscarded = true;
4928            return false;            return false;
4929          }          }
4930          nev = nevt;          nev = nevt;
# Line 4871  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4938  Bool_t PamLevel2::CheckLevel2File(TStrin
4938          if (nev && nevt != nev) {          if (nev && nevt != nev) {
4939            cout << "File: " << f->GetName() << " discarded ---- OrbitalInfo tree has " << nevt << " events instead of "            cout << "File: " << f->GetName() << " discarded ---- OrbitalInfo tree has " << nevt << " events instead of "
4940                << nev << endl;                << nev << endl;
4941              fDiscarded = true;
4942            return false;            return false;
4943          }          }
4944          nev = nevt;          nev = nevt;
# Line 4884  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4952  Bool_t PamLevel2::CheckLevel2File(TStrin
4952          if (nev && nevt != nev) {          if (nev && nevt != nev) {
4953            cout << "File: " << f->GetName() << " discarded ---- Tracker tree has " << nevt << " events instead of "            cout << "File: " << f->GetName() << " discarded ---- Tracker tree has " << nevt << " events instead of "
4954                << nev << endl;                << nev << endl;
4955              fDiscarded = true;
4956            return false;            return false;
4957          }          }
4958          nev = nevt;          nev = nevt;
# Line 4907  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4976  Bool_t PamLevel2::CheckLevel2File(TStrin
4976          if (nev && nevt != nev) {          if (nev && nevt != nev) {
4977            cout << "File: " << f->GetName() << " discarded ---- Calorimeter tree has " << nevt << " events instead of "            cout << "File: " << f->GetName() << " discarded ---- Calorimeter tree has " << nevt << " events instead of "
4978                << nev << endl;                << nev << endl;
4979              fDiscarded = true;
4980            return false;            return false;
4981          }          }
4982          nev = nevt;          nev = nevt;
# Line 4929  Bool_t PamLevel2::CheckLevel2File(TStrin Line 4999  Bool_t PamLevel2::CheckLevel2File(TStrin
4999          if (nev && nevt != nev) {          if (nev && nevt != nev) {
5000            cout << "File: " << f->GetName() << " discarded ---- h20 tree has " << nevt << " events instead of " << nev            cout << "File: " << f->GetName() << " discarded ---- h20 tree has " << nevt << " events instead of " << nev
5001                << endl;                << endl;
5002              fDiscarded = true;
5003            return false;            return false;
5004          }          }
5005          nev = nevt;          nev = nevt;
# Line 4941  Bool_t PamLevel2::CheckLevel2File(TStrin Line 5012  Bool_t PamLevel2::CheckLevel2File(TStrin
5012      SELLI = (Int_t) SELLI__ok;      SELLI = (Int_t) SELLI__ok;
5013    if (SELLI == 0 && SELLI__ok) {    if (SELLI == 0 && SELLI__ok) {
5014      cout << "File: " << f->GetName() << " discarded ---- found SelectionList (it is not a full-event file)" << endl;      cout << "File: " << f->GetName() << " discarded ---- found SelectionList (it is not a full-event file)" << endl;
5015        fDiscarded = true;
5016      return false;      return false;
5017    }    }
5018    if (SELLI == 1 && !SELLI__ok) {    if (SELLI == 1 && !SELLI__ok) {
5019      cout << "File: " << f->GetName() << " discarded ---- SelectionList missing" << endl;      cout << "File: " << f->GetName() << " discarded ---- SelectionList missing" << endl;
5020        fDiscarded = true;
5021      return false;      return false;
5022    }    }
5023    
# Line 4983  Bool_t PamLevel2::CheckLevel2File(TStrin Line 5056  Bool_t PamLevel2::CheckLevel2File(TStrin
5056    
5057    if (CAL1 && !CAL1__ok) {    if (CAL1 && !CAL1__ok) {
5058      cout << "File: " << f->GetName() << " discarded ---- Missing CaloLevel1 branch" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing CaloLevel1 branch" << endl;
5059        fDiscarded = true;
5060      return false;      return false;
5061    };    };
5062    if (CAL2 && !CAL2__ok) {    if (CAL2 && !CAL2__ok) {
5063      cout << "File: " << f->GetName() << " discarded ---- Missing CaloLevel2 branch" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing CaloLevel2 branch" << endl;
5064        fDiscarded = true;
5065      return false;      return false;
5066    };    };
5067    if (TRK2 && !TRK2__ok) {    if (TRK2 && !TRK2__ok) {
5068      cout << "File: " << f->GetName() << " discarded ---- Missing TrkLevel2 branch" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing TrkLevel2 branch" << endl;
5069        fDiscarded = true;
5070      return false;      return false;
5071    };    };
5072    if (TRK1 && !TRK1__ok) {    if (TRK1 && !TRK1__ok) {
5073      cout << "File: " << f->GetName() << " discarded ---- Missing TrkLevel1 branch" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing TrkLevel1 branch" << endl;
5074        fDiscarded = true;
5075      return false;      return false;
5076    };    };
5077    if (TRKh && !TRKh__ok) {    if (TRKh && !TRKh__ok) {
5078      cout << "File: " << f->GetName() << " discarded ---- Missing TrkHough branch" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing TrkHough branch" << endl;
5079        fDiscarded = true;
5080      return false;      return false;
5081    };    };
5082    if (ORB && !ORB__ok) {    if (ORB && !ORB__ok) {
5083      cout << "File: " << f->GetName() << " discarded ---- Missing ORB tree" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing ORB tree" << endl;
5084        fDiscarded = true;
5085      return false;      return false;
5086    };    };
5087    if (AC && !AC__ok) {    if (AC && !AC__ok) {
5088      cout << "File: " << f->GetName() << " discarded ---- Missing AC tree" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing AC tree" << endl;
5089        fDiscarded = true;
5090      return false;      return false;
5091    };    };
5092    if (S4 && !S4__ok) {    if (S4 && !S4__ok) {
5093      cout << "File: " << f->GetName() << " discarded ---- Missing S4 tree" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing S4 tree" << endl;
5094        fDiscarded = true;
5095      return false;      return false;
5096    };    };
5097    if (TOF && !TOF__ok) {    if (TOF && !TOF__ok) {
5098      cout << "File: " << f->GetName() << " discarded ---- Missing ToF tree" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing ToF tree" << endl;
5099        fDiscarded = true;
5100      return false;      return false;
5101    };    };
5102    
5103    if (ND && !ND__ok) {    if (ND && !ND__ok) {
5104      cout << "File: " << f->GetName() << " discarded ---- Missing ND tree" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing ND tree" << endl;
5105        fDiscarded = true;
5106      return false;      return false;
5107    };    };
5108    if (TRG && !TRG__ok) {    if (TRG && !TRG__ok) {
5109      cout << "File: " << f->GetName() << " discarded ---- Missing Trigger tree" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing Trigger tree" << endl;
5110        fDiscarded = true;
5111      return false;      return false;
5112    };    };
5113    if (GP && !GP__ok) {    if (GP && !GP__ok) {
5114      cout << "File: " << f->GetName() << " discarded ---- Missing h20 tree" << endl;      cout << "File: " << f->GetName() << " discarded ---- Missing h20 tree" << endl;
5115        fDiscarded = true;
5116      return false;      return false;
5117    };    };
5118    
# Line 5127  void PamLevel2::CreateCloneTrees(TFile * Line 5212  void PamLevel2::CreateCloneTrees(TFile *
5212      //  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;
5213      run_tree_clone->Fill();      run_tree_clone->Fill();
5214    }    }
5215    cout << "----------------------------------------------------" << endl;    //
5216      // replicate processinginfo tree
5217      //
5218      if ( PROC ){ // EMEMEM
5219        proc_tree_clone = new TTree("ProcessingInfo","Log of data processing");
5220        proc_tree_clone->Branch("ProcInfo", "ProcInfo", GetPointerTo("ProcInfo"));
5221        cout << "ProcessingInfo: branch ProcessingInfo" << endl;  
5222        // ------------------
5223        // replicate processinginfo tree
5224        // ------------------
5225        //    cout << "----------------------------------------------------"<<endl;
5226        //    cout << "irun\t | RUN\t NEVENTS\t absolute time"<<endl;
5227        for (Int_t i = 0; i < proc_tree->GetEntries(); i++) {
5228          proc_tree->GetEntry(i);
5229          //      cout << i<< "\t | "<<endl;
5230          proc_tree_clone->Fill();
5231        }
5232        if ( SELLI != 2 ){      
5233          proc_obj->runID = 0;
5234          TTimeStamp *dt = new TTimeStamp();
5235          proc_obj->date = dt->AsString();
5236          delete dt;
5237          proc_obj->commandLine = Form("PamelaLevel2 was called");
5238          proc_obj->outputFilename = "";
5239          proc_obj->localDir = gSystem->WorkingDirectory();
5240          proc_obj->uname = gSystem->GetFromPipe("uname -a");
5241          //      proc_obj->DB = Form("mysql://%s/%s",dbc->GetHost(),dbc->GetDB());
5242          proc_obj->DB = "";
5243          proc_tree_clone->Fill();
5244        }
5245        cout << "----------------------------------------------------" << endl;
5246      }
5247    // ------------------------------------    // ------------------------------------
5248    // add branch with dead and live times    // add branch with dead and live times
5249    // ------------------------------------    // ------------------------------------
# Line 5319  TTree* PamLevel2::GetCloneTree(TString n Line 5434  TTree* PamLevel2::GetCloneTree(TString n
5434      if (!name.CompareTo(na))      if (!name.CompareTo(na))
5435        return sel_tree_clone;        return sel_tree_clone;
5436    }    }
5437      if (proc_tree_clone && PROC) {
5438        TString na = proc_tree_clone->GetName();
5439        if (!name.CompareTo(na))
5440          return proc_tree_clone;
5441      }
5442    return NULL;    return NULL;
5443    
5444  }  }
# Line 5338  void PamLevel2::WriteCloneTrees() { Line 5458  void PamLevel2::WriteCloneTrees() {
5458    for (Int_t i = 0; i < NCLONES; i++) {    for (Int_t i = 0; i < NCLONES; i++) {
5459      if (pam_tree_clone[i]) {      if (pam_tree_clone[i]) {
5460        cout << pam_tree_clone[i]->GetName() << endl;        cout << pam_tree_clone[i]->GetName() << endl;
5461        pam_tree_clone[i]->Write();        pam_tree_clone[i]->Write(pam_tree_clone[i]->GetName(),TObject::kOverwrite);
5462      };      };
5463    }    }
5464      
5465      if ( PROC ){//EMEMEMEM
5466        proc_tree_clone->Write("ProcessingInfo",TObject::kOverwrite);
5467      }
5468    cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" << endl;    cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" << endl;
5469    
5470  }  }
Legend:
Removed from v.1.96  
changed lines
  Added in v.1.107
  ViewVC Help
Powered by ViewVC 1.1.23