[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::kSIM&&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::kSIM)
      rv[specie] = CheckSIM(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::CheckSIM(Int_t specie, TObjArray* list)
{
  // This method checks the QA histograms associated with simulation
  //
  // For TPC this is:
  // Digits : 
  // The digit histogram gives the ADC distribution for all sigbnals
  // above threshold. The check is just that there are digits.
  // Hits : The hit histograms are checked to see that they are not
  // empty. They contain a lot of detailed information on the
  // energyloss model (they were used to debug the AliRoot TPC use of
  // FLUKA).
  //
  // The check methods are simple:
  // We do not know if it is bad that histograms are missing because
  // this will always be the case for summable digits. So this check
  // is not possible here.
  // If digit histogram is empty (set error)
  // If one of the hit histograms are empty (set error)
  if(fDebug>0)
    AliInfo("In AliTPCQAChecker::CheckSIM");
  
  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* hHitsNhits = static_cast<TH1*> 
    (list->FindObject(Form("%s_hHitsNhits",specieName)));
  TH1* hHitsElectrons         = static_cast<TH1*> 
    (list->FindObject(Form("%s_hHitsElectrons",specieName)));
  TH1* hHitsRadius        = static_cast<TH1*> 
    (list->FindObject(Form("%s_hHitsRadius",specieName)));
  TH1* histHitsPrimPerCm          = static_cast<TH1*> 
    (list->FindObject(Form("%s_histHitsPrimPerCm",specieName)));
  TH1* histHitsElectronsPerCm          = static_cast<TH1*> 
    (list->FindObject(Form("%s_histHitsElectronsPerCm",specieName)));
  
//   if (!(hDigits) ||                                             // digit hists
//       !(hHitsNhits && hHitsElectrons && hHitsRadius && histHitsPrimPerCm && histHitsElectronsPerCm)) // hit hists
//     return -0.5; // fatal
  
  if (hDigits) {
    if(hDigits->GetEntries()==0) 
      return 0.25; // error
  }

  if (hHitsNhits && hHitsElectrons && hHitsRadius && 
      histHitsPrimPerCm && histHitsElectronsPerCm) {
    if (hHitsNhits->GetEntries()==0 || hHitsElectrons->GetEntries()==0 ||
	hHitsRadius->GetEntries()==0 || histHitsPrimPerCm->GetEntries()==0 ||
	histHitsElectronsPerCm->GetEntries()==0) 
      return 0.25; // error
  }

  return 1; // 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
db4f08c1	72	if (index!=AliQAv1::kRAW&&index!=AliQAv1::kSIM&&index!=AliQAv1::kREC&&index!=AliQAv1::kESD) {
ce0175fa	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]);
db4f08c1	89	if (index == AliQAv1::kSIM)
db4f08c1	90	rv[specie] = CheckSIM(specie, list[specie]);
ce0175fa	91	if (index == AliQAv1::kREC)
	92	rv[specie] = CheckREC(specie, list[specie]);
	93	if (index == AliQAv1::kESD)
	94	rv[specie] = CheckESD(specie, list[specie]);
	95
	96	if(fDebug>3)
	97	AliInfo(Form("Specie: %s. Task: %s. Value: %f",
	98	AliRecoParam::GetEventSpecieName(specie),
	99	AliQAv1::GetAliTaskName(index),
	100	rv[specie]));
	101	}
	102
	103	return rv;
	104	}
	105
	106	//______________________________________________________________________________
	107	Double_t AliTPCQAChecker::CheckRAW(Int_t specie, TObjArray* list)
	108	{
	109	/// Check ESD
	110	if(fDebug>0)
	111	AliInfo("In AliTPCQAChecker::CheckRAW");
	112
	113	if(fDebug>2)
	114	list->Print();
	115
	116	Char_t specieName[256];
	117	sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));
	118
	119	TH1* hRawsOccupancyVsSector = static_cast<TH1*>
	120	(list->FindObject(Form("%s_hRawsOccupancyVsSector",specieName)));
	121	TH1* hRawsQmaxVsSector = static_cast<TH1*>
	122	(list->FindObject(Form("%s_hRawsQmaxVsSector",specieName)));
	123	TH1* hRawsOccupancyVsEvent = static_cast<TH1*>
	124	(list->FindObject(Form("%s_hRawsOccupancyVsEvent",specieName)));
	125	TH1* hRawsNclustersVsEvent = static_cast<TH1*>
	126	(list->FindObject(Form("%s_hRawsNclustersVsEvent",specieName)));
	127
	128	if (!hRawsOccupancyVsSector \|\| !hRawsQmaxVsSector \|\|
	129	!hRawsOccupancyVsEvent \|\| !hRawsNclustersVsEvent)
	130	return -0.5; // fatal
	131
	132	if(hRawsOccupancyVsSector->GetEntries()==0) {
	133	return 0.25; // error - No TPC data!
	134	}
	135
	136	Int_t nBinsX = hRawsOccupancyVsSector->GetNbinsX();
	137	for(Int_t i = 1; i <= nBinsX; i++) {
	138
	139	if(hRawsOccupancyVsSector->GetBinContent(i)==0)
	140	return 0.75; // warning - no TPC data for at least one sector
	141	}
	142
	143	return 1.0; // ok
	144	}
	145
	146	//______________________________________________________________________________
