[u/mrichter/AliRoot.git] / T0 / AliT0QAChecker.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.                  *
 **************************************************************************/
//...
//  Checks the quality assurance. 
//  By comparing with reference data
//  Skeleton for T0
//---------------------------------------------
//checkig without reference data:
//for RAW QA all histograms should have approximatly the same 
//number of entries as RefPoint
//for Rec Points checks 
//  - amplitude measured by 2 methos
//  - online and offline T0 measurements
// for ESD quality of reconstruction ( and measurements):
// RMS of vertex and T0 less than 75ps
//
// Alla.Maevskaya@cern.ch   
//...

// --- ROOT system ---
#include <Riostream.h>
#include <TClass.h>
#include <TH1F.h> 
#include <TF1.h> 
#include <TFitResultPtr.h>
#include <TH2.h> 
#include <TIterator.h> 
#include <TKey.h> 
#include <TFile.h> 
#include <TMath.h>
#include <TString.h>
#include <TPaveText.h>

// --- Standard library ---

// --- AliRoot header files ---
#include "AliLog.h"
#include "AliQAv1.h"
#include "AliQAChecker.h"
#include "AliCDBEntry.h"
#include "AliQAManager.h"
#include "AliT0QAChecker.h"

ClassImp(AliT0QAChecker)
//____________________________________________________________________________
AliT0QAChecker::AliT0QAChecker() :
AliQACheckerBase("T0","T0 Quality Assurance Checker")

{
  // Standard constructor

}

//____________________________________________________________________________
AliT0QAChecker::AliT0QAChecker(const AliT0QAChecker& qac):
  AliQACheckerBase(qac.GetName(), qac.GetTitle()) 
{
  // copy constructor
  AliError("Copy should not be used with this class\n");
}
//____________________________________________________________________________
AliT0QAChecker& AliT0QAChecker::operator=(const AliT0QAChecker& qac){
  // assignment operator
  this->~AliT0QAChecker();
  new(this)AliT0QAChecker(qac);
  return *this;
}


//____________________________________________________________________________
AliT0QAChecker::~AliT0QAChecker(){
  // destructor

}

//__________________________________________________________________
void AliT0QAChecker::Check(Double_t *  test, AliQAv1::ALITASK_t index, TObjArray ** list, const AliDetectorRecoParam * /*recoParam*/)
{

  // Super-basic check on the QA histograms on the input list:
  // look whether they are empty!
    
 char * detOCDBDir = Form("T0/%s/%s", AliQAv1::GetRefOCDBDirName(), AliQAv1::GetRefDataDirName()) ; 

  AliCDBEntry *QARefRec = AliQAManager::QAManager()->Get(detOCDBDir);
  //  QARefRec->Dump();
  if( !QARefRec){
    AliInfo("QA reference data NOT retrieved for Reconstruction check. No T0 reference distribution");
  }
  
  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) 
    test[specie]    = 10.0 ; 

  for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
    //  TString dataType = AliQAv1::GetAliTaskName(index);
    if (list[specie]->GetEntries() == 0){
      test[specie] = 1. ; // nothing to check
    }
    else {
      if (index == AliQAv1::kRAW && AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCalib)
       //      if (index == AliQAv1::kRAW )
	{
	  test[specie] = CheckRaw(list[specie]);
	}
      
       if (index == AliQAv1::kESD && AliRecoParam::Convert(specie) != AliRecoParam::kCalib)
	  test[specie] = CheckESD(list[specie]);
    }
  
  }
}

