[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 (!(AliQAv1::Instance()->IsEventSpecieSet(specie) && list[specie]) || list[specie]->GetEntries() == 0) {
      test[specie] = 1. ; // nothing to check
      continue;
    }
    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};
   Int_t nnn[4] = { 169, 207, 208, 209};
   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);	

     }
    delete linelowyellow;
    delete linelowred;
    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 %d mean %f  #sigma %f", (int)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);
92664bc8	108	if (!(AliQAv1::Instance()->IsEventSpecieSet(specie) && list[specie]) \|\| list[specie]->GetEntries() == 0) {
57acd2d2	109	test[specie] = 1. ; // nothing to check
92664bc8	110	continue;
57acd2d2	111	}
92664bc8	112	if (index == AliQAv1::kRAW && AliRecoParam::ConvertIndex(specie) == AliRecoParam::kCalib)
	113	// if (index == AliQAv1::kRAW )
	114	{
	115	test[specie] = CheckRaw(list[specie]);
	116	}
	117	if (index == AliQAv1::kESD && AliRecoParam::Convert(specie) != AliRecoParam::kCalib)
	118	test[specie] = CheckESD(list[specie]);
	119	//
be3ada9f	120	}
be3ada9f	121	}
394c1a6d	122
be3ada9f	123	//--------------------------------------------------------------------------
	124	Double_t AliT0QAChecker::CheckRaw(TObjArray *listrec) const
	125	{
	126	Float_t checkr = 0;
	127
	128	TString hname[10];
	129	TH1*hdata;
	130	const char *cname;
	131	TH1 *fhRawEff[10];
	132	Int_t nh=0;
	133
394c1a6d	134
f88549d9	135	// Int_t nnn[4] = { 420, 458, 459, 460};
13579355	136	Int_t nnn[4] = { 169, 207, 208, 209};
be3ada9f	137	for (Int_t ir=0; ir<4; ir++)
	138	{
	139	hdata = (TH1*) listrec->UncheckedAt(nnn[ir]);
	140	if(hdata) {
	141	cname = hdata->GetName();
	142	hname[ir] = cname;
	143	fhRawEff[nh] = hdata;
	144	nh++;
	145	}
	146	}
	147
	148	TLine *linelowyellow = new TLine(0, 0.7, 24, 0.7);
	149	linelowyellow->SetLineColor(5);
	150	linelowyellow->SetLineStyle(3);
	151	linelowyellow->SetLineWidth(4);
	152	TLine *linelowred = new TLine(0, 0.2, 24, 0.2);
	153	linelowred->SetLineColor(2);
	154	linelowred->SetLineStyle(3);
	155	linelowred->SetLineWidth(4);
8825d145	156
be3ada9f	157	Float_t thryell = 0.7;
	158	// TPaveText* text = new TPaveText(0.30,0.50,0.99,0.99,"NDC");
	159	Float_t thrred = 0.2;
	160	Float_t chcont =0;
	161	for (Int_t ih= 0; ih<4; ih++)
	162	{
	163	fhRawEff[ih]->SetLineWidth(2);
	164	fhRawEff[ih]->SetMaximum(2.);
	165	fhRawEff[ih]->SetMinimum(0.);
	166	fhRawEff[ih]->GetListOfFunctions()->Add(linelowyellow);
	167	fhRawEff[ih]->GetListOfFunctions()->Add(linelowred);
8825d145	168
be3ada9f	169	Int_t nbins= fhRawEff[ih]->GetNbinsX();
	170	Bool_t yell = kFALSE;
	171	Bool_t red = kFALSE;
	172	for (Int_t in=1; in<nbins-1; in++)
	173	{
	174	if(ih==0 && in==5) continue;
	175	chcont=fhRawEff[ih]->GetBinContent(in);
	176	if (chcont < thryell ) yell = kTRUE;
	177	if (chcont < thrred ) red = kTRUE;
	178	}
	179
	180	if (! yell && !red) {
	181	AliDebug(AliQAv1::GetQADebugLevel(), Form(" efficiency in all channes %s is good", fhRawEff[ih]->GetName() ));
	182	checkr=1.;
	183	// text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	184	// text->AddText(Form(" No problems "));
	185	// text->SetFillColor(3);
	186	}
	187
	188	if(red ) {
	189	checkr = 0.;
	190	AliDebug(AliQAv1::GetQADebugLevel(), Form(" efficiency in all channes %s is not so good", fhRawEff[ih]->GetName() ));
	191	// text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	192	// text->AddText(Form("Very serious problem, call expert "));
	193	// text->SetFillColor(2);
	194	}
	195
	196	if ( yell && !red) {
	197	AliDebug(AliQAv1::GetQADebugLevel(), Form(" efficiency in all channes %s is not so good", fhRawEff[ih]->GetName() ));
	198	checkr=0.75;
	199	// text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	200	// text->AddText(Form("Some problems "));
	201	// text->SetFillColor(5);
57acd2d2	202	}
be3ada9f	203	// fhRawEff[ih]->GetListOfFunctions()->Add(text);
	204
	205	}
1d6bd73f	206	delete linelowyellow;
	207	delete linelowred;
	208	return checkr;
be3ada9f	209
3685c21c	210	}
5ed41460	211
be3ada9f	212
394c1a6d	213	//--------------------------------------------------------------------------
be3ada9f	214	Double_t AliT0QAChecker::CheckESD(TObjArray *listrec ) const
394c1a6d	215	{
394c1a6d	216	Float_t checkr = 0;
be3ada9f	217	TH1 *fhESD;
	218
	219	fhESD = (TH1*) listrec->UncheckedAt(2);
	220	if(fhESD){
	221	AliDebug(AliQAv1::GetQADebugLevel(), Form("count %s ", fhESD->GetName()) );
	222	TF1 *f1 = new TF1("f1","gaus",-1,1);
	223	fhESD->Fit("f1","R","Q", -1,1);
	224	Double_t par[3];
	225	f1->GetParameters(&par[0]);
	226
	227	TPaveText* text = new TPaveText(0.30,0.50,0.99,0.99,"NDC");
	228
	229	text->AddText(Form("T0 RUN %d ",AliCDBManager::Instance()->GetRun()));
	230
db33b25a	231	AliDebug(AliQAv1::GetQADebugLevel(), Form("numentries %d mean %f #sigma %f", (int)fhESD->GetEntries(),par[1], par[2]));
be3ada9f	232
	233
	234	if (par[2] > 0.07 && par[2] < 1.) {
	235	checkr=0.5;
	236	text->AddText(Form("not good resolution :\n %f ns\n", par[2] ));
	237	text->SetFillColor(5);
	238	printf("T0 detector resolution is not good enouph: %f ns\n",par[2] );
	239	}
	240	if(TMath::Abs(par[1])>0.05) {
	241	checkr = 0.5;
	242	text->AddText(Form(" Check clock shift on %f ns", par[1]));
	243	text->SetFillColor(5);
	244	}
	245	if (par[2] > 1. \|\| TMath::Abs(par[1])>0.1) {
	246	checkr = 0.25;
	247	text->AddText(Form(" Bad resolution:\n mean %f ns sigma %f ns", par[1], par[2]));
	248	text->SetFillColor(2);
8825d145	249
be3ada9f	250	fhESD->GetListOfFunctions()->Add(text);
be3ada9f	251	AliDebug(AliQAv1::GetQADebugLevel(),
	252	Form("Please, check calibration: shift= %f resolution %f test=%f\n",
	253	par[1], par[2], checkr) ) ;
	254	}
394c1a6d	255	}
be3ada9f	256	else
	257	{
	258	AliDebug(AliQAv1::GetQADebugLevel(),
	259	Form("No ESD QA histogram found, nothing to check"));
	260	checkr=0;
	261	}
	262
332b9234	263
394c1a6d	264	return checkr;
394c1a6d	265	}