[u/mrichter/AliRoot.git] / TOF / AliTOFPreprocessor.cxx

/**************************************************************************
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 *                                                                        *
 * Author: The ALICE Off-line Project.                                    *
 * Contributors are mentioned in the code where appropriate.              *
 *                                                                        *
 * Permission to use, copy, modify and distribute this software and its   *
 * documentation strictly for non-commercial purposes is hereby granted   *
 * without fee, provided that the above copyright notice appears in all   *
 * copies and that both the copyright notice and this permission notice   *
 * appear in the supporting documentation. The authors make no claims     *
 * about the suitability of this software for any purpose. It is          *
 * provided "as is" without express or implied warranty.                  *
 **************************************************************************/

/* $Id$ */

#include <Riostream.h>
#include <stdio.h>
#include <stdlib.h>

#include <TFile.h>
#include <TH1.h>
#include <TH1F.h>
#include <TH1S.h>
#include <TH2S.h>
#include <TMath.h>
#include <TObjArray.h>
#include <TObjString.h>
#include <TTimeStamp.h>

#include "AliCDBMetaData.h"
#include "AliLog.h"
#include "AliTOFChannelOnline.h"
#include "AliTOFDataDCS.h"
#include "AliTOFGeometry.h"
#include "AliTOFPreprocessor.h"

class TF1;
class AliDCSValue;
class AliTOFGeometry;

// TOF preprocessor class.
// It takes data from DCS and passes them to the class AliTOFDataDCS, which
// processes them. The result is then written to the CDB.
// analogously, it takes data form DAQ (both at Run level and inclusive - 
// of all the runs - level, processes them, and stores both Reference Data
// and Online Calibration files in the CDB. 


ClassImp(AliTOFPreprocessor)

const Int_t AliTOFPreprocessor::fgkBinRangeAve    = 13; // number of bins where to calculate the mean 
const Double_t AliTOFPreprocessor::fgkThrPar    = 0.013; // parameter used to trigger the calculation of the delay

//_____________________________________________________________________________

AliTOFPreprocessor::AliTOFPreprocessor(AliShuttleInterface* shuttle) :
  AliPreprocessor("TOF", shuttle),
  fData(0),
  fh2(0),
  fCal(0),
  fNChannels(0),
  fStoreRefData(kTRUE)
{
  // constructor

}

//_____________________________________________________________________________

AliTOFPreprocessor::~AliTOFPreprocessor()
{
  // destructor
  if (fData){
    delete fData;
    fData = 0;
  }
  if (fh2){
    delete fh2;
    fh2 = 0;
  }
  if (fCal){
    //    fCal->Clear();
    delete fCal;
    fCal = 0;
  }
}

//______________________________________________________________________________
void AliTOFPreprocessor::Initialize(Int_t run, UInt_t startTime,
	UInt_t endTime)
{
  // Creates AliTOFDataDCS object

  AliPreprocessor::Initialize(run, startTime, endTime);

	AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run,
		TTimeStamp(startTime).AsString(),
		TTimeStamp(endTime).AsString()));

	fData = new AliTOFDataDCS(fRun, fStartTime, fEndTime);
	fh2 = 0x0;
	fNChannels = AliTOFGeometry::NSectors()*(2*(AliTOFGeometry::NStripC()+AliTOFGeometry::NStripB())+AliTOFGeometry::NStripA())*AliTOFGeometry::NpadZ()*AliTOFGeometry::NpadX();
	fCal = new TObjArray(fNChannels);
	fCal->SetOwner();
	for (Int_t ich = 0; ich<fNChannels; ich ++){
	  AliTOFChannelOnline * calChOnline = new AliTOFChannelOnline();
	  fCal->AddAt(calChOnline,ich);
	}
}

//_____________________________________________________________________________

