[u/mrichter/AliRoot.git] / HMPID / AliHMPIDQualAssDataMaker.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 "AliHMPIDQualAssDataMaker.h"

ClassImp(AliHMPIDQualAssDataMaker)
           
//____________________________________________________________________________ 
  AliHMPIDQualAssDataMaker::AliHMPIDQualAssDataMaker() : 
  AliQualAssDataMaker(AliQualAss::GetDetName(AliQualAss::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()) ;  
}

//____________________________________________________________________________ 
AliHMPIDQualAssDataMaker::AliHMPIDQualAssDataMaker(const AliHMPIDQualAssDataMaker& qadm) :
  AliQualAssDataMaker(), 
  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()); 
}

//__________________________________________________________________
AliHMPIDQualAssDataMaker& AliHMPIDQualAssDataMaker::operator = (const AliHMPIDQualAssDataMaker& qadm )
{
  // Equal operator.
  this->~AliHMPIDQualAssDataMaker();
  new(this) AliHMPIDQualAssDataMaker(qadm);
  return *this;
}
 
//____________________________________________________________________________ 
void AliHMPIDQualAssDataMaker::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 AliHMPIDQualAssDataMaker::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 AliHMPIDQualAssDataMaker::InitSDigits()
{
  // create SDigits histograms in SDigits subdir
      fhSDigits     = new TH1F("hHmpidSDigits",    "SDigits Q  distribution in HMPID",  500, 0., 5000.) ; 
}

//____________________________________________________________________________ 

void AliHMPIDQualAssDataMaker::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 AliHMPIDQualAssDataMaker::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 AliHMPIDQualAssDataMaker::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 AliHMPIDQualAssDataMaker::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 AliHMPIDQualAssDataMaker::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 AliHMPIDQualAssDataMaker::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 AliHMPIDQualAssDataMaker::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"
	42	#include "AliHMPIDQualAssDataMaker.h"
	43
	44	ClassImp(AliHMPIDQualAssDataMaker)
	45
	46	//____________________________________________________________________________
	47	AliHMPIDQualAssDataMaker::AliHMPIDQualAssDataMaker() :
	48	AliQualAssDataMaker(AliQualAss::GetDetName(AliQualAss::kHMPID), "HMPID Quality Assurance Data Maker"),
	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;
2943bde2	68	// fDetectorDir = fOutput->GetDirectory(GetName()) ;
	69	// if (!fDetectorDir)
	70	// fDetectorDir = fOutput->mkdir(GetName()) ;
b39da4dc	71	}
	72
	73	//____________________________________________________________________________
	74	AliHMPIDQualAssDataMaker::AliHMPIDQualAssDataMaker(const AliHMPIDQualAssDataMaker& qadm) :
	75	AliQualAssDataMaker(),
	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	//__________________________________________________________________
	101	AliHMPIDQualAssDataMaker& AliHMPIDQualAssDataMaker::operator = (const AliHMPIDQualAssDataMaker& qadm )
	102	{
	103	// Equal operator.
	104	this->~AliHMPIDQualAssDataMaker();
	105	new(this) AliHMPIDQualAssDataMaker(qadm);
	106	return *this;
	107	}
	108
	109	//____________________________________________________________________________
	110	void AliHMPIDQualAssDataMaker::InitHits()
	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	//____________________________________________________________________________
	118	void AliHMPIDQualAssDataMaker::InitDigits()
	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	//____________________________________________________________________________
	127	void AliHMPIDQualAssDataMaker::InitSDigits()
	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
135	void AliHMPIDQualAssDataMaker::InitRecPoints()
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	//____________________________________________________________________________
146	void AliHMPIDQualAssDataMaker::InitESDs()
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	//____________________________________________________________________________
d76c31f4	161	void AliHMPIDQualAssDataMaker::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	//____________________________________________________________________________
d76c31f4	178	void AliHMPIDQualAssDataMaker::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	//____________________________________________________________________________
d76c31f4	200	void AliHMPIDQualAssDataMaker::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	//____________________________________________________________________________
d76c31f4	220	void AliHMPIDQualAssDataMaker::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	//____________________________________________________________________________
d76c31f4	249	void AliHMPIDQualAssDataMaker::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