[u/mrichter/AliRoot.git] / HMPID / AliHMPIDQADataMaker.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.                  *
 **************************************************************************/


/* $Id$ */

//---
//  Produces the data needed to calculate the quality assurance. 
//  All data must be mergeable objects.
//  A. Mastroserio
//---

// --- ROOT system ---
#include <TClonesArray.h>
#include <TFile.h> 
#include <TH1F.h> 
#include <TH2F.h>
#include <TH1I.h> 
#include <TDirectory.h>
#include <Riostream.h>
// --- Standard library ---

// --- AliRoot header files ---
#include "AliESDCaloCluster.h"
#include "AliESDEvent.h"
#include "AliLog.h"
#include "AliHMPIDDigit.h"
#include "AliHMPIDHit.h"
#include "AliHMPIDCluster.h"
#include "AliHMPIDQADataMaker.h"

ClassImp(AliHMPIDQADataMaker)
           
//____________________________________________________________________________ 
  AliHMPIDQADataMaker::AliHMPIDQADataMaker() : 
  AliQADataMaker(AliQA::GetDetName(AliQA::kHMPID), "HMPID Quality Assurance Data Maker"),
  fhHitQdc(0x0), 
  fhSDigits(0x0),
  fhDigPcEvt(0x0),
  fhDigChEvt(0x0),
  fhDigQ(0x0),
  fhCluEvt(0x0),
  fhCluChi2(0x0),
  fhCluQ(0x0),
  fhCluFlg(0x0), 
  fhCluSize(0x0),  
  fhMipCluSize(0x0),
  fhCkovP(0x0),
  fhSigP(0x0),
  fhMipXY(0x0),
  fhDifXY(0x0)
{
  // ctor
  for(Int_t i=0; i<7; i++) fhHitMap[i]=0x0;
  for(Int_t j=0; j<5; j++) fhPid[j]=0x0;
//   fDetectorDir = fOutput->GetDirectory(GetName()) ;  
//   if (!fDetectorDir) 
//     fDetectorDir = fOutput->mkdir(GetName()) ;  
}

//____________________________________________________________________________ 
AliHMPIDQADataMaker::AliHMPIDQADataMaker(const AliHMPIDQADataMaker& qadm) :
  AliQADataMaker(), 
  fhHitQdc(qadm.fhHitQdc), 
  fhSDigits(qadm.fhSDigits),
  fhDigPcEvt(qadm.fhDigPcEvt),
  fhDigChEvt(qadm.fhDigChEvt),
  fhDigQ(qadm.fhDigQ),
  fhCluEvt(qadm.fhCluEvt),
  fhCluChi2(qadm.fhCluChi2),
  fhCluQ(qadm.fhCluQ),
  fhCluFlg(qadm.fhCluFlg),
  fhCluSize(qadm.fhCluSize),
  fhMipCluSize(qadm.fhMipCluSize),
  fhCkovP(qadm.fhCkovP),
  fhSigP(qadm.fhSigP),
  fhMipXY(qadm.fhMipXY),
  fhDifXY(qadm.fhDifXY)
{
  //copy ctor 
  for(Int_t i=0; i<7; i++) fhHitMap[i]=qadm.fhHitMap[i];
  for(Int_t j=0; j<5; j++) fhPid[j]=qadm.fhPid[j];

  SetName((const char*)qadm.GetName()) ; 
  SetTitle((const char*)qadm.GetTitle()); 
}

//__________________________________________________________________
AliHMPIDQADataMaker& AliHMPIDQADataMaker::operator = (const AliHMPIDQADataMaker& qadm )
{
  // Equal operator.
  this->~AliHMPIDQADataMaker();
  new(this) AliHMPIDQADataMaker(qadm);
  return *this;
}
 
//____________________________________________________________________________ 
void AliHMPIDQADataMaker::InitHits()
{
  // create Hits histograms in Hits subdir
     fhHitQdc=new TH1F("HitQdc","HMPID Hit Qdc all chamber;QDC",500,0,4000);
     for(Int_t iCh=0;iCh<7;iCh++) fhHitMap[iCh]=new TH2F(Form("HMPID HitMap%i",iCh),Form("Ch%i;x_{Hit};y_{Hit}",iCh),162,-1,161,146,-1,145);   
}