UInt_t AliTOFPreprocessor::Process(TMap* dcsAliasMap)
{
  // Fills data into a AliTOFDataDCS object
  // return codes:
  // return=0 : all ok
  // return=1 : no DCS input data Map
  // return=2 : no DCS input data processing
  // return=3 : no DCS processed data was stored in Ref Data
  // return=4 : no DAQ input for Ref Data
  // return=5 : failed to store Ref Data
  // return=6 : failed to retrieve DAQ data for calibration 
  // return=7 : problems in histos in the input DAQ file 
  // return=8 : failed to store Online Delays

  TH1::AddDirectory(0);

  Bool_t resultDCSMap=kFALSE;
  Bool_t resultDCSStore=kFALSE;
  Bool_t resultDAQ=kFALSE;
  Bool_t resultDAQRef=kFALSE;

  // processing DCS

  if (!dcsAliasMap){
    Log("No DCS map found: TOF exiting from Shuttle");
    return 1;// return error Code for DCS input data not found 
  }
  else {
  // The processing of the DCS input data is forwarded to AliTOFDataDCS
    resultDCSMap=fData->ProcessData(*dcsAliasMap);
    if(!resultDCSMap){
      Log("Some problems occurred while processing DCS data, TOF exiting from Shuttle");
      return 2;// return error Code for processed DCS data not stored 
    }
    else{
      AliCDBMetaData metaDataDCS;
      metaDataDCS.SetBeamPeriod(0);
      metaDataDCS.SetResponsible("Chiara Zampolli");
      metaDataDCS.SetComment("This preprocessor fills an AliTOFDataDCS object.");
      AliInfo("Storing DCS Data");
      resultDCSStore = StoreReferenceData("Calib","DCSData",fData, &metaDataDCS);
      if (!resultDCSStore){
	Log("Some problems occurred while storing DCS data results in Reference Data, TOF exiting from Shuttle");
	return 3;// return error Code for processed DCS data not stored 
	         // in reference data
	
      }
    }
  }
  
  // processing DAQ
  
  TFile * daqFile=0x0;
  
  if(fStoreRefData){
    //retrieving data at Run level
    TList* list = GetFileSources(kDAQ, "RUNLevel");
    if (list)
      {
	AliInfo("The following sources produced files with the id RUNLevel");
	list->Print();
	for (Int_t jj=0;jj<list->GetEntries();jj++){
	  TObjString * str = dynamic_cast<TObjString*> (list->At(jj));
	  AliInfo(Form("found source %s", str->String().Data()));
	  // file to be stored run per run
	  TString fileNameRun = GetFile(kDAQ, "RUNLevel", str->GetName());
	  if (fileNameRun.Length()>0){
	    AliInfo(Form("Got the file %s, now we can store the Reference Data for the current Run.", fileNameRun.Data()));
	    daqFile = new TFile(fileNameRun.Data(),"READ");
	    fh2 = (TH2S*) daqFile->Get("htof");
	    AliCDBMetaData metaDataHisto;
	    metaDataHisto.SetBeamPeriod(0);
	    metaDataHisto.SetResponsible("Chiara Zampolli");
	    metaDataHisto.SetComment("This preprocessor stores the array of histos object as Reference Data.");
	    AliInfo("Storing Reference Data");
	    resultDAQRef = StoreReferenceData("Calib","DAQData",fh2, &metaDataHisto);
	    if (!resultDAQRef){
	      Log("some problems occurred::No Reference Data stored, TOF exiting from Shuttle");
	      return 5;//return error code for failure in storing Ref Data 
	    }
	    daqFile->Close();
	    delete daqFile;
	  }
	  
	  else{
	    Log("The input data file from DAQ (run-level) was not found, TOF exiting from Shuttle "); 
	    return 4;//return error code for failure in retrieving Ref Data 
	  }
	}
      }
    else{
      Log("The input data file list from DAQ (run-level) was not found, TOF exiting from Shuttle "); 
      return 4;//return error code for failure in retrieving Ref Data 
    }	
  }


//Total files, with cumulative histos
  
  TList* listTot = GetFileSources(kDAQ, "DELAYS");
  if (listTot)
    {
      AliInfo("The following sources produced files with the id DELAYS");
      listTot->Print();
      for (Int_t jj=0;jj<listTot->GetEntries();jj++){
	TObjString * str = dynamic_cast<TObjString*> (listTot->At(jj));
	AliInfo(Form("found source %s", str->String().Data()));

	// file with summed histos, to extract calib params
	TString fileName = GetFile(kDAQ, "DELAYS", str->GetName());
	if (fileName.Length()>0){
	  AliInfo(Form("Got the file %s, now we can extract some values.", fileName.Data()));

	  daqFile = new TFile(fileName.Data(),"READ");
	  if (fh2) delete fh2;
	  fh2 = (TH2S*) daqFile->Get("htoftot");
	  if (!fh2){
	    Log("some problems occurred:: No histo retrieved, TOF exiting from Shuttle");
	    delete daqFile;
	    return 7; //return error code for histograms not existing/junky
	  }
	  else {
	    static const Int_t kSize=fh2->GetNbinsX();
	    static const Int_t kNBins=fh2->GetNbinsY();
	    static const Double_t kXBinmin=fh2->GetYaxis()->GetBinLowEdge(1);
	    if (kSize != fNChannels){
	      Log(" number of bins along x different from number of pads, found only a subset of the histograms, TOF exiting from Shuttle");
	      delete daqFile;
              return 7; //return error code for histograms not existing/junky
	    }
	    for (Int_t ich=0;ich<kSize;ich++){
	      TH1S *h1 = new TH1S("h1","h1",kNBins,kXBinmin-0.5,kNBins*1.+kXBinmin-0.5);
	      for (Int_t ibin=0;ibin<kNBins;ibin++){
		h1->SetBinContent(ibin+1,fh2->GetBinContent(ich+1,ibin+1));
	      }
	      Bool_t found=kFALSE; 
	      Float_t minContent=h1->Integral()*fgkThrPar; 
	      Int_t nbinsX = h1->GetNbinsX();
	      Int_t startBin=1;
	      for (Int_t j=1; j<=nbinsX; j++){
		if ((
		     h1->GetBinContent(j) +     
		     h1->GetBinContent(j+1)+
		     h1->GetBinContent(j+2)+ 
		     h1->GetBinContent(j+3))>minContent){
		  found=kTRUE;
		  startBin=j;
		  break;
		}
	      }
	      if(!found) AliInfo(Form("WARNING!!! no start of fit found for histo # %i",ich));
	      // Now calculate the mean over the interval. 
	      Double_t mean = 0;
	      Double_t sumw2 = 0;
	      Double_t nent = 0;
	      for(Int_t k=0;k<fgkBinRangeAve;k++){
		mean=mean+h1->GetBinCenter(startBin+k)*h1->GetBinContent(startBin+k);                 
		nent=nent+h1->GetBinContent(startBin+k);                 
		sumw2=sumw2+(h1->GetBinCenter(startBin+k))*(h1->GetBinCenter(startBin+k))*(h1->GetBinContent(startBin+k));
	      }
	      mean= mean/nent; //<x>
	      sumw2=sumw2/nent; //<x^2>
	      Double_t rmsmean= 0;
	      rmsmean = TMath::Sqrt((sumw2-mean*mean)/nent);
	      if (ich<fNChannels) {
		AliTOFChannelOnline * ch = (AliTOFChannelOnline *)fCal->At(ich);
		ch->SetDelay((Double_t)mean*AliTOFGeometry::TdcBinWidth()*1.E-3);  // delay in ns
	      }
	    delete h1;
	    h1=0x0;
	    }
	  }
	  daqFile->Close();
	  delete daqFile;
	  AliCDBMetaData metaData;
	  metaData.SetBeamPeriod(0);
	  metaData.SetResponsible("Chiara Zampolli");
	  metaData.SetComment("This preprocessor fills an AliTOFCal object.");
	  AliInfo("Storing Calibration Data");
	  resultDAQ = Store("Calib","ParOnline",fCal, &metaData);
          if(!resultDAQ){
	    Log("Some problems occurred while storing DAQ data processing results");
	    return 8;//return error code for problems in storing DAQ data 
	  }
	}
	else{
	  Log("The Cumulative data file from DAQ does not exist, TOF exiting from Shuttle"); 
          return 6;//return error code for problems in retrieving DAQ data 
	}
      }
    }
  else{
    Log("Problem: no list for Cumulative data file from DAQ was found, TOF exiting from Shuttle");
    return 6; //return error code for problems in retrieving DAQ data 
  }

  daqFile=0;

  return 0;
}
Commit	Line	Data
c9fe8530	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
9edefa04	16	/* $Id$ */
9d883ed9	17
9edefa04	18	#include <Riostream.h>
	19	#include <stdio.h>
	20	#include <stdlib.h>