//--------------------------------------------------------------------------
Double_t AliT0QAChecker::CheckRaw(TObjArray *listrec) const
{
   Float_t checkr = 0;
   
  TString hname[10];
   TH1*hdata; 
   const char *cname;
   TH1 *fhRawEff[10];
   Int_t nh=0;
 

   Int_t nnn[4] = { 420, 458, 459, 460};
   for (Int_t ir=0; ir<4; ir++)
     {
	 hdata = (TH1*) listrec->UncheckedAt(nnn[ir]);
	 if(hdata) {
	   cname = hdata->GetName();
	   hname[ir] = cname;
	   fhRawEff[nh] = hdata;
	   nh++;
	 }
    }
     
   TLine *linelowyellow = new TLine(0, 0.7, 24, 0.7);    
   linelowyellow->SetLineColor(5);
   linelowyellow->SetLineStyle(3);
   linelowyellow->SetLineWidth(4);
   TLine *linelowred = new TLine(0, 0.2, 24, 0.2);    
   linelowred->SetLineColor(2);
   linelowred->SetLineStyle(3);
   linelowred->SetLineWidth(4);

   Float_t thryell = 0.7;
   //   TPaveText* text = new TPaveText(0.30,0.50,0.99,0.99,"NDC");    
   Float_t thrred = 0.2;
   Float_t chcont =0;
    for (Int_t ih= 0; ih<4; ih++)
     { 
       fhRawEff[ih]->SetLineWidth(2);
       fhRawEff[ih]->SetMaximum(2.);
       fhRawEff[ih]->SetMinimum(0.);
       fhRawEff[ih]->GetListOfFunctions()->Add(linelowyellow);
       fhRawEff[ih]->GetListOfFunctions()->Add(linelowred);

       Int_t nbins= fhRawEff[ih]->GetNbinsX();
       Bool_t yell = kFALSE;
       Bool_t red = kFALSE;
       for (Int_t in=1; in<nbins-1; in++)
	 {
	   if(ih==0 && in==5) continue;
	   chcont=fhRawEff[ih]->GetBinContent(in);
	   if (chcont < thryell  ) yell = kTRUE;
	   if (chcont < thrred  ) red = kTRUE;
	 }
       
       if (! yell && !red) {
	 AliDebug(AliQAv1::GetQADebugLevel(), Form(" efficiency in all channes %s  is good", fhRawEff[ih]->GetName() ));
	 checkr=1.;
	 //	 text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	 //	 text->AddText(Form(" No problems "));
	 //	 text->SetFillColor(3);
       }
       
       if(red ) {
	 checkr = 0.;
	 AliDebug(AliQAv1::GetQADebugLevel(), Form(" efficiency in all channes %s  is not so good", fhRawEff[ih]->GetName() ));
	 //	 text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	 //	 text->AddText(Form("Very serious problem, call expert "));
	 //	 text->SetFillColor(2);
       }
       
       if ( yell && !red) {
	 AliDebug(AliQAv1::GetQADebugLevel(), Form(" efficiency in all channes %s  is not so good", fhRawEff[ih]->GetName() ));
	 checkr=0.75;
	 //	 text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	 //	 text->AddText(Form("Some problems "));
	 //	 text->SetFillColor(5);
      }
       // fhRawEff[ih]->GetListOfFunctions()->Add(text);	

     }
  return checkr;
  
}


//--------------------------------------------------------------------------
Double_t AliT0QAChecker::CheckESD(TObjArray *listrec ) const
{
  Float_t checkr = 0;
  TH1 *fhESD;
 
  fhESD  = (TH1*) listrec->UncheckedAt(2);
  if(fhESD){
    AliDebug(AliQAv1::GetQADebugLevel(), Form("count %s ", fhESD->GetName()) );
     TF1 *f1 = new TF1("f1","gaus",-1,1);
    fhESD->Fit("f1","R","Q", -1,1);
    Double_t par[3];
    f1->GetParameters(&par[0]);
    
    TPaveText* text = new TPaveText(0.30,0.50,0.99,0.99,"NDC");
    
    text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
    
    AliDebug(AliQAv1::GetQADebugLevel(), Form("numentries %f mean %f  #sigma %f", fhESD->GetEntries(),par[1], par[2]));
    
    
    if (par[2] > 0.07 && par[2] < 1.) {
      checkr=0.5;
       text->AddText(Form("not good resolution :\n %f ns\n", par[2] ));
       text->SetFillColor(5);
       printf("T0 detector resolution is not good enouph: %f ns\n",par[2] );
    }
    if(TMath::Abs(par[1])>0.05) {
      checkr = 0.5;
      text->AddText(Form(" Check clock shift on %f ns", par[1]));
      text->SetFillColor(5);
    }
    if (par[2] >  1. || TMath::Abs(par[1])>0.1) {
      checkr = 0.25;
      text->AddText(Form(" Bad resolution:\n mean %f ns sigma %f ns", par[1], par[2]));
      text->SetFillColor(2);
      
      fhESD->GetListOfFunctions()->Add(text);	
       AliDebug(AliQAv1::GetQADebugLevel(),
		Form("Please, check calibration: shift= %f resolution %f test=%f\n",
		     par[1], par[2], checkr) ) ; 
    }
  }
  else
    {
      AliDebug(AliQAv1::GetQADebugLevel(),
	       Form("No ESD QA histogram found, nothing to check"));
	checkr=0;
    }
  
  
  return checkr;
}
Commit	Line	Data
c03351ac	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	**************************************************************************/
c03351ac	15	//...
	16	// Checks the quality assurance.
	17	// By comparing with reference data
	18	// Skeleton for T0
