[u/mrichter/AliRoot.git] / TOF / AliTOFDataDCS.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.                  *
 **************************************************************************/

/*
$Log: AliTOFDataDCS.cxx,v $
Revision 1.9  2007/05/04 14:02:45  decaro
AliTOFDataDCS::Draw(Option_t *) method declared const: compiling warning suppression

Revision 1.8  2007/05/03 09:45:09  decaro
Coding convention: RN11 violation -> suppression

Revision 1.7  2007/05/02 14:09:39  arcelli
Retrieval of Env. Temperature removed (will get it from the GRP)

Revision 1.6  2007/04/04 17:19:19  arcelli
Moved some printout to debug level

Revision 1.5  2007/02/20 15:57:00  decaro
Raw data update: to read the TOF raw data defined in UNPACKED mode

Revision 1.4  2007/02/19 15:41:55  decaro
Coding convention: few corrections

Revision 1.3  2007/01/24 11:19:58  arcelli
Modify ProcessData to return a logical (CZ)

Revision 1.2  2006/12/18 18:17:38  arcelli
Updated Aliases for DCS TOF datapoints (C.Zampolli)

Revision 1.1  2006/10/26 09:10:52  arcelli
Class for handling the TOF DCS data in the Shuttle (C.Zampolli)

*/  

// AliTOFDataDCS class
// main aim to introduce the aliases for the TOF DCS
// data points to be then
// stored in the OCDB, and to process them. 
// Process() method called by TOFPrepr

#include "TString.h"
#include "TTimeStamp.h"
#include "TMap.h"
#include "TCanvas.h"

#include "AliDCSValue.h"
#include "AliLog.h"

#include "AliTOFDataDCS.h"
#include "AliTOFFormatDCS.h"

class TH2;
class AliCDBMetaData;
class TDatime;

ClassImp(AliTOFDataDCS)

//---------------------------------------------------------------
AliTOFDataDCS::AliTOFDataDCS():
	TObject(),
	fRun(0),
	fStartTime(0),
	fEndTime(0),
	fStartTimeDCSQuery(0),
	fEndTimeDCSQuery(0),
	fIsProcessed(kFALSE),
	fFDR(kFALSE)
{

  // main constructor 

  for(int i=0;i<kNHV;i++) {
    fHVvp[i]=0x0;
    fHVvn[i]=0x0;
    fHVip[i]=0x0;
    fHVin[i]=0x0;
  }
    
  
}

//---------------------------------------------------------------
AliTOFDataDCS::AliTOFDataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime, UInt_t startTimeDCSQuery, UInt_t endTimeDCSQuery):
	TObject(),
	fRun(nRun),
	fStartTime(startTime),
	fEndTime(endTime),
	fStartTimeDCSQuery(startTimeDCSQuery),
	fEndTimeDCSQuery(endTimeDCSQuery),
	fIsProcessed(kFALSE),
	fFDR(kFALSE)
{

  // constructor with arguments

	AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s \n\tStartTime DCS Query %s \n\tEndTime DCS Query %s", nRun,
	TTimeStamp(startTime).AsString(),
	TTimeStamp(endTime).AsString(), 
	TTimeStamp(startTimeDCSQuery).AsString(), 
        TTimeStamp(endTimeDCSQuery).AsString()));

	Init();

}

//---------------------------------------------------------------

AliTOFDataDCS::AliTOFDataDCS(const AliTOFDataDCS & data):
  TObject(data), 
  fRun(data.fRun),
  fStartTime(data.fStartTime),
  fEndTime(data.fEndTime),
  fStartTimeDCSQuery(data.fStartTimeDCSQuery),
  fEndTimeDCSQuery(data.fEndTimeDCSQuery),
  fIsProcessed(data.fIsProcessed),
  fFDR(data.fFDR)

{

  // copy constructor

  for(int i=0;i<kNAliases;i++) {
    fAliasNames[i]=data.fAliasNames[i];
  }
 
  for(int i=0;i<kNHV;i++) {
    fHVvp[i]=data.fHVvp[i];
    fHVvn[i]=data.fHVvn[i];
    fHVip[i]=data.fHVip[i];
    fHVin[i]=data.fHVin[i];
  }
  
    
}
//---------------------------------------------------------------