//____________________________________________________________________________ 
void AliHMPIDQADataMaker::InitDigits()
{
  // create Digits histograms in Digits subdir
      fhDigPcEvt=new TH1F("hDigPcEvt","PC occupancy",156,-1,77);
      fhDigChEvt=new TH1F("hDigChEvt","Chamber occupancy",32,-1,7);
      fhDigQ  =new TH1F("Q        "," digit charge               ",3000,0,3000);
}

//____________________________________________________________________________ 
void AliHMPIDQADataMaker::InitSDigits()
{
  // create SDigits histograms in SDigits subdir
      fhSDigits     = new TH1F("hHmpidSDigits",    "SDigits Q  distribution in HMPID",  500, 0., 5000.) ; 
}

//____________________________________________________________________________ 

void AliHMPIDQADataMaker::InitRecPoints()
{
  // create cluster histograms in RecPoint subdir
      fhCluEvt=new TH1F("CluPerEvt","# clusters per chamber",16,-1,7);
      fhCluChi2  =new TH1F("CluChi2"  ,"Chi2 "               ,1000,0,100);
      fhCluQ   =new TH1F("CluQ"   ,"Cluster charge"        ,3000,0,3000);
      fhCluFlg   =new TH1F("CluFlg"   ,"Cluster flag"        ,14,-1.5,12.5);
      fhCluSize  =new TH1F("CluSize"  ,"Raw cluster size    ",100,0,100);
      fhMipCluSize =new TH1F("MipCluSize"  ,"Mip cluster size    ",100,0,100);
}
//____________________________________________________________________________
void AliHMPIDQADataMaker::InitESDs()
{
  //create ESDs histograms in ESDs subdir
     fhCkovP  = new TH2F("CkovP" , "#theta_{c}, [rad];P, [GeV]"   , 150,   0,  7  ,100, 0, 1)   ;
     fhSigP   = new TH2F("SigP"  ,"#sigma_{#theta_c} [mrad];[GeV]", 150,   0,  7  ,100, 0, 1)   ;
     fhMipXY  = new TH2F("MipXY" ,"mip position"                  , 260,   0,130  ,252, 0,126)  ;
     fhDifXY  = new TH2F("DifXY" ,"diff"                          , 200, -10, 10  ,200,-10,10)  ;
     fhPid[0] = new TH1F("PidE" ,"PID: e yellow #mu magenta"  ,100,0,1)                         ;
     fhPid[1] = new TH1F("PidMu","pid of #mu"                 ,100,0,1)                         ;
     fhPid[2] = new TH1F("PidPi","PID: #pi red K green p blue",100,0,1)                         ;
     fhPid[3] = new TH1F("PidK" ,"pid of K"                   ,100,0,1)                         ;
     fhPid[4] = new TH1F("PidP" ,"pid of p"                   ,100,0,1)                         ;
}

//____________________________________________________________________________
void AliHMPIDQADataMaker::MakeHits(TObject * data)
{
  //fills QA histos for Hits
  TClonesArray * hits = dynamic_cast<TClonesArray *>(data) ; 
  if (!hits){
    AliError("Wrong type of hits container") ; 
  } else {
    TIter next(hits); 
    AliHMPIDHit * hit ; 
    while ( (hit = dynamic_cast<AliHMPIDHit *>(next())) ) {
      if(hit->Pid()<500000) fhHitQdc->Fill(hit->Q()) ;
      if(hit->Pid()<500000) fhHitMap[hit->Ch()]->Fill(hit->LorsX(),hit->LorsY());
    }
  } 
}

