[u/mrichter/AliRoot.git] / ACORDE / AliACORDEDataDCS.cxx

/****************************************************
	
  AliACORDEDataDCS class
  create to make a pointer to the
  ACORDE data DCS points

  Author: Pedro Gonzalez (CIEMAT, Madrid)
  ACORDE-DCS creator: Mario Ivan Martinez Hdez
			<mim@fcfm.buap.mx>

  Last update: June 13th 2014 from Mario Rodríguez Cahuantzi (CINVESTAV, mrodriguez@fis.cinvestav.mx)
	==> wrong ACORDE aliases for AMANDA fixed (https://alice.its.cern.ch/jira/browse/ALIROOT-5479)
		-- old alias: ACO_HV_MODULE[0..59]_VMON
		-- new alias: ACO_HV_MODULE[0..59]_INSIDE_VMON ()/ ACO_HV_MODULE[0..59]_OUTSIDE_VMON

  Last update: Fix of coding violations
  Mario Rodriguez C. (FCFM-BUAP)
  <mrodrigu@mail.cern.ch>

*****************************************************/
#include "AliACORDEDataDCS.h"

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

#include <TTimeStamp.h>
#include <TObjString.h>
#include <TH2F.h>
#include <TProfile.h>
#include <TGraph.h>
#include <TDatime.h>
#include <TStyle.h>
#include <TCanvas.h>

ClassImp(AliACORDEDataDCS)

//---------------------------------------------------------------
AliACORDEDataDCS::AliACORDEDataDCS():
	TObject(),
	fRun(0),
	fStartTime(0),
	fEndTime(0),
	fGraphs("TGraph",kNGraphs),
        fFunc(0),
	fIsProcessed(kFALSE)
{
	for(int i=0;i<kNHistos;i++) 
	{
		fHv[i]=0x0;
		fMean[i] = fWidth[i] = 0.0;
	}
        
}

//---------------------------------------------------------------
AliACORDEDataDCS::AliACORDEDataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime):
	TObject(),
	fRun(nRun),
	fStartTime(startTime),
	fEndTime(endTime),
	fGraphs("TGraph",kNGraphs),
        fFunc(0),
	fIsProcessed(kFALSE)
{
// Init of class AliACORDEDataDCS
// Gettin the TimeStamp an put it on a string

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

       
	Init();

}

//---------------------------------------------------------------
AliACORDEDataDCS::~AliACORDEDataDCS() {

	for(int i=0;i<kNHistos;i++) 
	{
		delete fHv[i]; fHv[i]=0;
	}
	fGraphs.Clear("C");
	fFunc=0;
}
//---------------------------------------------------------------

AliACORDEDataDCS::AliACORDEDataDCS(const AliACORDEDataDCS & data):
TObject(),
fRun(0),
fStartTime(0),
fEndTime(0),
fGraphs("TGraph",kNGraphs),
fFunc(0),
fIsProcessed(kFALSE)
{
// Setting the initial values
// fRUn, Start of Run, End of Run, IsProcessed

	fRun=data.fRun;
	fStartTime=data.fStartTime;
	fEndTime=data.fEndTime;
	fFunc=data.fFunc;
	fIsProcessed=data.fIsProcessed;


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

        for(int i=0;i<kNHistos;i++)
	{
		fHv[i]=data.fHv[i];
		fMean[i]=data.fMean[i];
		fWidth[i]=data.fWidth[i];
	}


}
//--------------------------------------------------------------
AliACORDEDataDCS& AliACORDEDataDCS:: operator=(const AliACORDEDataDCS & data) { 

	// use ctor copy

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

  	new (this) AliACORDEDataDCS(data);
  	return *this;
 
}
//---------------------------------------------------------------
void AliACORDEDataDCS::ProcessData(TMap& aliasMap)
{
// Process of the data from the aliases DCS-data points

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

	TObjArray *aliasArr;
	AliDCSValue* aValue;

	for(int j=0; j<kNAliases; j++)
        {
		aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[j].Data());
		if(!aliasArr)
                {
			AliError(Form("Alias %s not found!", fAliasNames[j].Data()));
			continue;
		}
		Introduce(j, aliasArr);

		if(aliasArr->GetEntries()<2)
                {
		AliError(Form("Alias %s has just %d entries!",
					fAliasNames[j].Data(),aliasArr->GetEntries()));
			continue;
		}

		TIter iterarray(aliasArr);

		Double_t *time = new Double_t[aliasArr->GetEntries()];
		Double_t *val  = new Double_t[aliasArr->GetEntries()];

		UInt_t ne=0;

		while ((aValue = (AliDCSValue*) iterarray.Next())) 
                {
                val[ne] = aValue->GetFloat();
		time[ne] = (Double_t) (aValue->GetTimeStamp());
		fHv[j]->Fill(val[ne]);
		ne++;
		}

               
		CreateGraph(j, aliasArr->GetEntries(), time, val);
		delete[] val;
		delete[] time;
	}


	// calculate mean and rms of the first two histos
	for(int i=0;i<kNHistos;i++)
        {
		fMean[i] = fHv[i]->GetMean();
		fWidth[i] = fHv[i]->GetRMS();
	}


	fIsProcessed=kTRUE;


}

