58 |
* @param debug debug flag. |
* @param debug debug flag. |
59 |
* @param tlefilename ascii file with TLE 3 line elements. |
* @param tlefilename ascii file with TLE 3 line elements. |
60 |
*/ |
*/ |
61 |
PamelaDBOperations::PamelaDBOperations(TString host, TString user, TString password, TString filerawname, TString filerootname, UInt_t boot, UInt_t tsync, UInt_t obt0, Bool_t debug, TString tlefilename, UInt_t dwinput, Bool_t staticp, Bool_t gpamela){ |
PamelaDBOperations::PamelaDBOperations(TString host, TString user, TString password, TString filerawname, TString filerootname, UInt_t boot, UInt_t tsync, UInt_t obt0, Bool_t debug, TString tlefilename, UInt_t dwinput, Bool_t staticp, Bool_t gpamela, Bool_t keepenv){ |
62 |
// |
// |
63 |
|
RELAXED = false; |
64 |
chewbacca = false; |
chewbacca = false; |
65 |
chminentry = 0; |
chminentry = 0; |
66 |
chID = 0; |
chID = 0; |
77 |
// |
// |
78 |
}; |
}; |
79 |
// |
// |
80 |
|
KEEPENV = false; |
81 |
|
if ( keepenv ) KEEPENV=true; |
82 |
STATIC=false; |
STATIC=false; |
83 |
if ( staticp ) STATIC=true; |
if ( staticp ) STATIC=true; |
84 |
// |
// |
88 |
// |
// |
89 |
glrun = new GL_RUN(); |
glrun = new GL_RUN(); |
90 |
// |
// |
91 |
if ( !chewbacca ){ |
// if ( !chewbacca ){ |
92 |
if ( !boot ) SetNOBOOT(false); |
// if ( !boot ) SetNOBOOT(false); |
93 |
SetTsync(tsync,gpamela); |
// SetTsync(tsync,gpamela); |
94 |
SetBOOTnumber(boot,gpamela); |
// SetBOOTnumber(boot,gpamela); |
95 |
SetObt0(obt0); |
// SetObt0(obt0); |
96 |
}; |
// }; |
97 |
// |
// |
98 |
SetTLEPath(tlefilename); |
SetTLEPath(tlefilename); |
99 |
// |
// |
102 |
if (INSERT_RAW) SetRawName(filerawname); |
if (INSERT_RAW) SetRawName(filerawname); |
103 |
// |
// |
104 |
INSERT_ROOT = !filerootname.IsNull(); |
INSERT_ROOT = !filerootname.IsNull(); |
105 |
|
if ( INSERT_ROOT ) this->SetRootName(filerootname); |
106 |
this->SetOrbitNo(dwinput); |
this->SetOrbitNo(dwinput); |
107 |
// |
// |
108 |
this->SetID_RAW(0); |
this->SetID_RAW(0); |
126 |
// |
// |
127 |
// SETTERS |
// SETTERS |
128 |
// |
// |
129 |
|
void PamelaDBOperations::NotChewbacca(UInt_t boot, UInt_t tsync, UInt_t obt0, Bool_t gpamela){ |
130 |
|
if ( !boot ) SetNOBOOT(false); |
131 |
|
SetTsync(tsync,gpamela); |
132 |
|
SetBOOTnumber(boot,gpamela); |
133 |
|
SetObt0(obt0); |
134 |
|
}; |
135 |
|
|
136 |
void PamelaDBOperations::OpenL0File(TString filerootname){ |
void PamelaDBOperations::OpenL0File(TString filerootname){ |
137 |
if( INSERT_ROOT ){ |
if( INSERT_ROOT ){ |
142 |
stringstream qu; |
stringstream qu; |
143 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
144 |
TSQLResult *result2 = 0; |
TSQLResult *result2 = 0; |
145 |
|
TSQLResult *result3 = 0; |
146 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
147 |
TSQLRow *row2 = 0; |
TSQLRow *row2 = 0; |
148 |
|
TSQLRow *row3 = 0; |
149 |
TString chpath; |
TString chpath; |
150 |
TString chfile; |
TString chfile; |
151 |
UInt_t ridn = 0; |
UInt_t ridn = 0; |
152 |
qu.str(""); |
qu.str(""); |
153 |
qu << "SELECT ROOT_ID_N,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL from ROOT_TABLE_MERGING where ID_N=" << chID << ";"; |
qu << "SELECT ROOT_ID_N,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL,INSERTED_BY,REAL_TIME_INIT from ROOT_TABLE_MERGING where ID_N=" << chID << ";"; |
154 |
if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str()); |
if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str()); |
155 |
result = conn->Query(qu.str().c_str()); |
result = conn->Query(qu.str().c_str()); |
156 |
if ( result ){ |
if ( result ){ |
161 |
chpktmax = (UInt_t)atoll(row->GetField(2)); |
chpktmax = (UInt_t)atoll(row->GetField(2)); |
162 |
chobtmin = (UInt_t)atoll(row->GetField(3)); |
chobtmin = (UInt_t)atoll(row->GetField(3)); |
163 |
chobtmax = (UInt_t)atoll(row->GetField(4)); |
chobtmax = (UInt_t)atoll(row->GetField(4)); |
164 |
|
chiby = (TString)(row->GetField(5)); |
165 |
|
chrtinit = (UInt_t)atoll(row->GetField(6)); |
166 |
} else { |
} else { |
167 |
throw -84; |
throw -84; |
168 |
}; |
}; |
171 |
}; |
}; |
172 |
delete result; |
delete result; |
173 |
qu.str(""); |
qu.str(""); |
174 |
qu << "SELECT FOLDER_NAME,FILE_NAME,OBT_TIME_SYNC,LAST_TIME_SYNC_INFO,TIME_OFFSET,BOOT_NUMBER,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL from ROOT_TABLE where ID_N=" << ridn << ";"; |
qu << "SELECT FOLDER_NAME,FILE_NAME,OBT_TIME_SYNC,LAST_TIME_SYNC_INFO,TIME_OFFSET,BOOT_NUMBER,PKT_NUMBER_INIT,PKT_NUMBER_FINAL,PKT_OBT_INIT,PKT_OBT_FINAL,REAL_TIME_INIT from ROOT_TABLE where ID_N=" << ridn << ";"; |
175 |
if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str()); |
if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str()); |
176 |
result2 = conn->Query(qu.str().c_str()); |
result2 = conn->Query(qu.str().c_str()); |
177 |
if ( result2 ){ |
if ( result2 ){ |
178 |
row2 = result2->Next(); |
row2 = result2->Next(); |
179 |
if ( row2 ){ |
if ( row2 ){ |
180 |
chpath = (TString)gSystem->ExpandPathName(row2->GetField(0))+'/'; |
if ( KEEPENV ){ |
181 |
|
chpath = (TString)(row2->GetField(0))+'/'; |
182 |
|
} else { |
183 |
|
chpath = (TString)gSystem->ExpandPathName(row2->GetField(0))+'/'; |
184 |
|
}; |
185 |
chfile = (TString)(row2->GetField(1)); |
chfile = (TString)(row2->GetField(1)); |
186 |
chobtts = (UInt_t)atoll(row2->GetField(2)); |
chobtts = (UInt_t)atoll(row2->GetField(2)); |
187 |
chlastts = (UInt_t)atoll(row2->GetField(3)); |
chlastts = (UInt_t)atoll(row2->GetField(3)); |
188 |
chresursts = (UInt_t)atoll(row2->GetField(4)); |
chresursts = (UInt_t)atoll(row2->GetField(4)); |
189 |
chboot = (UInt_t)atoll(row2->GetField(5)); |
chboot = (UInt_t)atoll(row2->GetField(5)); |
190 |
// |
// |
191 |
|
if ( !chboot ){ |
192 |
|
if ( debug ) printf(" mmm... no BOOT number determined by chewbacca, try to find it looking at DB \n"); |
193 |
|
// |
194 |
|
// new feature: check in ROOT_TABLE (timewise) the previous and the next BOOT number != 0 , if it is the same we know this one! |
195 |
|
// if it is different roll back to the old behaviour |
196 |
|
// |
197 |
|
qu.str(""); |
198 |
|
qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND REAL_TIME_INIT>"<< row2->GetField(10) <<" order by REAL_TIME_INIT asc limit 1;"; |
199 |
|
if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str()); |
200 |
|
UInt_t chbootA = 0; |
201 |
|
UInt_t chbootB = 1; |
202 |
|
result3 = conn->Query(qu.str().c_str()); |
203 |
|
if ( result3 ){ |
204 |
|
row3 = result3->Next(); |
205 |
|
if ( row3 ){ |
206 |
|
chbootA = (UInt_t)atoll(row3->GetField(0)); |
207 |
|
if ( debug ) printf(" Found boot_number A = %u \n",chbootA); |
208 |
|
}; |
209 |
|
}; |
210 |
|
delete result3; |
211 |
|
qu.str(""); |
212 |
|
qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND REAL_TIME_INIT<"<< row2->GetField(10) <<" order by REAL_TIME_INIT desc limit 1;"; |
213 |
|
if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str()); |
214 |
|
result3 = conn->Query(qu.str().c_str()); |
215 |
|
if ( result3 ){ |
216 |
|
row3 = result3->Next(); |
217 |
|
if ( row3 ){ |
218 |
|
chbootB = (UInt_t)atoll(row3->GetField(0)); |
219 |
|
if ( debug ) printf(" Found boot_number B = %u \n",chbootB); |
220 |
|
}; |
221 |
|
}; |
222 |
|
if ( chbootA == chbootB ){ |
223 |
|
chboot = chbootA; |
224 |
|
if ( debug ) printf(" Found boot_number! it is %u \n",chboot); |
225 |
|
} else { |
226 |
|
qu.str(""); |
227 |
|
// qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<10000 AND ABS(REAL_TIME_INIT-"<< row2->GetField(10) <<")<7200 group by BOOT_NUMBER order by BOOT_NUMBER asc;"; |
228 |
|
// qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<40000 AND ABS(REAL_TIME_INIT-"<< row2->GetField(10) <<")<40000 group by BOOT_NUMBER order by BOOT_NUMBER asc;"; // relax conditions... 090112 [8RED: error -29] |
229 |
|
qu << "SELECT BOOT_NUMBER from ROOT_TABLE where BOOT_NUMBER!=0 AND ABS(LAST_TIME_SYNC_INFO-"<< chlastts <<")<50000 AND ABS(REAL_TIME_INIT-"<< row2->GetField(10) <<")<50000 group by BOOT_NUMBER order by BOOT_NUMBER asc;"; // even more relaxed conditions... 091214 [9REDtest3: error -29] |
230 |
|
if ( debug ) printf(" chewbacca: query is %s \n",qu.str().c_str()); |
231 |
|
result3 = conn->Query(qu.str().c_str()); |
232 |
|
if ( result3 ){ |
233 |
|
row3 = result3->Next(); |
234 |
|
if ( row3 && result3->GetRowCount() == 1 ){ |
235 |
|
chboot = (UInt_t)atoll(row3->GetField(0)); |
236 |
|
if ( debug ) printf(" Found boot_number = %u \n",chboot); |
237 |
|
} else { |
238 |
|
if ( debug ) printf(" AGH CANNOT DETERMINE THE BOOT NUMBER... \n"); |
239 |
|
throw -29; |
240 |
|
}; |
241 |
|
}; |
242 |
|
}; |
243 |
|
}; |
244 |
|
// |
245 |
chpktinit = (UInt_t)atoll(row2->GetField(6)); |
chpktinit = (UInt_t)atoll(row2->GetField(6)); |
246 |
chpktfinal = (UInt_t)atoll(row2->GetField(7)); |
chpktfinal = (UInt_t)atoll(row2->GetField(7)); |
247 |
chobtinit = (UInt_t)atoll(row2->GetField(8)); |
chobtinit = (UInt_t)atoll(row2->GetField(8)); |
254 |
throw -85; |
throw -85; |
255 |
}; |
}; |
256 |
filerootname = chpath + chfile;// + ".root"; |
filerootname = chpath + chfile;// + ".root"; |
|
if ( debug ) printf(" chewbacca: filename is %s \n",filerootname.Data()); |
|
257 |
}; |
}; |
258 |
this->SetRootName(filerootname); |
this->SetRootName(filerootname); |
|
file = TFile::Open(this->GetRootName().Data()); |
|
|
} else { |
|
|
this->SetRootName(""); |
|
259 |
}; |
}; |
260 |
|
if ( debug ) printf(" Filename is %s \n",filerootname.Data()); |
261 |
|
file = TFile::Open(this->GetRootName().Data()); |
262 |
|
// } else { |
263 |
|
// this->SetRootName(""); |
264 |
|
// this->SetRootName(filerootname); |
265 |
|
// }; |
266 |
} |
} |
267 |
|
|
268 |
// |
// |
462 |
return; |
return; |
463 |
}; |
}; |
464 |
// |
// |
465 |
TString name = this->GetRootFile(); |
if ( !chewbacca ){ |
466 |
Int_t nlength = name.Length(); |
TString name = this->GetRootFile(); |
467 |
if ( nlength < 5 ) return; |
Int_t nlength = name.Length(); |
468 |
TString dwo = 0; |
if ( nlength < 5 ){ |
469 |
for (Int_t i = 0; i<5; i++){ |
if ( IsDebug() ) printf(" Agh problems determining the orbit number! name = %s \n",name.Data()); |
470 |
dwo.Append(name[i],1); |
return; |
471 |
}; |
}; |
472 |
if ( dwo.IsDigit() ){ |
TString dwo = ""; |
473 |
dworbit = (UInt_t)dwo.Atoi(); |
for (Int_t i = 0; i<5; i++){ |
|
} else { |
|
|
dwo=""; |
|
|
for (Int_t i = 8; i<13; i++){ |
|
474 |
dwo.Append(name[i],1); |
dwo.Append(name[i],1); |
475 |
}; |
}; |
476 |
if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi(); |
if ( dwo.IsDigit() ){ |
477 |
|
dworbit = (UInt_t)dwo.Atoi(); |
478 |
|
} else { |
479 |
|
dwo=""; |
480 |
|
for (Int_t i = 8; i<13; i++){ |
481 |
|
dwo.Append(name[i],1); |
482 |
|
}; |
483 |
|
if ( dwo.IsDigit() ) dworbit = (UInt_t)dwo.Atoi(); |
484 |
|
}; |
485 |
|
if ( IsDebug() ) printf(" Downlink orbit is %i (dwo = %s) \n",dworbit,dwo.Data()); |
486 |
}; |
}; |
|
if ( IsDebug() ) printf(" Downlink orbit is %i (dwo = %s) \n",dworbit,dwo.Data()); |
|
487 |
return; |
return; |
488 |
}; |
}; |
489 |
|
|
514 |
|
|
515 |
TString PamelaDBOperations::GetRootPath(){ |
TString PamelaDBOperations::GetRootPath(){ |
516 |
if ( STATIC ){ |
if ( STATIC ){ |
517 |
return((TString)gSystem->DirName(filerootname.Data())+'/'); |
return((TString)gSystem->DirName(this->GetRootName().Data())+'/'); |
518 |
} else { |
} else { |
519 |
return((TString)gSystem->ExpandPathName("$PAM_L0")+'/'); |
if ( KEEPENV ){ |
520 |
|
return((TString)gSystem->ExpandPathName(gSystem->DirName(filerootname.Data()))+'/'); |
521 |
|
} else { |
522 |
|
return((TString)gSystem->ExpandPathName("$PAM_L0")+'/'); |
523 |
|
}; |
524 |
}; |
}; |
525 |
}; |
}; |
526 |
|
|
732 |
UInt_t tjump = 50000; |
UInt_t tjump = 50000; |
733 |
//UInt_t tjump = 100; |
//UInt_t tjump = 100; |
734 |
while ( tjump > 0 ){ |
while ( tjump > 0 ){ |
735 |
pktlast = numeric_limits<UInt_t>::max(); |
// pktlast = numeric_limits<UInt_t>::max(); |
736 |
while ( pktlast > chpktmax && (Int_t)(nevent-1-it) >= 0 ){ |
pktlast = chpktmax + 1; |
737 |
|
while ( PKT(pktlast) > PKT(chpktmax) && (Int_t)(nevent-1-it) >= 0 ){ |
738 |
if ( (Int_t)(nevent-1-it) >= 0 ){ |
if ( (Int_t)(nevent-1-it) >= 0 ){ |
739 |
T->GetEntry(nevent-1-it); |
T->GetEntry(nevent-1-it); |
740 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
742 |
} else { |
} else { |
743 |
pktlast = chpktmax + 1; |
pktlast = chpktmax + 1; |
744 |
}; |
}; |
745 |
if ( (!(it%1000) || abs((int)pktlast - (int)chpktmax)<1000 ) && debug ) printf(" look for up %i %i %i nevent %u (nevent-1-it) %i \n",it,pktlast,chpktmax,nevent,(Int_t)(nevent-1-it)); |
if ( (!(it%1000) || abs((int)(PKT(pktlast) - PKT(chpktmax)))<1000 ) && debug ) printf(" look for up %i %i %i nevent %u (nevent-1-it) %i \n",it,pktlast,chpktmax,nevent,(Int_t)(nevent-1-it)); |
746 |
it += tjump; |
it += tjump; |
747 |
}; |
}; |
748 |
if ( tjump > 1 ) it -= 2*tjump; |
if ( tjump > 1 ) it -= 2*tjump; |
766 |
tjump = 50000; |
tjump = 50000; |
767 |
//tjump = 100; |
//tjump = 100; |
768 |
while ( tjump > 0 ){ |
while ( tjump > 0 ){ |
769 |
pktlast = 0; |
// pktlast = 0; |
770 |
while ( pktlast < chpktmin && it < (Int_t)nevent ){ |
pktlast = chpktmin - 1; |
771 |
|
if ( debug ) printf("LLlook for down %i %i %llu %llu \n",it,pktlast,PKT(pktlast),PKT(chpktmin)); |
772 |
|
while ( PKT(pktlast) < PKT(chpktmin) && it < (Int_t)nevent ){ |
773 |
if ( it < (Int_t)nevent ){ |
if ( it < (Int_t)nevent ){ |
774 |
T->GetEntry(it); |
T->GetEntry(it); |
775 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
863 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
864 |
pktlast = ph->GetCounter(); |
pktlast = ph->GetCounter(); |
865 |
nevent = upperentry - chminentry; |
nevent = upperentry - chminentry; |
866 |
|
if ( ((Int_t)upperentry - (Int_t)chminentry) < 0 ) nevent = 0; // nevent is UInt_t! 090112 [8RED: error -88] |
867 |
|
if ( IsDebug() ) printf(" Setting nevent to %u - upperentry %i chminentry %i \n",nevent,upperentry,chminentry); |
868 |
// |
// |
869 |
} else { |
} else { |
870 |
upperpkt = PKT(pktlast); |
upperpkt = PKT(pktlast); |
873 |
}; |
}; |
874 |
// |
// |
875 |
if ( chewbacca && nevent < 1 ) { |
if ( chewbacca && nevent < 1 ) { |
876 |
|
if ( IsDebug() ) printf(" chewbacca and no events in the file \n"); |
877 |
pktfirst = chpktmin; |
pktfirst = chpktmin; |
878 |
upperpkt = PKT(chpktmax); |
upperpkt = PKT(chpktmax); |
879 |
pktlast = chpktmax; |
pktlast = chpktmax; |
899 |
// if ( nevent < jump ) jump = int(nevent/10); |
// if ( nevent < jump ) jump = int(nevent/10); |
900 |
// if ( !jump ) jump = 1; |
// if ( !jump ) jump = 1; |
901 |
// |
// |
902 |
if ( (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt) && nevent > deltapkt ){ |
|
903 |
|
if ( ((PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) || (labs(PKT(pktlast)-PKT(pktfirst))<deltapkt && labs(OBT(obtlast)-OBT(obtfirst))<deltaobt)) && nevent > deltapkt ){ |
904 |
// |
// |
905 |
if ( IsDebug() ) printf(" starting jump %i \n",jump); |
if ( IsDebug() ) printf(" starting jump %i \n",jump); |
906 |
// if ( IsDebug() ) printf(" (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) %llu < %llu && %llu > %llu \n OR \n (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) %llu > %llu && %llu < %llu \n",PKT(pktlast),PKT(pktfirst),OBT(obtlast),OBT(obtfirst),PKT(pktlast),PKT(pktfirst),OBT(obtlast),OBT(obtfirst)); |
// if ( IsDebug() ) printf(" (PKT(pktlast) < PKT(pktfirst) && OBT(obtlast) > OBT(obtfirst)) %llu < %llu && %llu > %llu \n OR \n (PKT(pktlast) > PKT(pktfirst) && OBT(obtlast) < OBT(obtfirst)) %llu > %llu && %llu < %llu \n",PKT(pktlast),PKT(pktfirst),OBT(obtlast),OBT(obtfirst),PKT(pktlast),PKT(pktfirst),OBT(obtlast),OBT(obtfirst)); |
1215 |
stringstream oss; |
stringstream oss; |
1216 |
// |
// |
1217 |
oss.str(""); |
oss.str(""); |
1218 |
oss << "INSERT INTO _RUNID_GEN VALUES (NULL);"; |
if ( chewbacca ){// if chewbacca and tag=none then use chewbacca tag (chiby = chewbacca inserted by), if chewbacca and tag!=none then use tag, if not chewbacca use tag. |
1219 |
|
if ( !strcmp(tag.Data(),"NONE") ){ |
1220 |
|
oss << "INSERT INTO _RUNID_GEN VALUES (NULL,'"<< chiby.Data() <<"');"; |
1221 |
|
} else { |
1222 |
|
oss << "INSERT INTO _RUNID_GEN VALUES (NULL,'"<< tag.Data() <<"');"; |
1223 |
|
}; |
1224 |
|
} else { |
1225 |
|
oss << "INSERT INTO _RUNID_GEN VALUES (NULL,'"<< tag.Data() <<"');"; |
1226 |
|
}; |
1227 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
1228 |
if ( !result ) throw -10; |
if ( !result ) throw -10; |
1229 |
oss.str(""); |
oss.str(""); |
1428 |
// |
// |
1429 |
// if ( IsDebug() ) printf(" obt conversion: obt is %u obtfirst is %u (numeric_limits<UInt_t>::max()/2) is %u \n",obt,pobtfirst,(UInt_t)(numeric_limits<UInt_t>::max()/2)); |
// if ( IsDebug() ) printf(" obt conversion: obt is %u obtfirst is %u (numeric_limits<UInt_t>::max()/2) is %u \n",obt,pobtfirst,(UInt_t)(numeric_limits<UInt_t>::max()/2)); |
1430 |
// |
// |
1431 |
if ( obt < (pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ){ |
if ( obt < ((Long64_t)pobtfirst/2) && pobtfirst > (numeric_limits<UInt_t>::max()/2) ){ |
1432 |
// if ( IsDebug() ) printf(" rise up obt %lld \n",(Long64_t)obt+(Long64_t)numeric_limits<UInt_t>::max()); |
// if ( IsDebug() ) printf(" rise up obt %lld \n",(Long64_t)obt+(Long64_t)numeric_limits<UInt_t>::max()); |
1433 |
return((Long64_t)obt+(Long64_t)numeric_limits<UInt_t>::max()); |
return((Long64_t)obt+(Long64_t)numeric_limits<UInt_t>::max()); |
1434 |
}; |
}; |
1547 |
glrun->SetRUNTRAILER_PKT(lastPkt); |
glrun->SetRUNTRAILER_PKT(lastPkt); |
1548 |
// |
// |
1549 |
}; |
}; |
1550 |
|
// |
1551 |
|
// ad hoc solution for INF EV_TO [8RED] 090113 |
1552 |
|
// |
1553 |
|
if ( !firstev && lastev == -1 ){ |
1554 |
|
if ( IsDebug() ) printf(" EE firstev %i lastev %i nevents %i \n",firstev,lastev,lastev-firstev+1); |
1555 |
|
firstev = 1; |
1556 |
|
lastev = 0; |
1557 |
|
}; |
1558 |
glrun->SetEV_FROM((UInt_t)firstev); |
glrun->SetEV_FROM((UInt_t)firstev); |
1559 |
glrun->SetEV_TO((UInt_t)lastev); |
glrun->SetEV_TO((UInt_t)lastev); |
1560 |
glrun->SetNEVENTS((UInt_t)lastev-(UInt_t)firstev+1); |
glrun->SetNEVENTS((UInt_t)lastev-(UInt_t)firstev+1); |
1570 |
// PUBLIC FUNCTIONS |
// PUBLIC FUNCTIONS |
1571 |
// |
// |
1572 |
|
|
|
// /** |
|
|
// * Insert a new row into GL_RAW table. |
|
|
// */ |
|
|
// Int_t PamelaDBOperations::insertPamelaRawFile(){ |
|
|
// // |
|
|
// stringstream oss; |
|
|
// // |
|
|
// Bool_t idr = this->SetID_RAW(); |
|
|
// if ( idr ) return(1); |
|
|
// // |
|
|
// oss.str(""); |
|
|
// if ( STATIC ){ |
|
|
// oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('" |
|
|
// << this->GetRawPath().Data() << "', '" << this->GetRawFile().Data() << "')"; |
|
|
// } else { |
|
|
// oss << "INSERT INTO GL_RAW (PATH, NAME) VALUES ('$PAM_RAW', '" << this->GetRawFile().Data() << "')"; |
|
|
// }; |
|
|
// if ( conn->Query(oss.str().c_str()) == 0 ) throw -4; |
|
|
// // |
|
|
// idr = this->SetID_RAW(); |
|
|
// if ( !idr ) throw -11; |
|
|
// // |
|
|
// return(0); |
|
|
// } |
|
1573 |
/** |
/** |
1574 |
* Insert a new row into GL_RAW table. |
* Insert a new row into GL_RAW table. |
1575 |
*/ |
*/ |
1583 |
|
|
1584 |
glraw.PATH = GetRawPath(); |
glraw.PATH = GetRawPath(); |
1585 |
glraw.NAME = GetRawFile(); |
glraw.NAME = GetRawFile(); |
1586 |
glraw.BOOT_NUMBER = 0;//??? |
// glraw.BOOT_NUMBER = 0;//??? |
1587 |
|
glraw.BOOT_NUMBER = this->GetBOOTnumber(); |
1588 |
|
|
1589 |
if( insertPamelaRawFile(&glraw) )return(1); |
if( insertPamelaRawFile(&glraw) )return(1); |
1590 |
// |
// |
1598 |
*/ |
*/ |
1599 |
Int_t PamelaDBOperations::insertPamelaRawFile(GL_RAW *glraw){ |
Int_t PamelaDBOperations::insertPamelaRawFile(GL_RAW *glraw){ |
1600 |
// |
// |
1601 |
if(!glraw)return(1);//?? |
if(!glraw)return(1);//?? ok I think |
1602 |
// |
// |
1603 |
stringstream oss; |
stringstream oss; |
1604 |
// |
// |
1661 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
1662 |
if ( !result ) throw -10; |
if ( !result ) throw -10; |
1663 |
row = result->Next(); |
row = result->Next(); |
1664 |
|
if ( !row ) throw -92; |
1665 |
idresof = (UInt_t)atoll(row->GetField(0)); |
idresof = (UInt_t)atoll(row->GetField(0)); |
1666 |
existsts = true; |
existsts = true; |
1667 |
goto eout; |
goto eout; |
1683 |
if ( !row ){ |
if ( !row ){ |
1684 |
oss.str(""); |
oss.str(""); |
1685 |
oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE),ID FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< " |
oss << "SELECT YEAR(OFFSET_DATE),MONTH(OFFSET_DATE),DAY(OFFSET_DATE),HOUR(OFFSET_DATE),MINUTE(OFFSET_DATE),SECOND(OFFSET_DATE),ID FROM GL_RESURS_OFFSET WHERE FROM_ORBIT< " |
1686 |
<< dworbit << " order by FROM_ORBIT desc limit 1;"; |
<< dworbit << " AND SPECIAL_FILE='' order by FROM_ORBIT desc limit 1;"; |
1687 |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" %s \n",oss.str().c_str()); |
1688 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
1689 |
if ( !result ) throw -10; |
if ( !result ) throw -10; |
1730 |
TYPE = 0; |
TYPE = 0; |
1731 |
TSYNC = 0; |
TSYNC = 0; |
1732 |
// |
// |
1733 |
Double_t minimum = 0.; |
// Double_t minimum = 0.; |
1734 |
Double_t maximum = 0.; |
// Double_t maximum = 0.; |
1735 |
Double_t minimum2 = 0.; |
//Double_t minimum2 = 0.; |
1736 |
Double_t maximum2 = 0.; |
Double_t maximum2 = 0.; |
1737 |
// |
// |
1738 |
// |
// |
1740 |
pamela::McmdRecord *mcrc = 0; |
pamela::McmdRecord *mcrc = 0; |
1741 |
TArrayC *mcmddata = 0; |
TArrayC *mcmddata = 0; |
1742 |
// |
// |
1743 |
minimum = numeric_limits<Double_t>::max(); |
// minimum = numeric_limits<Double_t>::max(); |
1744 |
maximum = numeric_limits<Double_t>::min(); |
// maximum = numeric_limits<Double_t>::min(); |
1745 |
minimum2 = numeric_limits<Double_t>::max(); |
// minimum2 = numeric_limits<Double_t>::max(); |
1746 |
maximum2 = numeric_limits<Double_t>::min(); |
maximum2 = numeric_limits<Double_t>::min(); |
1747 |
// |
// |
1748 |
T = (TTree*)file->Get("Mcmd"); |
T = (TTree*)file->Get("Mcmd"); |
1936 |
return(signal); |
return(signal); |
1937 |
} |
} |
1938 |
|
|
|
// /** |
|
|
// * Insert all the new rows into GL_ROOT. |
|
|
// * The raw file indicates in the parameters should be already been stored in the database. |
|
|
// */ |
|
|
// Int_t PamelaDBOperations::insertPamelaRootFile(){ |
|
|
// stringstream oss; |
|
|
// TSQLResult *result = 0; |
|
|
// TSQLRow *row = 0; |
|
|
// UInt_t idtimesync = 0; |
|
|
// // |
|
|
// // |
|
|
// if ( chewbacca ){ |
|
|
// oss.str(""); |
|
|
// oss << " SELECT ID FROM GL_TIMESYNC where TIMESYNC="<<chlastts<<" AND OBT0="<<chobtts*1000<<" limit 1;"; |
|
|
// if ( debug ) printf(" query is %s \n",oss.str().c_str()); |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// // |
|
|
// if ( !result ) throw -3; |
|
|
// // |
|
|
// row = result->Next(); |
|
|
// // |
|
|
// if ( !row ) throw -3; |
|
|
// idtimesync = (UInt_t)atoll(row->GetField(0)); |
|
|
// } else { |
|
|
// oss.str(""); |
|
|
// if ( STATIC ){ |
|
|
// oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
|
|
// << " LEFT JOIN GL_ROOT " |
|
|
// << " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
|
|
// << " WHERE GL_RAW.PATH = '" << this->GetRawPath().Data() << "' AND " |
|
|
// << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
|
|
// } else { |
|
|
// oss << " SELECT COUNT(GL_ROOT.ID_RAW),GL_RAW.ID,GL_ROOT.ID FROM GL_RAW " |
|
|
// << " LEFT JOIN GL_ROOT " |
|
|
// << " ON GL_RAW.ID = GL_ROOT.ID_RAW " |
|
|
// << " WHERE GL_RAW.PATH = '$PAM_RAW' AND " |
|
|
// << " GL_RAW.NAME = '" << this->GetRawFile().Data() << "' GROUP BY GL_RAW.ID "; |
|
|
// }; |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// // |
|
|
// if ( !result ) throw -12; |
|
|
// // |
|
|
// row = result->Next(); |
|
|
// // |
|
|
// if ( !row ) throw -10; |
|
|
// if ( row != NULL && (UInt_t)atoll(row->GetField(0))>0 ){ |
|
|
// idroot = (UInt_t)atoll(row->GetField(2)); |
|
|
// return(1); |
|
|
// }; |
|
|
// // |
|
|
// // determine which timesync has to be used |
|
|
// // |
|
|
// oss.str(""); |
|
|
// oss << "SELECT GL_TIMESYNC.ID FROM GL_TIMESYNC LEFT JOIN GL_RAW ON GL_RAW.ID = GL_TIMESYNC.ID_RAW ORDER BY GL_TIMESYNC.ID DESC LIMIT 1;"; |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// // |
|
|
// if ( !result ) throw -3; |
|
|
// // |
|
|
// row = result->Next(); |
|
|
// // |
|
|
// if ( !row ) throw -3; |
|
|
// idtimesync = (UInt_t)atoll(row->GetField(0)); |
|
|
// }; |
|
|
// // |
|
|
// oss.str(""); |
|
|
// if ( STATIC ){ |
|
|
// oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
|
|
// << this->GetID_RAW() << "', '" << idtimesync << "', '" << this->GetRootPath().Data() << "', '" << this->GetRootFile().Data() << "')"; |
|
|
// } else { |
|
|
// oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
|
|
// << this->GetID_RAW() << "', '" << idtimesync << "', '$PAM_L0', '" << this->GetRootFile().Data() << "')"; |
|
|
// }; |
|
|
// // |
|
|
// if ( debug ) printf(" query is %s \n",oss.str().c_str()); |
|
|
// if (conn->Query(oss.str().c_str()) == 0) throw -4; |
|
|
// // |
|
|
// delete result; |
|
|
// // |
|
|
// oss.str(""); |
|
|
// // oss << "SELECT ID FROM GL_ROOT WHERE ID_RAW=" << this->GetID_RAW() << ";"; |
|
|
// if ( STATIC ){ |
|
|
// oss << "SELECT ID FROM GL_ROOT WHERE PATH='" << this->GetRootPath().Data() << "' and NAME='"<< this->GetRootFile().Data() <<"';"; |
|
|
// } else { |
|
|
// oss << "SELECT ID FROM GL_ROOT WHERE PATH='$PAM_L0' and NAME='"<< this->GetRootFile().Data() <<"';"; |
|
|
// }; |
|
|
// // |
|
|
// if ( debug ) printf(" query is %s \n",oss.str().c_str()); |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// if ( !result ) throw -12; |
|
|
// row = result->Next(); |
|
|
// if ( !row ) throw -3; |
|
|
// this->SetID_ROOT((UInt_t)atoll(row->GetField(0))); |
|
|
// if ( debug ) printf(" The ID of the ROOT file is %u \n",this->GetID_ROOT()); |
|
|
// // |
|
|
// delete result; |
|
|
// // |
|
|
// return(0); |
|
|
// } |
|
1939 |
/** |
/** |
1940 |
* Insert all the new rows into GL_ROOT. |
* Insert all the new rows into GL_ROOT. |
1941 |
* The raw file indicates in the parameters should be already been stored in the database. |
* The raw file indicates in the parameters should be already been stored in the database. |
1945 |
stringstream oss; |
stringstream oss; |
1946 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
1947 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
1948 |
// |
// |
1949 |
// ---------------------- |
// ---------------------- |
1950 |
// determine the timesync |
// determine the timesync |
1951 |
// ---------------------- |
// ---------------------- |
2017 |
|
|
2018 |
glroot.ID_RAW = GetID_RAW(); |
glroot.ID_RAW = GetID_RAW(); |
2019 |
glroot.ID_TIMESYNC = idtimesync; |
glroot.ID_TIMESYNC = idtimesync; |
2020 |
glroot.PATH = GetRootPath(); |
if ( STATIC ){ |
2021 |
|
glroot.PATH = GetRootPath(); |
2022 |
|
} else { |
2023 |
|
if ( KEEPENV ){ |
2024 |
|
glroot.PATH = gSystem->DirName(filerootname.Data()); |
2025 |
|
} else { |
2026 |
|
glroot.PATH = "$PAM_L0"; |
2027 |
|
}; |
2028 |
|
}; |
2029 |
glroot.NAME = GetRootFile(); |
glroot.NAME = GetRootFile(); |
2030 |
|
|
2031 |
if ( insertPamelaRootFile(&glroot) )return 1; |
if ( insertPamelaRootFile(&glroot) )return 1; |
2040 |
* The raw file indicates in the parameters should be already been stored in the database. |
* The raw file indicates in the parameters should be already been stored in the database. |
2041 |
*/ |
*/ |
2042 |
Int_t PamelaDBOperations::insertPamelaRootFile(GL_ROOT *glroot){ |
Int_t PamelaDBOperations::insertPamelaRootFile(GL_ROOT *glroot){ |
|
|
|
|
|
|
|
|
|
2043 |
stringstream oss; |
stringstream oss; |
2044 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
2045 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
2046 |
// |
// |
2047 |
// |
// |
2048 |
oss.str(""); |
oss.str(""); |
2049 |
if ( STATIC ){ |
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
2050 |
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
<< glroot->ID_RAW << "', '" << glroot->ID_TIMESYNC << "', '" << glroot->PATH << "', '" << glroot->NAME << "')"; |
|
<< glroot->ID_RAW << "', '" << glroot->ID_TIMESYNC << "', '" << glroot->PATH << "', '" << glroot->NAME << "')"; |
|
|
} else { |
|
|
oss << "INSERT INTO GL_ROOT (ID_RAW, ID_TIMESYNC,PATH, NAME) VALUES ('" |
|
|
<< glroot->ID_RAW << "', '" << glroot->ID_TIMESYNC << "', '$PAM_L0', '" << glroot->NAME << "')"; |
|
|
}; |
|
2051 |
// |
// |
2052 |
if ( debug ) printf("%s \n",oss.str().c_str()); |
if ( debug ) printf("%s \n",oss.str().c_str()); |
2053 |
if (conn->Query(oss.str().c_str()) == 0) throw -4; |
if (conn->Query(oss.str().c_str()) == 0) throw -4; |
2056 |
// |
// |
2057 |
oss.str(""); |
oss.str(""); |
2058 |
// oss << "SELECT ID FROM GL_ROOT WHERE ID_RAW=" << this->GetID_RAW() << ";"; |
// oss << "SELECT ID FROM GL_ROOT WHERE ID_RAW=" << this->GetID_RAW() << ";"; |
2059 |
if ( STATIC ){ |
oss << "SELECT ID FROM GL_ROOT WHERE PATH='" << glroot->PATH << "' and NAME='"<< glroot->NAME <<"';"; |
|
oss << "SELECT ID FROM GL_ROOT WHERE PATH='" << glroot->PATH << "' and NAME='"<< glroot->NAME <<"';"; |
|
|
} else { |
|
|
oss << "SELECT ID FROM GL_ROOT WHERE PATH='$PAM_L0' and NAME='"<< glroot->NAME <<"';"; |
|
|
}; |
|
2060 |
// |
// |
2061 |
if ( debug ) printf("%s \n",oss.str().c_str()); |
if ( debug ) printf("%s \n",oss.str().c_str()); |
2062 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
2270 |
// |
// |
2271 |
// no runtrailers in the file! |
// no runtrailers in the file! |
2272 |
// |
// |
2273 |
if ( !rtev ){ |
if ( IsDebug() ) printf(" Start rtev %i rhev %i nrtev %i nrhev %i \n",rtev,rhev,nrtev,nrhev); |
2274 |
if ( !rhev ){ |
if ( !rtev || (chewbacca && !nrtev) ){ |
2275 |
|
if ( IsDebug() ) printf(" No runtrailers \n"); |
2276 |
|
if ( !rhev || (chewbacca && !nrhev) ){ |
2277 |
if ( IsDebug() ) printf(" No runheaders nor runtrailers!! \n"); |
if ( IsDebug() ) printf(" No runheaders nor runtrailers!! \n"); |
2278 |
if ( !(upperentry-chminentry) ){ |
if ( !(upperentry-chminentry) ){ |
2279 |
if ( IsDebug() ) printf(" No physics events nor runs in the file \n"); // di nuovo potrebbe esserci un runtrailer senza eventi (riempimento MM) |
if ( IsDebug() ) printf(" No physics events nor runs in the file \n"); // di nuovo potrebbe esserci un runtrailer senza eventi (riempimento MM) |
2310 |
// |
// |
2311 |
} else { |
} else { |
2312 |
// |
// |
2313 |
|
if ( IsDebug() ) printf(" We have runtrailers \n"); |
2314 |
Int_t conptt = -1; |
Int_t conptt = -1; |
2315 |
for (Int_t ptt=0; ptt<rtev; ptt++){ |
for (Int_t ptt=0; ptt<rtev; ptt++){ |
2316 |
// |
// |
2327 |
ptht = cod->Get(pctp->RunHeader) - 1; |
ptht = cod->Get(pctp->RunHeader) - 1; |
2328 |
// evbeft = cod->Get(pctp->Physics); |
// evbeft = cod->Get(pctp->Physics); |
2329 |
tcod = (UInt_t)cod->Get(pctp->Physics); |
tcod = (UInt_t)cod->Get(pctp->Physics); |
2330 |
if ( !tcod ) tcod = 1; |
if ( !tcod ) tcod = 1; // IS A BUG SINCE RUNS WITH ZERO EVENTS WILL LOOK LIKE RUNS WITH ONE EVENT |
2331 |
evbeft = TMath::Min(upperentry,(tcod-1)); |
evbeft = TMath::Min(upperentry,(tcod-1)); |
2332 |
if ( debug ) printf(" Loop in runtrailers, evbeft is %u upperentry %u cod->getetc %u \n",evbeft,upperentry,cod->Get(pctp->Physics)); |
if ( debug ) printf(" Loop in runtrailers, evbeft is %u upperentry %u cod->getetc %u \n",evbeft,upperentry,cod->Get(pctp->Physics)); |
2333 |
// |
// |
2366 |
obth = phh->GetOrbitalTime(); |
obth = phh->GetOrbitalTime(); |
2367 |
cod = ehh->GetCounter(); |
cod = ehh->GetCounter(); |
2368 |
tcod = (UInt_t)cod->Get(pctp->Physics); |
tcod = (UInt_t)cod->Get(pctp->Physics); |
2369 |
if ( !tcod ) tcod = 1; |
if ( !tcod ) tcod = 1; // IS A BUG SINCE RUNS WITH ZERO EVENTS WILL LOOK LIKE RUNS WITH ONE EVENT |
2370 |
evbefh = TMath::Max(chminentry,(tcod-1)); |
evbefh = TMath::Max(chminentry,(tcod-1)); |
2371 |
// |
// |
2372 |
if ( PKT(pkth) >= PKT(pktfirst) && PKT(pkth) <= upperpkt ){ |
if ( PKT(pkth) >= PKT(pktfirst) && PKT(pkth) <= upperpkt ){ |
2405 |
// |
// |
2406 |
this->HandleRunFragments(true,true,chminentry,(evbefh)); |
this->HandleRunFragments(true,true,chminentry,(evbefh)); |
2407 |
// |
// |
2408 |
|
// what about pth=ptht-1 |
2409 |
}; |
}; |
2410 |
}; |
}; |
2411 |
// |
// |
2428 |
cod = ehh->GetCounter(); |
cod = ehh->GetCounter(); |
2429 |
// evbeft = cod->Get(pctp->Physics); |
// evbeft = cod->Get(pctp->Physics); |
2430 |
tcod = (UInt_t)cod->Get(pctp->Physics); |
tcod = (UInt_t)cod->Get(pctp->Physics); |
2431 |
if ( !tcod ) tcod = 1; |
if ( !tcod ) tcod = 1; // IS A BUG SINCE RUNS WITH ZERO EVENTS WILL LOOK LIKE RUNS WITH ONE EVENT |
2432 |
evbeft = TMath::Min(upperentry,(tcod-1)); |
evbeft = TMath::Min(upperentry,(tcod-1)); |
2433 |
rh->GetEntry(pth); |
rh->GetEntry(pth); |
2434 |
phh = ehh->GetPscuHeader(); |
phh = ehh->GetPscuHeader(); |
2468 |
// if ( ptt+1 == rtev){ |
// if ( ptt+1 == rtev){ |
2469 |
if ( conptt+1 == nrtev ){ |
if ( conptt+1 == nrtev ){ |
2470 |
// if ( conptt+1 == (nrtev+nrtbef )){ |
// if ( conptt+1 == (nrtev+nrtbef )){ |
2471 |
|
if ( IsDebug() ) printf(" >>>>>>>>>>> %i %u %i %u \n",ptht,rhev,nrtev,conptt); |
2472 |
ptht++; |
ptht++; |
2473 |
if ( ptht < rhev ){ |
if ( ptht < rhev ){ |
2474 |
rh->GetEntry(ptht); |
rh->GetEntry(ptht); |
2475 |
|
//pth++; |
2476 |
|
//if ( pth < rhev ){ |
2477 |
|
//rh->GetEntry(pth); |
2478 |
phh = ehh->GetPscuHeader(); |
phh = ehh->GetPscuHeader(); |
2479 |
pkth = phh->GetCounter(); |
pkth = phh->GetCounter(); |
2480 |
obth = phh->GetOrbitalTime(); |
obth = phh->GetOrbitalTime(); |
2498 |
if ( (UInt_t)evbeft < upperentry-1 && upperentry>0 ){ |
if ( (UInt_t)evbeft < upperentry-1 && upperentry>0 ){ |
2499 |
if ( IsDebug() ) printf(" Piece of run at the end of the file with NO RUNHEADER! evbeft %u upperentry-1 %u \n",(UInt_t)evbeft,upperentry-1); |
if ( IsDebug() ) printf(" Piece of run at the end of the file with NO RUNHEADER! evbeft %u upperentry-1 %u \n",(UInt_t)evbeft,upperentry-1); |
2500 |
// |
// |
2501 |
if ( (ptt-1) < 0 ) throw -15; // should never arrive here! |
// if ( (ptt-1) < 0 ) throw -15; // should never arrive here! |
2502 |
rt->GetEntry(ptt-1); |
if ( (rtev-1) < 0 || ptt < 0 ) throw -15; // should never arrive here! |
2503 |
|
// rt->GetEntry(ptt-1); |
2504 |
|
rt->GetEntry(rtev-1); |
2505 |
cod = eht->GetCounter(); |
cod = eht->GetCounter(); |
2506 |
tcod = (UInt_t)cod->Get(pctp->Physics); |
tcod = (UInt_t)cod->Get(pctp->Physics)+1; |
2507 |
evbefh = TMath::Max(chminentry,tcod); |
evbefh = TMath::Max(chminentry,tcod); |
2508 |
// evbefh = cod->Get(pctp->Physics); |
// evbefh = cod->Get(pctp->Physics); |
2509 |
rt->GetEntry(ptt); |
rt->GetEntry(ptt); |
2734 |
// |
// |
2735 |
// we have now the good first piece of a run, fill the glrun object |
// we have now the good first piece of a run, fill the glrun object |
2736 |
// |
// |
2737 |
if ( rhfirstev != firstev && !mishead ) mishead = true; |
if ( firstev != (lastev+1) ){ // could be a problem sometimes (?) |
2738 |
if ( rtlastev != lastev && !mistrail ) mistrail = true; |
if ( rhfirstev != firstev && !mishead ) mishead = true; |
2739 |
|
if ( rtlastev != lastev && !mistrail ) mistrail = true; |
2740 |
|
}; |
2741 |
// |
// |
2742 |
|
if ( IsDebug() ) printf(" bhere firstev is %i lastev is %i \n",firstev,lastev); |
2743 |
this->FillClass(mishead,mistrail,firstev,lastev); |
this->FillClass(mishead,mistrail,firstev,lastev); |
2744 |
|
if ( IsDebug() ) printf(" chere firstev is %i lastev is %i \n",firstev,lastev); |
2745 |
// |
// |
2746 |
if ( IsDebug() ) printf("The run is good, is it the other piece in the GL_RUN_FRAGMENTS table?\n"); |
if ( IsDebug() ) printf("The run is good, is it the other piece in the GL_RUN_FRAGMENTS table?\n"); |
2747 |
if ( IsDebug() ) printf(" C THIS RUN: RUNHEADER_OBT %u RUNTRAILER_OBT %u RUNHEADER_PKT %u RUNTRAILER_PKT %u \n", glrun->GetRUNHEADER_OBT(),glrun->GetRUNTRAILER_OBT(),glrun->GetRUNHEADER_PKT(),glrun->GetRUNTRAILER_PKT()); |
if ( IsDebug() ) printf(" C THIS RUN: RUNHEADER_OBT %u RUNTRAILER_OBT %u RUNHEADER_PKT %u RUNTRAILER_PKT %u \n", glrun->GetRUNHEADER_OBT(),glrun->GetRUNTRAILER_OBT(),glrun->GetRUNHEADER_PKT(),glrun->GetRUNTRAILER_PKT()); |
2783 |
glrun->Fill_GL_RUN_FRAGMENTS(conn); |
glrun->Fill_GL_RUN_FRAGMENTS(conn); |
2784 |
// |
// |
2785 |
} else { |
} else { |
2786 |
if ( IsDebug() ) printf(" The run is already present in the fragment table \n"); |
if ( IsDebug() ) printf(" The run is already present in the fragment table, relaxed %u \n",RELAXED); |
2787 |
if ( PEDANTIC ) throw -69; |
if ( !RELAXED ){ |
2788 |
return; |
if ( PEDANTIC ) throw -69; |
2789 |
|
return; |
2790 |
|
} else { |
2791 |
|
glrun->SetID((UInt_t)atoll(row->GetField(0))); |
2792 |
|
glrun->SetID_RUN_FRAG(0); |
2793 |
|
}; |
2794 |
}; |
}; |
2795 |
// |
// |
2796 |
if ( chewbacca && mishead && mistrail ) goto justcheck; |
if ( chewbacca && mishead && mistrail ) goto justcheck; |
2841 |
// |
// |
2842 |
found = false; // default value |
found = false; // default value |
2843 |
// |
// |
2844 |
if ( IsDebug() ) printf(" Found a possible candidate, checking if it is the good one... \n"); |
if ( IsDebug() ) printf(" A Found a possible candidate, checking if it is the good one... \n"); |
2845 |
// |
// |
2846 |
// if we have both runheader and runtrailer we can check with pkt_counter: |
// if we have both runheader and runtrailer we can check with pkt_counter: |
2847 |
// |
// |
2857 |
if ( IsDebug() ) printf(" FOUND!!! check %llu pktt %llu \n",chkpkt,pktt); |
if ( IsDebug() ) printf(" FOUND!!! check %llu pktt %llu \n",chkpkt,pktt); |
2858 |
// |
// |
2859 |
found = true; |
found = true; |
2860 |
|
if ( IsDebug() ) printf(" where firstev is %i lastev is %i \n",firstev,lastev); |
2861 |
// |
// |
2862 |
} else { |
} else { |
2863 |
// |
// |
2923 |
oss << " ID="<<row->GetField(0)<<";"; |
oss << " ID="<<row->GetField(0)<<";"; |
2924 |
// |
// |
2925 |
glrun1->Query_GL_RUN_FRAGMENTS(oss.str().c_str(),conn); // here we have runheader infos |
glrun1->Query_GL_RUN_FRAGMENTS(oss.str().c_str(),conn); // here we have runheader infos |
2926 |
|
if ( IsDebug() ) printf(" there firstev is %i lastev is %i \n",firstev,lastev); |
2927 |
// |
// |
2928 |
// merge infos |
// merge infos |
2929 |
// |
// |
2954 |
glrun1->SetRUNTRAILER_OBT(glrun->GetRUNTRAILER_OBT()); |
glrun1->SetRUNTRAILER_OBT(glrun->GetRUNTRAILER_OBT()); |
2955 |
glrun1->SetRUNTRAILER_PKT(glrun->GetRUNTRAILER_PKT()); |
glrun1->SetRUNTRAILER_PKT(glrun->GetRUNTRAILER_PKT()); |
2956 |
// |
// |
2957 |
|
if ( IsDebug() ) printf(" here firstev is %i lastev is %i \n",firstev,lastev); |
2958 |
|
// |
2959 |
glrun->SetEV_FROM(firstev); |
glrun->SetEV_FROM(firstev); |
2960 |
glrun->SetNEVENTS(lastev-firstev+1); |
glrun->SetNEVENTS(lastev-firstev+1); |
2961 |
// |
// |
3048 |
// |
// |
3049 |
found = false; // default value |
found = false; // default value |
3050 |
// |
// |
3051 |
if ( IsDebug() ) printf(" Found a possible candidate, checking if it is the good one... \n"); |
if ( IsDebug() ) printf(" B Found a possible candidate, checking if it is the good one... \n"); |
3052 |
// |
// |
3053 |
// if we have both runheader and runtrailer we can check with pkt_counter: |
// if we have both runheader and runtrailer we can check with pkt_counter: |
3054 |
// |
// |
3148 |
apkt = PKT(ph->GetCounter()); |
apkt = PKT(ph->GetCounter()); |
3149 |
aobt = OBT(ph->GetOrbitalTime()); |
aobt = OBT(ph->GetOrbitalTime()); |
3150 |
lastev--; |
lastev--; |
3151 |
if ( PEDANTIC ) throw -72; |
if ( PEDANTIC && !RELAXED) throw -72; |
3152 |
}; |
}; |
3153 |
if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev); |
if ( IsDebug() ) printf(" Check overlapping events done: %u %u %llu %llu lastev is %i\n",apkt,bpkt,aobt,bobt,lastev); |
3154 |
// |
// |
3216 |
// not found, has this run already inserted in the GL_RUN or in the GL_RUN_FRAGMENTS table? |
// not found, has this run already inserted in the GL_RUN or in the GL_RUN_FRAGMENTS table? |
3217 |
// |
// |
3218 |
oss.str(""); |
oss.str(""); |
3219 |
oss << " SELECT ID FROM GL_RUN WHERE " |
oss << " SELECT ID,ID_ROOT_L0 FROM GL_RUN WHERE " |
3220 |
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND (" |
<< " BOOT_NUMBER=" << this->GetBOOTnumber() << " AND (" |
3221 |
<< " (RUNHEADER_TIME>=" << (UInt_t)(glrun->GetRUNHEADER_TIME()-10) << " AND " |
<< " (RUNHEADER_TIME>=" << (UInt_t)(glrun->GetRUNHEADER_TIME()-10) << " AND " |
3222 |
<< " RUNTRAILER_TIME<=" << (UInt_t)(glrun->GetRUNTRAILER_TIME()+10) << " AND (" |
<< " RUNTRAILER_TIME<=" << (UInt_t)(glrun->GetRUNTRAILER_TIME()+10) << " AND (" |
3239 |
row = result->Next(); |
row = result->Next(); |
3240 |
// |
// |
3241 |
if ( row ){ |
if ( row ){ |
3242 |
if ( IsDebug() ) printf(" The run is already present in the GL_RUN table \n"); |
// |
3243 |
if ( PEDANTIC ) throw -70; |
// we end up here if chewbacca and we have RH---| small gap |---RT, in this case the two pieces are attached and moved to GL_RUN. We are now processing the |
3244 |
|
// "small gap" piece... we recognize this since: we have two entries from this query, the pkt number is consistent with our pkt number. |
3245 |
|
// |
3246 |
|
// 090112 [8RED (-70): RUN ALREADY INSERTED] |
3247 |
|
// 093012 [8RED (-70): RUN ALREADY INSERTED] try to be more general... the run has not been instered yet if the runheader-trailer PKT interval is not covered by already inserted runs, NOTICE: we must look in the GL_RUN_TRASH table! |
3248 |
|
// |
3249 |
|
Bool_t OK = false; |
3250 |
|
UInt_t IDRL2A = 0; |
3251 |
|
UInt_t IDRL2B = 0; |
3252 |
|
if ( chewbacca ){ |
3253 |
|
if ( result->GetRowCount() >= 2 ){ |
3254 |
|
// |
3255 |
|
if ( IsDebug() ) printf(" RH---| gap |---RT case, number of pieces in the GL_RUN table: %i \n",result->GetRowCount()); |
3256 |
|
OK = true; |
3257 |
|
IDRL2A = (UInt_t)atoll(row->GetField(0)); |
3258 |
|
// |
3259 |
|
while ( row ){ |
3260 |
|
// |
3261 |
|
TSQLResult *trresult = 0; |
3262 |
|
TSQLRow *trrow = 0; |
3263 |
|
// |
3264 |
|
stringstream tross; |
3265 |
|
tross.str(""); |
3266 |
|
// |
3267 |
|
tross << " SELECT RUNHEADER_PKT, RUNTRAILER_PKT FROM GL_RUN_TRASH where ID=" << row->GetField(0) << ";"; |
3268 |
|
if ( IsDebug() ) printf(" check in the gl_run_trash table for pkt intervals: query is \n %s \n",tross.str().c_str()); |
3269 |
|
trresult = conn->Query(tross.str().c_str()); |
3270 |
|
if ( !trresult ) throw -4; |
3271 |
|
trrow = trresult->Next(); |
3272 |
|
if ( !trrow || trresult->GetRowCount() != 1 ){ |
3273 |
|
OK = false; |
3274 |
|
if ( IsDebug() ) printf(" OPS! no such run (or multiple runs!) in GL_RUN_TRASH table, something wrong is going on! \n"); |
3275 |
|
break; |
3276 |
|
}; |
3277 |
|
// |
3278 |
|
UInt_t pktH = (UInt_t)atoll(row->GetField(0)); |
3279 |
|
UInt_t pktT = (UInt_t)atoll(row->GetField(1)); |
3280 |
|
delete trresult; |
3281 |
|
if ( |
3282 |
|
(PKT(pktH) >= PKT(glrun->GetRUNHEADER_PKT()) && PKT(pktH) <= PKT(glrun->GetRUNTRAILER_PKT())) || |
3283 |
|
(PKT(pktT) >= PKT(glrun->GetRUNHEADER_PKT()) && PKT(pktT) <= PKT(glrun->GetRUNTRAILER_PKT())) || |
3284 |
|
(PKT(pktH) <= PKT(glrun->GetRUNTRAILER_PKT()) && PKT(pktT) >= PKT(glrun->GetRUNTRAILER_PKT())) ){ |
3285 |
|
OK = false; |
3286 |
|
if ( IsDebug() ) printf(" The run is overlapping with already inserted ones!! \n"); |
3287 |
|
break; |
3288 |
|
}; |
3289 |
|
//PKT(glrun->GetRUNHEADER_PKT()) PKT(glrun->GetRUNTRAILER_PKT()) |
3290 |
|
row = result->Next(); |
3291 |
|
}; |
3292 |
|
|
3293 |
|
// if ( result->GetRowCount() == 2 ) { |
3294 |
|
// IDRL2A = (UInt_t)atoll(row->GetField(0)); |
3295 |
|
// UInt_t IDRL0A = (UInt_t)atoll(row->GetField(1)); |
3296 |
|
// row = result->Next(); |
3297 |
|
// IDRL2B = (UInt_t)atoll(row->GetField(0)); |
3298 |
|
// UInt_t IDRL0B = (UInt_t)atoll(row->GetField(1)); |
3299 |
|
// if ( IsDebug() ) printf(" IDRL0A %u B %u IDRL2A %u B %u \n",IDRL0A,IDRL0B,IDRL2A,IDRL2B); |
3300 |
|
// // if ( IDRL0A == IDRL0B ){ //091230 why this condition???? RH could be in a file and RT in another one... |
3301 |
|
// if ( true ){ |
3302 |
|
// TSQLResult *result2 = 0; |
3303 |
|
// TSQLRow *row2 = 0; |
3304 |
|
// oss.str(""); |
3305 |
|
// oss << "select PKT_NUMBER_FINAL from ROOT_TABLE_MERGING where REAL_TIME_INIT<" << chrtinit << " order by REAL_TIME_INIT desc limit 1;"; |
3306 |
|
// if ( IsDebug() ) printf(" Check if we are in the case RH---| small gap |---RT: query is \n %s \n",oss.str().c_str()); |
3307 |
|
// result2 = conn->Query(oss.str().c_str()); |
3308 |
|
// // |
3309 |
|
// if ( !result2 ) throw -4; |
3310 |
|
// // |
3311 |
|
// row2 = result2->Next(); |
3312 |
|
// // |
3313 |
|
// if ( row2 ){ |
3314 |
|
// UInt_t PKA = (UInt_t)atoll(row2->GetField(0)); |
3315 |
|
// delete result2; |
3316 |
|
// oss.str(""); |
3317 |
|
// oss << "select PKT_NUMBER_INIT from ROOT_TABLE_MERGING where REAL_TIME_INIT>" << chrtinit << " order by REAL_TIME_INIT asc limit 1;"; |
3318 |
|
// if ( IsDebug() ) printf(" Check if we are in the case RH---| small gap |---RT: query is \n %s \n",oss.str().c_str()); |
3319 |
|
// result2 = conn->Query(oss.str().c_str()); |
3320 |
|
// // |
3321 |
|
// if ( !result2 ) throw -4; |
3322 |
|
// // |
3323 |
|
// row2 = result2->Next(); |
3324 |
|
// // |
3325 |
|
// if ( row2 ){ |
3326 |
|
// UInt_t PKB = (UInt_t)atoll(row2->GetField(0)); |
3327 |
|
// // |
3328 |
|
// if ( IsDebug() ) printf(" PKT(PKA) + 1 %llu == runheaderpkt %llu && PKB = runtrailer %llu + 1 %llu \n",PKT(PKA)+1LL, PKT(glrun->GetRUNHEADER_PKT()),PKT(PKB), PKT(glrun->GetRUNTRAILER_PKT())+1LL); |
3329 |
|
// if ( PKT(PKA)+1LL == PKT(glrun->GetRUNHEADER_PKT()) && PKT(PKB) == PKT(glrun->GetRUNTRAILER_PKT())+1LL ){ |
3330 |
|
// if ( IsDebug() ) printf(" Ok, we are in the case: RH---| small gap |---RT \n"); |
3331 |
|
// OK = true; |
3332 |
|
// }; |
3333 |
|
// }; |
3334 |
|
// }; |
3335 |
|
// }; |
3336 |
|
// }; |
3337 |
|
}; |
3338 |
|
}; |
3339 |
|
if ( OK ){ |
3340 |
|
// |
3341 |
|
// this is the case in which we must insert a piece of run between two fragments in the GL_RUN table // we arrive here ONLY with internal pieces of runs, never runheader or runtrailer by definition |
3342 |
|
// |
3343 |
|
// we have to update with runheader/trailer infos our run and we have to change the ROOT_ID_FRAG of the header piece ( 1/2 2/1 must become 1/3 3/2 2/1 ) |
3344 |
|
// |
3345 |
|
if ( IsDebug() ) printf(" Ok, we are in the case: RH---| small gap |---RT \n"); |
3346 |
|
// |
3347 |
|
GL_RUN *glA = new GL_RUN(); |
3348 |
|
glA->Query_GL_RUN(IDRL2A,conn); |
3349 |
|
// |
3350 |
|
glrun->SetRUNHEADER_TIME(glA->GetRUNHEADER_TIME()); |
3351 |
|
glrun->SetRUNHEADER_OBT(glA->GetRUNHEADER_OBT()); |
3352 |
|
glrun->SetRUNHEADER_PKT(glA->GetRUNHEADER_PKT()); |
3353 |
|
// |
3354 |
|
glrun->SetRUNTRAILER_TIME(glA->GetRUNTRAILER_TIME()); |
3355 |
|
glrun->SetRUNTRAILER_OBT(glA->GetRUNTRAILER_OBT()); |
3356 |
|
glrun->SetRUNTRAILER_PKT(glA->GetRUNTRAILER_PKT()); |
3357 |
|
// |
3358 |
|
if ( glA->GetPHYSENDRUN_MASK_S3S2S12() ) glrun->SetPHYSENDRUN_MASK_S3S2S12(glA->GetPHYSENDRUN_MASK_S3S2S12()); |
3359 |
|
if ( glA->GetPHYSENDRUN_MASK_S11CRC() ) glrun->SetPHYSENDRUN_MASK_S11CRC(glA->GetPHYSENDRUN_MASK_S11CRC()); |
3360 |
|
// |
3361 |
|
if ( glA->GetACQ_BUILD_INFO() != 0 ){ |
3362 |
|
// |
3363 |
|
// the first piece contains a good runheader we can update all the other runs with correct infos! |
3364 |
|
// |
3365 |
|
glrun->SetCOMPILATIONTIMESTAMP(glA->GetCOMPILATIONTIMESTAMP()); |
3366 |
|
glrun->SetFAV_WRK_SCHEDULE(glA->GetFAV_WRK_SCHEDULE()); |
3367 |
|
glrun->SetEFF_WRK_SCHEDULE(glA->GetEFF_WRK_SCHEDULE()); |
3368 |
|
glrun->SetPRH_VAR_TRG_MODE_A(glA->GetPRH_VAR_TRG_MODE_A()); |
3369 |
|
glrun->SetPRH_VAR_TRG_MODE_B(glA->GetPRH_VAR_TRG_MODE_B()); |
3370 |
|
glrun->SetACQ_BUILD_INFO(glA->GetACQ_BUILD_INFO()); |
3371 |
|
glrun->SetACQ_VAR_INFO(glA->GetACQ_VAR_INFO()); |
3372 |
|
glrun->SetRM_ACQ_AFTER_CALIB(glA->GetRM_ACQ_AFTER_CALIB()); |
3373 |
|
glrun->SetRM_ACQ_SETTING_MODE(glA->GetRM_ACQ_SETTING_MODE()); |
3374 |
|
glrun->SetTRK_CALIB_USED(glA->GetTRK_CALIB_USED()); |
3375 |
|
glrun->SetCAL_DSP_MASK(glA->GetCAL_DSP_MASK()); |
3376 |
|
glrun->SetLAST_TIMESYNC(glA->GetLAST_TIMESYNC()); |
3377 |
|
glrun->SetOBT_TIMESYNC(glA->GetOBT_TIMESYNC()); |
3378 |
|
// |
3379 |
|
}; |
3380 |
|
// |
3381 |
|
if ( glA->GetPKT_READY_COUNTER() != 0 ){ |
3382 |
|
// |
3383 |
|
// the first piece contains a good runtrailer we can update all the other runs with correct infos! |
3384 |
|
// |
3385 |
|
glrun->SetPKT_COUNTER(glA->GetPKT_COUNTER()); |
3386 |
|
glrun->SetPKT_READY_COUNTER(glA->GetPKT_READY_COUNTER()); |
3387 |
|
}; |
3388 |
|
// |
3389 |
|
// update runheader ROOT_ID_FRAG |
3390 |
|
// |
3391 |
|
oss.str(""); |
3392 |
|
oss << "UPDATE GL_RUN SET ID_RUN_FRAG = " << glrun->GetID() << " where ID = " << IDRL2A << ";"; |
3393 |
|
if ( IsDebug() ) printf(" update gl_run to have cross indexing: %s\n",oss.str().c_str()); |
3394 |
|
conn->Query(oss.str().c_str()); |
3395 |
|
// |
3396 |
|
oss.str(""); |
3397 |
|
oss << "SELECT ID FROM GL_RUN WHERE ID_RUN_FRAG = " << IDRL2A << ";"; |
3398 |
|
if ( IsDebug() ) printf(" update gl_run to have cross indexing: %s\n",oss.str().c_str()); |
3399 |
|
TSQLResult *tsresult = 0; |
3400 |
|
TSQLRow *tsrow = 0; |
3401 |
|
tsresult=conn->Query(oss.str().c_str()); |
3402 |
|
if ( !tsresult ) throw -4; |
3403 |
|
tsrow = tsresult->Next(); |
3404 |
|
if ( !tsrow || tsresult->GetRowCount() != 1 ) throw -4; |
3405 |
|
IDRL2B = (UInt_t)atoll(tsrow->GetField(0)); |
3406 |
|
glrun->SetID_RUN_FRAG(IDRL2B); |
3407 |
|
// |
3408 |
|
// fill the new run in GL_RUN |
3409 |
|
// |
3410 |
|
glrun->Fill_GL_RUN(conn); |
3411 |
|
// |
3412 |
|
glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS"); |
3413 |
|
delete glA; |
3414 |
|
// |
3415 |
|
} else { |
3416 |
|
// // |
3417 |
|
// // this is the case in which we must insert a piece of run between two fragments in the GL_RUN table |
3418 |
|
// // |
3419 |
|
// // we have to update with runheader/trailer infos our run and we have to change the ROOT_ID_FRAG of the header piece ( 1/2 2/1 must become 1/3 3/2 2/1 ) |
3420 |
|
// // |
3421 |
|
// GL_RUN *glA = new GL_RUN(); |
3422 |
|
// glA->Query_GL_RUN(IDRL2A,conn); |
3423 |
|
// // |
3424 |
|
// if ( glA->GetACQ_BUILD_INFO() != 0 ){ |
3425 |
|
// // |
3426 |
|
// // the first piece contains a good runheader we can update all the other runs with correct infos! |
3427 |
|
// // |
3428 |
|
// oss.str(""); |
3429 |
|
// oss << "UPDATE GL_RUN_FRAGMENTS SET " |
3430 |
|
// << " RUNHEADER_TIME=" << glA->GetRUNHEADER_TIME()<< " , " |
3431 |
|
// << " RUNHEADER_OBT=" << glA->GetRUNHEADER_OBT()<< " , " |
3432 |
|
// << " RUNHEADER_PKT=" << glA->GetRUNHEADER_PKT()<< " , " |
3433 |
|
// << " COMPILATIONTIMESTAMP=" << glA->GetCOMPILATIONTIMESTAMP()<< " , " |
3434 |
|
// << " FAV_WRK_SCHEDULE=" << glA->GetFAV_WRK_SCHEDULE()<< " , " |
3435 |
|
// << " EFF_WRK_SCHEDULE=" << glA->GetEFF_WRK_SCHEDULE()<< " , " |
3436 |
|
// << " PRH_VAR_TRG_MODE_A=" << glA->GetPRH_VAR_TRG_MODE_A()<< " , " |
3437 |
|
// << " PRH_VAR_TRG_MODE_B=" << glA->GetPRH_VAR_TRG_MODE_B()<< " , " |
3438 |
|
// << " ACQ_BUILD_INFO=" << glA->GetACQ_BUILD_INFO()<< " , " |
3439 |
|
// << " ACQ_VAR_INFO=" << glA->GetACQ_VAR_INFO()<< " , " |
3440 |
|
// << " RM_ACQ_AFTER_CALIB=" << glA->GetRM_ACQ_AFTER_CALIB()<< " , " |
3441 |
|
// << " RM_ACQ_SETTING_MODE=" << glA->GetRM_ACQ_SETTING_MODE()<< " , " |
3442 |
|
// << " TRK_CALIB_USED=" << glA->GetTRK_CALIB_USED()<< " , " |
3443 |
|
// << " CAL_DSP_MASK=" << glA->GetCAL_DSP_MASK()<< " , " |
3444 |
|
// << " LAST_TIMESYNC=" << glA->GetLAST_TIMESYNC()<< " , "; |
3445 |
|
// // |
3446 |
|
// if ( glA->GetPHYSENDRUN_MASK_S3S2S12() ) |
3447 |
|
// oss << " PHYSENDRUN_MASK_S3S2S12=" << glA->GetPHYSENDRUN_MASK_S3S2S12() << " , "; |
3448 |
|
// if ( glA->GetPHYSENDRUN_MASK_S11CRC() ) |
3449 |
|
// oss << " PHYSENDRUN_MASK_S11CRC=" << glA->GetPHYSENDRUN_MASK_S11CRC() << " , "; |
3450 |
|
// // |
3451 |
|
// oss << " OBT_TIMESYNC=" << glA->GetOBT_TIMESYNC(); |
3452 |
|
// oss << " WHERE ID=" << glrun->GetID() << ";"; |
3453 |
|
// if ( IsDebug() ) printf(" update with correct infos: %s\n",oss.str().c_str()); |
3454 |
|
// conn->Query(oss.str().c_str()); |
3455 |
|
// // |
3456 |
|
// } else { |
3457 |
|
// // |
3458 |
|
// // sig no runheader, let set anyway what is possible... |
3459 |
|
// // |
3460 |
|
// oss.str(""); |
3461 |
|
// oss << "UPDATE GL_RUN_FRAGMENTS SET " |
3462 |
|
// << " RUNHEADER_TIME=" << glA->GetRUNHEADER_TIME()<< " , "; |
3463 |
|
// // |
3464 |
|
// if ( glA->GetPHYSENDRUN_MASK_S3S2S12() ) |
3465 |
|
// oss << " PHYSENDRUN_MASK_S3S2S12=" << glA->GetPHYSENDRUN_MASK_S3S2S12()<< " , "; |
3466 |
|
// if ( glA->GetPHYSENDRUN_MASK_S11CRC() ) |
3467 |
|
// oss << " PHYSENDRUN_MASK_S11CRC=" << glA->GetPHYSENDRUN_MASK_S11CRC()<< " , "; |
3468 |
|
// // |
3469 |
|
// oss << " RUNHEADER_OBT=" << glA->GetRUNHEADER_OBT()<< " , " |
3470 |
|
// << " RUNHEADER_PKT=" << glA->GetRUNHEADER_PKT(); |
3471 |
|
// oss << " WHERE ID=" << glrun->GetID() << ";"; |
3472 |
|
// if ( IsDebug() ) printf(" update with correct infos2: %s\n",oss.str().c_str()); |
3473 |
|
// conn->Query(oss.str().c_str()); |
3474 |
|
// }; |
3475 |
|
// // |
3476 |
|
// // update runheader ROOT_ID_FRAG |
3477 |
|
// // |
3478 |
|
// oss.str(""); |
3479 |
|
// oss << "UPDATE GL_RUN SET ID_RUN_FRAG = " << glrun->GetID() << " where ID = " << IDRL2A << ";"; |
3480 |
|
// if ( IsDebug() ) printf(" update gl_run to have cross indexing: %s\n",oss.str().c_str()); |
3481 |
|
// conn->Query(oss.str().c_str()); |
3482 |
|
// // |
3483 |
|
// // now let's look for runtrailer if any in the last run |
3484 |
|
// // |
3485 |
|
// glA->Query_GL_RUN(IDRL2B,conn); |
3486 |
|
// // |
3487 |
|
// if ( glA->GetPKT_READY_COUNTER() != 0 ){ |
3488 |
|
// // |
3489 |
|
// // the first piece contains a good runtrailer we can update all the other runs with correct infos! |
3490 |
|
// // |
3491 |
|
// oss.str(""); |
3492 |
|
// oss << "UPDATE GL_RUN_FRAGMENTS SET " |
3493 |
|
// << " RUNTRAILER_TIME=" << glA->GetRUNTRAILER_TIME()<< " , " |
3494 |
|
// << " RUNTRAILER_OBT=" << glA->GetRUNTRAILER_OBT()<< " , " |
3495 |
|
// << " RUNTRAILER_PKT=" << glA->GetRUNTRAILER_PKT()<< " , " |
3496 |
|
// << " PKT_COUNTER=" << glA->GetPKT_COUNTER()<< " , "; |
3497 |
|
// // |
3498 |
|
// if ( glA->GetPHYSENDRUN_MASK_S3S2S12() ){ |
3499 |
|
// oss << " PHYSENDRUN_MASK_S3S2S12=" << glA->GetPHYSENDRUN_MASK_S3S2S12()<< " , "; }; |
3500 |
|
// if ( glA->GetPHYSENDRUN_MASK_S11CRC() ) { |
3501 |
|
// oss << " PHYSENDRUN_MASK_S11CRC=" << glA->GetPHYSENDRUN_MASK_S11CRC()<< " , "; }; |
3502 |
|
// // |
3503 |
|
// oss << " PKT_READY_COUNTER=" << glA->GetPKT_READY_COUNTER() |
3504 |
|
// << " WHERE ID=" << glrun->GetID() << ";"; |
3505 |
|
// if ( IsDebug() ) printf(" update with correct trailer infos: %s\n",oss.str().c_str()); |
3506 |
|
// conn->Query(oss.str().c_str()); |
3507 |
|
// // |
3508 |
|
// } else { |
3509 |
|
// // |
3510 |
|
// // sig no runtrailer, let set anyway what is possible... |
3511 |
|
// // |
3512 |
|
// oss.str(""); |
3513 |
|
// oss << "UPDATE GL_RUN_FRAGMENTS SET " |
3514 |
|
// << " RUNTRAILER_TIME=" << glA->GetRUNTRAILER_TIME()<< " , " |
3515 |
|
// << " RUNTRAILER_OBT=" << glA->GetRUNTRAILER_OBT()<< " , "; |
3516 |
|
// // |
3517 |
|
// if ( glA->GetPHYSENDRUN_MASK_S3S2S12() ){ |
3518 |
|
// oss << " PHYSENDRUN_MASK_S3S2S12=" << glA->GetPHYSENDRUN_MASK_S3S2S12()<< " , "; }; |
3519 |
|
// if ( glA->GetPHYSENDRUN_MASK_S11CRC() ){ |
3520 |
|
// oss << " PHYSENDRUN_MASK_S11CRC=" << glA->GetPHYSENDRUN_MASK_S11CRC()<< " , "; }; |
3521 |
|
// // |
3522 |
|
// oss << " RUNTRAILER_PKT=" << glrun->GetRUNTRAILER_PKT() |
3523 |
|
// << " WHERE ID=" << glrun->GetID() << ";"; |
3524 |
|
// if ( IsDebug() ) printf(" update with correct trailer infos2: %s\n",oss.str().c_str()); |
3525 |
|
// conn->Query(oss.str().c_str()); |
3526 |
|
// }; |
3527 |
|
// // |
3528 |
|
// UInt_t myi = glrun->GetID(); |
3529 |
|
// oss.str(""); |
3530 |
|
// oss << " ID= "<< myi; |
3531 |
|
// // |
3532 |
|
// glrun->Query_GL_RUN_FRAGMENTS(oss.str().c_str(),conn); |
3533 |
|
// // |
3534 |
|
// // fill the new run in GL_RUN |
3535 |
|
// // |
3536 |
|
// glrun->SetID_RUN_FRAG(IDRL2B); |
3537 |
|
// glrun->Fill_GL_RUN(conn); |
3538 |
|
// glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS"); |
3539 |
|
// delete glA; |
3540 |
|
// // |
3541 |
|
// } else { |
3542 |
|
// |
3543 |
|
// is just a repetition |
3544 |
|
// |
3545 |
|
if ( IsDebug() ) printf(" The run is already present in the GL_RUN table \n"); |
3546 |
|
if ( PEDANTIC ) throw -70; |
3547 |
|
}; |
3548 |
|
// |
3549 |
} else { |
} else { |
3550 |
if ( NoFrag() ){ |
if ( NoFrag() ){ |
3551 |
glrun->SetID_RUN_FRAG(glrun->GetID()); |
glrun->SetID_RUN_FRAG(glrun->GetID()); |
3582 |
glrun->SetID(this->AssignRunID()); |
glrun->SetID(this->AssignRunID()); |
3583 |
glrun->SetID_RUN_FRAG(0); |
glrun->SetID_RUN_FRAG(0); |
3584 |
glrun->Fill_GL_RUN(conn); // it'ok we arrive here only inside a file hence in the middle of the runs... |
glrun->Fill_GL_RUN(conn); // it'ok we arrive here only inside a file hence in the middle of the runs... |
3585 |
}; |
}; |
3586 |
// |
// |
3587 |
}; |
}; |
3588 |
// |
// |
3598 |
// |
// |
3599 |
EventCounter *code=0; |
EventCounter *code=0; |
3600 |
// |
// |
3601 |
UInt_t nevent = 0; |
// UInt_t nevent = 0; |
3602 |
UInt_t checkfirst = 0; |
UInt_t checkfirst = 0; |
3603 |
UInt_t checklast = 0; |
UInt_t checklast = 0; |
3604 |
UInt_t firstentry = 0; |
UInt_t firstentry = 0; |
3620 |
EventHeader *eh = 0; |
EventHeader *eh = 0; |
3621 |
PscuHeader *ph = 0; |
PscuHeader *ph = 0; |
3622 |
T->SetBranchAddress("Header", &eh); |
T->SetBranchAddress("Header", &eh); |
3623 |
nevent = T->GetEntries(); |
// nevent = T->GetEntries(); |
3624 |
// |
// |
3625 |
// |
// |
3626 |
if ( firstev == lastev+1 || lastev == firstev ) { // no events inside the run! |
if ( firstev == lastev+1 || lastev == firstev ) { // no events inside the run! |
3627 |
//if ( firstev <= lastev+1 ) { // no events inside the run! |
//if ( firstev <= lastev+1 ) { // no events inside the run! |
3628 |
if ( IsDebug() ) printf(" Checking but no events in the run! \n"); |
if ( IsDebug() ) printf(" Checking but no events in the run! \n"); |
3629 |
|
firstev = lastev+1; // this is necessary for files with no Physics entries, should have no influence on other files |
3630 |
// return true is correct |
// return true is correct |
3631 |
return(true); |
return(true); |
3632 |
// |
// |
3656 |
// |
// |
3657 |
} else { |
} else { |
3658 |
// |
// |
3659 |
if ( IsDebug() ) printf(" There are no-phyics packets inside the run!\n"); |
if ( IsDebug() ) printf(" There are no-physics packets inside the run!\n"); |
3660 |
// |
// |
3661 |
// HERE WE MUST HANDLE THAT RUNS AND GO BACK |
// HERE WE MUST HANDLE THAT RUNS AND GO BACK |
3662 |
// |
// |
3666 |
UInt_t check = 0; |
UInt_t check = 0; |
3667 |
UInt_t lastevtemp = lastev; |
UInt_t lastevtemp = lastev; |
3668 |
UInt_t firstevno = firstev; |
UInt_t firstevno = firstev; |
3669 |
|
UInt_t rhchko=0; |
3670 |
|
UInt_t rhchk=0; |
3671 |
// |
// |
3672 |
for (UInt_t i=firstev; i<=lastev; i++){ |
for (UInt_t i=firstev; i<=lastev; i++){ |
3673 |
// |
// |
3676 |
// |
// |
3677 |
check = 0; |
check = 0; |
3678 |
// |
// |
3679 |
|
|
3680 |
|
// if we have a runheader set lastev then exit |
3681 |
|
// |
3682 |
for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){ |
for(Iter = packetsNames.begin(); Iter != packetsNames.end(); Iter++){ |
3683 |
if ( strcmp(*Iter,"Physics") ) check += code->Get(GetPacketType(*Iter)); |
if ( strcmp(*Iter,"Physics") ) check += code->Get(GetPacketType(*Iter)); |
3684 |
}; |
}; |
3685 |
|
// check here if we have a runheader |
3686 |
|
rhchko = rhchk; |
3687 |
|
rhchk = code->Get(GetPacketType("RunHeader")); |
3688 |
// |
// |
3689 |
if ( checkfirst < check || i == lastev ){ |
if ( checkfirst < check || i == lastev ){ |
3690 |
// |
// |
3768 |
// this->HandleRunFragments(true,mistrail,firstentry,lastentry); // cannot call it here since it enters a loop which will destroy the already stored variables if we arrive here from HandleRunFragments |
// this->HandleRunFragments(true,mistrail,firstentry,lastentry); // cannot call it here since it enters a loop which will destroy the already stored variables if we arrive here from HandleRunFragments |
3769 |
// |
// |
3770 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3771 |
mishead = true; |
mishead = true; |
3772 |
|
|
3773 |
|
|
3882 |
// |
// |
3883 |
found = false; // default value |
found = false; // default value |
3884 |
// |
// |
3885 |
if ( IsDebug() ) printf(" Found a possible candidate, checking if it is the good one... \n"); |
if ( IsDebug() ) printf(" C Found a possible candidate, checking if it is the good one... \n"); |
3886 |
// |
// |
3887 |
// if we have both runheader and runtrailer we can check with pkt_counter: |
// if we have both runheader and runtrailer we can check with pkt_counter: |
3888 |
// |
// |
4087 |
}; // EEE |
}; // EEE |
4088 |
|
|
4089 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4090 |
} else { |
} else { |
4091 |
if ( !IsRunAlreadyInserted() ){ |
if ( !IsRunAlreadyInserted() ){ |
4092 |
glrun->SetID(this->AssignRunID()); |
glrun->SetID(this->AssignRunID()); |
4101 |
// |
// |
4102 |
}; |
}; |
4103 |
// |
// |
4104 |
|
if ( i > firstev ){ |
4105 |
|
if ( rhchko != rhchk ){ |
4106 |
|
if ( IsDebug() ) printf("oh oh... we have a runheader! stop here and handle later the remaining piece\n"); |
4107 |
|
lastev = i; |
4108 |
|
return(false); |
4109 |
|
}; |
4110 |
|
}; |
4111 |
|
// |
4112 |
if ( check == checklast && i != lastev ){ |
if ( check == checklast && i != lastev ){ |
4113 |
lastevtemp = i - 1; |
lastevtemp = i - 1; |
4114 |
i = lastev - 1; |
i = lastev - 1; |
4141 |
EventCounter *code=0; |
EventCounter *code=0; |
4142 |
UInt_t firstev = 0; |
UInt_t firstev = 0; |
4143 |
UInt_t lastev = 0; |
UInt_t lastev = 0; |
4144 |
UInt_t nevent = 0; |
// UInt_t nevent = 0; |
4145 |
UInt_t checkfirst = 0; |
UInt_t checkfirst = 0; |
4146 |
UInt_t checklast = 0; |
UInt_t checklast = 0; |
4147 |
UInt_t firstentry = 0; |
UInt_t firstentry = 0; |
4166 |
EventHeader *eh = 0; |
EventHeader *eh = 0; |
4167 |
PscuHeader *ph = 0; |
PscuHeader *ph = 0; |
4168 |
T->SetBranchAddress("Header", &eh); |
T->SetBranchAddress("Header", &eh); |
4169 |
nevent = T->GetEntries(); |
// nevent = T->GetEntries(); |
4170 |
// |
// |
4171 |
codt = eht->GetCounter(); |
codt = eht->GetCounter(); |
4172 |
codh = ehh->GetCounter(); |
codh = ehh->GetCounter(); |
4400 |
UInt_t totime = 0; |
UInt_t totime = 0; |
4401 |
UInt_t obt = 0; |
UInt_t obt = 0; |
4402 |
UInt_t pkt = 0; |
UInt_t pkt = 0; |
4403 |
|
Float_t totped = 0.; |
4404 |
// |
// |
4405 |
tr = (TTree*)file->Get("CalibCalPed"); |
tr = (TTree*)file->Get("CalibCalPed"); |
4406 |
if ( !tr || tr->IsZombie() ) throw -21; |
if ( !tr || tr->IsZombie() ) throw -21; |
4418 |
if ( calibCalPed->cstwerr[section] ){ |
if ( calibCalPed->cstwerr[section] ){ |
4419 |
valid = 1; |
valid = 1; |
4420 |
if ( calibCalPed->cperror[section] ) valid = 0; |
if ( calibCalPed->cperror[section] ) valid = 0; |
4421 |
|
// |
4422 |
|
// check pedestal values for one plane, if all zeros calibration is not valid (calorimeter power problems) [8th data reduction bug, fixed on 25/11/2009 by E.M.] |
4423 |
|
// |
4424 |
|
totped = 0.; |
4425 |
|
Int_t ns = 0; |
4426 |
|
if ( section == 2 ) ns = 3; |
4427 |
|
if ( section == 3 ) ns = 1; |
4428 |
|
if ( section == 1 ) ns = 2; |
4429 |
|
for (UInt_t ss=0; ss<96; ss++){ |
4430 |
|
totped += fabs(calibCalPed->calped[ns][0][ss]); |
4431 |
|
} |
4432 |
|
if ( totped < 1. ){ |
4433 |
|
if ( IsDebug() ) printf(" Section %i totped %f - No calibration data! Calorimeter power problems? \n",section,totped); |
4434 |
|
valid = 0; |
4435 |
|
}; |
4436 |
|
// |
4437 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
4438 |
obt = ph->GetOrbitalTime(); |
obt = ph->GetOrbitalTime(); |
4439 |
pkt = ph->GetCounter(); |
pkt = ph->GetCounter(); |
4601 |
tr->SetBranchAddress("Header", &eh); |
tr->SetBranchAddress("Header", &eh); |
4602 |
nevents = tr->GetEntries(); |
nevents = tr->GetEntries(); |
4603 |
// |
// |
4604 |
|
Float_t totpul = 0.; |
4605 |
|
// |
4606 |
if ( nevents > 0 ){ |
if ( nevents > 0 ){ |
4607 |
// |
// |
4608 |
for (UInt_t i=0; i < nevents; i++){ |
for (UInt_t i=0; i < nevents; i++){ |
4612 |
if ( cp1->pstwerr[section] && cp1->unpackError == 0 ){ |
if ( cp1->pstwerr[section] && cp1->unpackError == 0 ){ |
4613 |
valid = 1; |
valid = 1; |
4614 |
if ( cp1->pperror[section] ) valid = 0; |
if ( cp1->pperror[section] ) valid = 0; |
4615 |
|
// |
4616 |
|
// check pulse values for one plane, if all zeros calibration is not valid (calorimeter power problems) [8th data reduction bug, fixed on 25/11/2009 by E.M.] |
4617 |
|
// |
4618 |
|
totpul = 0.; |
4619 |
|
Int_t ns = 0; |
4620 |
|
if ( section == 2 ) ns = 3; |
4621 |
|
if ( section == 3 ) ns = 1; |
4622 |
|
if ( section == 1 ) ns = 2; |
4623 |
|
for (UInt_t ss=0; ss<96; ss++){ |
4624 |
|
totpul += cp1->calpuls[ns][0][ss]; |
4625 |
|
} |
4626 |
|
if ( totpul >= 3145632. ){ |
4627 |
|
if ( IsDebug() ) printf(" PULSE1 Section %i totpul %f - No calibration data! Calorimeter power problems? \n",section,totpul); |
4628 |
|
valid = 0; |
4629 |
|
}; |
4630 |
|
// |
4631 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
4632 |
obt = ph->GetOrbitalTime(); |
obt = ph->GetOrbitalTime(); |
4633 |
pkt = ph->GetCounter(); |
pkt = ph->GetCounter(); |
4785 |
if ( cp2->pstwerr[section] && cp2->unpackError == 0 ){ |
if ( cp2->pstwerr[section] && cp2->unpackError == 0 ){ |
4786 |
valid = 1; |
valid = 1; |
4787 |
if ( cp2->pperror[section] ) valid = 0; |
if ( cp2->pperror[section] ) valid = 0; |
4788 |
|
// |
4789 |
|
// check pulse values for one plane, if all zeros calibration is not valid (calorimeter power problems) [8th data reduction bug, fixed on 25/11/2009 by E.M.] |
4790 |
|
// |
4791 |
|
totpul = 0.; |
4792 |
|
Int_t ns = 0; |
4793 |
|
if ( section == 2 ) ns = 3; |
4794 |
|
if ( section == 3 ) ns = 1; |
4795 |
|
if ( section == 1 ) ns = 2; |
4796 |
|
for (UInt_t ss=0; ss<96; ss++){ |
4797 |
|
totpul += cp2->calpuls[ns][0][ss]; |
4798 |
|
} |
4799 |
|
if ( totpul >= 3145632. ){ |
4800 |
|
if ( IsDebug() ) printf(" PULSE2 Section %i totpul %f - No calibration data! Calorimeter power problems? \n",section,totpul); |
4801 |
|
valid = 0; |
4802 |
|
}; |
4803 |
|
// |
4804 |
ph = eh->GetPscuHeader(); |
ph = eh->GetPscuHeader(); |
4805 |
obt = ph->GetOrbitalTime(); |
obt = ph->GetOrbitalTime(); |
4806 |
pkt = ph->GetCounter(); |
pkt = ph->GetCounter(); |
4976 |
return(0); |
return(0); |
4977 |
}; |
}; |
4978 |
|
|
|
// /** |
|
|
// * Fill the GL_TRK_CALIB table |
|
|
// */ |
|
|
// void PamelaDBOperations::HandleTRK_CALIB(Bool_t pk1, Bool_t pk2){ |
|
|
// // |
|
|
// TSQLResult *result = 0; |
|
|
// TSQLRow *row = 0; |
|
|
// // |
|
|
// stringstream oss; |
|
|
// oss.str(""); |
|
|
// // |
|
|
// UInt_t totime = 0; |
|
|
// // |
|
|
// if ( !pk1 && !pk2 ){ |
|
|
// if ( IsDebug() ) printf(" Cannot handle trk calibration with both packet missing!\n"); |
|
|
// return; |
|
|
// }; |
|
|
// // |
|
|
// // check if the calibration has already been inserted |
|
|
// // |
|
|
// oss.str(""); |
|
|
// oss << " SELECT ID FROM GL_TRK_CALIB WHERE " |
|
|
// << " BOOT_NUMBER = "<< this->GetBOOTnumber(); // |
|
|
// oss << " AND ( ( "; |
|
|
// if ( pk1 ){ |
|
|
// oss << " OBT1 = "<< obt1 << " AND " |
|
|
// << " PKT1 = "<< pkt1 |
|
|
// << " ) OR ( "; |
|
|
// } else { |
|
|
// oss << " PKT1 = "<< pkt2-1 |
|
|
// << " ) OR ( "; |
|
|
// }; |
|
|
// if ( pk2 ){ |
|
|
// oss << " OBT2 = "<< obt2 << " AND " |
|
|
// << " PKT2 = "<< pkt2; |
|
|
// } else { |
|
|
// oss << " PKT2 = "<< pkt1+1; |
|
|
// }; |
|
|
// oss << " ) );"; |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" Check if the trk calibration has already been inserted: query is \n %s \n",oss.str().c_str()); |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// // |
|
|
// if ( !result ) throw -4; |
|
|
// // |
|
|
// row = result->Next(); |
|
|
// // |
|
|
// if ( row ){ |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" Trk calibration already inserted in the DB\n"); |
|
|
// if ( PEDANTIC ) throw -80; |
|
|
// // |
|
|
// } else { |
|
|
// // |
|
|
// // we have to insert a new calibration, check where to place it |
|
|
// // |
|
|
// oss.str(""); |
|
|
// oss << " SELECT ID,TO_TIME FROM GL_TRK_CALIB WHERE " |
|
|
// << " FROM_TIME < "<< fromtime << " AND " |
|
|
// << " TO_TIME > "<< fromtime << ";"; |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" Check where to place the trk calibration: query is \n %s \n",oss.str().c_str()); |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// // |
|
|
// if ( !result ) throw -4; |
|
|
// // |
|
|
// row = result->Next(); |
|
|
// // |
|
|
// if ( !row ){ |
|
|
// // |
|
|
// // no calibrations in the db contain our calibration |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" Calibration with fromtime lower than others to be inserted in the DB\n"); |
|
|
// if ( fromtime < 1150871000 ) fromtime = 0; // the first flight calibration was taken at about 1150863300 s, this line allows to analyze first runs in raw mode |
|
|
// // |
|
|
// oss.str(""); |
|
|
// oss << " SELECT FROM_TIME FROM GL_TRK_CALIB WHERE " |
|
|
// << " FROM_TIME > "<< fromtime << " ORDER BY FROM_TIME ASC LIMIT 1;"; |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" Check the upper limit for calibration: query is \n %s \n",oss.str().c_str()); |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// // |
|
|
// if ( !result ) throw -4; |
|
|
// // |
|
|
// row = result->Next(); |
|
|
// if ( !row ){ |
|
|
// totime = numeric_limits<UInt_t>::max(); |
|
|
// } else { |
|
|
// totime = (UInt_t)atoll(row->GetField(0)); |
|
|
// }; |
|
|
// // |
|
|
// } else { |
|
|
// // |
|
|
// // determine upper and lower limits and make space for the new calibration |
|
|
// // |
|
|
// totime = (UInt_t)atoll(row->GetField(1)); |
|
|
// // |
|
|
// oss.str(""); |
|
|
// oss << " UPDATE GL_TRK_CALIB SET " |
|
|
// << " TO_TIME = "<< fromtime << " WHERE " // NOTICE: to_time is equal to from_time of the calibration before, so the interval is: [from_time,to_time[ |
|
|
// << " ID = "<< row->GetField(0) << ";"; |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" Make space for the new trk calibration: query is \n %s \n",oss.str().c_str()); |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// // |
|
|
// if ( !result ) throw -4; |
|
|
// // |
|
|
// }; |
|
|
// // |
|
|
// oss.str(""); |
|
|
// oss << " INSERT INTO GL_TRK_CALIB (ID,ID_ROOT_L0,EV_ROOT_CALIBTRK1,EV_ROOT_CALIBTRK2,FROM_TIME,TO_TIME,OBT1,PKT1,OBT2,PKT2,BOOT_NUMBER,VALIDATION) " |
|
|
// << " VALUES (NULL,' " |
|
|
// << idroot << "',"; |
|
|
// // |
|
|
// if ( !pk1 ){ |
|
|
// oss << "NULL,"; |
|
|
// } else { |
|
|
// oss << "'" |
|
|
// << t1 << "',"; |
|
|
// }; |
|
|
// // |
|
|
// if ( !pk2 ){ |
|
|
// oss << "NULL,'"; |
|
|
// } else { |
|
|
// oss << "'" |
|
|
// << t2 << "','"; |
|
|
// }; |
|
|
// // |
|
|
// oss << fromtime << "','" |
|
|
// << totime << "','" |
|
|
// << obt1 << "','" |
|
|
// << pkt1 << "','" |
|
|
// << obt2 << "','" |
|
|
// << pkt2 << "','" |
|
|
// << this->GetBOOTnumber() << "','" |
|
|
// << valid << "');"; |
|
|
// // |
|
|
// if ( IsDebug() ) printf(" Insert the new trk calibration: query is \n %s \n",oss.str().c_str()); |
|
|
// // |
|
|
// result = conn->Query(oss.str().c_str()); |
|
|
// // |
|
|
// if ( !result ) throw -4; |
|
|
// // |
|
|
// }; |
|
|
// // |
|
|
// }; |
|
4979 |
/** |
/** |
4980 |
* Fill the GL_TRK_CALIB table |
* Fill the GL_TRK_CALIB table |
4981 |
*/ |
*/ |
5166 |
if ( !result ) throw -4;; |
if ( !result ) throw -4;; |
5167 |
row = result->Next(); |
row = result->Next(); |
5168 |
if(row)glcal->ID = (UInt_t)atoll(row->GetField(0)); |
if(row)glcal->ID = (UInt_t)atoll(row->GetField(0)); |
5169 |
|
// |
5170 |
|
delete result; // mmm... "Error in <TMySQLResult::Next>: result set closed" solved? 090112 [8RED error messages] |
5171 |
// |
// |
5172 |
}; |
}; |
5173 |
|
|
5377 |
// |
// |
5378 |
}; |
}; |
5379 |
// |
// |
5380 |
|
|
5381 |
|
if( !(MISSING_pkt1&MISSING_pkt2) ){ |
5382 |
|
this->HandleTRK_CALIB(!MISSING_pkt1,!MISSING_pkt2); |
5383 |
|
ncalib++; |
5384 |
|
if( MISSING_pkt1||MISSING_pkt2||!valid )try_to_recover=true; |
5385 |
|
} |
5386 |
|
|
5387 |
|
|
5388 |
} else { |
} else { |
5389 |
// |
// |
5390 |
if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u OUTSIDE the considered time interval \n",fromtime,obt1,pkt1); |
if ( IsDebug() ) printf(" Trk calibration1 at time %u obt %u pkt %u OUTSIDE the considered time interval \n",fromtime,obt1,pkt1); |
5391 |
// if ( PEDANTIC ) throw -79; |
// if ( PEDANTIC ) throw -79; |
5392 |
// |
// |
5393 |
}; |
}; |
5394 |
// |
// |
5395 |
|
|
|
if( !(MISSING_pkt1&MISSING_pkt2) ){ |
|
|
this->HandleTRK_CALIB(!MISSING_pkt1,!MISSING_pkt2); |
|
|
ncalib++; |
|
|
if( MISSING_pkt1||MISSING_pkt2||!valid )try_to_recover=true; |
|
|
} |
|
|
|
|
5396 |
}; //end loop on pkt1 |
}; //end loop on pkt1 |
5397 |
|
|
5398 |
|
|
5449 |
|
|
5450 |
if ( IsDebug() ) cout << endl << ">>>> TRY TO RECOVER TRACKER CALIBRATIONS <<<<"<<endl; |
if ( IsDebug() ) cout << endl << ">>>> TRY TO RECOVER TRACKER CALIBRATIONS <<<<"<<endl; |
5451 |
|
|
5452 |
TSQLResult *result = 0; |
TSQLResult *rresult = 0; // 10RED bug, result is re-used inside result looop... |
5453 |
|
TSQLResult *result = 0; // 10RED bug, result is re-used inside result looop... |
5454 |
TSQLRow *row = 0; |
TSQLRow *row = 0; |
5455 |
// |
// |
5456 |
stringstream oss; |
stringstream oss; |
5480 |
oss << " ORDER BY ID DESC LIMIT "<<ncalib<<"; "; |
oss << " ORDER BY ID DESC LIMIT "<<ncalib<<"; "; |
5481 |
if ( IsDebug() ) cout << oss.str().c_str() << endl; |
if ( IsDebug() ) cout << oss.str().c_str() << endl; |
5482 |
|
|
5483 |
result = conn->Query(oss.str().c_str()); |
rresult = conn->Query(oss.str().c_str()); |
5484 |
if ( !result ) throw -4;; |
if ( !result ) throw -4;; |
5485 |
if ( IsDebug() ) cout <<"Rows: "<<result->GetRowCount()<<endl; |
if ( IsDebug() ) cout <<"Rows: "<<result->GetRowCount()<<endl; |
5486 |
|
|
5489 |
// ----------------------------------- |
// ----------------------------------- |
5490 |
UInt_t nn=0; |
UInt_t nn=0; |
5491 |
do { |
do { |
5492 |
row = result->Next(); |
row = rresult->Next(); |
5493 |
if(!row)break; |
if(!row)break; |
5494 |
|
|
5495 |
UInt_t id = (UInt_t)atoll(row->GetField(0)); |
UInt_t id = (UInt_t)atoll(row->GetField(0)); |
5528 |
TString path = ""; |
TString path = ""; |
5529 |
TString name = ""; |
TString name = ""; |
5530 |
TString raw = ""; |
TString raw = ""; |
5531 |
UInt_t obt0 = 0; |
// UInt_t obt0 = 0; |
5532 |
UInt_t timesync = 0; |
// UInt_t timesync = 0; |
5533 |
UInt_t boot_number = 0; |
// UInt_t boot_number = 0; |
5534 |
bool FOUND = false; |
bool FOUND = false; |
5535 |
|
|
5536 |
if ( IsDebug() ) cout << "------------------------------------------------------------" <<endl; |
if ( IsDebug() ) cout << "------------------------------------------------------------" <<endl; |
5603 |
path = (TString)row2->GetField(0); |
path = (TString)row2->GetField(0); |
5604 |
name = (TString)row2->GetField(1); |
name = (TString)row2->GetField(1); |
5605 |
raw = (TString)row2->GetField(5); |
raw = (TString)row2->GetField(5); |
5606 |
obt0 = (UInt_t)atoll(row2->GetField(2)); |
// obt0 = (UInt_t)atoll(row2->GetField(2)); |
5607 |
timesync = (UInt_t)atoll(row2->GetField(3)); |
// timesync = (UInt_t)atoll(row2->GetField(3)); |
5608 |
boot_number = (UInt_t)atoll(row2->GetField(4)); |
// boot_number = (UInt_t)atoll(row2->GetField(4)); |
5609 |
|
|
5610 |
if ( IsDebug() ) cout << "- - - - - - - - - - -" <<endl; |
if ( IsDebug() ) cout << "- - - - - - - - - - -" <<endl; |
5611 |
// cout << path <<endl; |
// cout << path <<endl; |
5612 |
cout << "File : " <<name <<endl; |
// cout << "File : " <<name <<endl; |
5613 |
// cout << obt0 <<endl; |
// cout << obt0 <<endl; |
5614 |
// cout << timesync <<endl; |
// cout << timesync <<endl; |
5615 |
cout << "boot n. : "<<boot_number <<endl; |
// cout << "boot n. : "<<boot_number <<endl; |
5616 |
// cout << raw <<endl; |
// cout << raw <<endl; |
5617 |
|
|
5618 |
// ------------------------------------------------------ |
// ------------------------------------------------------ |
5619 |
// ... retrieve the calibration packets. |
// ... retrieve the calibration packets. |
5620 |
// ------------------------------------------------------ |
// ------------------------------------------------------ |
5621 |
TFile *file = new TFile(path+"/"+name); |
if ( IsDebug() ) printf(" file is %s/%s \n",((TString)gSystem->ExpandPathName(path.Data())).Data(),name.Data()); |
5622 |
|
TFile *file = new TFile(((TString)gSystem->ExpandPathName(path.Data()))+"/"+name); // EM, path could be symbolic and we must expand it |
5623 |
if(!file)throw -100; |
if(!file)throw -100; |
5624 |
if(file->IsZombie())throw -100; |
if(file->IsZombie())throw -100; |
5625 |
// |
// |
5669 |
if(FOUND)break; |
if(FOUND)break; |
5670 |
|
|
5671 |
}while(1);//endl loop over root table entries |
}while(1);//endl loop over root table entries |
5672 |
|
|
5673 |
|
delete result2; // 10RED potential memory leak |
5674 |
if(FOUND)break; |
if(FOUND)break; |
5675 |
|
|
5676 |
}//end loop over tables |
}//end loop over tables |
5677 |
|
|
5678 |
if(FOUND){ |
if(FOUND){ |
5679 |
|
|
5680 |
if ( IsDebug() ) cout << " >>> REPETITION FOUND :-) <<<" <<endl; |
if ( IsDebug() ) cout << " >>> REPETITION FOUND :-) <<<" <<endl; |
|
|
|
|
//////////////////////////////////////////// |
|
|
// insert a new entry in GL_TRK_CALIB and |
|
|
// modify the time-tag of the previous one |
|
|
//////////////////////////////////////////// |
|
5681 |
|
|
5682 |
// --------------------------------------------------------------------- |
//////////////////////////////////////////// |
5683 |
// step 1: insert a new raw file in GL_RAW |
// insert a new entry in GL_TRK_CALIB and |
5684 |
// --------------------------------------------------------------------- |
// modify the time-tag of the previous one |
5685 |
// |
//////////////////////////////////////////// |
5686 |
// check if the raw file already exist |
|
5687 |
// |
// --------------------------------------------------------------------- |
5688 |
UInt_t id_raw = 0; |
// step 1: insert a new raw file in GL_RAW |
5689 |
oss.str(""); |
// --------------------------------------------------------------------- |
5690 |
oss << "SELECT ID FROM GL_RAW where NAME=\"" << gSystem->BaseName(raw.Data()) <<"\";"; |
// |
5691 |
if ( IsDebug() ) cout << oss.str().c_str() << endl; |
// check if the raw file already exist |
5692 |
|
// |
5693 |
|
UInt_t id_raw = 0; // EM GL_RAW is there only for backward compatibility so we do not need to fill it when in "chewbacca" mode |
5694 |
|
// oss.str(""); |
5695 |
|
// oss << "SELECT ID FROM GL_RAW where NAME=\"" << gSystem->BaseName(raw.Data()) <<"\";"; |
5696 |
|
// if ( IsDebug() ) cout << oss.str().c_str() << endl; |
5697 |
|
|
5698 |
result = conn->Query(oss.str().c_str()); |
// result = conn->Query(oss.str().c_str()); |
5699 |
if ( !result ) throw -4;; |
// if ( !result ) throw -4;; |
5700 |
if ( IsDebug() ) cout <<"Rows: "<<result->GetRowCount()<<endl; |
// if ( IsDebug() ) cout <<"Rows: "<<result->GetRowCount()<<endl; |
5701 |
if( result->GetRowCount() == 0){ |
// if( result->GetRowCount() == 0){ |
5702 |
if ( IsDebug() ) cout << " << Insert new RAW file >> "<<endl; |
// if ( IsDebug() ) cout << " << Insert new RAW file >> "<<endl; |
5703 |
// - - - - - - - - - - - |
// // - - - - - - - - - - - |
5704 |
// insert new raw file |
// // insert new raw file |
5705 |
// - - - - - - - - - - - |
// // - - - - - - - - - - - |
5706 |
GL_RAW glraw = GL_RAW(); |
// GL_RAW glraw = GL_RAW(); |
5707 |
glraw.PATH = gSystem->DirName(raw.Data()); |
// glraw.PATH = gSystem->DirName(raw.Data()); |
5708 |
glraw.NAME = gSystem->BaseName(raw.Data()); |
// glraw.NAME = gSystem->BaseName(raw.Data()); |
5709 |
glraw.BOOT_NUMBER = boot_number; |
// glraw.BOOT_NUMBER = boot_number; |
5710 |
// |
// // |
5711 |
insertPamelaRawFile( &glraw ); |
// insertPamelaRawFile( &glraw ); |
5712 |
// |
// // |
5713 |
id_raw = glraw.ID; |
// id_raw = glraw.ID; |
5714 |
}else{ |
// }else{ |
5715 |
row = result->Next(); |
// row = result->Next(); |
5716 |
id_raw = (UInt_t)atoll(row->GetField(0)); |
// id_raw = (UInt_t)atoll(row->GetField(0)); |
5717 |
} |
// } |
5718 |
if ( IsDebug() ) cout << "ID_RAW = "<<id_raw<<endl; |
// if ( IsDebug() ) cout << "ID_RAW = "<<id_raw<<endl; |
5719 |
|
|
5720 |
// --------------------------------------------------------------------- |
// --------------------------------------------------------------------- |
5721 |
// step 1(bis): retrieve the timesync id associated to the file |
// step 1(bis): retrieve the timesync id associated to the file |
5722 |
// (NB, uso lo stesso associato al file iniziale) |
// (NB, uso lo stesso associato al file iniziale) |
5723 |
// --------------------------------------------------------------------- |
// --------------------------------------------------------------------- |
5724 |
UInt_t idtimesync = 0; |
UInt_t idtimesync = 0; |
5725 |
oss.str(""); |
oss.str(""); |
5726 |
oss << " SELECT ID FROM GL_TIMESYNC where TIMESYNC="<<chlastts<<" AND OBT0="<<chobtts*1000<<" limit 1;"; |
oss << " SELECT ID FROM GL_TIMESYNC where TIMESYNC="<<chlastts<<" AND OBT0="<<chobtts*1000<<" limit 1;"; |
5727 |
if ( debug ) printf(" %s \n",oss.str().c_str()); |
if ( debug ) printf(" %s \n",oss.str().c_str()); |
5728 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
5729 |
if ( !result ) throw -3; |
if ( !result ) throw -3; |
5730 |
row = result->Next(); |
row = result->Next(); |
5731 |
if ( !row ) throw -3; |
if ( !row ) throw -3; |
5732 |
idtimesync = (UInt_t)atoll(row->GetField(0)); |
idtimesync = (UInt_t)atoll(row->GetField(0)); |
5733 |
if ( IsDebug() ) cout << "ID_TIMESYNC = "<<idtimesync<<endl; |
if ( IsDebug() ) cout << "ID_TIMESYNC = "<<idtimesync<<endl; |
5734 |
|
|
5735 |
|
delete result;// mmm... "Error in <TMySQLResult::Next>: result set closed" solved? 090112 [8RED error messages] |
5736 |
// --------------------------------------------------------------------- |
// --------------------------------------------------------------------- |
5737 |
// step 2: insert a new root file in GL_ROOT |
// step 2: insert a new root file in GL_ROOT |
5738 |
// --------------------------------------------------------------------- |
// --------------------------------------------------------------------- |
5739 |
// |
// |
5740 |
// check if the raw file already exist |
// check if the root file already exist |
5741 |
// |
// |
5742 |
UInt_t id_root = 0; |
UInt_t id_root = 0; |
5743 |
oss.str(""); |
oss.str(""); |
5755 |
GL_ROOT glroot = GL_ROOT(); |
GL_ROOT glroot = GL_ROOT(); |
5756 |
glroot.ID_RAW = id_raw; |
glroot.ID_RAW = id_raw; |
5757 |
glroot.ID_TIMESYNC = idtimesync; |
glroot.ID_TIMESYNC = idtimesync; |
5758 |
glroot.PATH = path; |
// |
5759 |
|
// EM STATIC = the full expanded path must be put in the DB, KEEPENV = the path given as input (or found in ROOT_TABLE) must be used, |
5760 |
|
// NOT STATIC NOT KEEPENV = $PAM_L0 must be used in the DB |
5761 |
|
// |
5762 |
|
if ( STATIC ){ |
5763 |
|
glroot.PATH = (TString)gSystem->ExpandPathName(path); |
5764 |
|
} else { |
5765 |
|
if ( KEEPENV ){ |
5766 |
|
glroot.PATH = path; |
5767 |
|
} else { |
5768 |
|
glroot.PATH = "$PAM_L0"; |
5769 |
|
}; |
5770 |
|
}; |
5771 |
|
// glroot.PATH = path; |
5772 |
glroot.NAME = name; |
glroot.NAME = name; |
5773 |
// |
// |
5774 |
insertPamelaRootFile( &glroot ); |
insertPamelaRootFile( &glroot ); |
5780 |
} |
} |
5781 |
if ( IsDebug() ) cout << "ID_ROOT = "<<id_root<<endl; |
if ( IsDebug() ) cout << "ID_ROOT = "<<id_root<<endl; |
5782 |
|
|
5783 |
|
delete result;// mmm... "Error in <TMySQLResult::Next>: result set closed" solved? 090112 [8RED error messages] |
5784 |
// --------------------------------------------------------------------- |
// --------------------------------------------------------------------- |
5785 |
// step 3: modify time-tag of corrupted GL_TRK_CALIB entry |
// step 3: modify time-tag of corrupted GL_TRK_CALIB entry |
5786 |
// --------------------------------------------------------------------- |
// --------------------------------------------------------------------- |
5793 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
5794 |
if ( !result ) throw -4;; |
if ( !result ) throw -4;; |
5795 |
|
|
5796 |
|
delete result; // mmm... "Error in <TMySQLResult::Next>: result set closed" solved? 090112 [8RED error messages] |
5797 |
|
|
5798 |
// --------------------------------------------------------------------- |
// --------------------------------------------------------------------- |
5799 |
// step 4: insert the new calibration: |
// step 4: insert the new calibration: |
5800 |
// --------------------------------------------------------------------- |
// --------------------------------------------------------------------- |
5824 |
}while(1);//end loop over calibrations |
}while(1);//end loop over calibrations |
5825 |
|
|
5826 |
|
|
5827 |
if( result )delete result; |
if( rresult )delete rresult; // 10RED bug |
5828 |
if( row )delete row; |
if( row )delete row; |
5829 |
|
|
5830 |
|
|
6125 |
memset(myid,0,500*sizeof(UInt_t)); |
memset(myid,0,500*sizeof(UInt_t)); |
6126 |
// |
// |
6127 |
oss.str(""); |
oss.str(""); |
6128 |
oss << "SELECT ID,RUNTRAILER_TIME,RUNTRAILER_PKT,BOOT_NUMBER FROM GL_RUN_FRAGMENTS WHERE INSERT_TIME <= '" << clean_time->AsSQLString() << "'order BY RUNHEADER_TIME asc;"; |
oss << "SELECT ID,RUNTRAILER_TIME,RUNTRAILER_PKT,BOOT_NUMBER FROM GL_RUN_FRAGMENTS WHERE INSERT_TIME <= '" << clean_time->AsSQLString() << "' order BY RUNHEADER_TIME asc;"; |
6129 |
if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS runs older than %s : query is \n %s \n",clean_time->AsSQLString(),oss.str().c_str()); |
if ( IsDebug() ) printf(" Select from GL_RUN_FRAGMENTS runs older than %s : query is \n %s \n",clean_time->AsSQLString(),oss.str().c_str()); |
6130 |
nresult = conn->Query(oss.str().c_str()); |
nresult = conn->Query(oss.str().c_str()); |
6131 |
// |
// |
6246 |
oss << " PHYSENDRUN_MASK_S11CRC=" << glrun->GetPHYSENDRUN_MASK_S11CRC()<< " , "; |
oss << " PHYSENDRUN_MASK_S11CRC=" << glrun->GetPHYSENDRUN_MASK_S11CRC()<< " , "; |
6247 |
// |
// |
6248 |
oss << " RUNHEADER_OBT=" << glrun->GetRUNHEADER_OBT()<< " , " |
oss << " RUNHEADER_OBT=" << glrun->GetRUNHEADER_OBT()<< " , " |
6249 |
<< " RUNHEADER_PKT=" << glrun->GetRUNHEADER_PKT()<< ";"; |
<< " RUNHEADER_PKT=" << glrun->GetRUNHEADER_PKT(); |
6250 |
|
oss << " WHERE ID=" << myid[u] << ";"; // BUG!!! 090112 |
6251 |
conn->Query(oss.str().c_str()); |
conn->Query(oss.str().c_str()); |
6252 |
}; |
}; |
6253 |
}; |
}; |
6296 |
if ( glrun->GetPHYSENDRUN_MASK_S11CRC() ){ |
if ( glrun->GetPHYSENDRUN_MASK_S11CRC() ){ |
6297 |
oss << " PHYSENDRUN_MASK_S11CRC=" << glrun->GetPHYSENDRUN_MASK_S11CRC()<< " , "; }; |
oss << " PHYSENDRUN_MASK_S11CRC=" << glrun->GetPHYSENDRUN_MASK_S11CRC()<< " , "; }; |
6298 |
// |
// |
6299 |
oss << " RUNTRAILER_PKT=" << glrun->GetRUNTRAILER_PKT()<< ";"; |
oss << " RUNTRAILER_PKT=" << glrun->GetRUNTRAILER_PKT(); |
6300 |
|
oss << " WHERE ID=" << myid[u] << ";"; // BUG!!! 090112 |
6301 |
conn->Query(oss.str().c_str()); |
conn->Query(oss.str().c_str()); |
6302 |
}; |
}; |
6303 |
}; |
}; |
6348 |
row2 = result2->Next(); |
row2 = result2->Next(); |
6349 |
// |
// |
6350 |
if ( row2 ){ |
if ( row2 ){ |
6351 |
if ( IsDebug() ) printf(" The already exist in the GL_RUN table! \n"); |
if ( IsDebug() ) printf(" The run already exists in the GL_RUN table! \n"); |
6352 |
if ( PEDANTIC ) throw -83; |
if ( PEDANTIC ) throw -83; |
6353 |
row2->Close(); |
row2->Close(); |
6354 |
}; |
}; |
6425 |
glrun->Query_GL_RUN_FRAGMENTS(oss.str().c_str(),conn); |
glrun->Query_GL_RUN_FRAGMENTS(oss.str().c_str(),conn); |
6426 |
// |
// |
6427 |
oss.str(""); |
oss.str(""); |
6428 |
oss << " SELECT ID,NEVENTS,TRK_CALIB_USED,PKT_COUNTER FROM GL_RUN WHERE " |
oss << " SELECT ID,NEVENTS,TRK_CALIB_USED,PKT_COUNTER,ID_ROOT_L0 FROM GL_RUN WHERE " |
6429 |
<< " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND (" |
<< " BOOT_NUMBER=" << glrun->GetBOOT_NUMBER() << " AND (" |
6430 |
<< " (RUNHEADER_TIME>=" << (UInt_t)(glrun->GetRUNHEADER_TIME()-10) << " AND " |
<< " (RUNHEADER_TIME>=" << (UInt_t)(glrun->GetRUNHEADER_TIME()-10) << " AND " |
6431 |
<< " RUNTRAILER_TIME<=" << (UInt_t)(glrun->GetRUNTRAILER_TIME()+10) << " AND (" |
<< " RUNTRAILER_TIME<=" << (UInt_t)(glrun->GetRUNTRAILER_TIME()+10) << " AND (" |
6483 |
moved++; |
moved++; |
6484 |
// |
// |
6485 |
} else { |
} else { |
6486 |
if ( IsDebug() ) printf(" The already exist in the GL_RUN table! \n"); |
|
6487 |
if ( PEDANTIC ) throw -83; |
// questa parte mah mah mah... da controllare |
6488 |
|
|
6489 |
|
// |
6490 |
|
// we end up here if chewbacca and we have RH---| small gap |---RT, in this case the two pieces are attached and moved to GL_RUN. We are now processing the |
6491 |
|
// "small gap" piece... we recognize this since: we have two entries from this query, the pkt number is consistent with our pkt number. |
6492 |
|
// |
6493 |
|
// 090112 [8RED (-70): RUN ALREADY INSERTED] |
6494 |
|
// |
6495 |
|
Bool_t OK = false; |
6496 |
|
UInt_t IDRL2A = 0; |
6497 |
|
UInt_t IDRL2B = 0; |
6498 |
|
if ( result2->GetRowCount() == 2 ){ |
6499 |
|
IDRL2A = (UInt_t)atoll(row2->GetField(0)); |
6500 |
|
UInt_t IDRL0A = (UInt_t)atoll(row2->GetField(4)); |
6501 |
|
row2 = result2->Next(); |
6502 |
|
IDRL2B = (UInt_t)atoll(row2->GetField(0)); |
6503 |
|
UInt_t IDRL0B = (UInt_t)atoll(row2->GetField(4)); |
6504 |
|
if ( IsDebug() ) printf(" IDRL0A %u B %u IDRL2A %u B %u \n",IDRL0A,IDRL0B,IDRL2A,IDRL2B); |
6505 |
|
// if ( IDRL0A == IDRL0B ){ |
6506 |
|
// TSQLResult *result2a = 0; |
6507 |
|
// TSQLRow *row2a = 0; |
6508 |
|
// oss.str(""); |
6509 |
|
// oss << "select PKT_NUMBER_FINAL from ROOT_TABLE_MERGING where REAL_TIME_INIT<" << chrtinit << " order by REAL_TIME_INIT desc limit 1;"; |
6510 |
|
// if ( IsDebug() ) printf(" Check if we are in the case RH---| small gap |---RT: query is \n %s \n",oss.str().c_str()); |
6511 |
|
// result2a = conn->Query(oss.str().c_str()); |
6512 |
|
// // |
6513 |
|
// if ( !result2a ) throw -4; |
6514 |
|
// // |
6515 |
|
// row2a = result2a->Next(); |
6516 |
|
// // |
6517 |
|
// if ( row2a ){ |
6518 |
|
// UInt_t PKA = (UInt_t)atoll(row2a->GetField(0)); |
6519 |
|
// delete result2a; |
6520 |
|
// oss.str(""); |
6521 |
|
// oss << "select PKT_NUMBER_INIT from ROOT_TABLE_MERGING where REAL_TIME_INIT>" << chrtinit << " order by REAL_TIME_INIT asc limit 1;"; |
6522 |
|
// if ( IsDebug() ) printf(" Check if we are in the case RH---| small gap |---RT: query is \n %s \n",oss.str().c_str()); |
6523 |
|
// result2a = conn->Query(oss.str().c_str()); |
6524 |
|
// // |
6525 |
|
// if ( !result2a ) throw -4; |
6526 |
|
// // |
6527 |
|
// row2a = result2a->Next(); |
6528 |
|
// // |
6529 |
|
// if ( row2a ){ |
6530 |
|
// UInt_t PKB = (UInt_t)atoll(row2a->GetField(0)); |
6531 |
|
// // |
6532 |
|
// if ( IsDebug() ) printf(" PKT(PKA) + 1 %llu == runheaderpkt %llu && PKB = runtrailer %llu + 1 %llu \n",PKT(PKA)+1LL, PKT(glrun->GetRUNHEADER_PKT()),PKT(PKB), PKT(glrun->GetRUNTRAILER_PKT())+1LL); |
6533 |
|
// if ( PKT(PKA)+1LL == PKT(glrun->GetRUNHEADER_PKT()) && PKT(PKB) == PKT(glrun->GetRUNTRAILER_PKT())+1LL ){ |
6534 |
|
// if ( IsDebug() ) printf(" Ok, we are in the case: RH---| small gap |---RT \n"); |
6535 |
|
OK = true; |
6536 |
|
// }; |
6537 |
|
// }; |
6538 |
|
// }; |
6539 |
|
// }; |
6540 |
|
}; |
6541 |
|
if ( OK ){ |
6542 |
|
// |
6543 |
|
// this is the case in which we must insert a piece of run between two fragments in the GL_RUN table |
6544 |
|
// |
6545 |
|
// we have to update with runheader/trailer infos our run and we have to change the ROOT_ID_FRAG of the header piece ( 1/2 2/1 must become 1/3 3/2 2/1 ) |
6546 |
|
// |
6547 |
|
GL_RUN *glA = new GL_RUN(); |
6548 |
|
glA->Query_GL_RUN(IDRL2A,conn); |
6549 |
|
// |
6550 |
|
if ( glA->GetACQ_BUILD_INFO() != 0 ){ |
6551 |
|
// |
6552 |
|
// the first piece contains a good runheader we can update all the other runs with correct infos! |
6553 |
|
// |
6554 |
|
oss.str(""); |
6555 |
|
oss << "UPDATE GL_RUN_FRAGMENTS SET " |
6556 |
|
<< " RUNHEADER_TIME=" << glA->GetRUNHEADER_TIME()<< " , " |
6557 |
|
<< " RUNHEADER_OBT=" << glA->GetRUNHEADER_OBT()<< " , " |
6558 |
|
<< " RUNHEADER_PKT=" << glA->GetRUNHEADER_PKT()<< " , " |
6559 |
|
<< " COMPILATIONTIMESTAMP=" << glA->GetCOMPILATIONTIMESTAMP()<< " , " |
6560 |
|
<< " FAV_WRK_SCHEDULE=" << glA->GetFAV_WRK_SCHEDULE()<< " , " |
6561 |
|
<< " EFF_WRK_SCHEDULE=" << glA->GetEFF_WRK_SCHEDULE()<< " , " |
6562 |
|
<< " PRH_VAR_TRG_MODE_A=" << glA->GetPRH_VAR_TRG_MODE_A()<< " , " |
6563 |
|
<< " PRH_VAR_TRG_MODE_B=" << glA->GetPRH_VAR_TRG_MODE_B()<< " , " |
6564 |
|
<< " ACQ_BUILD_INFO=" << glA->GetACQ_BUILD_INFO()<< " , " |
6565 |
|
<< " ACQ_VAR_INFO=" << glA->GetACQ_VAR_INFO()<< " , " |
6566 |
|
<< " RM_ACQ_AFTER_CALIB=" << glA->GetRM_ACQ_AFTER_CALIB()<< " , " |
6567 |
|
<< " RM_ACQ_SETTING_MODE=" << glA->GetRM_ACQ_SETTING_MODE()<< " , " |
6568 |
|
<< " TRK_CALIB_USED=" << glA->GetTRK_CALIB_USED()<< " , " |
6569 |
|
<< " CAL_DSP_MASK=" << glA->GetCAL_DSP_MASK()<< " , " |
6570 |
|
<< " LAST_TIMESYNC=" << glA->GetLAST_TIMESYNC()<< " , "; |
6571 |
|
// |
6572 |
|
if ( glA->GetPHYSENDRUN_MASK_S3S2S12() ) |
6573 |
|
oss << " PHYSENDRUN_MASK_S3S2S12=" << glA->GetPHYSENDRUN_MASK_S3S2S12() << " , "; |
6574 |
|
if ( glA->GetPHYSENDRUN_MASK_S11CRC() ) |
6575 |
|
oss << " PHYSENDRUN_MASK_S11CRC=" << glA->GetPHYSENDRUN_MASK_S11CRC() << " , "; |
6576 |
|
// |
6577 |
|
oss << " OBT_TIMESYNC=" << glA->GetOBT_TIMESYNC(); |
6578 |
|
oss << " WHERE ID=" << glrun->GetID() << ";"; |
6579 |
|
if ( IsDebug() ) printf(" update with correct infos: %s\n",oss.str().c_str()); |
6580 |
|
conn->Query(oss.str().c_str()); |
6581 |
|
// |
6582 |
|
} else { |
6583 |
|
// |
6584 |
|
// sig no runheader, let set anyway what is possible... |
6585 |
|
// |
6586 |
|
oss.str(""); |
6587 |
|
oss << "UPDATE GL_RUN_FRAGMENTS SET " |
6588 |
|
<< " RUNHEADER_TIME=" << glA->GetRUNHEADER_TIME()<< " , "; |
6589 |
|
// |
6590 |
|
if ( glA->GetPHYSENDRUN_MASK_S3S2S12() ) |
6591 |
|
oss << " PHYSENDRUN_MASK_S3S2S12=" << glA->GetPHYSENDRUN_MASK_S3S2S12()<< " , "; |
6592 |
|
if ( glA->GetPHYSENDRUN_MASK_S11CRC() ) |
6593 |
|
oss << " PHYSENDRUN_MASK_S11CRC=" << glA->GetPHYSENDRUN_MASK_S11CRC()<< " , "; |
6594 |
|
// |
6595 |
|
oss << " RUNHEADER_OBT=" << glA->GetRUNHEADER_OBT()<< " , " |
6596 |
|
<< " RUNHEADER_PKT=" << glA->GetRUNHEADER_PKT(); |
6597 |
|
oss << " WHERE ID=" << glrun->GetID() << ";"; |
6598 |
|
if ( IsDebug() ) printf(" update with correct infos2: %s\n",oss.str().c_str()); |
6599 |
|
conn->Query(oss.str().c_str()); |
6600 |
|
}; |
6601 |
|
// |
6602 |
|
// update runheader ROOT_ID_FRAG |
6603 |
|
// |
6604 |
|
oss.str(""); |
6605 |
|
oss << "UPDATE GL_RUN SET ID_RUN_FRAG = " << glrun->GetID() << " where ID = " << IDRL2A << ";"; |
6606 |
|
if ( IsDebug() ) printf(" update gl_run to have cross indexing: %s\n",oss.str().c_str()); |
6607 |
|
conn->Query(oss.str().c_str()); |
6608 |
|
// |
6609 |
|
// now let's look for runtrailer if any in the last run |
6610 |
|
// |
6611 |
|
glA->Query_GL_RUN(IDRL2B,conn); |
6612 |
|
// |
6613 |
|
if ( glA->GetPKT_READY_COUNTER() != 0 ){ |
6614 |
|
// |
6615 |
|
// the first piece contains a good runtrailer we can update all the other runs with correct infos! |
6616 |
|
// |
6617 |
|
oss.str(""); |
6618 |
|
oss << "UPDATE GL_RUN_FRAGMENTS SET " |
6619 |
|
<< " RUNTRAILER_TIME=" << glA->GetRUNTRAILER_TIME()<< " , " |
6620 |
|
<< " RUNTRAILER_OBT=" << glA->GetRUNTRAILER_OBT()<< " , " |
6621 |
|
<< " RUNTRAILER_PKT=" << glA->GetRUNTRAILER_PKT()<< " , " |
6622 |
|
<< " PKT_COUNTER=" << glA->GetPKT_COUNTER()<< " , "; |
6623 |
|
// |
6624 |
|
if ( glA->GetPHYSENDRUN_MASK_S3S2S12() ){ |
6625 |
|
oss << " PHYSENDRUN_MASK_S3S2S12=" << glA->GetPHYSENDRUN_MASK_S3S2S12()<< " , "; }; |
6626 |
|
if ( glA->GetPHYSENDRUN_MASK_S11CRC() ) { |
6627 |
|
oss << " PHYSENDRUN_MASK_S11CRC=" << glA->GetPHYSENDRUN_MASK_S11CRC()<< " , "; }; |
6628 |
|
// |
6629 |
|
oss << " PKT_READY_COUNTER=" << glA->GetPKT_READY_COUNTER() |
6630 |
|
<< " WHERE ID=" << glrun->GetID() << ";"; |
6631 |
|
if ( IsDebug() ) printf(" update with correct trailer infos: %s\n",oss.str().c_str()); |
6632 |
|
conn->Query(oss.str().c_str()); |
6633 |
|
// |
6634 |
|
} else { |
6635 |
|
// |
6636 |
|
// sig no runtrailer, let set anyway what is possible... |
6637 |
|
// |
6638 |
|
oss.str(""); |
6639 |
|
oss << "UPDATE GL_RUN_FRAGMENTS SET " |
6640 |
|
<< " RUNTRAILER_TIME=" << glA->GetRUNTRAILER_TIME()<< " , " |
6641 |
|
<< " RUNTRAILER_OBT=" << glA->GetRUNTRAILER_OBT()<< " , "; |
6642 |
|
// |
6643 |
|
if ( glA->GetPHYSENDRUN_MASK_S3S2S12() ){ |
6644 |
|
oss << " PHYSENDRUN_MASK_S3S2S12=" << glA->GetPHYSENDRUN_MASK_S3S2S12()<< " , "; }; |
6645 |
|
if ( glA->GetPHYSENDRUN_MASK_S11CRC() ){ |
6646 |
|
oss << " PHYSENDRUN_MASK_S11CRC=" << glA->GetPHYSENDRUN_MASK_S11CRC()<< " , "; }; |
6647 |
|
// |
6648 |
|
oss << " RUNTRAILER_PKT=" << glrun->GetRUNTRAILER_PKT() |
6649 |
|
<< " WHERE ID=" << glrun->GetID() << ";"; |
6650 |
|
if ( IsDebug() ) printf(" update with correct trailer infos2: %s\n",oss.str().c_str()); |
6651 |
|
conn->Query(oss.str().c_str()); |
6652 |
|
}; |
6653 |
|
// |
6654 |
|
UInt_t myi = glrun->GetID(); |
6655 |
|
oss.str(""); |
6656 |
|
oss << " ID= "<< myi; |
6657 |
|
// |
6658 |
|
glrun->Query_GL_RUN_FRAGMENTS(oss.str().c_str(),conn); |
6659 |
|
// |
6660 |
|
// fill the new run in GL_RUN |
6661 |
|
// |
6662 |
|
glrun->SetID_RUN_FRAG(IDRL2B); |
6663 |
|
glrun->Fill_GL_RUN(conn); |
6664 |
|
glrun->DeleteRun(conn,0,"GL_RUN_FRAGMENTS"); |
6665 |
|
delete glA; |
6666 |
|
moved++; |
6667 |
|
// |
6668 |
|
} else { |
6669 |
|
// |
6670 |
|
// is just a repetition |
6671 |
|
// |
6672 |
|
if ( IsDebug() ) printf(" The run %u is already present in the GL_RUN table...\n",glrun->GetID()); |
6673 |
|
// printf(" CCCCCCCCICCCCCCCCCCCCIOOOOOOOOOO si muove Ciccio! %u \n",glrun->GetID()); |
6674 |
|
// if ( PEDANTIC && glrun->GetID() != 61563 ) throw -83; |
6675 |
|
if ( PEDANTIC ) throw -83; |
6676 |
|
}; |
6677 |
|
// if ( IsDebug() ) printf(" The run already exists in the GL_RUN table! ...\n"); |
6678 |
|
// if ( PEDANTIC ) throw -83; |
6679 |
}; |
}; |
6680 |
if ( IsDebug() ) printf(" Delete run %s from the GL_RUN_FRAGMENTS table \n",row->GetField(0)); |
if ( IsDebug() ) printf(" Delete run %s from the GL_RUN_FRAGMENTS table \n",row->GetField(0)); |
6681 |
// |
// |
6741 |
// 2) get the OBT of the last validated run |
// 2) get the OBT of the last validated run |
6742 |
// -------------------------------------------------------------- |
// -------------------------------------------------------------- |
6743 |
oss.str(""); |
oss.str(""); |
6744 |
oss << " SELECT * FROM GL_RUN WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start |
// oss << " SELECT * FROM GL_RUN WHERE VALIDATION=1 AND RUNHEADER_TIME<="<< t_start |
6745 |
|
oss << " SELECT * FROM GL_RUN WHERE VALIDATION>0 AND RUNHEADER_TIME<="<< t_start |
6746 |
<<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; |
<<" ORDER BY RUNHEADER_TIME DESC LIMIT 1;"; |
6747 |
if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str()); |
if ( IsDebug() ) printf(" Get stop validation-time: query is \n %s \n",oss.str().c_str()); |
6748 |
result = conn->Query(oss.str().c_str()); |
result = conn->Query(oss.str().c_str()); |
6885 |
if( interval >= calibtime )CHECK = true; //more than calibtime s => there might be a calibration |
if( interval >= calibtime )CHECK = true; //more than calibtime s => there might be a calibration |
6886 |
|
|
6887 |
if( !CHECK && this_run->VALIDATION ){ |
if( !CHECK && this_run->VALIDATION ){ |
6888 |
for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true); |
// for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],true); |
6889 |
|
for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],this_run->VALIDATION); |
6890 |
nseq=0; |
nseq=0; |
6891 |
} |
} |
6892 |
|
|
6907 |
if( CHECK ){ |
if( CHECK ){ |
6908 |
// check if calibration exists |
// check if calibration exists |
6909 |
if ( IsDebug() )printf("DT %i ===> CHECK Missing calibration\n",interval); |
if ( IsDebug() )printf("DT %i ===> CHECK Missing calibration\n",interval); |
6910 |
Bool_t MISSING = MissingTRK_CALIB(t1,t2); |
// Bool_t MISSING = MissingTRK_CALIB(t1,t2); |
6911 |
for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],!MISSING); |
UInt_t MISSING = MissingTRK_CALIB(t1,t2); |
6912 |
|
UInt_t val = 0; |
6913 |
|
if ( MISSING == 1 ) val = 0; |
6914 |
|
if ( MISSING == 0 ) val = 1; |
6915 |
|
if ( MISSING == 2 ) val = 2; |
6916 |
|
for (Int_t irun = 0; irun < nseq; irun++)assignVALIDATION(sequence[irun],val); |
6917 |
nseq=0; |
nseq=0; |
6918 |
}; |
}; |
6919 |
//-------------- |
//-------------- |
6943 |
* @param t2 To absolute time |
* @param t2 To absolute time |
6944 |
* @return true if there might be a missing calibration |
* @return true if there might be a missing calibration |
6945 |
*/ |
*/ |
6946 |
Bool_t PamelaDBOperations::MissingTRK_CALIB(UInt_t t1,UInt_t t2){ |
//Bool_t PamelaDBOperations::MissingTRK_CALIB(UInt_t t1,UInt_t t2){ |
6947 |
|
UInt_t PamelaDBOperations::MissingTRK_CALIB(UInt_t t1,UInt_t t2){ |
6948 |
|
|
6949 |
GL_TRK_CALIB* trkcalib = new GL_TRK_CALIB(); |
GL_TRK_CALIB* trkcalib = new GL_TRK_CALIB(); |
6950 |
|
|
6951 |
// get the closest calibration before the run start (t2) |
// get the closest calibration before the run start (t2) |
6952 |
if ( trkcalib->Query_GL_TRK_CALIB(t2, conn) )return(true); //>>> missing |
// if ( trkcalib->Query_GL_TRK_CALIB(t2, conn) )return(true); //>>> missing |
6953 |
|
if ( trkcalib->Query_GL_TRK_CALIB(t2, conn) )return(1); //>>> missing |
6954 |
|
|
6955 |
if ( trkcalib->TO_TIME < t2 ) return(true); //>>> missing |
// if ( trkcalib->TO_TIME < t2 ) return(true); //>>> missing |
6956 |
|
if ( trkcalib->TO_TIME < t2 ) return(1); //>>> missing |
6957 |
|
|
6958 |
//============================================================== |
//============================================================== |
6959 |
// Check is done first on the basis of time between calibration, |
// Check is done first on the basis of time between calibration, |
6970 |
//============================================================== |
//============================================================== |
6971 |
Bool_t DOWNLOAD = false; |
Bool_t DOWNLOAD = false; |
6972 |
// check if the calib was skipped becouse of download .... DA FARE!! |
// check if the calib was skipped becouse of download .... DA FARE!! |
6973 |
if(DOWNLOAD)return(false); |
// if(DOWNLOAD)return(false); |
6974 |
|
if(DOWNLOAD)return(0); |
6975 |
|
|
6976 |
return(true); //>>> missing |
// return(true); //>>> missing |
6977 |
|
return(1); //>>> missing |
6978 |
|
|
6979 |
}; |
}; |
6980 |
|
|
6984 |
//============================================================== |
//============================================================== |
6985 |
// the long time interval bewteen runs might be due to download |
// the long time interval bewteen runs might be due to download |
6986 |
if ( IsDebug() )printf("Short time between calib and run start %u :-) ==> OK! \n",t2 - trkcalib->FROM_TIME); |
if ( IsDebug() )printf("Short time between calib and run start %u :-) ==> OK! \n",t2 - trkcalib->FROM_TIME); |
6987 |
return(false); |
if ( trkcalib->VALIDATION ) return(0); |
6988 |
|
if ( IsDebug() )printf("Calibration is not validated... :-/ ==> OK but with VALIDATION=2! \n"); |
6989 |
|
return(2); |
6990 |
|
|
6991 |
}; |
}; |
6992 |
/** |
/** |
6995 |
* @param validation true/false |
* @param validation true/false |
6996 |
*/ |
*/ |
6997 |
Int_t PamelaDBOperations::assignVALIDATION(UInt_t idrun, Bool_t validation){ |
Int_t PamelaDBOperations::assignVALIDATION(UInt_t idrun, Bool_t validation){ |
6998 |
|
return(this->assignVALIDATION(idrun,(UInt_t)validation)); |
6999 |
|
} |
7000 |
|
|
7001 |
|
Int_t PamelaDBOperations::assignVALIDATION(UInt_t idrun, UInt_t validation){ |
7002 |
TSQLResult *result = 0; |
TSQLResult *result = 0; |
7003 |
stringstream oss; |
stringstream oss; |
7004 |
oss.str(""); |
oss.str(""); |
7195 |
stringstream myquery; |
stringstream myquery; |
7196 |
// |
// |
7197 |
myquery.str(""); |
myquery.str(""); |
7198 |
myquery << " SELECT ID, ID_RAW FROM GL_ROOT where NAME='"<<remfile.Data() <<"';"; |
myquery << " SELECT ID,ID_RAW,ID_TIMESYNC FROM GL_ROOT where NAME='"<<remfile.Data() <<"';"; |
7199 |
// |
// |
7200 |
pResult = conn->Query(myquery.str().c_str()); |
pResult = conn->Query(myquery.str().c_str()); |
7201 |
// |
// |
7211 |
// |
// |
7212 |
this->SetID_ROOT((UInt_t)atoll(Row->GetField(0))); |
this->SetID_ROOT((UInt_t)atoll(Row->GetField(0))); |
7213 |
this->SetID_RAW((UInt_t)atoll(Row->GetField(1))); |
this->SetID_RAW((UInt_t)atoll(Row->GetField(1))); |
7214 |
|
UInt_t idtsy=(UInt_t)atoll(Row->GetField(2)); |
7215 |
// |
// |
7216 |
this->ValidationOFF(); |
this->ValidationOFF(); |
7217 |
// |
// |
7219 |
// |
// |
7220 |
this->RemoveRUNS(); |
this->RemoveRUNS(); |
7221 |
// |
// |
7222 |
this->RemoveFILES(); |
this->RemoveFILES(idtsy); |
7223 |
// |
// |
7224 |
this->SetID_ROOT(0); |
this->SetID_ROOT(0); |
7225 |
this->SetID_RAW(0); |
this->SetID_RAW(0); |
7439 |
* Rearrange calibration tables |
* Rearrange calibration tables |
7440 |
* |
* |
7441 |
**/ |
**/ |
7442 |
void PamelaDBOperations::RemoveFILES(){ |
void PamelaDBOperations::RemoveFILES(UInt_t idtsy){ |
7443 |
stringstream myquery; |
stringstream myquery; |
7444 |
// |
// |
7445 |
myquery.str(""); |
myquery.str(""); |
7449 |
// |
// |
7450 |
conn->Query(myquery.str().c_str()); |
conn->Query(myquery.str().c_str()); |
7451 |
// |
// |
7452 |
|
myquery.str(""); |
7453 |
|
myquery << " DELETE FROM GL_ROOT WHERE ID=" <<this->GetID_ROOT() <<";"; |
7454 |
|
// |
7455 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
7456 |
|
// |
7457 |
|
conn->Query(myquery.str().c_str()); // |
7458 |
|
// |
7459 |
|
if ( !chewbacca ){ |
7460 |
|
myquery.str(""); |
7461 |
|
myquery << " DELETE FROM GL_TIMESYNC WHERE ID=" << idtsy <<";"; |
7462 |
|
// |
7463 |
|
if ( IsDebug() ) printf(" query is \n %s \n",myquery.str().c_str()); |
7464 |
|
// |
7465 |
|
conn->Query(myquery.str().c_str()); |
7466 |
|
}; |
7467 |
|
// |
7468 |
}; |
}; |
7469 |
|
|
7470 |
/** |
/** |
7814 |
|
|
7815 |
|
|
7816 |
} |
} |
|
// UInt_t PamelaDBOperations::ValidateTrkCalib( CalibTrk1Event* caltrk, EventHeader *eh ){ |
|
|
|
|
|
// if(!caltrk || !eh) return 0; |
|
|
|
|
|
// Int_t vorder[]={5,5,3,3,4,4,2,2,1,1,0,0}; |
|
|
// UInt_t timeaftercalib=120000; //2000; |
|
|
// // ---------- |
|
|
// // Check CRCs |
|
|
// // ---------- |
|
|
// for(Int_t ipkt=0; ipkt<6; ipkt++){ |
|
|
// if( caltrk->crc_hcal[ipkt] )return 0; // :-( |
|
|
// for(Int_t ilad=0; ilad<3; ilad++)if( caltrk->crc_cal[ipkt][ilad] )return 0; // :-( |
|
|
// } |
|
|
// // ----------------------- |
|
|
// // Check missing packets: |
|
|
// // ----------------------- |
|
|
// // Readout order: |
|
|
// // ------------------ |
|
|
// // DSP packet board |
|
|
// // ------------------ |
|
|
// // 12 0 1 |
|
|
// // 10 1 1 |
|
|
// // 8 2 1 |
|
|
// // 4 3 1 |
|
|
// // 6 4 1 |
|
|
// // 2 5 1 |
|
|
// // ------------------ |
|
|
// // 11 0 2 |
|
|
// // 9 1 2 |
|
|
// // 7 2 2 |
|
|
// // 3 3 2 |
|
|
// // 5 4 2 |
|
|
// // 1 5 2 |
|
|
// // ------------------ |
|
|
// // ------------------------------------------------- |
|
|
// // Check if it is first or second calibration packet |
|
|
// // ------------------------------------------------- |
|
|
// UInt_t build=0; |
|
|
// TString classname = caltrk->GetName(); |
|
|
// UInt_t base=0; |
|
|
// UInt_t mask=0; |
|
|
// if(classname.Contains("CalibTrk1Event")){ |
|
|
// base=12; |
|
|
// mask=0x03F000; |
|
|
// } |
|
|
// if(classname.Contains("CalibTrk2Event")){ |
|
|
// base=18; |
|
|
// mask=0xFC0000; |
|
|
// } |
|
|
// // ------------------------------------------------- |
|
|
// // Count number of packets and set build variable |
|
|
// // ------------------------------------------------- |
|
|
// Int_t npkts=0; |
|
|
// for(Int_t ipkt=0; ipkt<6; ipkt++){ |
|
|
// if(caltrk->DSPnumber[ipkt]>0 && caltrk->DSPnumber[ipkt]<=12){ |
|
|
// npkts++; |
|
|
// build = build | ( 1<<(base+vorder[caltrk->DSPnumber[ipkt]-1]) ); |
|
|
// } |
|
|
// } |
|
|
// // if( npkts==6 )return 1; // :-) |
|
|
|
|
|
// // cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime()<<endl; |
|
|
|
|
|
// // ----------------------------------------------- |
|
|
// // If missing packets: check the acq configuration |
|
|
// // (some DSPs might be excluded from acquisition) |
|
|
// // ----------------------------------------------- |
|
|
|
|
|
// // ----------------------------------------------- |
|
|
// // retrieve the first run header after calib |
|
|
// // ----------------------------------------------- |
|
|
// PacketType *pctp; |
|
|
// EventCounter *cod; |
|
|
// cod = eh->GetCounter(); |
|
|
// Int_t irun = cod->Get(pctp->RunHeader); |
|
|
// TTree *rh=(TTree*)file->Get("RunHeader"); |
|
|
// if ( !rh || rh->IsZombie() ) throw -17; |
|
|
// if( rh->GetEntries() == irun ){ |
|
|
// if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (1) -- cannot validate :-( "<<endl; |
|
|
// return 0; // :-( |
|
|
// } |
|
|
|
|
|
// RunHeaderEvent *run = 0; |
|
|
// EventHeader *hrun = 0; |
|
|
// rh->SetBranchAddress("RunHeader", &run); |
|
|
// rh->SetBranchAddress("Header", &hrun); |
|
|
// rh->GetEntry(irun); |
|
|
// // cout << classname << " "<<eh->GetPscuHeader()->GetOrbitalTime() << " Run " << hrun->GetPscuHeader()->GetOrbitalTime() <<endl; |
|
|
|
|
|
// if( OBT(hrun->GetPscuHeader()->GetOrbitalTime()) < OBT(eh->GetPscuHeader()->GetOrbitalTime())){ |
|
|
// if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) no runs after calib (2) -- cannot validate :-( "<<endl; |
|
|
// return 0; // :-( |
|
|
// } |
|
|
|
|
|
// if( !run->RM_ACQ_AFTER_CALIB ){ |
|
|
// if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) RM_ACQ_AFTER_CALIB=0 -- cannot validate :-( "<<endl; |
|
|
// return 0; // :-( |
|
|
// } |
|
|
|
|
|
// UInt_t dtime = OBT(hrun->GetPscuHeader()->GetOrbitalTime()) - OBT(eh->GetPscuHeader()->GetOrbitalTime()); |
|
|
// if( dtime > timeaftercalib ){ |
|
|
// if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) run after calib too far ( "<<dtime<<"ms ) -- cannot validate :-( "<<endl; |
|
|
// return 0; // :-( |
|
|
// } |
|
|
|
|
|
|
|
|
|
|
|
// if( (run->ACQ_BUILD_INFO & mask) != build ){ |
|
|
// if ( IsDebug() ) cout << "ValidateTrkCalib: (MISSING VIEW) ACQ_BUILD_INFO= >>> "<<hex << (run->ACQ_BUILD_INFO&mask) << " != "<< build << dec<<endl; |
|
|
// return 0; // :-( |
|
|
// } |
|
|
// return 1; // :-) |
|
|
|
|
|
// } |
|
7817 |
|
|
7818 |
/** |
/** |
7819 |
* |
* |