//____________________________________________________________________________
void AliHMPIDQADataMaker::MakeDigits( TObject * data)
{
  //fills QA histos for Digits
  TObjArray *chambers = dynamic_cast<TObjArray*>(data);
  if ( !chambers) {
    AliError("Wrong type of digits container") ; 
  } else {
    for(Int_t i =0; i< chambers->GetEntries(); i++)
      {
	TClonesArray * digits = dynamic_cast<TClonesArray*>(chambers->At(i)); 
	fhDigChEvt->Fill(i,digits->GetEntriesFast()/(48.*80.*6.));
	TIter next(digits); 
	AliHMPIDDigit * digit; 
	while ( (digit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
	  fhDigPcEvt->Fill(10.*i+digit->Pc(),1./(48.*80.));
	  fhDigQ->Fill(digit->Q());
	}  
      }
  }
}

//____________________________________________________________________________
void AliHMPIDQADataMaker::MakeSDigits( TObject * data)
{
  //fills QA histos for SDigits
  TClonesArray * sdigits = dynamic_cast<TClonesArray *>(data) ; 
  if (!sdigits) {
    AliError("Wrong type of sdigits container") ; 
  } else {
    AliHMPIDDigit *ref = (AliHMPIDDigit *)sdigits->At(0);
    Float_t zero = ref->GetTrack(0); 
    TIter next(sdigits) ; 
    AliHMPIDDigit * sdigit ; 
    while ( (sdigit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
      fhSDigits->Fill(sdigit->Q()) ;
      if(zero == sdigit->GetTrack(0)) continue;
      else zero = sdigit->GetTrack(0);
    } 
  }
}

//____________________________________________________________________________
void AliHMPIDQADataMaker::MakeRecPoints(TTree * clustersTree)
{
  //fills QA histos for clusters

  TClonesArray *clusters = new TClonesArray("AliHMPIDCluster");
  for(int i=AliHMPIDParam::kMinCh;i<=AliHMPIDParam::kMaxCh;i++){
    TBranch *branch = clustersTree->GetBranch(Form("HMPID%d",i));
    branch->SetAddress(&clusters);
    branch->GetEntry(0);

    fhCluEvt->Fill(i,clusters->GetEntries());
    TIter next(clusters);
    AliHMPIDCluster *clu;
    while ( (clu = dynamic_cast<AliHMPIDCluster *>(next())) ) {;
      fhCluFlg->Fill(clu->Status());  fhCluChi2->Fill(clu->Chi2());  fhCluSize->Fill(clu->Size());
      fhCluQ->Fill(clu->Q()); 
      Int_t qCut=100;
      if(clu->Q()>qCut) {
	fhMipCluSize->SetTitle(Form("Mip cluster size at a Qcut = %i ADC",qCut));
	fhMipCluSize->Fill(clu->Size());
      }
    }
  }

  clusters->Delete();
  delete clusters;
}

//____________________________________________________________________________
void AliHMPIDQADataMaker::MakeESDs(AliESDEvent * esd)
{
  //fills QA histos for ESD
  for(Int_t iTrk = 0 ; iTrk < esd->GetNumberOfTracks() ; iTrk++){
    AliESDtrack *pTrk = esd->GetTrack(iTrk) ;
    fhCkovP->Fill(pTrk->GetP(),pTrk->GetHMPIDsignal());
    fhSigP->Fill( pTrk->GetP(),TMath::Sqrt(pTrk->GetHMPIDchi2()));
    Float_t xm,ym; Int_t q,np;  
    pTrk->GetHMPIDmip(xm,ym,q,np);                       //mip info
    fhMipXY->Fill(xm,ym);
    Float_t xRad,yRad,th,ph;        
    pTrk->GetHMPIDtrk(xRad,yRad,th,ph);              //track info at the middle of the radiator
    Float_t xPc = xRad+9.25*TMath::Tan(th)*TMath::Cos(ph); // temporar: linear extrapol (B=0!)
    Float_t yPc = yRad+9.25*TMath::Tan(th)*TMath::Sin(ph); // temporar:          "
    fhDifXY->Fill(xm-xPc,ym-yPc); //track info
    Double_t pid[5] ;      pTrk->GetHMPIDpid(pid) ;
    for(Int_t i = 0 ; i < 5 ; i++) fhPid[i]->Fill(pid[i]) ;
  }
}
Commit	Line	Data
b39da4dc	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	**************************************************************************/
	15
	16
	17	/* $Id$ */
	18
	19	//---
	20	// Produces the data needed to calculate the quality assurance.
	21	// All data must be mergeable objects.
	22	// A. Mastroserio
	23	//---
	24
	25	// --- ROOT system ---
	26	#include <TClonesArray.h>
	27	#include <TFile.h>
	28	#include <TH1F.h>
	29	#include <TH2F.h>
	30	#include <TH1I.h>
	31	#include <TDirectory.h>
	32	#include <Riostream.h>
	33	// --- Standard library ---
	34
	35	// --- AliRoot header files ---
	36	#include "AliESDCaloCluster.h"
	37	#include "AliESDEvent.h"
	38	#include "AliLog.h"
	39	#include "AliHMPIDDigit.h"
	40	#include "AliHMPIDHit.h"
	41	#include "AliHMPIDCluster.h"
b990acb3	42	#include "AliHMPIDQADataMaker.h"
b39da4dc	43
b990acb3	44	ClassImp(AliHMPIDQADataMaker)
b39da4dc	45
b39da4dc	46	//____________________________________________________________________________
b990acb3	47	AliHMPIDQADataMaker::AliHMPIDQADataMaker() :
b990acb3	48	AliQADataMaker(AliQA::GetDetName(AliQA::kHMPID), "HMPID Quality Assurance Data Maker"),
b39da4dc	49	fhHitQdc(0x0),
	50	fhSDigits(0x0),
	51	fhDigPcEvt(0x0),
	52	fhDigChEvt(0x0),
	53	fhDigQ(0x0),
	54	fhCluEvt(0x0),
	55	fhCluChi2(0x0),
	56	fhCluQ(0x0),
	57	fhCluFlg(0x0),
	58	fhCluSize(0x0),
	59	fhMipCluSize(0x0),
	60	fhCkovP(0x0),
	61	fhSigP(0x0),
	62	fhMipXY(0x0),
	63	fhDifXY(0x0)
	64	{
	65	// ctor
	66	for(Int_t i=0; i<7; i++) fhHitMap[i]=0x0;
	67	for(Int_t j=0; j<5; j++) fhPid[j]=0x0;
a01035e3	68	// fDetectorDir = fOutput->GetDirectory(GetName()) ;
	69	// if (!fDetectorDir)
	70	// fDetectorDir = fOutput->mkdir(GetName()) ;
b39da4dc	71	}
	72
	73	//____________________________________________________________________________
b990acb3	74	AliHMPIDQADataMaker::AliHMPIDQADataMaker(const AliHMPIDQADataMaker& qadm) :
b990acb3	75	AliQADataMaker(),
b39da4dc	76	fhHitQdc(qadm.fhHitQdc),
	77	fhSDigits(qadm.fhSDigits),
	78	fhDigPcEvt(qadm.fhDigPcEvt),
	79	fhDigChEvt(qadm.fhDigChEvt),
	80	fhDigQ(qadm.fhDigQ),
	81	fhCluEvt(qadm.fhCluEvt),
	82	fhCluChi2(qadm.fhCluChi2),
	83	fhCluQ(qadm.fhCluQ),
	84	fhCluFlg(qadm.fhCluFlg),
	85	fhCluSize(qadm.fhCluSize),
	86	fhMipCluSize(qadm.fhMipCluSize),
	87	fhCkovP(qadm.fhCkovP),
	88	fhSigP(qadm.fhSigP),
	89	fhMipXY(qadm.fhMipXY),
	90	fhDifXY(qadm.fhDifXY)
	91	{
	92	//copy ctor
	93	for(Int_t i=0; i<7; i++) fhHitMap[i]=qadm.fhHitMap[i];
	94	for(Int_t j=0; j<5; j++) fhPid[j]=qadm.fhPid[j];
	95
	96	SetName((const char*)qadm.GetName()) ;
	97	SetTitle((const char*)qadm.GetTitle());
	98	}
	99
	100	//__________________________________________________________________
b990acb3	101	AliHMPIDQADataMaker& AliHMPIDQADataMaker::operator = (const AliHMPIDQADataMaker& qadm )
b39da4dc	102	{
b39da4dc	103	// Equal operator.
b990acb3	104	this->~AliHMPIDQADataMaker();
b990acb3	105	new(this) AliHMPIDQADataMaker(qadm);
b39da4dc	106	return *this;
	107	}
	108
	109	//____________________________________________________________________________
b990acb3	110	void AliHMPIDQADataMaker::InitHits()
b39da4dc	111	{
	112	// create Hits histograms in Hits subdir
	113	fhHitQdc=new TH1F("HitQdc","HMPID Hit Qdc all chamber;QDC",500,0,4000);
	114	for(Int_t iCh=0;iCh<7;iCh++) fhHitMap[iCh]=new TH2F(Form("HMPID HitMap%i",iCh),Form("Ch%i;x_{Hit};y_{Hit}",iCh),162,-1,161,146,-1,145);
	115	}
	116
	117	//____________________________________________________________________________
b990acb3	118	void AliHMPIDQADataMaker::InitDigits()
b39da4dc	119	{
	120	// create Digits histograms in Digits subdir
	121	fhDigPcEvt=new TH1F("hDigPcEvt","PC occupancy",156,-1,77);
	122	fhDigChEvt=new TH1F("hDigChEvt","Chamber occupancy",32,-1,7);
	123	fhDigQ =new TH1F("Q "," digit charge ",3000,0,3000);
	124	}
	125
	126	//____________________________________________________________________________
b990acb3	127	void AliHMPIDQADataMaker::InitSDigits()
b39da4dc	128	{
	129	// create SDigits histograms in SDigits subdir
	130	fhSDigits = new TH1F("hHmpidSDigits", "SDigits Q distribution in HMPID", 500, 0., 5000.) ;
	131	}
	132
	133	//____________________________________________________________________________
	134
b990acb3	135	void AliHMPIDQADataMaker::InitRecPoints()
b39da4dc	136	{
	137	// create cluster histograms in RecPoint subdir
	138	fhCluEvt=new TH1F("CluPerEvt","# clusters per chamber",16,-1,7);
	139	fhCluChi2 =new TH1F("CluChi2" ,"Chi2 " ,1000,0,100);
	140	fhCluQ =new TH1F("CluQ" ,"Cluster charge" ,3000,0,3000);
	141	fhCluFlg =new TH1F("CluFlg" ,"Cluster flag" ,14,-1.5,12.5);
	142	fhCluSize =new TH1F("CluSize" ,"Raw cluster size ",100,0,100);
	143	fhMipCluSize =new TH1F("MipCluSize" ,"Mip cluster size ",100,0,100);
	144	}
	145	//____________________________________________________________________________
b990acb3	146	void AliHMPIDQADataMaker::InitESDs()
b39da4dc	147	{
	148	//create ESDs histograms in ESDs subdir
	149	fhCkovP = new TH2F("CkovP" , "#theta_{c}, [rad];P, [GeV]" , 150, 0, 7 ,100, 0, 1) ;
	150	fhSigP = new TH2F("SigP" ,"#sigma_{#theta_c} [mrad];[GeV]", 150, 0, 7 ,100, 0, 1) ;
	151	fhMipXY = new TH2F("MipXY" ,"mip position" , 260, 0,130 ,252, 0,126) ;
	152	fhDifXY = new TH2F("DifXY" ,"diff" , 200, -10, 10 ,200,-10,10) ;
	153	fhPid[0] = new TH1F("PidE" ,"PID: e yellow #mu magenta" ,100,0,1) ;
	154	fhPid[1] = new TH1F("PidMu","pid of #mu" ,100,0,1) ;
	155	fhPid[2] = new TH1F("PidPi","PID: #pi red K green p blue",100,0,1) ;
	156	fhPid[3] = new TH1F("PidK" ,"pid of K" ,100,0,1) ;
	157	fhPid[4] = new TH1F("PidP" ,"pid of p" ,100,0,1) ;
	158	}
	159
	160	//____________________________________________________________________________
b990acb3	161	void AliHMPIDQADataMaker::MakeHits(TObject * data)
b39da4dc	162	{
b39da4dc	163	//fills QA histos for Hits
d76c31f4	164	TClonesArray * hits = dynamic_cast<TClonesArray *>(data) ;
	165	if (!hits){
	166	AliError("Wrong type of hits container") ;
	167	} else {
	168	TIter next(hits);
	169	AliHMPIDHit * hit ;
	170	while ( (hit = dynamic_cast<AliHMPIDHit *>(next())) ) {
	171	if(hit->Pid()<500000) fhHitQdc->Fill(hit->Q()) ;
	172	if(hit->Pid()<500000) fhHitMap[hit->Ch()]->Fill(hit->LorsX(),hit->LorsY());
	173	}
	174	}
	175	}
	176
b39da4dc	177	//____________________________________________________________________________
b990acb3	178	void AliHMPIDQADataMaker::MakeDigits( TObject * data)
b39da4dc	179	{
b39da4dc	180	//fills QA histos for Digits
d76c31f4	181	TObjArray chambers = dynamic_cast<TObjArray>(data);
	182	if ( !chambers) {
	183	AliError("Wrong type of digits container") ;
	184	} else {
	185	for(Int_t i =0; i< chambers->GetEntries(); i++)
	186	{
	187	TClonesArray * digits = dynamic_cast<TClonesArray*>(chambers->At(i));
	188	fhDigChEvt->Fill(i,digits->GetEntriesFast()/(48.80.6.));
	189	TIter next(digits);
	190	AliHMPIDDigit * digit;
	191	while ( (digit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
	192	fhDigPcEvt->Fill(10.i+digit->Pc(),1./(48.80.));
	193	fhDigQ->Fill(digit->Q());
	194	}
	195	}
	196	}
b39da4dc	197	}
	198
	199	//____________________________________________________________________________
b990acb3	200	void AliHMPIDQADataMaker::MakeSDigits( TObject * data)
b39da4dc	201	{
b39da4dc	202	//fills QA histos for SDigits
d76c31f4	203	TClonesArray * sdigits = dynamic_cast<TClonesArray *>(data) ;
	204	if (!sdigits) {
	205	AliError("Wrong type of sdigits container") ;
	206	} else {
	207	AliHMPIDDigit ref = (AliHMPIDDigit )sdigits->At(0);
	208	Float_t zero = ref->GetTrack(0);
	209	TIter next(sdigits) ;
	210	AliHMPIDDigit * sdigit ;
	211	while ( (sdigit = dynamic_cast<AliHMPIDDigit *>(next())) ) {
	212	fhSDigits->Fill(sdigit->Q()) ;
	213	if(zero == sdigit->GetTrack(0)) continue;
	214	else zero = sdigit->GetTrack(0);
	215	}
	216	}
b39da4dc	217	}
b39da4dc	218
d76c31f4	219	//____________________________________________________________________________
b990acb3	220	void AliHMPIDQADataMaker::MakeRecPoints(TTree * clustersTree)
b39da4dc	221	{
b39da4dc	222	//fills QA histos for clusters
b39da4dc	223
d76c31f4	224	TClonesArray *clusters = new TClonesArray("AliHMPIDCluster");
	225	for(int i=AliHMPIDParam::kMinCh;i<=AliHMPIDParam::kMaxCh;i++){
	226	TBranch *branch = clustersTree->GetBranch(Form("HMPID%d",i));
	227	branch->SetAddress(&clusters);
	228	branch->GetEntry(0);
	229
b39da4dc	230	fhCluEvt->Fill(i,clusters->GetEntries());
	231	TIter next(clusters);
	232	AliHMPIDCluster *clu;
	233	while ( (clu = dynamic_cast<AliHMPIDCluster *>(next())) ) {;
d76c31f4	234	fhCluFlg->Fill(clu->Status()); fhCluChi2->Fill(clu->Chi2()); fhCluSize->Fill(clu->Size());
	235	fhCluQ->Fill(clu->Q());
	236	Int_t qCut=100;
	237	if(clu->Q()>qCut) {
	238	fhMipCluSize->SetTitle(Form("Mip cluster size at a Qcut = %i ADC",qCut));
	239	fhMipCluSize->Fill(clu->Size());
	240	}
b39da4dc	241	}
b39da4dc	242	}
d76c31f4	243
	244	clusters->Delete();
	245	delete clusters;
b39da4dc	246	}
	247
	248	//____________________________________________________________________________
b990acb3	249	void AliHMPIDQADataMaker::MakeESDs(AliESDEvent * esd)
b39da4dc	250	{
b39da4dc	251	//fills QA histos for ESD
b39da4dc	252	for(Int_t iTrk = 0 ; iTrk < esd->GetNumberOfTracks() ; iTrk++){
	253	AliESDtrack *pTrk = esd->GetTrack(iTrk) ;
	254	fhCkovP->Fill(pTrk->GetP(),pTrk->GetHMPIDsignal());
	255	fhSigP->Fill( pTrk->GetP(),TMath::Sqrt(pTrk->GetHMPIDchi2()));
	256	Float_t xm,ym; Int_t q,np;
	257	pTrk->GetHMPIDmip(xm,ym,q,np); //mip info
	258	fhMipXY->Fill(xm,ym);
	259	Float_t xRad,yRad,th,ph;
	260	pTrk->GetHMPIDtrk(xRad,yRad,th,ph); //track info at the middle of the radiator
	261	Float_t xPc = xRad+9.25TMath::Tan(th)TMath::Cos(ph); // temporar: linear extrapol (B=0!)
	262	Float_t yPc = yRad+9.25TMath::Tan(th)TMath::Sin(ph); // temporar: "
	263	fhDifXY->Fill(xm-xPc,ym-yPc); //track info
	264	Double_t pid[5] ; pTrk->GetHMPIDpid(pid) ;
	265	for(Int_t i = 0 ; i < 5 ; i++) fhPid[i]->Fill(pid[i]) ;
	266	}
	267	}
	268