c9fe8530	21
9edefa04	22	#include <TFile.h>
	23	#include <TH1.h>
	24	#include <TH1F.h>
	25	#include <TH1S.h>
	26	#include <TH2S.h>
	27	#include <TMath.h>
	28	#include <TObjArray.h>
	29	#include <TObjString.h>
	30	#include <TTimeStamp.h>
c9fe8530	31
	32	#include "AliCDBMetaData.h"
	33	#include "AliLog.h"
c9fe8530	34	#include "AliTOFChannelOnline.h"
9edefa04	35	#include "AliTOFDataDCS.h"
10056aa6	36	#include "AliTOFGeometry.h"
9edefa04	37	#include "AliTOFPreprocessor.h"
c9fe8530	38
	39	class TF1;
	40	class AliDCSValue;
	41	class AliTOFGeometry;
	42
	43	// TOF preprocessor class.
	44	// It takes data from DCS and passes them to the class AliTOFDataDCS, which
	45	// processes them. The result is then written to the CDB.
	46	// analogously, it takes data form DAQ (both at Run level and inclusive -
	47	// of all the runs - level, processes them, and stores both Reference Data
	48	// and Online Calibration files in the CDB.
	49
	50
	51	ClassImp(AliTOFPreprocessor)
	52
	53	const Int_t AliTOFPreprocessor::fgkBinRangeAve = 13; // number of bins where to calculate the mean