AliTOFDataDCS& AliTOFDataDCS:: operator=(const AliTOFDataDCS & data) { 

  // assignment operator

  if (this == &data)
    return *this;

  TObject::operator=(data);
  fRun=data.GetRun();
  fStartTime=data.GetStartTime();
  fEndTime=data.GetEndTime();
  fStartTimeDCSQuery=data.GetStartTimeDCSQuery();
  fEndTimeDCSQuery=data.GetEndTimeDCSQuery();

  for(int i=0;i<kNAliases;i++) {
    fAliasNames[i]=data.GetAliasName(i);
  }

  for(int i=0;i<kNHV;i++) {
    fHVvp[i]=data.GetHVvp(i);
    fHVvn[i]=data.GetHVvn(i);
    fHVip[i]=data.GetHVip(i);
    fHVin[i]=data.GetHVin(i);
  }


  return *this;
}
//---------------------------------------------------------------
AliTOFDataDCS::~AliTOFDataDCS() {

  // destructor

  for(int i=0;i<kNHV;i++) {
    delete fHVvp[i];
    fHVvp[i]=0;
    delete fHVvn[i];
    fHVvn[i]=0;
    delete fHVip[i];
    fHVip[i]=0;
    delete fHVin[i];
    fHVin[i]=0;
  }
  
}