db4f08c1	147	Double_t AliTPCQAChecker::CheckSIM(Int_t specie, TObjArray* list)
	148	{
	149	// This method checks the QA histograms associated with simulation
	150	//
	151	// For TPC this is:
	152	// Digits :
	153	// The digit histogram gives the ADC distribution for all sigbnals
	154	// above threshold. The check is just that there are digits.
	155	// Hits : The hit histograms are checked to see that they are not
	156	// empty. They contain a lot of detailed information on the
	157	// energyloss model (they were used to debug the AliRoot TPC use of
	158	// FLUKA).
	159	//
	160	// The check methods are simple:
	161	// We do not know if it is bad that histograms are missing because
	162	// this will always be the case for summable digits. So this check
	163	// is not possible here.
	164	// If digit histogram is empty (set error)
	165	// If one of the hit histograms are empty (set error)
	166	if(fDebug>0)
	167	AliInfo("In AliTPCQAChecker::CheckSIM");
	168
	169	if(fDebug>2)
	170	list->Print();
	171
	172	Char_t specieName[256];
	173	sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));
	174
	175	TH1* hDigits = static_cast<TH1*>
	176	(list->FindObject(Form("%s_hDigitsADC",specieName)));
	177	TH1* hHitsNhits = static_cast<TH1*>
	178	(list->FindObject(Form("%s_hHitsNhits",specieName)));
	179	TH1* hHitsElectrons = static_cast<TH1*>
	180	(list->FindObject(Form("%s_hHitsElectrons",specieName)));
	181	TH1* hHitsRadius = static_cast<TH1*>
	182	(list->FindObject(Form("%s_hHitsRadius",specieName)));
	183	TH1* histHitsPrimPerCm = static_cast<TH1*>
	184	(list->FindObject(Form("%s_histHitsPrimPerCm",specieName)));
	185	TH1* histHitsElectronsPerCm = static_cast<TH1*>
	186	(list->FindObject(Form("%s_histHitsElectronsPerCm",specieName)));
	187
	188	// if (!(hDigits) \|\| // digit hists
	189	// !(hHitsNhits && hHitsElectrons && hHitsRadius && histHitsPrimPerCm && histHitsElectronsPerCm)) // hit hists
	190	// return -0.5; // fatal
	191
	192	if (hDigits) {
	193	if(hDigits->GetEntries()==0)
	194	return 0.25; // error
	195	}
	196
	197	if (hHitsNhits && hHitsElectrons && hHitsRadius &&
	198	histHitsPrimPerCm && histHitsElectronsPerCm) {
	199	if (hHitsNhits->GetEntries()==0 \|\| hHitsElectrons->GetEntries()==0 \|\|
	200	hHitsRadius->GetEntries()==0 \|\| histHitsPrimPerCm->GetEntries()==0 \|\|
	201	histHitsElectronsPerCm->GetEntries()==0)
	202	return 0.25; // error
	203	}
	204
	205	return 1; // ok
	206	}
	207
	208	//______________________________________________________________________________