9d883ed9	54	const Double_t AliTOFPreprocessor::fgkThrPar = 0.013; // parameter used to trigger the calculation of the delay
c9fe8530	55
	56	//_____________________________________________________________________________
	57
708db10b	58	AliTOFPreprocessor::AliTOFPreprocessor(AliShuttleInterface* shuttle) :
708db10b	59	AliPreprocessor("TOF", shuttle),
c9fe8530	60	fData(0),
708db10b	61	fh2(0),
c9fe8530	62	fCal(0),
10056aa6	63	fNChannels(0),
5936ab02	64	fStoreRefData(kTRUE)
c9fe8530	65	{
c9fe8530	66	// constructor
10056aa6	67
c9fe8530	68	}
	69
	70	//_____________________________________________________________________________
	71
	72	AliTOFPreprocessor::~AliTOFPreprocessor()
	73	{
	74	// destructor
7981612d	75	if (fData){
	76	delete fData;
	77	fData = 0;
	78	}
	79	if (fh2){
	80	delete fh2;
	81	fh2 = 0;
	82	}
	83	if (fCal){
64f2148a	84	// fCal->Clear();
7981612d	85	delete fCal;
	86	fCal = 0;
	87	}
c9fe8530	88	}
	89
	90	//______________________________________________________________________________
	91	void AliTOFPreprocessor::Initialize(Int_t run, UInt_t startTime,
	92	UInt_t endTime)
	93	{
	94	// Creates AliTOFDataDCS object
	95
	96	AliPreprocessor::Initialize(run, startTime, endTime);
	97
	98	AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", run,
	99	TTimeStamp(startTime).AsString(),
	100	TTimeStamp(endTime).AsString()));
	101
	102	fData = new AliTOFDataDCS(fRun, fStartTime, fEndTime);