f6c2d5e2	19	//---------------------------------------------
	20	//checkig without reference data:
	21	//for RAW QA all histograms should have approximatly the same
	22	//number of entries as RefPoint
	23	//for Rec Points checks
	24	// - amplitude measured by 2 methos
	25	// - online and offline T0 measurements
	26	// for ESD quality of reconstruction ( and measurements):
	27	// RMS of vertex and T0 less than 75ps
	28	//
	29	// Alla.Maevskaya@cern.ch
c03351ac	30	//...
	31
	32	// --- ROOT system ---
b09247a2	33	#include <Riostream.h>
c03351ac	34	#include <TClass.h>
c03351ac	35	#include <TH1F.h>
be3ada9f	36	#include <TF1.h>
be3ada9f	37	#include <TFitResultPtr.h>
0d849919	38	#include <TH2.h>
c03351ac	39	#include <TIterator.h>
	40	#include <TKey.h>
	41	#include <TFile.h>
446d6ec4	42	#include <TMath.h>
0d849919	43	#include <TString.h>
be3ada9f	44	#include <TPaveText.h>
c03351ac	45
	46	// --- Standard library ---
	47
	48	// --- AliRoot header files ---
	49	#include "AliLog.h"
4e25ac79	50	#include "AliQAv1.h"
c03351ac	51	#include "AliQAChecker.h"
394c1a6d	52	#include "AliCDBEntry.h"
394c1a6d	53	#include "AliQAManager.h"
c03351ac	54	#include "AliT0QAChecker.h"
	55
	56	ClassImp(AliT0QAChecker)
be3ada9f	57	//____________________________________________________________________________
	58	AliT0QAChecker::AliT0QAChecker() :
	59	AliQACheckerBase("T0","T0 Quality Assurance Checker")
	60
	61	{
	62	// Standard constructor
	63
	64	}
	65
	66	//____________________________________________________________________________
	67	AliT0QAChecker::AliT0QAChecker(const AliT0QAChecker& qac):
	68	AliQACheckerBase(qac.GetName(), qac.GetTitle())
	69	{
	70	// copy constructor
	71	AliError("Copy should not be used with this class\n");
	72	}
	73	//____________________________________________________________________________
	74	AliT0QAChecker& AliT0QAChecker::operator=(const AliT0QAChecker& qac){
	75	// assignment operator
	76	this->~AliT0QAChecker();
	77	new(this)AliT0QAChecker(qac);
	78	return *this;
	79	}
	80
	81
	82	//____________________________________________________________________________
	83	AliT0QAChecker::~AliT0QAChecker(){
	84	// destructor
	85
	86	}
