Bug fix in return code
[u/mrichter/AliRoot.git] / STEER / AliGRPPreprocessor.cxx
1 /**************************************************************************
2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
3  *                                                                        *
4  * Author: The ALICE Off-line Project.                                    *
5  * Contributors are mentioned in the code where appropriate.              *
6  *                                                                        *
7  * Permission to use, copy, modify and distribute this software and its   *
8  * documentation strictly for non-commercial purposes is hereby granted   *
9  * without fee, provided that the above copyright notice appears in all   *
10  * copies and that both the copyright notice and this permission notice   *
11  * appear in the supporting documentation. The authors make no claims     *
12  * about the suitability of this software for any purpose. It is          *
13  * provided "as is" without express or implied warranty.                  *
14  **************************************************************************/
15
16 /* $Id$ */
17
18 //-------------------------------------------------------------------------
19 //                          Class AliGRPPreprocessor
20 //                  Global Run Parameters (GRP) preprocessor
21 //    Origin: Panos Christakoglou, UOA-CERN, Panos.Christakoglou@cern.ch
22 //-------------------------------------------------------------------------
23
24 #include <TChain.h>
25 #include <TList.h>
26 #include <TMap.h>
27 #include <TObjString.h>
28 #include <TTimeStamp.h>
29 #include <TSystem.h>
30 #include <TFile.h>
31
32 #include "AliGRPPreprocessor.h"
33 #include "AliGRPDCS.h"
34 #include "AliDCSSensorArray.h"
35
36 #include "AliCDBMetaData.h"
37 #include "AliLog.h"
38
39 class AliDCSValue;
40 class AliShuttleInterface;
41
42 #include <TH1.h>
43
44 const Double_t kFitFraction = 0.7;                 // Fraction of DCS sensor fits required
45
46 ClassImp(AliGRPPreprocessor)
47
48 //_______________________________________________________________
49   const char* AliGRPPreprocessor::fgkDCSDataPoints[12] = {"LHCState","LHCPeriod","LHCLuminosity","BeamIntensity","L3Current","L3Polarity","DipoleCurrent","DipolePolarity","CavernTemperature","CavernAtmosPressure","gva_cr5AtmosphericPressure","gva_meyrinAtmosphericPressure"};
50
51 //_______________________________________________________________
52 AliGRPPreprocessor::AliGRPPreprocessor():
53   AliPreprocessor("GRP",0), fPressure(0) {
54   // default constructor - Don't use this!
55   
56 }
57
58 //_______________________________________________________________
59 AliGRPPreprocessor::AliGRPPreprocessor(AliShuttleInterface* shuttle):
60   AliPreprocessor("GRP",shuttle), fPressure(0) {
61   // constructor - shuttle must be instantiated!
62   
63 }
64
65 //_______________________________________________________________
66 AliGRPPreprocessor::~AliGRPPreprocessor() {
67   //destructor
68   delete fPressure;
69 }
70
71 //_______________________________________________________________
72 void AliGRPPreprocessor::Initialize(Int_t run, UInt_t startTime, UInt_t endTime) {
73   // Initialize preprocessor
74   
75   AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run, TTimeStamp(startTime).AsString(), TTimeStamp(endTime).AsString()));
76   
77   fRun = run;
78   fStartTime = startTime;
79   fEndTime = endTime;
80   AliInfo("Initialization of the GRP preprocessor.");
81
82   TClonesArray * array = new TClonesArray("AliDCSSensor",2);
83   for(Int_t j = 0; j < 2; j++) {
84     AliDCSSensor * sens = new ((*array)[j])AliDCSSensor;
85     sens->SetStringID(fgkDCSDataPoints[j+10]);
86   }
87   AliInfo(Form("Pressure Entries: %d",array->GetEntries()));
88
89   fPressure = new AliDCSSensorArray(fStartTime, fEndTime, array);
90 }
91
92 //_______________________________________________________________
93 UInt_t AliGRPPreprocessor::Process(TMap* valueMap) {
94   // process data retrieved by the Shuttle
95   
96   //=================//
97   // DAQ logbook     //
98   //=================//
99   
100   TList *daqlblist = ProcessDaqLB();
101   if(!daqlblist) {
102     Log(Form("Problem with the DAQ logbook parameters!!!"));
103     return 1;
104   }
105   TMap *m1 = (TMap *)daqlblist->At(0);
106   TObjString *s1 = (TObjString *)m1->GetValue("fAliceStartTime");
107   UInt_t iStartTime = atoi(s1->String().Data());
108   TMap *m2 = (TMap *)daqlblist->At(1);
109   TObjString *s2 = (TObjString *)m2->GetValue("fAliceStopTime");
110   UInt_t iStopTime = atoi(s2->String().Data());
111
112   //=================//
113   // DAQ FXS         //
114   //=================//
115   UInt_t iDaqFxs = ProcessDaqFxs();
116   if(iDaqFxs == 0) {
117         Log(Form("ProcessDaqFxs successful!"));
118   } else {
119         Log(Form("Could not store run raw tag file!"));
120         return 1;
121   }
122   
123   
124   //=================//
125   // DCS data points //
126   //=================//
127   TList *dcsdplist = ProcessDcsDPs(valueMap, iStartTime, iStopTime);
128   if(!dcsdplist) {
129     Log(Form("Problem with the DCS data points!!!"));
130     return 1; 
131   }    
132   if(dcsdplist->GetEntries() != 10) {
133     Log(Form("Problem with the DCS data points!!!"));
134     // return 1; // TODO:COMMENTED FOR TESTING PURPOSES!
135   }
136   //NEEDS TO BE REVISED - BREAKS!!!
137 //   AliDCSSensorArray *dcsSensorArray = GetPressureMap(valueMap,fPressure);
138 //   if(!dcsSensorArray) {
139 //     Log(Form("Problem with the pressure sensor values!!!"));
140 //     return 0;
141 //   }
142
143   TList * list = new TList();
144   list = GetGlobalList(daqlblist,dcsdplist);
145   list->SetOwner(1);
146   AliInfo(Form("Final list entries: %d",list->GetEntries()));
147   
148   AliCDBMetaData md;
149   md.SetResponsible("Panos Christakoglou");
150   md.SetComment("Output parameters from the GRP preprocessor.");
151   
152   Bool_t result = Store("GRP", "Data", list, &md);
153   
154   delete list;
155   
156   if (result)
157     return 0;
158   else
159     return 1;
160 }
161
162 //_______________________________________________________________
163 TList *AliGRPPreprocessor::GetGlobalList(TList *l1, TList *l2) {
164   //Getting the global output TList
165   TList *list = new TList();
166   TMap *map = new TMap();
167   for(Int_t i = 0; i < l1->GetEntries(); i++) 
168     list->AddLast(map = (TMap *)l1->At(i));
169   for(Int_t i = 0; i < l2->GetEntries(); i++) 
170     list->AddLast(map = (TMap *)l2->At(i));
171
172   return list;
173 }
174
175 //_______________________________________________________________
176 TList *AliGRPPreprocessor::ProcessDaqLB() {
177   //Getting the DAQ lb informnation
178   const char* timeStart = GetRunParameter("time_start");
179   const char* timeEnd = GetRunParameter("time_end");
180   const char* beamEnergy = GetRunParameter("beamEnergy");
181   const char* beamType = GetRunParameter("beamType");
182   const char* numberOfDetectors = GetRunParameter("numberOfDetectors");
183   const char* detectorMask = GetRunParameter("detectorMask");
184   const char* lhcPeriod = GetRunParameter("LHCperiod");
185
186   if (timeStart) {
187     Log(Form("Start time for run %d: %s",fRun, timeStart));
188   } else {
189     Log(Form("Start time not put in logbook!"));
190   }
191   TMap *mapDAQ1 = new TMap();
192   mapDAQ1->Add(new TObjString("fAliceStartTime"),new TObjString(timeStart));
193
194   if (timeEnd) {
195     Log(Form("End time for run %d: %s",fRun, timeEnd));
196   } else {
197     Log(Form("End time not put in logbook!"));
198   }
199   TMap *mapDAQ2 = new TMap();
200   mapDAQ2->Add(new TObjString("fAliceStopTime"),new TObjString(timeEnd));
201
202   if (beamEnergy) {
203     Log(Form("Beam energy for run %d: %s",fRun, beamEnergy));
204   } else {
205     Log(Form("Beam energy not put in logbook!"));
206   }
207   TMap *mapDAQ3 = new TMap();
208   mapDAQ3->Add(new TObjString("fAliceBeamEnergy"),new TObjString(beamEnergy));
209
210   if (beamType) {
211     Log(Form("Beam type for run %d: %s",fRun, beamType));
212   } else {
213     Log(Form("Beam type not put in logbook!"));
214   }
215   TMap *mapDAQ4 = new TMap();
216   mapDAQ4->Add(new TObjString("fAliceBeamType"),new TObjString(beamType));
217
218   if (numberOfDetectors) {
219     Log(Form("Number of active detectors for run %d: %s",fRun, numberOfDetectors));
220   } else {
221     Log(Form("Number of active detectors not put in logbook!"));
222   }
223   TMap *mapDAQ5 = new TMap();
224   mapDAQ5->Add(new TObjString("fNumberOfDetectors"),new TObjString(numberOfDetectors));
225
226   if (detectorMask) {
227     Log(Form("Detector mask for run %d: %s",fRun, detectorMask));
228   } else {
229     Log(Form("Detector mask not put in logbook!"));
230   }
231   TMap *mapDAQ6 = new TMap();
232   mapDAQ6->Add(new TObjString("fDetectorMask"),new TObjString(detectorMask));
233
234   if (lhcPeriod) {
235     Log(Form("LHC period (DAQ) for run %d: %s",fRun, lhcPeriod));
236   } else {
237     Log(Form("LHCperiod not put in logbook!"));
238   }
239   TMap *mapDAQ7 = new TMap();
240   mapDAQ7->Add(new TObjString("fLHCPeriod"),new TObjString(lhcPeriod));
241
242   TList *list = new TList();
243   list->Add(mapDAQ1); list->Add(mapDAQ2);
244   list->Add(mapDAQ3); list->Add(mapDAQ4);
245   list->Add(mapDAQ5); list->Add(mapDAQ6);
246   list->Add(mapDAQ7);
247
248   return list;
249 }
250
251 //_______________________________________________________________
252 UInt_t AliGRPPreprocessor::ProcessDaqFxs() {
253   //======DAQ FXS======//
254
255   TList* list = GetFileSources(kDAQ);  
256   if (!list) {
257     Log("No raw data tag list: connection problems with DAQ FXS logbook!");
258     return 1;
259   }
260   
261   if (list->GetEntries() == 0)
262   {
263         Log("no raw data tags in this run: nothing to merge!");
264         delete  list; list=0;
265         return 0;
266   }
267
268   TChain *fRawTagChain = new TChain("T");
269   Int_t nFiles=0;
270   TIterator* iter = list->MakeIterator();
271   TObject* obj = 0;
272   while ((obj = iter->Next())) {
273     TObjString* objStr = dynamic_cast<TObjString*> (obj);
274     if (objStr) {
275       Log(Form("Found source %s", objStr->String().Data()));
276       TList* list2 = GetFileIDs(kDAQ, objStr->String());
277       if (!list2) {
278         Log("No list with ids from DAQ was found: connection problems with DAQ FXS logbook!");
279         delete fRawTagChain; fRawTagChain=0;
280         return 1;
281       }
282       Log(Form("Number of ids: %d",list2->GetEntries()));
283       for(Int_t i = 0; i < list2->GetEntries(); i++) {
284         TObjString *idStr = (TObjString *)list2->At(i);
285         TString fileName = GetFile(kDAQ,idStr->String().Data(),objStr->String().Data());
286         if (fileName.Length() > 0) 
287         {      
288                 Log(Form("Adding file in the chain: %s",fileName.Data()));
289                 fRawTagChain->Add(fileName.Data());
290                 nFiles++;
291         } else {
292                 Log(Form("Could not retrieve file with id %s from source %s: "
293                         "connection problems with DAQ FXS!",
294                                 idStr->String().Data(),objStr->String().Data()));
295                 delete list; list=0;
296                 delete list2; list2=0;
297                 delete fRawTagChain; fRawTagChain=0;
298                 return 2;
299         }
300       }
301       delete list2;
302     }
303   }
304   
305   TString fRawDataFileName = "GRP_Merged.tag.root";
306   Log(Form("Merging %d raw data tags into file: %s", nFiles, fRawDataFileName.Data()));
307   fRawTagChain->Merge(fRawDataFileName);
308   
309   TString outputfile = Form("Run%d.Merged.RAW.tag.root", fRun);
310   Bool_t result = StoreRunMetadataFile(fRawDataFileName.Data(),outputfile.Data());
311   
312   if (!result)
313   {
314         Log("Problem storing raw data tags in local file!!");
315   } else {
316         Log("Raw data tags merged successfully!!");  
317   }
318   
319   delete iter;
320   delete list;
321   delete fRawTagChain; fRawTagChain=0;
322   
323   if (result == kFALSE)
324   {
325         return 3;
326   }
327   
328   return 3;
329   
330 }
331
332 //_______________________________________________________________
333 TList *AliGRPPreprocessor::ProcessDcsDPs(TMap* valueMap, UInt_t iStartTime, UInt_t iStopTime) {
334   //Getting the DCS dps
335   //===========//
336   
337   TList *list = new TList();
338
339   //DCS data points
340   //===========//
341   AliInfo(Form("==========LHCState==========="));
342   TObjArray *aliasLHCState = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[0]);
343   if(!aliasLHCState) {
344     Log(Form("LHCState not found!!!"));
345     return list;
346   }
347   AliGRPDCS *dcs1 = new AliGRPDCS(aliasLHCState,iStartTime,iStopTime);
348   TString sLHCState = dcs1->ProcessDCS(3);  
349   if (sLHCState) {
350     Log(Form("<LHCState> for run %d: %s",fRun, sLHCState.Data()));
351   } else {
352     Log(Form("LHCState not put in TMap!"));
353   }
354   TMap *mapDCS1 = new TMap();
355   mapDCS1->Add(new TObjString("fLHCState"),new TObjString(sLHCState));
356   list->Add(mapDCS1);
357
358   AliInfo(Form("==========LHCPeriod==========="));
359   TObjArray *aliasLHCPeriod = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[1]);
360   if(!aliasLHCPeriod) {
361     Log(Form("LHCPeriod not found!!!"));
362     return list;
363   }
364   AliGRPDCS *dcs2 = new AliGRPDCS(aliasLHCPeriod,iStartTime,iStopTime);
365   TString sLHCPeriod = dcs2->ProcessDCS(3);  
366   if (sLHCPeriod) {
367     Log(Form("<LHCPeriod> for run %d: %s",fRun, sLHCPeriod.Data()));
368   } else {
369     Log(Form("LHCPeriod not put in TMap!"));
370   }
371   TMap *mapDCS2 = new TMap();
372   mapDCS2->Add(new TObjString("fLHCCondition"),new TObjString(sLHCPeriod));
373   list->Add(mapDCS2);
374
375   AliInfo(Form("==========LHCLuminosity==========="));
376   TObjArray *aliasLHCLuminosity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[2]);
377   if(!aliasLHCLuminosity) {
378     Log(Form("LHCLuminosity not found!!!"));
379     return list;
380   }
381   AliGRPDCS *dcs3 = new AliGRPDCS(aliasLHCLuminosity,iStartTime,iStopTime);
382   TString sMeanLHCLuminosity = dcs3->ProcessDCS(2);  
383   if (sMeanLHCLuminosity) {
384     Log(Form("<LHCLuminosity> for run %d: %s",fRun, sMeanLHCLuminosity.Data()));
385   } else {
386     Log(Form("LHCLuminosity not put in TMap!"));
387   }
388   TMap *mapDCS3 = new TMap();
389   mapDCS3->Add(new TObjString("fLHCLuminosity"),new TObjString(sMeanLHCLuminosity));
390   list->Add(mapDCS3);
391
392   AliInfo(Form("==========BeamIntensity==========="));
393   TObjArray *aliasBeamIntensity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[3]);
394   if(!aliasBeamIntensity) {
395     Log(Form("BeamIntensity not found!!!"));
396     return list;
397   }
398   AliGRPDCS *dcs4 = new AliGRPDCS(aliasBeamIntensity,iStartTime,iStopTime);
399   TString sMeanBeamIntensity = dcs4->ProcessDCS(2);  
400   if (sMeanBeamIntensity) {
401     Log(Form("<BeamIntensity> for run %d: %s",fRun, sMeanBeamIntensity.Data()));
402   } else {
403     Log(Form("BeamIntensity not put in TMap!"));
404   }
405   TMap *mapDCS4 = new TMap();
406   mapDCS4->Add(new TObjString("fBeamIntensity"),new TObjString(sMeanBeamIntensity));
407   list->Add(mapDCS4);
408
409   AliInfo(Form("==========L3Current==========="));
410   TObjArray *aliasL3Current = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[4]);
411   if(!aliasL3Current) {
412     Log(Form("L3Current not found!!!"));
413     return list;
414   }
415   AliGRPDCS *dcs5 = new AliGRPDCS(aliasL3Current,iStartTime,iStopTime);
416   TString sMeanL3Current = dcs5->ProcessDCS(2);  
417   if (sMeanL3Current) {
418     Log(Form("<L3Current> for run %d: %s",fRun, sMeanL3Current.Data()));
419   } else {
420     Log(Form("L3Current not put in TMap!"));
421   }
422   TMap *mapDCS5 = new TMap();
423   mapDCS5->Add(new TObjString("fL3Current"),new TObjString(sMeanL3Current));
424   list->Add(mapDCS5);
425
426   AliInfo(Form("==========L3Polarity==========="));
427   TObjArray *aliasL3Polarity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[5]);
428   if(!aliasL3Polarity) {
429     Log(Form("L3Polarity not found!!!"));
430     return list;
431   }
432   AliGRPDCS *dcs6 = new AliGRPDCS(aliasL3Polarity,iStartTime,iStopTime);
433   TString sL3Polarity = dcs6->ProcessDCS(4);  
434   if (sL3Polarity) {
435     Log(Form("<L3Polarity> for run %d: %s",fRun, sL3Polarity.Data()));
436   } else {
437     Log(Form("L3Polarity not put in TMap!"));
438   }
439   TMap *mapDCS6 = new TMap();
440   mapDCS6->Add(new TObjString("fL3Polarity"),new TObjString(sL3Polarity));
441   list->Add(mapDCS6);
442
443   AliInfo(Form("==========DipoleCurrent==========="));
444   TObjArray *aliasDipoleCurrent = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[6]);
445   if(!aliasDipoleCurrent) {
446     Log(Form("DipoleCurrent not found!!!"));
447     return list;
448   }
449   AliGRPDCS *dcs7 = new AliGRPDCS(aliasDipoleCurrent,iStartTime,iStopTime);
450   TString sMeanDipoleCurrent = dcs7->ProcessDCS(2);  
451   if (sMeanDipoleCurrent) {
452     Log(Form("<DipoleCurrent> for run %d: %s",fRun, sMeanDipoleCurrent.Data()));
453   } else {
454     Log(Form("DipoleCurrent not put in TMap!"));
455   }
456   TMap *mapDCS7 = new TMap();
457   mapDCS7->Add(new TObjString("fDipoleCurrent"),new TObjString(sMeanDipoleCurrent));
458   list->Add(mapDCS7);
459
460   AliInfo(Form("==========DipolePolarity==========="));
461   TObjArray *aliasDipolePolarity = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[7]);
462   if(!aliasDipolePolarity) {
463     Log(Form("DipolePolarity not found!!!"));
464     return list;
465   }
466   AliGRPDCS *dcs8 = new AliGRPDCS(aliasDipolePolarity,iStartTime,iStopTime);
467   TString sDipolePolarity = dcs8->ProcessDCS(4);  
468   if (sDipolePolarity) {
469     Log(Form("<DipolePolarity> for run %d: %s",fRun, sDipolePolarity.Data()));
470   } else {
471     Log(Form("DipolePolarity not put in TMap!"));
472   }
473   TMap *mapDCS8 = new TMap();
474   mapDCS8->Add(new TObjString("fDipolePolarity"),new TObjString(sDipolePolarity));
475   list->Add(mapDCS8);
476
477   AliInfo(Form("==========CavernTemperature==========="));
478   TObjArray *aliasCavernTemperature = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[8]);
479   if(!aliasCavernTemperature) {
480     Log(Form("CavernTemperature not found!!!"));
481     return list;
482   }
483   AliGRPDCS *dcs9 = new AliGRPDCS(aliasCavernTemperature,iStartTime,iStopTime);
484   TString sMeanCavernTemperature = dcs9->ProcessDCS(2);  
485   if (sMeanCavernTemperature) {
486     Log(Form("<CavernTemperature> for run %d: %s",fRun, sMeanCavernTemperature.Data()));
487   } else {
488     Log(Form("CavernTemperature not put in TMap!"));
489   }
490   TMap *mapDCS9 = new TMap();
491   mapDCS9->Add(new TObjString("fCavernTemperature"),new TObjString(sMeanCavernTemperature));
492   list->Add(mapDCS9);
493
494   AliInfo(Form("==========CavernPressure==========="));
495   TObjArray *aliasCavernPressure = (TObjArray *)valueMap->GetValue(fgkDCSDataPoints[9]);
496   if(!aliasCavernPressure) {
497     Log(Form("CavernPressure not found!!!"));
498     return list;
499   }
500   AliGRPDCS *dcs10 = new AliGRPDCS(aliasCavernPressure,iStartTime,iStopTime);
501   TString sMeanCavernPressure = dcs10->ProcessDCS(2);  
502   if (sMeanCavernPressure) {
503     Log(Form("<CavernPressure> for run %d: %s",fRun, sMeanCavernPressure.Data()));
504   } else {
505     Log(Form("CavernPressure not put in TMap!"));
506   }
507   TMap *mapDCS10 = new TMap();
508   mapDCS10->Add(new TObjString("fCavernPressure"),new TObjString(sMeanCavernPressure));
509   list->Add(mapDCS10);
510
511   return list;
512 }
513
514 //_______________________________________________________________
515 AliDCSSensorArray *AliGRPPreprocessor::GetPressureMap(TMap* dcsAliasMap, AliDCSSensorArray *fPressure) {
516   // extract DCS pressure maps. Perform fits to save space
517   
518   TMap *map = fPressure->ExtractDCS(dcsAliasMap);
519   if (map) {
520     fPressure->MakeSplineFit(map);
521     Double_t fitFraction = fPressure->NumFits()/fPressure->NumSensors(); 
522     if (fitFraction > kFitFraction ) {
523       AliInfo(Form("Pressure values extracted, fits performed.\n"));
524     } else { 
525       AliInfo("Too few pressure maps fitted. \n");
526     }
527   } else {
528     AliInfo("AliGRPDCS: no atmospheric pressure map extracted. \n");
529   }
530   delete map;
531  
532   return fPressure;
533 }
534
535 //_______________________________________________________________
536 /*UInt_t AliGRPPreprocessor::MapPressure(TMap* dcsAliasMap) {
537   // extract DCS pressure maps. Perform fits to save space
538   
539   UInt_t result=0;
540   TMap *map = fPressure->ExtractDCS(dcsAliasMap);
541   if (map) {
542     fPressure->MakeSplineFit(map);
543     Double_t fitFraction = fPressure->NumFits()/fPressure->NumSensors(); 
544     if (fitFraction > kFitFraction ) {
545       AliInfo(Form("Pressure values extracted, fits performed.\n"));
546     } else { 
547       Log ("Too few pressure maps fitted. \n");
548       result = 9;
549     }
550   } else {
551     Log("AliTPCPreprocsessor: no atmospheric pressure map extracted. \n");
552     result=9;
553   }
554   delete map;
555   // Now store the final CDB file
556   
557   if ( result == 0 ) {
558     AliCDBMetaData metaData;
559     metaData.SetBeamPeriod(0);
560     metaData.SetResponsible("Panos Christakoglou");
561     metaData.SetComment("Preprocessor AliGRP data base pressure entries.");
562     
563     Bool_t storeOK = Store("Calib", "Pressure", fPressure, &metaData, 0, 0);
564     if ( !storeOK ) result=1; 
565   }
566   
567   return result; 
568   }*/