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; |
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; |
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 |
|
|
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; |
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) |
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(); |
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 |
// |
// |
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; |
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; |
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 |
; |
; |
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"); |
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 |
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 |
}; |
}; |
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 |
|
|
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 |
|
|
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; |
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; |
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; |
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; |
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; |
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; |
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; |
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; |
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; |
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; |
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; |
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 |
|
|
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 |
|
|
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 || true ){ |
5233 |
|
cout << "|| "<<endl; |
5234 |
|
proc_obj->runID = 0; |
5235 |
|
cout << "|||d "<<endl; |
5236 |
|
TTimeStamp *dt = new TTimeStamp(); |
5237 |
|
proc_obj->date = dt->AsString(); |
5238 |
|
delete dt; |
5239 |
|
cout << "|||f "<<endl; |
5240 |
|
proc_obj->commandLine = Form("PamelaLevel2 was called"); |
5241 |
|
cout << "|||g "<<endl; |
5242 |
|
proc_obj->outputFilename = ""; |
5243 |
|
cout << "|||h "<<endl; |
5244 |
|
proc_obj->localDir = gSystem->WorkingDirectory(); |
5245 |
|
cout << "|||j "<<endl; |
5246 |
|
proc_obj->uname = gSystem->GetFromPipe("uname -a"); |
5247 |
|
cout << "|||s "<<endl; |
5248 |
|
proc_obj->DB = Form("mysql://%s/%s",dbc->GetHost(),dbc->GetDB()); |
5249 |
|
cout << "||| "<<endl; |
5250 |
|
proc_tree_clone->Fill(); |
5251 |
|
} |
5252 |
|
cout << "----------------------------------------------------" << endl; |
5253 |
|
} |
5254 |
// ------------------------------------ |
// ------------------------------------ |
5255 |
// add branch with dead and live times |
// add branch with dead and live times |
5256 |
// ------------------------------------ |
// ------------------------------------ |
5441 |
if (!name.CompareTo(na)) |
if (!name.CompareTo(na)) |
5442 |
return sel_tree_clone; |
return sel_tree_clone; |
5443 |
} |
} |
5444 |
|
if (proc_tree_clone && PROC) { |
5445 |
|
TString na = proc_tree_clone->GetName(); |
5446 |
|
if (!name.CompareTo(na)) |
5447 |
|
return proc_tree_clone; |
5448 |
|
} |
5449 |
return NULL; |
return NULL; |
5450 |
|
|
5451 |
} |
} |
5465 |
for (Int_t i = 0; i < NCLONES; i++) { |
for (Int_t i = 0; i < NCLONES; i++) { |
5466 |
if (pam_tree_clone[i]) { |
if (pam_tree_clone[i]) { |
5467 |
cout << pam_tree_clone[i]->GetName() << endl; |
cout << pam_tree_clone[i]->GetName() << endl; |
5468 |
pam_tree_clone[i]->Write(); |
pam_tree_clone[i]->Write(pam_tree_clone[i]->GetName(),TObject::kOverwrite); |
5469 |
}; |
}; |
5470 |
} |
} |
5471 |
|
|
5472 |
|
if ( PROC && false ){//EMEMEMEM |
5473 |
|
proc_tree_clone->Write("ProcessingInfo",TObject::kOverwrite); |
5474 |
|
} |
5475 |
cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" << endl; |
cout << "+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+--+" << endl; |
5476 |
|
|
5477 |
} |
} |