c03351ac	87
4e67e173	88	//__________________________________________________________________
a42ceb0e	89	void AliT0QAChecker::Check(Double_t * test, AliQAv1::ALITASK_t index, TObjArray ** list, const AliDetectorRecoParam * /recoParam/)
3685c21c	90	{
	91
	92	// Super-basic check on the QA histograms on the input list:
	93	// look whether they are empty!
a42ceb0e	94
be3ada9f	95	char * detOCDBDir = Form("T0/%s/%s", AliQAv1::GetRefOCDBDirName(), AliQAv1::GetRefDataDirName()) ;
8825d145	96
394c1a6d	97	AliCDBEntry *QARefRec = AliQAManager::QAManager()->Get(detOCDBDir);
	98	// QARefRec->Dump();
	99	if( !QARefRec){
	100	AliInfo("QA reference data NOT retrieved for Reconstruction check. No T0 reference distribution");
	101	}
8825d145	102
be3ada9f	103	for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++)
57acd2d2	104	test[specie] = 10.0 ;
c724dd64	105
57acd2d2	106	for (Int_t specie = 0 ; specie < AliRecoParam::kNSpecies ; specie++) {
4e25ac79	107	// TString dataType = AliQAv1::GetAliTaskName(index);
57acd2d2	108	if (list[specie]->GetEntries() == 0){
	109	test[specie] = 1. ; // nothing to check
	110	}
	111	else {
be3ada9f	112	if (index == AliQAv1::kRAW && AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCalib)
	113	// if (index == AliQAv1::kRAW )
	114	{
	115	test[specie] = CheckRaw(list[specie]);
8825d145	116	}
be3ada9f	117
	118	if (index == AliQAv1::kESD && AliRecoParam::Convert(specie) != AliRecoParam::kCalib)
	119	test[specie] = CheckESD(list[specie]);
57acd2d2	120	}
be3ada9f	121
	122	}
	123	}