//---------------------------------------------------------------
Bool_t AliTOFDataDCS::ProcessData(TMap& aliasMap){

  // method to process the data

  if(!(fAliasNames[0])) Init();

  Float_t val=0;
  Float_t val0=0;
  Float_t val1=0;
  Float_t time=0; 
  Float_t time0=0; 
  Float_t time1=0; 
  Float_t delta[2];
  Float_t timedelta[2];

  AliInfo(Form(" Start Time = %i",fStartTime));
  AliInfo(Form(" End Time = %i",fEndTime));
  AliInfo(Form(" Start Time DCS Query= %i",fStartTimeDCSQuery));
  AliInfo(Form(" End Time DCS Query= %i",fEndTimeDCSQuery));

  if (fEndTime==fStartTime){
    AliError(Form(" Run with null time length: start time = %i = end time = %i",fStartTime,fEndTime));
    return kFALSE;
  }

  TObjArray *aliasArr;
  AliDCSValue* aValue;
  AliDCSValue* aValue0;
  AliDCSValue* aValue1;

  // starting loop on aliases
  for(int j=0; j<kNAliases; j++){
    for (Int_t k=0;k<2;k++) {
      delta[k]=0;
      timedelta[k]=0;
    }
    //AliInfo(Form("j = %i, with alias = %s",j,fAliasNames[j].Data()));
    aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[j].Data());
    if(!aliasArr){
      AliError(Form("Alias %s not found!", fAliasNames[j].Data()));
      if (!fFDR){
	return kFALSE;    // returning only in case we are not in a FDR run
      }
      else {
	continue;
      }
    }

    Introduce(j, aliasArr);
    
    if(aliasArr->GetEntries()==0){
      AliError(Form("Alias %s has no entries! Nothing will be stored",
		    fAliasNames[j].Data()));
      continue;
    }
    if(aliasArr->GetEntries() == 1){
      AliWarning(Form("Alias %s has just %d entries! Only first value will be stored",
		    fAliasNames[j].Data(),aliasArr->GetEntries()));
      aValue0 = (AliDCSValue*) aliasArr->At(0);
      val0 = aValue0->GetFloat();
      time0 = (Float_t) (aValue0->GetTimeStamp());
      if (j<kNHV){
	fHVvp[j]->SetFloat(0,val0);
	fHVvp[j]->SetTimeStampFloat(0,time0);
      }
      else if (j<kNHV*2){
	fHVvn[j-kNHV]->SetFloat(0,val0);
	fHVvn[j-kNHV]->SetTimeStampFloat(0,time0);
      }
      else if (j<kNHV*3){
	fHVip[j-2*kNHV]->SetFloat(0,val0);
	fHVip[j-2*kNHV]->SetTimeStampFloat(0,time0);
      }
      else if (j<kNHV*4){
	fHVin[j-3*kNHV]->SetFloat(0,val0);
	fHVin[j-3*kNHV]->SetTimeStampFloat(0,time0);
      }
      continue;
    }
    if(aliasArr->GetEntries()==2){
      AliWarning(Form("Alias %s has just %d entries! Storing first and second entry only and the delta between these two",
		    fAliasNames[j].Data(),aliasArr->GetEntries()));
      aValue0 = (AliDCSValue*) aliasArr->At(0);
      val0 = aValue0->GetFloat();
      time0 = (Float_t) (aValue0->GetTimeStamp());
      aValue1 = (AliDCSValue*) aliasArr->At(1);
      val1 = aValue1->GetFloat();
      time1 = (Float_t) (aValue1->GetTimeStamp());
      if (j<kNHV){
	fHVvp[j]->SetFloat(0,val0);
	fHVvp[j]->SetTimeStampFloat(0,time0);
	fHVvp[j]->SetFloat(1,val1);
	fHVvp[j]->SetTimeStampFloat(1,time1);
	fHVvp[j]->SetDelta(0,TMath::Abs(val1-val0));
	fHVvp[j]->SetTimeStampDelta(0,time1-time0);
      }
      else if (j<kNHV*2){
	fHVvn[j-kNHV]->SetFloat(0,val0);
	fHVvn[j-kNHV]->SetTimeStampFloat(0,time0);
	fHVvn[j-kNHV]->SetFloat(1,val1);
	fHVvn[j-kNHV]->SetTimeStampFloat(1,time1);
	fHVvn[j-kNHV]->SetDelta(0,TMath::Abs(val1-val0));
	fHVvn[j-kNHV]->SetTimeStampDelta(0,time1-time0);
      }
      else if (j<kNHV*3){
	fHVip[j-2*kNHV]->SetFloat(0,val0);
	fHVip[j-2*kNHV]->SetTimeStampFloat(0,time0);
	fHVip[j-2*kNHV]->SetFloat(1,val1);
	fHVip[j-2*kNHV]->SetTimeStampFloat(1,time1);
	fHVip[j-2*kNHV]->SetDelta(0,TMath::Abs(val1-val0));
	fHVip[j-2*kNHV]->SetTimeStampDelta(0,time1-time0);
      }
      else if (j<kNHV*4){
	fHVin[j-3*kNHV]->SetFloat(0,val0);
	fHVin[j-3*kNHV]->SetTimeStampFloat(0,time0);
	fHVin[j-3*kNHV]->SetFloat(1,val1);
	fHVin[j-3*kNHV]->SetTimeStampFloat(1,time1);
	fHVin[j-3*kNHV]->SetDelta(0,TMath::Abs(val1-val0));
	fHVin[j-3*kNHV]->SetTimeStampDelta(0,time1-time0);
      }
      continue;
    }
    
    TIter iterarray(aliasArr);
    
    Int_t nentries = aliasArr->GetEntries();
    Int_t deltaTimeStamp = (Int_t) nentries/3;
    Int_t deltaTimeStamp1 = (Int_t) nentries/2;

    // filling aliases with 10 floats+1 Usign
    Int_t index = 0;
    for (Int_t k=0;k<3;k++){
      index = deltaTimeStamp*k;
      if (k==0) {
	index=0;
      }
      else if (k==1) {
	index=deltaTimeStamp1;
      } 
      else if (k==2) {
	index=nentries-1; 
      }
      aValue = (AliDCSValue*) aliasArr->At(index);
      val = aValue->GetFloat();
      time = (Float_t) (aValue->GetTimeStamp());
      if (j<kNHV){
	fHVvp[j]->SetFloat(k,val);
	fHVvp[j]->SetTimeStampFloat(k,time);
      }
      else if (j<kNHV*2){
	fHVvn[j-kNHV]->SetFloat(k,val);
	fHVvn[j-kNHV]->SetTimeStampFloat(k,time);
      }
      else if (j<kNHV*3){
	fHVip[j-2*kNHV]->SetFloat(k,val);
	fHVip[j-2*kNHV]->SetTimeStampFloat(k,time);
      }
      else if (j<kNHV*4){
	fHVin[j-3*kNHV]->SetFloat(k,val);
	fHVin[j-3*kNHV]->SetTimeStampFloat(k,time);
      }
    }
  
    //computing the most significant variations

    //Float_t timeDiff = (Float_t)(fEndTime-fStartTime);
    Float_t timeDiff = (Float_t)(fEndTimeDCSQuery-fStartTimeDCSQuery);
    Int_t deltamin = (Int_t)(60/timeDiff*nentries); //sampling every minute
    Int_t klast = nentries-deltamin;
      
    for (Int_t k=0;k<klast;k++){
      aValue = (AliDCSValue*) aliasArr->At(k);
      aValue1 = (AliDCSValue*) aliasArr->At(k+deltamin);
      val = aValue->GetFloat();
      val1 = aValue1->GetFloat();
      if (delta[0]<=TMath::Abs(val1-val)) {
	delta[0]=TMath::Abs(val1-val);
	timedelta[0] = (Float_t)k;
      }
      if (delta[1]<=delta[0]) {
	Float_t temp = delta[1];
	Float_t timetemp = timedelta[1];
	delta[1]=delta[0];
	delta[0]=temp;
	timedelta[1]=timedelta[0];
	timedelta[0]=timetemp;
      }
    }
    
    for (Int_t kk=0;kk<2;kk++){
      if (j<kNHV){
	fHVvp[j]->SetDelta(kk,delta[kk]);
	fHVvp[j]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
      }
      else if (j<kNHV*2){
	fHVvn[j-kNHV]->SetDelta(kk,delta[kk]);
	fHVvn[j-kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
      }
      else if (j<kNHV*3){
	fHVip[j-2*kNHV]->SetDelta(kk,delta[kk]);
	fHVip[j-2*kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
      }
      else if (j<kNHV*4){
	fHVin[j-3*kNHV]->SetDelta(kk,delta[kk]);
	fHVin[j-3*kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
      }
    }
  }
  
  fIsProcessed=kTRUE;

  return kTRUE;
}

//---------------------------------------------------------------
void AliTOFDataDCS::Init(){

  // initialization of aliases and DCS data

  TString sindex;
  for(int i=0;i<kNAliases;i++){
    //HV, v
    if (i<kNHV){
	fAliasNames[i] = "tof_hv_vp_";
	sindex.Form("%02i",i);
	fAliasNames[i] += sindex;
	fHVvp[i] = new AliTOFFormatDCS();
    }
    else if (i<kNHV*2){
	fAliasNames[i] = "tof_hv_vn_";
	sindex.Form("%02i",i-kNHV);
	fAliasNames[i] += sindex;
	fHVvn[i-kNHV] = new AliTOFFormatDCS();
    }
    //HV, i
    else if (i<kNHV*3){
	fAliasNames[i] = "tof_hv_ip_";
	sindex.Form("%02i",i-2*kNHV);
	fAliasNames[i] += sindex;
	fHVip[i-2*kNHV] = new AliTOFFormatDCS();
    }
    else if (i<kNHV*4){
	fAliasNames[i] = "tof_hv_in_";
	sindex.Form("%02i",i-3*kNHV);
	fAliasNames[i] += sindex;
	fHVin[i-3*kNHV] = new AliTOFFormatDCS();
    }
  }
}


//---------------------------------------------------------------
void AliTOFDataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)const
{

  // method to introduce new aliases

  int entries=0;
  entries = aliasArr->GetEntries();
  int nal=0;
  nal=numAlias;
  AliDebug(2,Form("************ Alias: %s **********",fAliasNames[numAlias].Data()));
  AliDebug(2,Form("    	%d DP values collected",entries));

}

//---------------------------------------------------------------
void AliTOFDataDCS::Draw(const Option_t* /*option*/)
{
  // Draw all histos and graphs

  if(!fIsProcessed) return;

  TCanvas *ch;
  TString canvasHistoName="Histos";
  ch=new TCanvas(canvasHistoName,canvasHistoName,20,20,600,600);
  ch->cd();

  // to be implemented

}
Commit	Line	Data
345932e4	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	/*
7fffa85b	17	$Log: AliTOFDataDCS.cxx,v $
	18	Revision 1.9 2007/05/04 14:02:45 decaro
	19	AliTOFDataDCS::Draw(Option_t *) method declared const: compiling warning suppression
	20
b52b97c2	21	Revision 1.8 2007/05/03 09:45:09 decaro
	22	Coding convention: RN11 violation -> suppression
	23
9d5fd6fe	24	Revision 1.7 2007/05/02 14:09:39 arcelli
	25	Retrieval of Env. Temperature removed (will get it from the GRP)
	26
38e3bbf1	27	Revision 1.6 2007/04/04 17:19:19 arcelli
	28	Moved some printout to debug level
	29
72bf252d	30	Revision 1.5 2007/02/20 15:57:00 decaro
	31	Raw data update: to read the TOF raw data defined in UNPACKED mode
	32
15ec34b9	33	Revision 1.4 2007/02/19 15:41:55 decaro
	34	Coding convention: few corrections
	35
3a3ece53	36	Revision 1.3 2007/01/24 11:19:58 arcelli
	37	Modify ProcessData to return a logical (CZ)
	38
ebab534a	39	Revision 1.2 2006/12/18 18:17:38 arcelli
	40	Updated Aliases for DCS TOF datapoints (C.Zampolli)
	41
8259e6f8	42	Revision 1.1 2006/10/26 09:10:52 arcelli
	43	Class for handling the TOF DCS data in the Shuttle (C.Zampolli)
	44
345932e4	45	*/
345932e4	46
051a7fe6	47	// AliTOFDataDCS class
	48	// main aim to introduce the aliases for the TOF DCS
	49	// data points to be then
	50	// stored in the OCDB, and to process them.
	51	// Process() method called by TOFPrepr
	52
5c7c93fa	53	#include "TString.h"
5c7c93fa	54	#include "TTimeStamp.h"
	55	#include "TMap.h"
	56	#include "TCanvas.h"
345932e4	57
	58	#include "AliDCSValue.h"
	59	#include "AliLog.h"
	60
5c7c93fa	61	#include "AliTOFDataDCS.h"
345932e4	62	#include "AliTOFFormatDCS.h"
345932e4	63
	64	class TH2;
	65	class AliCDBMetaData;
	66	class TDatime;
	67
345932e4	68	ClassImp(AliTOFDataDCS)
	69
	70	//---------------------------------------------------------------
	71	AliTOFDataDCS::AliTOFDataDCS():
	72	TObject(),
	73	fRun(0),
	74	fStartTime(0),
	75	fEndTime(0),
ff326f84	76	fStartTimeDCSQuery(0),
ff326f84	77	fEndTimeDCSQuery(0),
7fffa85b	78	fIsProcessed(kFALSE),
7fffa85b	79	fFDR(kFALSE)
345932e4	80	{
	81
	82	// main constructor
	83
8259e6f8	84	for(int i=0;i<kNHV;i++) {
	85	fHVvp[i]=0x0;
	86	fHVvn[i]=0x0;
	87	fHVip[i]=0x0;
	88	fHVin[i]=0x0;
	89	}
d202b31a	90
8259e6f8	91
345932e4	92	}
	93
	94	//---------------------------------------------------------------
ff326f84	95	AliTOFDataDCS::AliTOFDataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime, UInt_t startTimeDCSQuery, UInt_t endTimeDCSQuery):
345932e4	96	TObject(),
	97	fRun(nRun),
	98	fStartTime(startTime),
	99	fEndTime(endTime),