//---------------------------------------------------------------
void AliACORDEDataDCS::Init()
{
// Init of AliACORDEDatDCS procedure
// Loop over the aliases

	TH1::AddDirectory(kFALSE);

	fGraphs.SetOwner(1);

        TString aliasName;

	for(int i=0;i<kNAliases;i++){
		if (i<kNAliases/2)
		{
                	aliasName.Form("ACO_HV_MODULE%02d_INSIDE_VMON",i);  //!::: first 60 entries for inside PMT's
		}else
		{
	               	aliasName.Form("ACO_HV_MODULE%02d_OUTSIDE_VMON",i-60); //!::: last 60 entries for outside PMT's
		}
		fAliasNames[i] = aliasName;
	}

	for(int i=0;i<kNHistos;i++)
        {
		fHv[i] = new TH1F(fAliasNames[i].Data(),fAliasNames[i].Data(), 20, kHvMin, kHvMax);
		fHv[i]->GetXaxis()->SetTitle("Hv");
	}
}

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

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

}

//---------------------------------------------------------------
void AliACORDEDataDCS::CreateGraph(int i, int dim, const Double_t *x, const Double_t *y)
{
// Create the plots for the ACORDE DCS 

	TGraph *gr = new(fGraphs[fGraphs.GetEntriesFast()]) TGraph(dim, x, y);

	gr->GetXaxis()->SetTimeDisplay(1);
	gr->SetTitle(fAliasNames[i].Data());

	AliInfo(Form("Array entries: %d",fGraphs.GetEntriesFast()));


}

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

  if(!fIsProcessed) return;

  TString canvasHistoName;
  TCanvas *ch[10];
  
  for (int i=0;i<10;i++)
  {
  canvasHistoName.Form("ACO_HV_MODULE");
  ch[i]=new TCanvas(canvasHistoName,canvasHistoName,20,20,600,600);
  ch[i]->Divide(2,3);

    for(int j=0;j<6;j++)
    { 
    ch[i]->cd(j+1);
    ((TGraph*) fGraphs.UncheckedAt(i*6+j))->SetMarkerStyle(20);
    ((TGraph*) fGraphs.UncheckedAt(i*6+j))->Draw("alp");
    }

  }

 
}
Commit	Line	Data
3b181530	1	/****************************************************
	2
	3	AliACORDEDataDCS class
	4	create to make a pointer to the
	5	ACORDE data DCS points
	6
	7	Author: Pedro Gonzalez (CIEMAT, Madrid)
	8	ACORDE-DCS creator: Mario Ivan Martinez Hdez
	9	<mim@fcfm.buap.mx>
	10
b3b7be1c	11	Last update: June 13th 2014 from Mario Rodríguez Cahuantzi (CINVESTAV, mrodriguez@fis.cinvestav.mx)
	12	==> wrong ACORDE aliases for AMANDA fixed (https://alice.its.cern.ch/jira/browse/ALIROOT-5479)
	13	-- old alias: ACO_HV_MODULE[0..59]_VMON
	14	-- new alias: ACO_HV_MODULE[0..59]_INSIDE_VMON ()/ ACO_HV_MODULE[0..59]_OUTSIDE_VMON
	15
3b181530	16	Last update: Fix of coding violations
	17	Mario Rodriguez C. (FCFM-BUAP)
	18	<mrodrigu@mail.cern.ch>
	19
	20	*****************************************************/
749347ef	21	#include "AliACORDEDataDCS.h"
	22
	23	#include "AliCDBMetaData.h"
	24	#include "AliDCSValue.h"
	25	#include "AliLog.h"
	26
	27	#include <TTimeStamp.h>
	28	#include <TObjString.h>
	29	#include <TH2F.h>
	30	#include <TProfile.h>
	31	#include <TGraph.h>
	32	#include <TDatime.h>
	33	#include <TStyle.h>
	34	#include <TCanvas.h>
	35
	36	ClassImp(AliACORDEDataDCS)
	37
	38	//---------------------------------------------------------------
	39	AliACORDEDataDCS::AliACORDEDataDCS():
	40	TObject(),
	41	fRun(0),
	42	fStartTime(0),
	43	fEndTime(0),
	44	fGraphs("TGraph",kNGraphs),