708db10b	103	fh2 = 0x0;
10056aa6	104	fNChannels = AliTOFGeometry::NSectors()(2(AliTOFGeometry::NStripC()+AliTOFGeometry::NStripB())+AliTOFGeometry::NStripA())AliTOFGeometry::NpadZ()AliTOFGeometry::NpadX();
	105	fCal = new TObjArray(fNChannels);
	106	fCal->SetOwner();
	107	for (Int_t ich = 0; ich<fNChannels; ich ++){
	108	AliTOFChannelOnline * calChOnline = new AliTOFChannelOnline();
	109	fCal->AddAt(calChOnline,ich);
	110	}
c9fe8530	111	}
	112
	113	//_____________________________________________________________________________
	114
	115	UInt_t AliTOFPreprocessor::Process(TMap* dcsAliasMap)
	116	{
	117	// Fills data into a AliTOFDataDCS object
5936ab02	118	// return codes:
5936ab02	119	// return=0 : all ok
7885c291	120	// return=1 : no DCS input data Map
7885c291	121	// return=2 : no DCS input data processing
b0d8f3fc	122	// return=3 : no DCS processed data was stored in Ref Data
7885c291	123	// return=4 : no DAQ input for Ref Data
b0d8f3fc	124	// return=5 : failed to store Ref Data
7885c291	125	// return=6 : failed to retrieve DAQ data for calibration
	126	// return=7 : problems in histos in the input DAQ file
	127	// return=8 : failed to store Online Delays
c9fe8530	128
708db10b	129	TH1::AddDirectory(0);
5936ab02	130
7885c291	131	Bool_t resultDCSMap=kFALSE;
	132	Bool_t resultDCSStore=kFALSE;
	133	Bool_t resultDAQ=kFALSE;
	134	Bool_t resultDAQRef=kFALSE;
c9fe8530	135
	136	// processing DCS
	137
	138	if (!dcsAliasMap){
5936ab02	139	Log("No DCS map found: TOF exiting from Shuttle");
5936ab02	140	return 1;// return error Code for DCS input data not found
c9fe8530	141	}
	142	else {
	143	// The processing of the DCS input data is forwarded to AliTOFDataDCS
7885c291	144	resultDCSMap=fData->ProcessData(*dcsAliasMap);
	145	if(!resultDCSMap){
	146	Log("Some problems occurred while processing DCS data, TOF exiting from Shuttle");
5936ab02	147	return 2;// return error Code for processed DCS data not stored
c9fe8530	148	}
7885c291	149	else{
	150	AliCDBMetaData metaDataDCS;
	151	metaDataDCS.SetBeamPeriod(0);
	152	metaDataDCS.SetResponsible("Chiara Zampolli");
	153	metaDataDCS.SetComment("This preprocessor fills an AliTOFDataDCS object.");
	154	AliInfo("Storing DCS Data");
b0d8f3fc	155	resultDCSStore = StoreReferenceData("Calib","DCSData",fData, &metaDataDCS);
7885c291	156	if (!resultDCSStore){
b0d8f3fc	157	Log("Some problems occurred while storing DCS data results in Reference Data, TOF exiting from Shuttle");
7885c291	158	return 3;// return error Code for processed DCS data not stored
b0d8f3fc	159	// in reference data
7885c291	160
	161	}
	162	}
c9fe8530	163	}
7885c291	164
c9fe8530	165	// processing DAQ
7885c291	166
c9fe8530	167	TFile * daqFile=0x0;
7885c291	168
5936ab02	169	if(fStoreRefData){
	170	//retrieving data at Run level
	171	TList* list = GetFileSources(kDAQ, "RUNLevel");
	172	if (list)
	173	{
	174	AliInfo("The following sources produced files with the id RUNLevel");
	175	list->Print();
	176	for (Int_t jj=0;jj<list->GetEntries();jj++){
	177	TObjString * str = dynamic_cast<TObjString*> (list->At(jj));
	178	AliInfo(Form("found source %s", str->String().Data()));
	179	// file to be stored run per run
	180	TString fileNameRun = GetFile(kDAQ, "RUNLevel", str->GetName());
	181	if (fileNameRun.Length()>0){
	182	AliInfo(Form("Got the file %s, now we can store the Reference Data for the current Run.", fileNameRun.Data()));
	183	daqFile = new TFile(fileNameRun.Data(),"READ");
	184	fh2 = (TH2S*) daqFile->Get("htof");
	185	AliCDBMetaData metaDataHisto;
	186	metaDataHisto.SetBeamPeriod(0);
	187	metaDataHisto.SetResponsible("Chiara Zampolli");
	188	metaDataHisto.SetComment("This preprocessor stores the array of histos object as Reference Data.");
	189	AliInfo("Storing Reference Data");
	190	resultDAQRef = StoreReferenceData("Calib","DAQData",fh2, &metaDataHisto);
	191	if (!resultDAQRef){
	192	Log("some problems occurred::No Reference Data stored, TOF exiting from Shuttle");
7885c291	193	return 5;//return error code for failure in storing Ref Data
5936ab02	194	}
	195	daqFile->Close();
	196	delete daqFile;
	197	}
	198
	199	else{
	200	Log("The input data file from DAQ (run-level) was not found, TOF exiting from Shuttle ");
7885c291	201	return 4;//return error code for failure in retrieving Ref Data
c9fe8530	202	}
c9fe8530	203	}
9bc469d1	204	}
5936ab02	205	else{
5936ab02	206	Log("The input data file list from DAQ (run-level) was not found, TOF exiting from Shuttle ");
7885c291	207	return 4;//return error code for failure in retrieving Ref Data
5936ab02	208	}
9bc469d1	209	}
9bc469d1	210
9bc469d1	211
5936ab02	212	//Total files, with cumulative histos
5936ab02	213
9bc469d1	214	TList* listTot = GetFileSources(kDAQ, "DELAYS");
	215	if (listTot)
	216	{
	217	AliInfo("The following sources produced files with the id DELAYS");
	218	listTot->Print();
	219	for (Int_t jj=0;jj<listTot->GetEntries();jj++){
	220	TObjString * str = dynamic_cast<TObjString*> (listTot->At(jj));
	221	AliInfo(Form("found source %s", str->String().Data()));
	222
c9fe8530	223	// file with summed histos, to extract calib params
9bc469d1	224	TString fileName = GetFile(kDAQ, "DELAYS", str->GetName());
	225	if (fileName.Length()>0){
	226	AliInfo(Form("Got the file %s, now we can extract some values.", fileName.Data()));
c9fe8530	227
9bc469d1	228	daqFile = new TFile(fileName.Data(),"READ");
64f2148a	229	if (fh2) delete fh2;
708db10b	230	fh2 = (TH2S*) daqFile->Get("htoftot");
708db10b	231	if (!fh2){
9d883ed9	232	Log("some problems occurred:: No histo retrieved, TOF exiting from Shuttle");
9d883ed9	233	delete daqFile;
7885c291	234	return 7; //return error code for histograms not existing/junky
c9fe8530	235	}
c9fe8530	236	else {
3a3ece53	237	static const Int_t kSize=fh2->GetNbinsX();
	238	static const Int_t kNBins=fh2->GetNbinsY();
	239	static const Double_t kXBinmin=fh2->GetYaxis()->GetBinLowEdge(1);
10056aa6	240	if (kSize != fNChannels){
9d883ed9	241	Log(" number of bins along x different from number of pads, found only a subset of the histograms, TOF exiting from Shuttle");
9d883ed9	242	delete daqFile;
7885c291	243	return 7; //return error code for histograms not existing/junky
9bc469d1	244	}
3a3ece53	245	for (Int_t ich=0;ich<kSize;ich++){
	246	TH1S h1 = new TH1S("h1","h1",kNBins,kXBinmin-0.5,kNBins1.+kXBinmin-0.5);
	247	for (Int_t ibin=0;ibin<kNBins;ibin++){
708db10b	248	h1->SetBinContent(ibin+1,fh2->GetBinContent(ich+1,ibin+1));
708db10b	249	}
c9fe8530	250	Bool_t found=kFALSE;
9bc469d1	251	Float_t minContent=h1->Integral()*fgkThrPar;
c9fe8530	252	Int_t nbinsX = h1->GetNbinsX();
	253	Int_t startBin=1;
	254	for (Int_t j=1; j<=nbinsX; j++){
	255	if ((
	256	h1->GetBinContent(j) +
	257	h1->GetBinContent(j+1)+
	258	h1->GetBinContent(j+2)+
	259	h1->GetBinContent(j+3))>minContent){
	260	found=kTRUE;
	261	startBin=j;
	262	break;
	263	}
	264	}
708db10b	265	if(!found) AliInfo(Form("WARNING!!! no start of fit found for histo # %i",ich));
c9fe8530	266	// Now calculate the mean over the interval.
	267	Double_t mean = 0;
	268	Double_t sumw2 = 0;
	269	Double_t nent = 0;
	270	for(Int_t k=0;k<fgkBinRangeAve;k++){
	271	mean=mean+h1->GetBinCenter(startBin+k)*h1->GetBinContent(startBin+k);
	272	nent=nent+h1->GetBinContent(startBin+k);
	273	sumw2=sumw2+(h1->GetBinCenter(startBin+k))(h1->GetBinCenter(startBin+k))(h1->GetBinContent(startBin+k));
	274	}
c9fe8530	275	mean= mean/nent; //<x>
	276	sumw2=sumw2/nent; //<x^2>
	277	Double_t rmsmean= 0;
	278	rmsmean = TMath::Sqrt((sumw2-mean*mean)/nent);
10056aa6	279	if (ich<fNChannels) {
	280	AliTOFChannelOnline * ch = (AliTOFChannelOnline *)fCal->At(ich);
	281	ch->SetDelay((Double_t)meanAliTOFGeometry::TdcBinWidth()1.E-3); // delay in ns
c9fe8530	282	}
708db10b	283	delete h1;
708db10b	284	h1=0x0;
c9fe8530	285	}
c9fe8530	286	}
	287	daqFile->Close();
	288	delete daqFile;
	289	AliCDBMetaData metaData;
	290	metaData.SetBeamPeriod(0);
	291	metaData.SetResponsible("Chiara Zampolli");
	292	metaData.SetComment("This preprocessor fills an AliTOFCal object.");
5936ab02	293	AliInfo("Storing Calibration Data");
a7e9db2a	294	resultDAQ = Store("Calib","ParOnline",fCal, &metaData);
9d883ed9	295	if(!resultDAQ){
9d883ed9	296	Log("Some problems occurred while storing DAQ data processing results");
7885c291	297	return 8;//return error code for problems in storing DAQ data
9d883ed9	298	}
c9fe8530	299	}
c9fe8530	300	else{
9d883ed9	301	Log("The Cumulative data file from DAQ does not exist, TOF exiting from Shuttle");
7885c291	302	return 6;//return error code for problems in retrieving DAQ data
c9fe8530	303	}
	304	}
	305	}
	306	else{
9d883ed9	307	Log("Problem: no list for Cumulative data file from DAQ was found, TOF exiting from Shuttle");
7885c291	308	return 6; //return error code for problems in retrieving DAQ data
c9fe8530	309	}
c9fe8530	310
c9fe8530	311	daqFile=0;
9d883ed9	312
5936ab02	313	return 0;
c9fe8530	314	}
	315
	316