[u/mrichter/AliRoot.git] / TPC / AliTPCQAChecker.cxx

/**************************************************************************
 * Copyright(c) 1998-2007, 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: $ */

/*
  Checks implemented a la AliMUONQAChecker.
  Checks the quality assurance by realzed checks on histogram content.
  P. Christiansen, Lund, September 2009.

  Based on AliPHOSQAChecker.
  Checks the quality assurance by comparing with reference data.
  P. Christiansen, Lund, January 2008.
*/

// --- ROOT header files ---
#include <TH1.h>

// --- AliRoot header files ---
#include "AliTPCQAChecker.h"
#include "AliTPCQADataMakerRec.h"

ClassImp(AliTPCQAChecker)

//__________________________________________________________________
Double_t *
AliTPCQAChecker::Check(AliQAv1::ALITASK_t index, TObjArray ** list, 
		       const AliDetectorRecoParam * /*recoParam*/)
{
  /* It is important to understand the destinction between indexed tasks (AliQAv1::TASKINDEX_t) which are used in the DataMaker classes and indexed tasks (AliQAv1::ALITASK_t) whihc are used in the checker class.

     From the AliQAChecker::Run() methods we have:
     AliQAv1::kRAW
     - AliQAv1::kRAWS 
     
     AliQAv1::kSIM 
     - AliQAv1::kHITS
     - AliQAv1::kSDIGITS
     - AliQAv1::kDIGITS
     
     AliQAv1::kREC
     - AliQAv1::kDIGITSR 
     - AliQAv1::kRECPOINTS
     - AliQAv1::kTRACKSEGMENTS 
     - AliQAv1::kRECPARTICLES
     
     AliQAv1::kESD ; 
     -AliQAv1::kESDS

     This means that for each group of tasks the Check will be called
     one or more times.  This also mean that we cannot know what
     histograms will be or not be there in a single call... And we
     also do not know the position in the list of the histogram.
  */
  
  /// Check objects in list
  if(fDebug>0)
    AliInfo("In AliTPCQAChecker::Check");
  
  if (index!=AliQAv1::kRAW&&index!=AliQAv1::kREC&&index!=AliQAv1::kESD) {
    
    AliWarning(Form("Checker for task %d not implement for the moment",index));
    return NULL;
  }
  
  Double_t * rv = new Double_t[AliRecoParam::kNSpecies] ; 
  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) 
    rv[specie] = 1.0; // All is fine 
  
  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
    
    if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) ) 
      continue ; 
    
    if (index == AliQAv1::kRAW)
      rv[specie] = CheckRAW(specie, list[specie]);
    if (index == AliQAv1::kREC)
      rv[specie] = CheckREC(specie, list[specie]);
    if (index == AliQAv1::kESD)
      rv[specie] = CheckESD(specie, list[specie]);

    if(fDebug>3)
      AliInfo(Form("Specie: %s. Task: %s. Value: %f",  
		   AliRecoParam::GetEventSpecieName(specie),
		   AliQAv1::GetAliTaskName(index),
		   rv[specie]));
  }
  
  return rv;
}

//______________________________________________________________________________
Double_t AliTPCQAChecker::CheckRAW(Int_t specie, TObjArray* list)
{
  /// Check ESD
  if(fDebug>0)
    AliInfo("In AliTPCQAChecker::CheckRAW");
  
  if(fDebug>2)
    list->Print();

  Char_t specieName[256];
  sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));

  TH1* hRawsOccupancyVsSector = static_cast<TH1*>
    (list->FindObject(Form("%s_hRawsOccupancyVsSector",specieName)));
  TH1* hRawsQmaxVsSector = static_cast<TH1*>
    (list->FindObject(Form("%s_hRawsQmaxVsSector",specieName)));
  TH1* hRawsOccupancyVsEvent = static_cast<TH1*>
    (list->FindObject(Form("%s_hRawsOccupancyVsEvent",specieName)));
  TH1* hRawsNclustersVsEvent = static_cast<TH1*>
    (list->FindObject(Form("%s_hRawsNclustersVsEvent",specieName)));
  
  if (!hRawsOccupancyVsSector || !hRawsQmaxVsSector ||
      !hRawsOccupancyVsEvent  || !hRawsNclustersVsEvent) 
    return -0.5; // fatal
  
  if(hRawsOccupancyVsSector->GetEntries()==0) {
    return 0.25; // error - No TPC data!
  }
  
  Int_t nBinsX = hRawsOccupancyVsSector->GetNbinsX();
  for(Int_t i = 1; i <= nBinsX; i++) {
    
    if(hRawsOccupancyVsSector->GetBinContent(i)==0)
      return 0.75; // warning - no TPC data for at least one sector
  }
  
  return 1.0; // ok
}