ce0175fa	209	Double_t AliTPCQAChecker::CheckREC(Int_t specie, TObjArray* list)
	210	{
	211	// This method checks the QA histograms associated with reconstruction
	212	//
	213	// For TPC this is:
	214	// DigitsR :
	215	// The digit histogram gives the ADC distribution for all sigbnals
	216	// above threshold. The check is just that there are digits.
	217	// RecPoints :
	218	// The cluster histograms are meant to give an idea about the gain
	219	// from the cluster charge and to indicate iof there are rows with
	220	// noise clusters, i.e., they are very visual.
	221	//
	222	// The check methods are simple:
	223	// If there are no histogram at all (set fatal)
	224	// If digit histogram is there, but there are no digits (set error)
	225	// If cluster histogram is there but there are less than 1000
	226	// clusters (set warning)
	227	// If there are more than 1000 clusters but no clusters for either short,
	228	// medium, or long pads (set error)
	229	if(fDebug>0)
	230	AliInfo("In AliTPCQAChecker::CheckREC");
	231
	232	if(fDebug>2)
	233	list->Print();
	234
	235	Char_t specieName[256];
	236	sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));
	237
	238	TH1* hDigits = static_cast<TH1*>
	239	(list->FindObject(Form("%s_hDigitsADC",specieName)));
	240	TH1* hNclustersVsRow = static_cast<TH1*>
	241	(list->FindObject(Form("%s_hRecPointsRow",specieName)));
	242	TH1* hQshort = static_cast<TH1*>
	243	(list->FindObject(Form("%s_hRecPointsQShort",specieName)));
	244	TH1* hQmedium = static_cast<TH1*>
	245	(list->FindObject(Form("%s_hRecPointsQMedium",specieName)));
	246	TH1* hQlong = static_cast<TH1*>
	247	(list->FindObject(Form("%s_hRecPointsQLong",specieName)));
	248	// The Qmax histograms are for now ignored
	249
	250	if (!hDigits && // digits missing
	251	(!hNclustersVsRow \|\| !hQshort \|\| !hQmedium \|\| !hQlong)) // 1 recpoint hist missing
	252	return -0.5; // fatal
	253
	254	if (hDigits && hDigits->GetEntries()==0)
	255	return 0.25; // error
	256
	257	if (hNclustersVsRow && hNclustersVsRow->GetEntries() < 1000) {
	258	return 0.75; // warning
	259	} else {
	260	if (!hQshort \|\| !hQlong \|\| !hQlong)
	261	return -0.5;// fatal - they should be there if the cluster vs row hist is there
	262	if (hQshort->GetEntries()==0 \|\| hQmedium->GetEntries()==0 \|\|
	263	hQlong->GetEntries()==0)
	264	return 0.25; // error
	265	}
	266	return 1; // ok
	267	}
	268
	269	//______________________________________________________________________________
	270	Double_t AliTPCQAChecker::CheckESD(Int_t specie, TObjArray* list)
	271	{
	272	// This method checks the QA histograms associated with ESDs
273	// (Note that there is aslo a globalQA which is running on all
274	// the ESD information so for now this is just a few basic
275	// histograms)
276	//
277	// The check methods are simple:
278	// If there are no histogram at all (set fatal)
279	//
280	if(fDebug>0)
281	AliInfo("In AliTPCQAChecker::CheckESD");
282
283	if(fDebug>2)
284	list->Print();
285
286	Char_t specieName[256];
287	sprintf(specieName, AliRecoParam::GetEventSpecieName(specie));
288
289	TH1* hESDclusters = static_cast<TH1*>
290	(list->FindObject(Form("%s_hESDclusters",specieName)));
291	TH1* hESDratio = static_cast<TH1*>
292	(list->FindObject(Form("%s_hESDratio",specieName)));
293	TH1* hESDpt = static_cast<TH1*>
294	(list->FindObject(Form("%s_hESDpt",specieName)));
295
296	if (!hESDclusters \|\| !hESDratio \|\| !hESDpt)
297	return -0.5; // fatal
298
299	return 1.0; // ok
300	}
301
302	//______________________________________________________________________________
303	void AliTPCQAChecker::Init(const AliQAv1::DETECTORINDEX_t det)
304	{
305	/// intialises QA and QA checker settings
306	if(fDebug>0)
307	AliInfo("In AliTPCQAChecker::Init");
308	AliQAv1::Instance(det) ;
309	Float_t hiValue[AliQAv1::kNBIT] ;
310	Float_t lowValue[AliQAv1::kNBIT] ;
311	hiValue[AliQAv1::kINFO] = 1.00;
312	lowValue[AliQAv1::kINFO] = 0.99;
313	hiValue[AliQAv1::kWARNING] = 0.99;
314	lowValue[AliQAv1::kWARNING] = 0.50;
315	hiValue[AliQAv1::kERROR] = 0.50;
316	lowValue[AliQAv1::kERROR] = 0.00;
317	hiValue[AliQAv1::kFATAL] = 0.00;
318	lowValue[AliQAv1::kFATAL] =-1.00;
319	// SetHiLo(&hiValue[0], &lowValue[0]) ;
320	SetHiLo(hiValue, lowValue) ;
321	}
322
323	//______________________________________________________________________________
324	void
325	AliTPCQAChecker::SetQA(AliQAv1::ALITASK_t index, Double_t * value) const
326	{
327	/// sets the QA according the return value of the Check
328
329	if(fDebug>0)
330	AliInfo("In AliTPCQAChecker::SetQA");
331
332	AliQAv1 * qa = AliQAv1::Instance(index);
333
334	for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
44f32dd2	335
ce0175fa	336	if (value==NULL) { // No checker is implemented, set all QA to Fatal
	337
	338	if(fDebug>1)
	339	AliInfo(Form("Fatal QA. Task: %s. Specie: %s",
	340	AliQAv1::GetAliTaskName(index),
	341	AliRecoParam::GetEventSpecieName(specie)));
	342	qa->Set(AliQAv1::kFATAL, specie) ;
	343	} else {
	344
	345	if ( value[specie] >= fLowTestValue[AliQAv1::kFATAL] &&
	346	value[specie] < fUpTestValue[AliQAv1::kFATAL] ) {
	347
	348	if(fDebug>1)
	349	AliInfo(Form("QA-Fatal. Task: %s. Specie: %s",
	350	AliQAv1::GetAliTaskName(index),
	351	AliRecoParam::GetEventSpecieName(specie)));
	352	qa->Set(AliQAv1::kFATAL, specie) ;
	353	} else if ( value[specie] > fLowTestValue[AliQAv1::kERROR] &&
	354	value[specie] <= fUpTestValue[AliQAv1::kERROR] ) {
	355
	356	if(fDebug>1)
	357	AliInfo(Form("QA-Error. Task: %s. Specie: %s",
	358	AliQAv1::GetAliTaskName(index),
	359	AliRecoParam::GetEventSpecieName(specie)));
	360	qa->Set(AliQAv1::kERROR, specie) ;
	361	} else if (value[specie] > fLowTestValue[AliQAv1::kWARNING] &&
	362	value[specie] <= fUpTestValue[AliQAv1::kWARNING]) {
	363
	364	if(fDebug>1)
	365	AliInfo(Form("QA-Warning. Task: %s. Specie: %s",
	366	AliQAv1::GetAliTaskName(index),
	367	AliRecoParam::GetEventSpecieName(specie)));
	368	qa->Set(AliQAv1::kWARNING, specie) ;
	369	} else if (value[specie] > fLowTestValue[AliQAv1::kINFO] &&
	370	value[specie] <= fUpTestValue[AliQAv1::kINFO] ) {
	371
	372	if(fDebug>1)
	373	AliInfo(Form("QA-Info. Task: %s. Specie: %s",
	374	AliQAv1::GetAliTaskName(index),
	375	AliRecoParam::GetEventSpecieName(specie)));
	376	qa->Set(AliQAv1::kINFO, specie) ;
	377	}
	378	}
	379	}
	380	}