ff326f84	100	fStartTimeDCSQuery(startTimeDCSQuery),
ff326f84	101	fEndTimeDCSQuery(endTimeDCSQuery),
7fffa85b	102	fIsProcessed(kFALSE),
7fffa85b	103	fFDR(kFALSE)
345932e4	104	{
	105
	106	// constructor with arguments
	107
ff326f84	108	AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s \n\tStartTime DCS Query %s \n\tEndTime DCS Query %s", nRun,
345932e4	109	TTimeStamp(startTime).AsString(),
ff326f84	110	TTimeStamp(endTime).AsString(),
	111	TTimeStamp(startTimeDCSQuery).AsString(),
	112	TTimeStamp(endTimeDCSQuery).AsString()));
345932e4	113
	114	Init();
	115
	116	}
	117
	118	//---------------------------------------------------------------
	119
	120	AliTOFDataDCS::AliTOFDataDCS(const AliTOFDataDCS & data):
8a190ba2	121	TObject(data),
	122	fRun(data.fRun),
	123	fStartTime(data.fStartTime),
	124	fEndTime(data.fEndTime),
	125	fStartTimeDCSQuery(data.fStartTimeDCSQuery),
	126	fEndTimeDCSQuery(data.fEndTimeDCSQuery),
	127	fIsProcessed(data.fIsProcessed),
	128	fFDR(data.fFDR)