//______________________________________________________________________________
Double_t AliTPCQAChecker::CheckREC(Int_t specie, TObjArray* list)
{
  // This method checks the QA histograms associated with reconstruction
  //
  // For TPC this is:
  // DigitsR : 
  // The digit histogram gives the ADC distribution for all sigbnals
  // above threshold. The check is just that there are digits.
  // RecPoints : 
  // The cluster histograms are meant to give an idea about the gain
  // from the cluster charge and to indicate iof there are rows with
  // noise clusters, i.e., they are very visual.
  //
  // The check methods are simple:
  // If there are no histogram at all (set fatal)
  // If digit histogram is there, but there are no digits (set error)
  // If cluster histogram is there but there are less than 1000 
  //    clusters (set warning)
  // If there are more than 1000 clusters but no clusters for either short, 
  // medium, or long pads (set error)
  if(fDebug>0)
    AliInfo("In AliTPCQAChecker::CheckREC");
  
  if(fDebug>2)
    list->Print();

  Char_t specieName[256];
  sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));

  TH1* hDigits = static_cast<TH1*>
    (list->FindObject(Form("%s_hDigitsADC",specieName)));
  TH1* hNclustersVsRow = static_cast<TH1*> 
    (list->FindObject(Form("%s_hRecPointsRow",specieName)));
  TH1* hQshort         = static_cast<TH1*> 
    (list->FindObject(Form("%s_hRecPointsQShort",specieName)));
  TH1* hQmedium        = static_cast<TH1*> 
    (list->FindObject(Form("%s_hRecPointsQMedium",specieName)));
  TH1* hQlong          = static_cast<TH1*> 
    (list->FindObject(Form("%s_hRecPointsQLong",specieName)));
  // The Qmax histograms are for now ignored

  if (!hDigits &&                                             // digits missing
      (!hNclustersVsRow || !hQshort || !hQmedium || !hQlong)) // 1 recpoint hist missing
    return -0.5; // fatal

  if (hDigits && hDigits->GetEntries()==0) 
    return 0.25; // error
  
  if (hNclustersVsRow && hNclustersVsRow->GetEntries() < 1000) {
    return 0.75; // warning
  } else { 
    if (!hQshort || !hQlong || !hQlong)
      return -0.5;// fatal - they should be there if the cluster vs row hist is there
    if (hQshort->GetEntries()==0 || hQmedium->GetEntries()==0 ||
	hQlong->GetEntries()==0) 
      return 0.25; // error
  }
  return 1; // ok
}

//______________________________________________________________________________
Double_t AliTPCQAChecker::CheckESD(Int_t specie, TObjArray* list)
{
  // This method checks the QA histograms associated with ESDs
  // (Note that there is aslo a globalQA which is running on all
  //  the ESD information so for now this is just a few basic
  //  histograms)
  //
  // The check methods are simple:
  // If there are no histogram at all (set fatal)
  // 
  if(fDebug>0)
    AliInfo("In AliTPCQAChecker::CheckESD");

  if(fDebug>2)
    list->Print();

  Char_t specieName[256];
  sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));

  TH1* hESDclusters = static_cast<TH1*>
    (list->FindObject(Form("%s_hESDclusters",specieName)));
  TH1* hESDratio = static_cast<TH1*>
    (list->FindObject(Form("%s_hESDratio",specieName)));
  TH1* hESDpt = static_cast<TH1*>
    (list->FindObject(Form("%s_hESDpt",specieName)));
  
  if (!hESDclusters || !hESDratio || !hESDpt) 
    return -0.5; // fatal

  return 1.0; // ok
}