5d8a7578	45	fFunc(0),
749347ef	46	fIsProcessed(kFALSE)
749347ef	47	{
ea3dd546	48	for(int i=0;i<kNHistos;i++)
	49	{
	50	fHv[i]=0x0;
	51	fMean[i] = fWidth[i] = 0.0;
	52	}
5d8a7578	53
749347ef	54	}
	55
	56	//---------------------------------------------------------------
	57	AliACORDEDataDCS::AliACORDEDataDCS(Int_t nRun, UInt_t startTime, UInt_t endTime):
	58	TObject(),
	59	fRun(nRun),
	60	fStartTime(startTime),
	61	fEndTime(endTime),
	62	fGraphs("TGraph",kNGraphs),
5d8a7578	63	fFunc(0),
749347ef	64	fIsProcessed(kFALSE)
749347ef	65	{
3b181530	66	// Init of class AliACORDEDataDCS
	67	// Gettin the TimeStamp an put it on a string
	68
749347ef	69	AliInfo(Form("\n\tRun %d \n\tStartTime %s \n\tEndTime %s", nRun,
	70	TTimeStamp(startTime).AsString(),
	71	TTimeStamp(endTime).AsString()));
	72
5d8a7578	73
749347ef	74	Init();
	75
	76	}
	77
	78	//---------------------------------------------------------------
	79	AliACORDEDataDCS::~AliACORDEDataDCS() {
	80
392813bf	81	for(int i=0;i<kNHistos;i++)
	82	{
	83	delete fHv[i]; fHv[i]=0;
	84	}
749347ef	85	fGraphs.Clear("C");
	86	fFunc=0;
	87	}
5d8a7578	88	//---------------------------------------------------------------
	89
	90	AliACORDEDataDCS::AliACORDEDataDCS(const AliACORDEDataDCS & data):
	91	TObject(),
	92	fRun(0),
	93	fStartTime(0),
	94	fEndTime(0),
	95	fGraphs("TGraph",kNGraphs),
	96	fFunc(0),
	97	fIsProcessed(kFALSE)
	98	{
3b181530	99	// Setting the initial values
3b181530	100	// fRUn, Start of Run, End of Run, IsProcessed
5d8a7578	101
	102	fRun=data.fRun;
	103	fStartTime=data.fStartTime;
	104	fEndTime=data.fEndTime;
	105	fFunc=data.fFunc;
	106	fIsProcessed=data.fIsProcessed;
	107
	108
	109	for(int i=0;i<kNAliases;i++){fAliasNames[i] = data.fAliasNames[i];}
	110
392813bf	111	for(int i=0;i<kNHistos;i++)
	112	{
	113	fHv[i]=data.fHv[i];
	114	fMean[i]=data.fMean[i];
	115	fWidth[i]=data.fWidth[i];
	116	}
5d8a7578	117
749347ef	118
5d8a7578	119
	120
	121	}
	122	//--------------------------------------------------------------
	123	AliACORDEDataDCS& AliACORDEDataDCS:: operator=(const AliACORDEDataDCS & data) {
	124
81f0f392	125	// use ctor copy
5d8a7578	126
81f0f392	127	if (&data == this) return *this;
5d8a7578	128
81f0f392	129	new (this) AliACORDEDataDCS(data);
81f0f392	130	return *this;
5d8a7578	131
5d8a7578	132	}
749347ef	133	//---------------------------------------------------------------
3b181530	134	void AliACORDEDataDCS::ProcessData(TMap& aliasMap)
	135	{
	136	// Process of the data from the aliases DCS-data points
749347ef	137
	138	if(!(fHv[0])) Init();
	139
	140	TObjArray *aliasArr;
	141	AliDCSValue* aValue;
	142
	143	for(int j=0; j<kNAliases; j++)
	144	{
	145	aliasArr = (TObjArray*) aliasMap.GetValue(fAliasNames[j].Data());
	146	if(!aliasArr)
	147	{
	148	AliError(Form("Alias %s not found!", fAliasNames[j].Data()));
	149	continue;
	150	}
	151	Introduce(j, aliasArr);
	152
	153	if(aliasArr->GetEntries()<2)
	154	{
ea3dd546	155	AliError(Form("Alias %s has just %d entries!",
749347ef	156	fAliasNames[j].Data(),aliasArr->GetEntries()));
	157	continue;
	158	}
	159
	160	TIter iterarray(aliasArr);
	161
	162	Double_t *time = new Double_t[aliasArr->GetEntries()];
	163	Double_t *val = new Double_t[aliasArr->GetEntries()];
	164
	165	UInt_t ne=0;
	166
	167	while ((aValue = (AliDCSValue*) iterarray.Next()))
	168	{
	169	val[ne] = aValue->GetFloat();
	170	time[ne] = (Double_t) (aValue->GetTimeStamp());
	171	fHv[j]->Fill(val[ne]);
	172	ne++;
	173	}
	174
	175
	176
	177	CreateGraph(j, aliasArr->GetEntries(), time, val);
	178	delete[] val;
	179	delete[] time;
	180	}
	181
	182
	183	// calculate mean and rms of the first two histos
	184	for(int i=0;i<kNHistos;i++)
	185	{
	186	fMean[i] = fHv[i]->GetMean();
	187	fWidth[i] = fHv[i]->GetRMS();
	188	}
	189
	190
	191	fIsProcessed=kTRUE;
	192
	193
	194	}
	195
	196	//---------------------------------------------------------------