345932e4	129
	130	{
	131
051a7fe6	132	// copy constructor
345932e4	133
345932e4	134	for(int i=0;i<kNAliases;i++) {
	135	fAliasNames[i]=data.fAliasNames[i];
	136	}
	137
	138	for(int i=0;i<kNHV;i++) {
8259e6f8	139	fHVvp[i]=data.fHVvp[i];
	140	fHVvn[i]=data.fHVvn[i];
	141	fHVip[i]=data.fHVip[i];
	142	fHVin[i]=data.fHVin[i];
345932e4	143	}
345932e4	144
38e3bbf1	145
345932e4	146	}
	147	//---------------------------------------------------------------
	148
	149	AliTOFDataDCS& AliTOFDataDCS:: operator=(const AliTOFDataDCS & data) {
	150
051a7fe6	151	// assignment operator
345932e4	152
8a190ba2	153	if (this == &data)
	154	return *this;
	155
	156	TObject::operator=(data);
	157	fRun=data.GetRun();
	158	fStartTime=data.GetStartTime();
	159	fEndTime=data.GetEndTime();
	160	fStartTimeDCSQuery=data.GetStartTimeDCSQuery();
	161	fEndTimeDCSQuery=data.GetEndTimeDCSQuery();
345932e4	162
345932e4	163	for(int i=0;i<kNAliases;i++) {
8a190ba2	164	fAliasNames[i]=data.GetAliasName(i);
345932e4	165	}
345932e4	166
345932e4	167	for(int i=0;i<kNHV;i++) {
8a190ba2	168	fHVvp[i]=data.GetHVvp(i);
	169	fHVvn[i]=data.GetHVvn(i);
	170	fHVip[i]=data.GetHVip(i);
	171	fHVin[i]=data.GetHVin(i);
345932e4	172	}
345932e4	173
345932e4	174
	175	return *this;
	176	}
	177	//---------------------------------------------------------------
	178	AliTOFDataDCS::~AliTOFDataDCS() {
	179
	180	// destructor
	181
8259e6f8	182	for(int i=0;i<kNHV;i++) {
	183	delete fHVvp[i];
	184	fHVvp[i]=0;
	185	delete fHVvn[i];
	186	fHVvn[i]=0;
	187	delete fHVip[i];
	188	fHVip[i]=0;
	189	delete fHVin[i];
	190	fHVin[i]=0;
	191	}
	192
345932e4	193	}
345932e4	194
345932e4	195	//---------------------------------------------------------------
ebab534a	196	Bool_t AliTOFDataDCS::ProcessData(TMap& aliasMap){
345932e4	197
3a3ece53	198	// method to process the data
3a3ece53	199
345932e4	200	if(!(fAliasNames[0])) Init();
345932e4	201
345932e4	202	Float_t val=0;
5f2d4eea	203	Float_t val0=0;
345932e4	204	Float_t val1=0;
345932e4	205	Float_t time=0;
5f2d4eea	206	Float_t time0=0;
5f2d4eea	207	Float_t time1=0;
345932e4	208	Float_t delta[2];
	209	Float_t timedelta[2];
	210
e83a0dd0	211	AliInfo(Form(" Start Time = %i",fStartTime));
e83a0dd0	212	AliInfo(Form(" End Time = %i",fEndTime));
ff326f84	213	AliInfo(Form(" Start Time DCS Query= %i",fStartTimeDCSQuery));
ff326f84	214	AliInfo(Form(" End Time DCS Query= %i",fEndTimeDCSQuery));
e83a0dd0	215
	216	if (fEndTime==fStartTime){
	217	AliError(Form(" Run with null time length: start time = %i = end time = %i",fStartTime,fEndTime));
	218	return kFALSE;
96fec15b	219	}
96fec15b	220
345932e4	221	TObjArray *aliasArr;
345932e4	222	AliDCSValue* aValue;
5f2d4eea	223	AliDCSValue* aValue0;
345932e4	224	AliDCSValue* aValue1;
345932e4	225
	226	// starting loop on aliases
	227	for(int j=0; j<kNAliases; j++){
	228	for (Int_t k=0;k<2;k++) {
	229	delta[k]=0;
	230	timedelta[k]=0;
	231	}
	232	//AliInfo(Form("j = %i, with alias = %s",j,fAliasNames[j].Data()));
	233	aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[j].Data());
	234	if(!aliasArr){
	235	AliError(Form("Alias %s not found!", fAliasNames[j].Data()));
7fffa85b	236	if (!fFDR){
	237	return kFALSE; // returning only in case we are not in a FDR run
	238	}
	239	else {
	240	continue;
	241	}
345932e4	242	}
	243
	244	Introduce(j, aliasArr);
	245
5f2d4eea	246	if(aliasArr->GetEntries()==0){
	247	AliError(Form("Alias %s has no entries! Nothing will be stored",
	248	fAliasNames[j].Data()));
	249	continue;
	250	}
	251	if(aliasArr->GetEntries() == 1){
	252	AliWarning(Form("Alias %s has just %d entries! Only first value will be stored",
345932e4	253	fAliasNames[j].Data(),aliasArr->GetEntries()));
5f2d4eea	254	aValue0 = (AliDCSValue*) aliasArr->At(0);
	255	val0 = aValue0->GetFloat();
	256	time0 = (Float_t) (aValue0->GetTimeStamp());
	257	if (j<kNHV){
	258	fHVvp[j]->SetFloat(0,val0);
	259	fHVvp[j]->SetTimeStampFloat(0,time0);
	260	}
	261	else if (j<kNHV*2){
	262	fHVvn[j-kNHV]->SetFloat(0,val0);
	263	fHVvn[j-kNHV]->SetTimeStampFloat(0,time0);
	264	}
	265	else if (j<kNHV*3){
	266	fHVip[j-2*kNHV]->SetFloat(0,val0);
	267	fHVip[j-2*kNHV]->SetTimeStampFloat(0,time0);
	268	}
	269	else if (j<kNHV*4){
	270	fHVin[j-3*kNHV]->SetFloat(0,val0);
	271	fHVin[j-3*kNHV]->SetTimeStampFloat(0,time0);
	272	}
	273	continue;
	274	}
	275	if(aliasArr->GetEntries()==2){
	276	AliWarning(Form("Alias %s has just %d entries! Storing first and second entry only and the delta between these two",
	277	fAliasNames[j].Data(),aliasArr->GetEntries()));
	278	aValue0 = (AliDCSValue*) aliasArr->At(0);
	279	val0 = aValue0->GetFloat();
	280	time0 = (Float_t) (aValue0->GetTimeStamp());
	281	aValue1 = (AliDCSValue*) aliasArr->At(1);
	282	val1 = aValue1->GetFloat();
	283	time1 = (Float_t) (aValue1->GetTimeStamp());
	284	if (j<kNHV){
	285	fHVvp[j]->SetFloat(0,val0);
	286	fHVvp[j]->SetTimeStampFloat(0,time0);
	287	fHVvp[j]->SetFloat(1,val1);
	288	fHVvp[j]->SetTimeStampFloat(1,time1);
	289	fHVvp[j]->SetDelta(0,TMath::Abs(val1-val0));
	290	fHVvp[j]->SetTimeStampDelta(0,time1-time0);
	291	}
	292	else if (j<kNHV*2){
	293	fHVvn[j-kNHV]->SetFloat(0,val0);
	294	fHVvn[j-kNHV]->SetTimeStampFloat(0,time0);
	295	fHVvn[j-kNHV]->SetFloat(1,val1);
	296	fHVvn[j-kNHV]->SetTimeStampFloat(1,time1);
	297	fHVvn[j-kNHV]->SetDelta(0,TMath::Abs(val1-val0));
	298	fHVvn[j-kNHV]->SetTimeStampDelta(0,time1-time0);
	299	}
	300	else if (j<kNHV*3){
	301	fHVip[j-2*kNHV]->SetFloat(0,val0);
	302	fHVip[j-2*kNHV]->SetTimeStampFloat(0,time0);
	303	fHVip[j-2*kNHV]->SetFloat(1,val1);
	304	fHVip[j-2*kNHV]->SetTimeStampFloat(1,time1);
	305	fHVip[j-2*kNHV]->SetDelta(0,TMath::Abs(val1-val0));
	306	fHVip[j-2*kNHV]->SetTimeStampDelta(0,time1-time0);
	307	}
	308	else if (j<kNHV*4){
	309	fHVin[j-3*kNHV]->SetFloat(0,val0);
	310	fHVin[j-3*kNHV]->SetTimeStampFloat(0,time0);
	311	fHVin[j-3*kNHV]->SetFloat(1,val1);
	312	fHVin[j-3*kNHV]->SetTimeStampFloat(1,time1);
	313	fHVin[j-3*kNHV]->SetDelta(0,TMath::Abs(val1-val0));
	314	fHVin[j-3*kNHV]->SetTimeStampDelta(0,time1-time0);
	315	}
345932e4	316	continue;
	317	}
	318
	319	TIter iterarray(aliasArr);
	320
	321	Int_t nentries = aliasArr->GetEntries();
345932e4	322	Int_t deltaTimeStamp = (Int_t) nentries/3;
345932e4	323	Int_t deltaTimeStamp1 = (Int_t) nentries/2;
345932e4	324
345932e4	325	// filling aliases with 10 floats+1 Usign
38e3bbf1	326	Int_t index = 0;
	327	for (Int_t k=0;k<3;k++){
	328	index = deltaTimeStamp*k;
	329	if (k==0) {
	330	index=0;
345932e4	331	}
38e3bbf1	332	else if (k==1) {
	333	index=deltaTimeStamp1;
	334	}
	335	else if (k==2) {
	336	index=nentries-1;
345932e4	337	}
38e3bbf1	338	aValue = (AliDCSValue*) aliasArr->At(index);
	339	val = aValue->GetFloat();
	340	time = (Float_t) (aValue->GetTimeStamp());
	341	if (j<kNHV){
	342	fHVvp[j]->SetFloat(k,val);
	343	fHVvp[j]->SetTimeStampFloat(k,time);
	344	}
	345	else if (j<kNHV*2){
	346	fHVvn[j-kNHV]->SetFloat(k,val);
	347	fHVvn[j-kNHV]->SetTimeStampFloat(k,time);
	348	}
	349	else if (j<kNHV*3){
	350	fHVip[j-2*kNHV]->SetFloat(k,val);
	351	fHVip[j-2*kNHV]->SetTimeStampFloat(k,time);
	352	}
	353	else if (j<kNHV*4){
	354	fHVin[j-3*kNHV]->SetFloat(k,val);
	355	fHVin[j-3*kNHV]->SetTimeStampFloat(k,time);
	356	}
345932e4	357	}
38e3bbf1	358
345932e4	359	//computing the most significant variations
345932e4	360
ff326f84	361	//Float_t timeDiff = (Float_t)(fEndTime-fStartTime);
ff326f84	362	Float_t timeDiff = (Float_t)(fEndTimeDCSQuery-fStartTimeDCSQuery);
e83a0dd0	363	Int_t deltamin = (Int_t)(60/timeDiff*nentries); //sampling every minute
345932e4	364	Int_t klast = nentries-deltamin;
38e3bbf1	365
345932e4	366	for (Int_t k=0;k<klast;k++){
	367	aValue = (AliDCSValue*) aliasArr->At(k);
	368	aValue1 = (AliDCSValue*) aliasArr->At(k+deltamin);
	369	val = aValue->GetFloat();
	370	val1 = aValue1->GetFloat();
	371	if (delta[0]<=TMath::Abs(val1-val)) {
	372	delta[0]=TMath::Abs(val1-val);
	373	timedelta[0] = (Float_t)k;
	374	}
	375	if (delta[1]<=delta[0]) {
	376	Float_t temp = delta[1];
	377	Float_t timetemp = timedelta[1];
	378	delta[1]=delta[0];
	379	delta[0]=temp;
	380	timedelta[1]=timedelta[0];
	381	timedelta[0]=timetemp;
	382	}
	383	}
	384
	385	for (Int_t kk=0;kk<2;kk++){
38e3bbf1	386	if (j<kNHV){
	387	fHVvp[j]->SetDelta(kk,delta[kk]);
	388	fHVvp[j]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
345932e4	389	}
38e3bbf1	390	else if (j<kNHV*2){
	391	fHVvn[j-kNHV]->SetDelta(kk,delta[kk]);
	392	fHVvn[j-kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
	393	}
	394	else if (j<kNHV*3){
	395	fHVip[j-2*kNHV]->SetDelta(kk,delta[kk]);
	396	fHVip[j-2*kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
	397	}
	398	else if (j<kNHV*4){
	399	fHVin[j-3*kNHV]->SetDelta(kk,delta[kk]);
	400	fHVin[j-3*kNHV]->SetTimeStampDelta(kk,(Float_t)timedelta[kk]);
	401	}
8259e6f8	402	}
345932e4	403	}
38e3bbf1	404
345932e4	405	fIsProcessed=kTRUE;
345932e4	406
ebab534a	407	return kTRUE;
345932e4	408	}
	409
	410	//---------------------------------------------------------------
	411	void AliTOFDataDCS::Init(){
	412
	413	// initialization of aliases and DCS data
	414
8259e6f8	415	TString sindex;
345932e4	416	for(int i=0;i<kNAliases;i++){
8259e6f8	417	//HV, v
345932e4	418	if (i<kNHV){
8259e6f8	419	fAliasNames[i] = "tof_hv_vp_";
	420	sindex.Form("%02i",i);
	421	fAliasNames[i] += sindex;
	422	fHVvp[i] = new AliTOFFormatDCS();
345932e4	423	}
345932e4	424	else if (i<kNHV*2){
8259e6f8	425	fAliasNames[i] = "tof_hv_vn_";
	426	sindex.Form("%02i",i-kNHV);
	427	fAliasNames[i] += sindex;
	428	fHVvn[i-kNHV] = new AliTOFFormatDCS();
345932e4	429	}
8259e6f8	430	//HV, i
345932e4	431	else if (i<kNHV*3){
8259e6f8	432	fAliasNames[i] = "tof_hv_ip_";
	433	sindex.Form("%02i",i-2*kNHV);
	434	fAliasNames[i] += sindex;
	435	fHVip[i-2*kNHV] = new AliTOFFormatDCS();
345932e4	436	}
345932e4	437	else if (i<kNHV*4){
8259e6f8	438	fAliasNames[i] = "tof_hv_in_";
	439	sindex.Form("%02i",i-3*kNHV);
	440	fAliasNames[i] += sindex;
	441	fHVin[i-3*kNHV] = new AliTOFFormatDCS();
345932e4	442	}
345932e4	443	}
	444	}
	445
	446
	447	//---------------------------------------------------------------
	448	void AliTOFDataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)const
	449	{
	450
	451	// method to introduce new aliases
	452
	453	int entries=0;
	454	entries = aliasArr->GetEntries();
	455	int nal=0;
	456	nal=numAlias;
72bf252d	457	AliDebug(2,Form("********** Alias: %s ********",fAliasNames[numAlias].Data()));
72bf252d	458	AliDebug(2,Form(" %d DP values collected",entries));
345932e4	459
	460	}
	461
	462	//---------------------------------------------------------------
b52b97c2	463	void AliTOFDataDCS::Draw(const Option_t* /option/)
345932e4	464	{
051a7fe6	465	// Draw all histos and graphs
345932e4	466
	467	if(!fIsProcessed) return;
	468
	469	TCanvas *ch;
	470	TString canvasHistoName="Histos";
	471	ch=new TCanvas(canvasHistoName,canvasHistoName,20,20,600,600);
	472	ch->cd();
	473
	474	// to be implemented
	475
	476	}
	477