//______________________________________________________________________________
void AliTPCQAChecker::Init(const AliQAv1::DETECTORINDEX_t det) 
{
  /// intialises QA and QA checker settings
  if(fDebug>0)
    AliInfo("In AliTPCQAChecker::Init");
  AliQAv1::Instance(det) ; 
  Float_t hiValue[AliQAv1::kNBIT] ; 
  Float_t lowValue[AliQAv1::kNBIT] ;
  hiValue[AliQAv1::kINFO]       = 1.00; 
  lowValue[AliQAv1::kINFO]      = 0.99; 
  hiValue[AliQAv1::kWARNING]    = 0.99; 
  lowValue[AliQAv1::kWARNING]   = 0.50; 
  hiValue[AliQAv1::kERROR]      = 0.50; 
  lowValue[AliQAv1::kERROR]     = 0.00; 
  hiValue[AliQAv1::kFATAL]      = 0.00; 
  lowValue[AliQAv1::kFATAL]     =-1.00; 
  //  SetHiLo(&hiValue[0], &lowValue[0]) ; 
  SetHiLo(hiValue, lowValue) ; 
}

//______________________________________________________________________________
void 
AliTPCQAChecker::SetQA(AliQAv1::ALITASK_t index, Double_t * value) const
{
  /// sets the QA according the return value of the Check

  if(fDebug>0)
    AliInfo("In AliTPCQAChecker::SetQA");

  AliQAv1 * qa = AliQAv1::Instance(index);
  
  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {

    if (value==NULL) { // No checker is implemented, set all QA to Fatal
      
      if(fDebug>1)
	AliInfo(Form("Fatal QA. Task: %s. Specie: %s", 
		     AliQAv1::GetAliTaskName(index), 
		     AliRecoParam::GetEventSpecieName(specie)));
      qa->Set(AliQAv1::kFATAL, specie) ; 
    } else {
      
      if ( value[specie] >= fLowTestValue[AliQAv1::kFATAL] && 
	   value[specie] < fUpTestValue[AliQAv1::kFATAL] ) {
	
	if(fDebug>1)
	  AliInfo(Form("QA-Fatal. Task: %s. Specie: %s", 
		       AliQAv1::GetAliTaskName(index), 
		       AliRecoParam::GetEventSpecieName(specie)));
        qa->Set(AliQAv1::kFATAL, specie) ; 
      } else if ( value[specie] > fLowTestValue[AliQAv1::kERROR] && 
		  value[specie] <= fUpTestValue[AliQAv1::kERROR]  ) {
	
	if(fDebug>1)
	  AliInfo(Form("QA-Error. Task: %s. Specie: %s", 
		       AliQAv1::GetAliTaskName(index), 
		       AliRecoParam::GetEventSpecieName(specie)));
        qa->Set(AliQAv1::kERROR, specie) ; 
      } else if (value[specie] > fLowTestValue[AliQAv1::kWARNING] && 
		 value[specie] <= fUpTestValue[AliQAv1::kWARNING]) {
	
	if(fDebug>1)
	  AliInfo(Form("QA-Warning. Task: %s. Specie: %s", 
		       AliQAv1::GetAliTaskName(index), 
		       AliRecoParam::GetEventSpecieName(specie)));
	qa->Set(AliQAv1::kWARNING, specie) ;
      } else if (value[specie] > fLowTestValue[AliQAv1::kINFO] && 
		 value[specie] <= fUpTestValue[AliQAv1::kINFO] ) { 
	
	if(fDebug>1)
	  AliInfo(Form("QA-Info. Task: %s. Specie: %s", 
		       AliQAv1::GetAliTaskName(index), 
		       AliRecoParam::GetEventSpecieName(specie)));
	qa->Set(AliQAv1::kINFO, specie) ; 	
      }
    }
  }
}
Commit	Line	Data
44f32dd2	1	/**************************************************************************
	2	* Copyright(c) 1998-2007, 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	/*
ce0175fa	19	Checks implemented a la AliMUONQAChecker.
	20	Checks the quality assurance by realzed checks on histogram content.
	21	P. Christiansen, Lund, September 2009.
	22
44f32dd2	23	Based on AliPHOSQAChecker.
	24	Checks the quality assurance by comparing with reference data.
	25	P. Christiansen, Lund, January 2008.
	26	*/
	27
ce0175fa	28	// --- ROOT header files ---
	29	#include <TH1.h>
	30