394c1a6d	124
be3ada9f	125	//--------------------------------------------------------------------------
	126	Double_t AliT0QAChecker::CheckRaw(TObjArray *listrec) const
	127	{
	128	Float_t checkr = 0;
	129
	130	TString hname[10];
	131	TH1*hdata;
	132	const char *cname;
	133	TH1 *fhRawEff[10];
	134	Int_t nh=0;
	135
394c1a6d	136
be3ada9f	137	Int_t nnn[4] = { 420, 458, 459, 460};
	138	for (Int_t ir=0; ir<4; ir++)
	139	{
	140	hdata = (TH1*) listrec->UncheckedAt(nnn[ir]);
	141	if(hdata) {
	142	cname = hdata->GetName();
	143	hname[ir] = cname;
	144	fhRawEff[nh] = hdata;
	145	nh++;
	146	}
	147	}
	148
	149	TLine *linelowyellow = new TLine(0, 0.7, 24, 0.7);
	150	linelowyellow->SetLineColor(5);
	151	linelowyellow->SetLineStyle(3);
	152	linelowyellow->SetLineWidth(4);
	153	TLine *linelowred = new TLine(0, 0.2, 24, 0.2);
	154	linelowred->SetLineColor(2);
	155	linelowred->SetLineStyle(3);
	156	linelowred->SetLineWidth(4);
8825d145	157
be3ada9f	158	Float_t thryell = 0.7;
	159	// TPaveText* text = new TPaveText(0.30,0.50,0.99,0.99,"NDC");
	160	Float_t thrred = 0.2;
	161	Float_t chcont =0;
	162	for (Int_t ih= 0; ih<4; ih++)
	163	{
	164	fhRawEff[ih]->SetLineWidth(2);
	165	fhRawEff[ih]->SetMaximum(2.);
	166	fhRawEff[ih]->SetMinimum(0.);
	167	fhRawEff[ih]->GetListOfFunctions()->Add(linelowyellow);
	168	fhRawEff[ih]->GetListOfFunctions()->Add(linelowred);
8825d145	169
be3ada9f	170	Int_t nbins= fhRawEff[ih]->GetNbinsX();
	171	Bool_t yell = kFALSE;
	172	Bool_t red = kFALSE;
	173	for (Int_t in=1; in<nbins-1; in++)
	174	{
	175	if(ih==0 && in==5) continue;
	176	chcont=fhRawEff[ih]->GetBinContent(in);
	177	if (chcont < thryell ) yell = kTRUE;
	178	if (chcont < thrred ) red = kTRUE;
	179	}
	180
	181	if (! yell && !red) {
	182	AliDebug(AliQAv1::GetQADebugLevel(), Form(" efficiency in all channes %s is good", fhRawEff[ih]->GetName() ));
	183	checkr=1.;
	184	// text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	185	// text->AddText(Form(" No problems "));
	186	// text->SetFillColor(3);
	187	}
	188
	189	if(red ) {
	190	checkr = 0.;
	191	AliDebug(AliQAv1::GetQADebugLevel(), Form(" efficiency in all channes %s is not so good", fhRawEff[ih]->GetName() ));
	192	// text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	193	// text->AddText(Form("Very serious problem, call expert "));
	194	// text->SetFillColor(2);
	195	}
	196
	197	if ( yell && !red) {
	198	AliDebug(AliQAv1::GetQADebugLevel(), Form(" efficiency in all channes %s is not so good", fhRawEff[ih]->GetName() ));
	199	checkr=0.75;
	200	// text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	201	// text->AddText(Form("Some problems "));
	202	// text->SetFillColor(5);
57acd2d2	203	}
be3ada9f	204	// fhRawEff[ih]->GetListOfFunctions()->Add(text);
	205
	206	}
	207	return checkr;
	208
3685c21c	209	}
5ed41460	210
be3ada9f	211
394c1a6d	212	//--------------------------------------------------------------------------
be3ada9f	213	Double_t AliT0QAChecker::CheckESD(TObjArray *listrec ) const
394c1a6d	214	{
394c1a6d	215	Float_t checkr = 0;
be3ada9f	216	TH1 *fhESD;
	217
	218	fhESD = (TH1*) listrec->UncheckedAt(2);
	219	if(fhESD){
	220	AliDebug(AliQAv1::GetQADebugLevel(), Form("count %s ", fhESD->GetName()) );
	221	TF1 *f1 = new TF1("f1","gaus",-1,1);
	222	fhESD->Fit("f1","R","Q", -1,1);
	223	Double_t par[3];
	224	f1->GetParameters(&par[0]);
	225
	226	TPaveText* text = new TPaveText(0.30,0.50,0.99,0.99,"NDC");
	227
	228	text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	229
b0fa4a0f	230	AliDebug(AliQAv1::GetQADebugLevel(), Form("numentries %f mean %f #sigma %f", fhESD->GetEntries(),par[1], par[2]));
be3ada9f	231
	232
	233	if (par[2] > 0.07 && par[2] < 1.) {
	234	checkr=0.5;
	235	text->AddText(Form("not good resolution :\n %f ns\n", par[2] ));
	236	text->SetFillColor(5);
	237	printf("T0 detector resolution is not good enouph: %f ns\n",par[2] );
	238	}
	239	if(TMath::Abs(par[1])>0.05) {
	240	checkr = 0.5;
	241	text->AddText(Form(" Check clock shift on %f ns", par[1]));
	242	text->SetFillColor(5);
	243	}
	244	if (par[2] > 1. \|\| TMath::Abs(par[1])>0.1) {
	245	checkr = 0.25;
	246	text->AddText(Form(" Bad resolution:\n mean %f ns sigma %f ns", par[1], par[2]));
	247	text->SetFillColor(2);
8825d145	248
be3ada9f	249	fhESD->GetListOfFunctions()->Add(text);
be3ada9f	250	AliDebug(AliQAv1::GetQADebugLevel(),
	251	Form("Please, check calibration: shift= %f resolution %f test=%f\n",
	252	par[1], par[2], checkr) ) ;
	253	}
394c1a6d	254	}
be3ada9f	255	else
	256	{
	257	AliDebug(AliQAv1::GetQADebugLevel(),
	258	Form("No ESD QA histogram found, nothing to check"));
	259	checkr=0;
	260	}
	261
332b9234	262
394c1a6d	263	return checkr;
394c1a6d	264	}