3b181530	197	void AliACORDEDataDCS::Init()
	198	{
	199	// Init of AliACORDEDatDCS procedure
	200	// Loop over the aliases
749347ef	201
	202	TH1::AddDirectory(kFALSE);
	203
	204	fGraphs.SetOwner(1);
	205
	206	TString aliasName;
	207
	208	for(int i=0;i<kNAliases;i++){
b3b7be1c	209	if (i<kNAliases/2)
	210	{
	211	aliasName.Form("ACO_HV_MODULE%02d_INSIDE_VMON",i); //!::: first 60 entries for inside PMT's
	212	}else
	213	{
	214	aliasName.Form("ACO_HV_MODULE%02d_OUTSIDE_VMON",i-60); //!::: last 60 entries for outside PMT's
	215	}
749347ef	216	fAliasNames[i] = aliasName;
	217	}
	218
	219	for(int i=0;i<kNHistos;i++)
	220	{
	221	fHv[i] = new TH1F(fAliasNames[i].Data(),fAliasNames[i].Data(), 20, kHvMin, kHvMax);
	222	fHv[i]->GetXaxis()->SetTitle("Hv");
	223	}
	224	}
	225
	226	//---------------------------------------------------------------
	227	void AliACORDEDataDCS::Introduce(UInt_t numAlias, const TObjArray* aliasArr)
	228	{
	229
	230	int entries=aliasArr->GetEntries();
	231	AliInfo(Form("********** Alias: %s ********",fAliasNames[numAlias].Data()));
	232	AliInfo(Form(" %d DP values collected",entries));
	233
	234	}
	235
	236	//---------------------------------------------------------------
	237	void AliACORDEDataDCS::CreateGraph(int i, int dim, const Double_t x, const Double_t y)
	238	{
3b181530	239	// Create the plots for the ACORDE DCS
749347ef	240
	241	TGraph *gr = new(fGraphs[fGraphs.GetEntriesFast()]) TGraph(dim, x, y);
	242
	243	gr->GetXaxis()->SetTimeDisplay(1);
	244	gr->SetTitle(fAliasNames[i].Data());
	245
	246	AliInfo(Form("Array entries: %d",fGraphs.GetEntriesFast()));
	247
	248
	249	}
	250
	251	//---------------------------------------------------------------
	252	void AliACORDEDataDCS::Draw(const Option_t* /option/)
	253	{
	254	// Draw all histos and graphs
	255
	256	if(!fIsProcessed) return;
	257
	258	TString canvasHistoName;
	259	TCanvas *ch[10];
	260
	261	for (int i=0;i<10;i++)
	262	{
	263	canvasHistoName.Form("ACO_HV_MODULE");
	264	ch[i]=new TCanvas(canvasHistoName,canvasHistoName,20,20,600,600);
	265	ch[i]->Divide(2,3);
	266
	267	for(int j=0;j<6;j++)
	268	{
	269	ch[i]->cd(j+1);
	270	((TGraph) fGraphs.UncheckedAt(i6+j))->SetMarkerStyle(20);
	271	((TGraph) fGraphs.UncheckedAt(i6+j))->Draw("alp");
	272	}
	273
	274	}
	275
	276
	277	}
	278