44f32dd2	31	// --- AliRoot header files ---
44f32dd2	32	#include "AliTPCQAChecker.h"
ce0175fa	33	#include "AliTPCQADataMakerRec.h"
44f32dd2	34
	35	ClassImp(AliTPCQAChecker)
	36
	37	//__________________________________________________________________
ce0175fa	38	Double_t *
ce0175fa	39	AliTPCQAChecker::Check(AliQAv1::ALITASK_t index, TObjArray ** list,
486788fc	40	const AliDetectorRecoParam * /recoParam/)
ce0175fa	41	{
	42	/* It is important to understand the destinction between indexed tasks (AliQAv1::TASKINDEX_t) which are used in the DataMaker classes and indexed tasks (AliQAv1::ALITASK_t) whihc are used in the checker class.
	43
	44	From the AliQAChecker::Run() methods we have:
	45	AliQAv1::kRAW
	46	- AliQAv1::kRAWS
	47
	48	AliQAv1::kSIM
	49	- AliQAv1::kHITS
	50	- AliQAv1::kSDIGITS
	51	- AliQAv1::kDIGITS
	52
	53	AliQAv1::kREC
	54	- AliQAv1::kDIGITSR
	55	- AliQAv1::kRECPOINTS
	56	- AliQAv1::kTRACKSEGMENTS
	57	- AliQAv1::kRECPARTICLES
	58
	59	AliQAv1::kESD ;
	60	-AliQAv1::kESDS
	61
	62	This means that for each group of tasks the Check will be called
	63	one or more times. This also mean that we cannot know what
	64	histograms will be or not be there in a single call... And we
	65	also do not know the position in the list of the histogram.
	66	*/
	67
	68	/// Check objects in list
	69	if(fDebug>0)
	70	AliInfo("In AliTPCQAChecker::Check");
	71
	72	if (index!=AliQAv1::kRAW&&index!=AliQAv1::kREC&&index!=AliQAv1::kESD) {
	73
	74	AliWarning(Form("Checker for task %d not implement for the moment",index));
	75	return NULL;
	76	}
	77
	78	Double_t * rv = new Double_t[AliRecoParam::kNSpecies] ;
	79	for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
	80	rv[specie] = 1.0; // All is fine
	81
	82	for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
	83
	84	if ( !AliQAv1::Instance()->IsEventSpecieSet(specie) )
	85	continue ;
	86
	87	if (index == AliQAv1::kRAW)
	88	rv[specie] = CheckRAW(specie, list[specie]);
	89	if (index == AliQAv1::kREC)
	90	rv[specie] = CheckREC(specie, list[specie]);
	91	if (index == AliQAv1::kESD)
	92	rv[specie] = CheckESD(specie, list[specie]);
	93
	94	if(fDebug>3)
	95	AliInfo(Form("Specie: %s. Task: %s. Value: %f",
	96	AliRecoParam::GetEventSpecieName(specie),
	97	AliQAv1::GetAliTaskName(index),
	98	rv[specie]));
	99	}
	100
	101	return rv;
	102	}
	103
	104	//______________________________________________________________________________
105	Double_t AliTPCQAChecker::CheckRAW(Int_t specie, TObjArray* list)
106	{
107	/// Check ESD
108	if(fDebug>0)
109	AliInfo("In AliTPCQAChecker::CheckRAW");
110
111	if(fDebug>2)
112	list->Print();
113
114	Char_t specieName[256];
115	sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));
116
117	TH1* hRawsOccupancyVsSector = static_cast<TH1*>
118	(list->FindObject(Form("%s_hRawsOccupancyVsSector",specieName)));
119	TH1* hRawsQmaxVsSector = static_cast<TH1*>
120	(list->FindObject(Form("%s_hRawsQmaxVsSector",specieName)));
121	TH1* hRawsOccupancyVsEvent = static_cast<TH1*>
122	(list->FindObject(Form("%s_hRawsOccupancyVsEvent",specieName)));
123	TH1* hRawsNclustersVsEvent = static_cast<TH1*>
124	(list->FindObject(Form("%s_hRawsNclustersVsEvent",specieName)));
125
126	if (!hRawsOccupancyVsSector \|\| !hRawsQmaxVsSector \|\|
127	!hRawsOccupancyVsEvent \|\| !hRawsNclustersVsEvent)
128	return -0.5; // fatal
129
130	if(hRawsOccupancyVsSector->GetEntries()==0) {
131	return 0.25; // error - No TPC data!
132	}
133
134	Int_t nBinsX = hRawsOccupancyVsSector->GetNbinsX();
135	for(Int_t i = 1; i <= nBinsX; i++) {
136
137	if(hRawsOccupancyVsSector->GetBinContent(i)==0)
138	return 0.75; // warning - no TPC data for at least one sector
139	}
140
141	return 1.0; // ok
142	}
143
144	//______________________________________________________________________________
145	Double_t AliTPCQAChecker::CheckREC(Int_t specie, TObjArray* list)
146	{
147	// This method checks the QA histograms associated with reconstruction
148	//
149	// For TPC this is:
150	// DigitsR :
151	// The digit histogram gives the ADC distribution for all sigbnals
152	// above threshold. The check is just that there are digits.
153	// RecPoints :
154	// The cluster histograms are meant to give an idea about the gain
155	// from the cluster charge and to indicate iof there are rows with
156	// noise clusters, i.e., they are very visual.
157	//
158	// The check methods are simple:
159	// If there are no histogram at all (set fatal)
160	// If digit histogram is there, but there are no digits (set error)
161	// If cluster histogram is there but there are less than 1000
162	// clusters (set warning)
163	// If there are more than 1000 clusters but no clusters for either short,
164	// medium, or long pads (set error)
165	if(fDebug>0)
166	AliInfo("In AliTPCQAChecker::CheckREC");
167
168	if(fDebug>2)
169	list->Print();
170
171	Char_t specieName[256];
172	sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));
173
174	TH1* hDigits = static_cast<TH1*>
175	(list->FindObject(Form("%s_hDigitsADC",specieName)));
176	TH1* hNclustersVsRow = static_cast<TH1*>
177	(list->FindObject(Form("%s_hRecPointsRow",specieName)));
178	TH1* hQshort = static_cast<TH1*>
179	(list->FindObject(Form("%s_hRecPointsQShort",specieName)));
180	TH1* hQmedium = static_cast<TH1*>
181	(list->FindObject(Form("%s_hRecPointsQMedium",specieName)));
182	TH1* hQlong = static_cast<TH1*>
183	(list->FindObject(Form("%s_hRecPointsQLong",specieName)));
184	// The Qmax histograms are for now ignored
185
186	if (!hDigits && // digits missing
187	(!hNclustersVsRow \|\| !hQshort \|\| !hQmedium \|\| !hQlong)) // 1 recpoint hist missing
188	return -0.5; // fatal
189
190	if (hDigits && hDigits->GetEntries()==0)
191	return 0.25; // error
192
193	if (hNclustersVsRow && hNclustersVsRow->GetEntries() < 1000) {
194	return 0.75; // warning
195	} else {
196	if (!hQshort \|\| !hQlong \|\| !hQlong)
197	return -0.5;// fatal - they should be there if the cluster vs row hist is there
198	if (hQshort->GetEntries()==0 \|\| hQmedium->GetEntries()==0 \|\|
199	hQlong->GetEntries()==0)
200	return 0.25; // error
201	}
202	return 1; // ok
203	}
204
205	//______________________________________________________________________________
206	Double_t AliTPCQAChecker::CheckESD(Int_t specie, TObjArray* list)
207	{
208	// This method checks the QA histograms associated with ESDs
209	// (Note that there is aslo a globalQA which is running on all
210	// the ESD information so for now this is just a few basic
211	// histograms)
212	//
213	// The check methods are simple:
214	// If there are no histogram at all (set fatal)
215	//
216	if(fDebug>0)
217	AliInfo("In AliTPCQAChecker::CheckESD");
218
219	if(fDebug>2)
220	list->Print();
221
222	Char_t specieName[256];
223	sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));
224
225	TH1* hESDclusters = static_cast<TH1*>
226	(list->FindObject(Form("%s_hESDclusters",specieName)));
227	TH1* hESDratio = static_cast<TH1*>
228	(list->FindObject(Form("%s_hESDratio",specieName)));
229	TH1* hESDpt = static_cast<TH1*>
230	(list->FindObject(Form("%s_hESDpt",specieName)));
231
232	if (!hESDclusters \|\| !hESDratio \|\| !hESDpt)
233	return -0.5; // fatal
234
235	return 1.0; // ok
236	}
237
238	//______________________________________________________________________________
239	void AliTPCQAChecker::Init(const AliQAv1::DETECTORINDEX_t det)
240	{
241	/// intialises QA and QA checker settings
242	if(fDebug>0)
243	AliInfo("In AliTPCQAChecker::Init");
244	AliQAv1::Instance(det) ;
245	Float_t hiValue[AliQAv1::kNBIT] ;
246	Float_t lowValue[AliQAv1::kNBIT] ;
247	hiValue[AliQAv1::kINFO] = 1.00;
248	lowValue[AliQAv1::kINFO] = 0.99;
249	hiValue[AliQAv1::kWARNING] = 0.99;
250	lowValue[AliQAv1::kWARNING] = 0.50;
251	hiValue[AliQAv1::kERROR] = 0.50;
252	lowValue[AliQAv1::kERROR] = 0.00;
253	hiValue[AliQAv1::kFATAL] = 0.00;
254	lowValue[AliQAv1::kFATAL] =-1.00;
255	// SetHiLo(&hiValue[0], &lowValue[0]) ;
256	SetHiLo(hiValue, lowValue) ;
257	}
258
259	//______________________________________________________________________________
260	void
261	AliTPCQAChecker::SetQA(AliQAv1::ALITASK_t index, Double_t * value) const
262	{
263	/// sets the QA according the return value of the Check
264
265	if(fDebug>0)
266	AliInfo("In AliTPCQAChecker::SetQA");
267
268	AliQAv1 * qa = AliQAv1::Instance(index);
269
270	for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
44f32dd2	271
ce0175fa	272	if (value==NULL) { // No checker is implemented, set all QA to Fatal
	273
	274	if(fDebug>1)
	275	AliInfo(Form("Fatal QA. Task: %s. Specie: %s",
	276	AliQAv1::GetAliTaskName(index),
	277	AliRecoParam::GetEventSpecieName(specie)));
	278	qa->Set(AliQAv1::kFATAL, specie) ;
	279	} else {
	280
	281	if ( value[specie] >= fLowTestValue[AliQAv1::kFATAL] &&
	282	value[specie] < fUpTestValue[AliQAv1::kFATAL] ) {
	283
	284	if(fDebug>1)
	285	AliInfo(Form("QA-Fatal. Task: %s. Specie: %s",
	286	AliQAv1::GetAliTaskName(index),
	287	AliRecoParam::GetEventSpecieName(specie)));
	288	qa->Set(AliQAv1::kFATAL, specie) ;
	289	} else if ( value[specie] > fLowTestValue[AliQAv1::kERROR] &&
	290	value[specie] <= fUpTestValue[AliQAv1::kERROR] ) {
	291
	292	if(fDebug>1)
	293	AliInfo(Form("QA-Error. Task: %s. Specie: %s",
	294	AliQAv1::GetAliTaskName(index),
	295	AliRecoParam::GetEventSpecieName(specie)));
	296	qa->Set(AliQAv1::kERROR, specie) ;
	297	} else if (value[specie] > fLowTestValue[AliQAv1::kWARNING] &&
	298	value[specie] <= fUpTestValue[AliQAv1::kWARNING]) {
	299
	300	if(fDebug>1)
	301	AliInfo(Form("QA-Warning. Task: %s. Specie: %s",
	302	AliQAv1::GetAliTaskName(index),
	303	AliRecoParam::GetEventSpecieName(specie)));
	304	qa->Set(AliQAv1::kWARNING, specie) ;
	305	} else if (value[specie] > fLowTestValue[AliQAv1::kINFO] &&
	306	value[specie] <= fUpTestValue[AliQAv1::kINFO] ) {
	307
	308	if(fDebug>1)
	309	AliInfo(Form("QA-Info. Task: %s. Specie: %s",
	310	AliQAv1::GetAliTaskName(index),
	311	AliRecoParam::GetEventSpecieName(specie)));
	312	qa->Set(AliQAv1::kINFO, specie) ;
	313	}
	314	}
	315	}
	316	}