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

#include "AliHMPID.h"
#include "AliHMPIDHit.h"   //OccupancyPrint(), HitQa()
#include "AliHMPIDDigit.h" //OccupancyPrint()
#include <TParticle.h>  //SummaryOfEvent(), HitQa()
#include <TBenchmark.h>  //HitQA()
#include <TPDGCode.h>    //HitQA()
#include <AliStack.h>   //OccupancyPrint(), SummaryOfEvent(), HitQa()
#include <AliRun.h>     //HitQa() 
#include <AliMC.h>       //ctor
#include <AliHeader.h>
#include <AliGenEventHeader.h>
#include <AliGenHijingEventHeader.h>
#include <TH1F.h>        //HitQA()
#include <AliLog.h>      //in many methods to print AliInfo 
ClassImp(AliHMPID)    
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
AliHMPID::AliHMPID(const char *name, const char *title):AliDetector(name,title),fSdi(0),fDig(0),fClu(0)
{
//Named ctor
  AliDebug(1,"Start.");
//AliDetector ctor deals with Hits and Digits (reset them to 0, does not create them)
  HitCreate();          gAlice->GetMCApp()->AddHitList(fHits);
  AliDebug(1,"Stop.");
}//AliHMPID::AliHMPID(const char *name, const char *title)
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
AliHMPID::~AliHMPID()
{
//dtor
  AliDebug(1,"Start.");

  
  if(fHits)      delete fHits;
  if(fDigits)    delete fDigits;
  if(fSdi)       delete fSdi;
  if(fDig)      {fDig->Delete();   delete fDig;}
  if(fClu)      {fClu->Delete();   delete fClu;}
  AliDebug(1,"Stop.");    
}//AliHMPID::~AliHMPID()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPID::MakeBranch(Option_t* option)
{
//Create Tree branches for the HMPID.
  AliDebug(1,Form("Start with option= %s.",option));
    
  const Int_t kBufSize = 4000;
      
  const char *cH = strstr(option,"H");
  const char *cD = strstr(option,"D");
  const char *cR = strstr(option,"R");
  const char *cS = strstr(option,"S");

  if(cH&&         TreeH()){HitCreate();                       MakeBranchInTree(         TreeH(),     "HMPID"     ,&fHits       ,kBufSize,0);}
  if(cS&&fLoader->TreeS()){SdiCreate();                       MakeBranchInTree(fLoader->TreeS(),     "HMPID"     ,&fSdi        ,kBufSize,0);}
  if(cD&&fLoader->TreeD()){DigCreate();for(Int_t i=0;i<7;i++) MakeBranchInTree(fLoader->TreeD(),Form("HMPID%d",i),&((*fDig)[i]),kBufSize,0);}
  if(cR&&fLoader->TreeR()){CluCreate();for(Int_t i=0;i<7;i++) MakeBranchInTree(fLoader->TreeR(),Form("HMPID%d",i),&((*fClu)[i]),kBufSize,0);}   
  
  AliDebug(1,"Stop.");   
}//void AliHMPID::MakeBranch(Option_t* option)
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPID::SetTreeAddress()
{
//Set branch address for the Hits and Digits Tree.
  AliDebug(1,"Start.");
  if(fLoader->TreeH() && fLoader->TreeH()->GetBranch("HMPID" )){HitCreate();                      fLoader->TreeH()->SetBranchAddress(     "HMPID"     ,&fHits       );}
  if(fLoader->TreeS() && fLoader->TreeS()->GetBranch("HMPID" )){SdiCreate();                      fLoader->TreeS()->SetBranchAddress(     "HMPID"     ,&fSdi        );}
  if(fLoader->TreeD() && fLoader->TreeD()->GetBranch("HMPID0")){DigCreate(); for(int i=0;i<7;i++) fLoader->TreeD()->SetBranchAddress(Form("HMPID%d",i),&((*fDig)[i]));}
  if(fLoader->TreeR() && fLoader->TreeR()->GetBranch("HMPID0")){CluCreate(); for(int i=0;i<7;i++) fLoader->TreeR()->SetBranchAddress(Form("HMPID%d",i),&((*fClu)[i]));}
  AliDebug(1,"Stop.");
}//void AliHMPID::SetTreeAddress()
//__________________________________________________________________________________________________
// AliHMPIDHit* AliHMPID::Hit(Int_t tid)const
// {
// // Search for the first HMPID hit belonging to the given tid
//   GetLoader()->LoadHits();
//   for(Int_t iPrimN=0;iPrimN<GetLoader()->TreeH()->GetEntries();iPrimN++){//prims loop      
//     GetLoader()->TreeH()->GetEntry(iPrimN);
//     for(Int_t iHitN=0;iHitN<Hits()->GetEntries();iHitN++){
//       AliHMPIDHit *pHit=(AliHMPIDHit*)Hits()->At(iHitN);
//       if(tid==pHit->Track()) {GetLoader()->UnloadHits();return pHit;}
//     }//hits
//   }//prims loop
//   GetLoader()->UnloadHits();
//   return 0;
// }
//__________________________________________________________________________________________________
void AliHMPID::HitPrint(Int_t iEvtN)const
{
//Prints a list of HMPID hits for a given event. Default is event number 0.
  if(GetLoader()->GetRunLoader()->GetEvent(iEvtN)) return;    
  AliInfo(Form("List of HMPID hits for event %i",iEvtN));
  if(GetLoader()->LoadHits()) return;
  
  Int_t iTotalHits=0;
  for(Int_t iPrimN=0;iPrimN<GetLoader()->TreeH()->GetEntries();iPrimN++){//prims loop
    GetLoader()->TreeH()->GetEntry(iPrimN);      
    Hits()->Print();
    iTotalHits+=Hits()->GetEntries();
  }
  GetLoader()->UnloadHits();
  AliInfo(Form("totally %i hits",iTotalHits));
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPID::SdiPrint(Int_t iEvt)const
{
//prints a list of HMPID sdigits  for a given event
  if(GetLoader()->GetRunLoader()->GetEvent(iEvt)) return;    
  Info("PrintSDigits","List of HMPID sdigits for event %i",iEvt);
  if(GetLoader()->LoadSDigits()) return;
  
  GetLoader()->TreeS()->GetEntry(0);
  SdiLst()->Print();
  GetLoader()->UnloadSDigits();
  Printf("totally %i sdigits",SdiLst()->GetEntries());
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPID::DigPrint(Int_t iEvt)const
{
//prints a list of HMPID digits  for a given event
  if(GetLoader()->GetRunLoader()->GetEvent(iEvt)) return;    
  Printf("List of HMPID digits for event %i",iEvt);
  if(GetLoader()->LoadDigits()) return;
  
  GetLoader()->TreeD()->GetEntry(0);
  DigLst()->Print();
  GetLoader()->UnloadDigits();
  Printf("totally %i Digits",DigLst()->GetEntries());
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPID::OccupancyPrint(Int_t iEvtNreq)
{
//prints occupancy for each chamber in a given event
  Int_t iEvtNmin,iEvtNmax;
  if(iEvtNreq==-1){
    iEvtNmin=0;
    iEvtNmax=gAlice->GetEventsPerRun();
  } else { 
    iEvtNmin=iEvtNreq;iEvtNmax=iEvtNreq+1;
  }
    
  if(GetLoader()->GetRunLoader()->LoadHeader()) return;    
  if(GetLoader()->GetRunLoader()->LoadKinematics()) return;    
  
//  Info("Occupancy","for event %i",iEvtN);
  if(GetLoader()->LoadHits()) return;
  if(GetLoader()->LoadDigits()) return;

  
  for(Int_t iEvtN=iEvtNmin;iEvtN<iEvtNmax;iEvtN++){    
    Int_t nDigCh[7]={0,0,0,0,0,0,0};  
    Int_t iChHits[7]={0,0,0,0,0,0,0};
    Int_t nPrim[7]={0,0,0,0,0,0,0};
    Int_t nSec[7]={0,0,0,0,0,0,0};
    AliInfo(Form("events processed %i",iEvtN));
    if(GetLoader()->GetRunLoader()->GetEvent(iEvtN)) return;    
    AliStack *pStack = GetLoader()->GetRunLoader()->Stack();
    for(Int_t iPrimN=0;iPrimN<GetLoader()->TreeH()->GetEntries();iPrimN++){//prims loop
      GetLoader()->TreeH()->GetEntry(iPrimN);      
      for(Int_t iHitN=0;iHitN<Hits()->GetEntries();iHitN++){
        AliHMPIDHit *pHit = (AliHMPIDHit*)Hits()->At(iHitN);
        if(pHit->E()>0){
          iChHits[pHit->Ch()]++;
          if(pStack->Particle(pHit->GetTrack())->Rho()<0.01) nPrim[pHit->Ch()]++;else nSec[pHit->Ch()]++;
        }
      }
    }
    
    GetLoader()->TreeD()->GetEntry(0);
    for(Int_t iCh=0;iCh<7;iCh++){
      for(Int_t iDig=0;iDig<DigLst(iCh)->GetEntries();iDig++){
        AliHMPIDDigit *pDig=(AliHMPIDDigit*)DigLst(iCh)->At(iDig);
        nDigCh[pDig->Ch()]++;
      }
      Printf("Occupancy for chamber %i = %4.2f %% and charged prim tracks %i and sec. tracks %i with total %i",
        iCh,Float_t(nDigCh[iCh])*100/AliHMPIDDigit::kPadAll,nPrim[iCh],nSec[iCh],iChHits[iCh]);
    }      
    
    
  }//events loop
  GetLoader()->UnloadHits();
  GetLoader()->UnloadDigits();
  GetLoader()->GetRunLoader()->UnloadHeader();    
  GetLoader()->GetRunLoader()->UnloadKinematics();    
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPID::CluPrint(Int_t iEvtN)const
{
//prints a list of HMPID clusters  for a given event
  Printf("List of HMPID clusters for event %i",iEvtN);
  GetLoader()->GetRunLoader()->GetEvent(iEvtN);    
  if(GetLoader()->LoadRecPoints()) return;
  
  Int_t iCluCnt=0;
  GetLoader()->TreeR()->GetEntry(0);
  for(Int_t iCh=0;iCh<7;iCh++){
    TClonesArray *pCluLst=(TClonesArray*)fClu->At(iCh);    iCluCnt+=pCluLst->GetEntries();    pCluLst->Print();
  }
  GetLoader()->UnloadRecPoints();
  Printf("totally %i clusters for event %i",iCluCnt,iEvtN);
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void AliHMPID::SummaryOfEvent(Int_t iEvtN) const
{
//prints a summary for a given event
  AliInfo(Form("Summary of event %i",iEvtN));
  GetLoader()->GetRunLoader()->GetEvent(iEvtN);    
  if(GetLoader()->GetRunLoader()->LoadHeader()) return;
  if(GetLoader()->GetRunLoader()->LoadKinematics()) return;
  AliStack *pStack=GetLoader()->GetRunLoader()->Stack();
  
  AliGenEventHeader* pGenHeader =  gAlice->GetHeader()->GenEventHeader();
  if(pGenHeader->InheritsFrom("AliGenHijingEventHeader")) {
    AliInfo(Form(" Hijing event with impact parameter b = %.2f (fm)",((AliGenHijingEventHeader*) pGenHeader)->ImpactParameter()));
  }
  Int_t nChargedPrimaries=0;
  for(Int_t i=0;i<pStack->GetNtrack();i++) {
    TParticle *pParticle = pStack->Particle(i);
    if(pParticle->IsPrimary()&&pParticle->GetPDG()->Charge()!=0) nChargedPrimaries++;
    }
  AliInfo(Form("Total number of         primaries %i",pStack->GetNprimary()));
  AliInfo(Form("Total number of charged primaries %i",nChargedPrimaries));
  AliInfo(Form("Total n. of tracks in stack(+sec) %i",pStack->GetNtrack()));
  GetLoader()->GetRunLoader()->UnloadHeader();
  GetLoader()->GetRunLoader()->UnloadKinematics();
}
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
Commit	Line	Data
d3da6dc4	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	#include "AliHMPID.h"
	17	#include "AliHMPIDHit.h" //OccupancyPrint(), HitQa()
	18	#include "AliHMPIDDigit.h" //OccupancyPrint()
	19	#include <TParticle.h> //SummaryOfEvent(), HitQa()
	20	#include <TBenchmark.h> //HitQA()
	21	#include <TPDGCode.h> //HitQA()
	22	#include <AliStack.h> //OccupancyPrint(), SummaryOfEvent(), HitQa()
	23	#include <AliRun.h> //HitQa()
	24	#include <AliMC.h> //ctor
	25	#include <AliHeader.h>
	26	#include <AliGenEventHeader.h>
	27	#include <AliGenHijingEventHeader.h>
	28	#include <TH1F.h> //HitQA()
	29	#include <AliLog.h> //in many methods to print AliInfo
	30	ClassImp(AliHMPID)
	31	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	32	AliHMPID::AliHMPID(const char name, const char title):AliDetector(name,title),fSdi(0),fDig(0),fClu(0)
	33	{
	34	//Named ctor
	35	AliDebug(1,"Start.");
	36	//AliDetector ctor deals with Hits and Digits (reset them to 0, does not create them)
	37	HitCreate(); gAlice->GetMCApp()->AddHitList(fHits);
	38	AliDebug(1,"Stop.");
	39	}//AliHMPID::AliHMPID(const char name, const char title)
	40	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	41	AliHMPID::~AliHMPID()
	42	{
	43	//dtor
	44	AliDebug(1,"Start.");
	45
	46
	47	if(fHits) delete fHits;
	48	if(fDigits) delete fDigits;
	49	if(fSdi) delete fSdi;
	50	if(fDig) {fDig->Delete(); delete fDig;}
	51	if(fClu) {fClu->Delete(); delete fClu;}
	52	AliDebug(1,"Stop.");
	53	}//AliHMPID::~AliHMPID()
	54	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
	55	void AliHMPID::MakeBranch(Option_t* option)
	56	{
	57	//Create Tree branches for the HMPID.
	58	AliDebug(1,Form("Start with option= %s.",option));
	59
	60	const Int_t kBufSize = 4000;
	61
	62	const char *cH = strstr(option,"H");
	63	const char *cD = strstr(option,"D");
	64	const char *cR = strstr(option,"R");
65	const char *cS = strstr(option,"S");
66
67	if(cH&& TreeH()){HitCreate(); MakeBranchInTree( TreeH(), "HMPID" ,&fHits ,kBufSize,0);}
68	if(cS&&fLoader->TreeS()){SdiCreate(); MakeBranchInTree(fLoader->TreeS(), "HMPID" ,&fSdi ,kBufSize,0);}
69	if(cD&&fLoader->TreeD()){DigCreate();for(Int_t i=0;i<7;i++) MakeBranchInTree(fLoader->TreeD(),Form("HMPID%d",i),&((*fDig)[i]),kBufSize,0);}
70	if(cR&&fLoader->TreeR()){CluCreate();for(Int_t i=0;i<7;i++) MakeBranchInTree(fLoader->TreeR(),Form("HMPID%d",i),&((*fClu)[i]),kBufSize,0);}
71
72	AliDebug(1,"Stop.");
73	}//void AliHMPID::MakeBranch(Option_t* option)
74	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
75	void AliHMPID::SetTreeAddress()
76	{
77	//Set branch address for the Hits and Digits Tree.
78	AliDebug(1,"Start.");
79	if(fLoader->TreeH() && fLoader->TreeH()->GetBranch("HMPID" )){HitCreate(); fLoader->TreeH()->SetBranchAddress( "HMPID" ,&fHits );}
80	if(fLoader->TreeS() && fLoader->TreeS()->GetBranch("HMPID" )){SdiCreate(); fLoader->TreeS()->SetBranchAddress( "HMPID" ,&fSdi );}
81	if(fLoader->TreeD() && fLoader->TreeD()->GetBranch("HMPID0")){DigCreate(); for(int i=0;i<7;i++) fLoader->TreeD()->SetBranchAddress(Form("HMPID%d",i),&((*fDig)[i]));}
82	if(fLoader->TreeR() && fLoader->TreeR()->GetBranch("HMPID0")){CluCreate(); for(int i=0;i<7;i++) fLoader->TreeR()->SetBranchAddress(Form("HMPID%d",i),&((*fClu)[i]));}
83	AliDebug(1,"Stop.");
84	}//void AliHMPID::SetTreeAddress()
85	//__________________________________________________________________________________________________
86	// AliHMPIDHit* AliHMPID::Hit(Int_t tid)const
87	// {
88	// // Search for the first HMPID hit belonging to the given tid
89	// GetLoader()->LoadHits();
90	// for(Int_t iPrimN=0;iPrimN<GetLoader()->TreeH()->GetEntries();iPrimN++){//prims loop
91	// GetLoader()->TreeH()->GetEntry(iPrimN);
92	// for(Int_t iHitN=0;iHitN<Hits()->GetEntries();iHitN++){
93	// AliHMPIDHit pHit=(AliHMPIDHit)Hits()->At(iHitN);
94	// if(tid==pHit->Track()) {GetLoader()->UnloadHits();return pHit;}
95	// }//hits
96	// }//prims loop
97	// GetLoader()->UnloadHits();
98	// return 0;
99	// }
100	//__________________________________________________________________________________________________
101	void AliHMPID::HitPrint(Int_t iEvtN)const
102	{
103	//Prints a list of HMPID hits for a given event. Default is event number 0.
104	if(GetLoader()->GetRunLoader()->GetEvent(iEvtN)) return;
105	AliInfo(Form("List of HMPID hits for event %i",iEvtN));
106	if(GetLoader()->LoadHits()) return;
107
108	Int_t iTotalHits=0;
109	for(Int_t iPrimN=0;iPrimN<GetLoader()->TreeH()->GetEntries();iPrimN++){//prims loop
110	GetLoader()->TreeH()->GetEntry(iPrimN);
111	Hits()->Print();
112	iTotalHits+=Hits()->GetEntries();
113	}
114	GetLoader()->UnloadHits();
115	AliInfo(Form("totally %i hits",iTotalHits));
116	}
117	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
118	void AliHMPID::SdiPrint(Int_t iEvt)const
119	{
120	//prints a list of HMPID sdigits for a given event
121	if(GetLoader()->GetRunLoader()->GetEvent(iEvt)) return;
122	Info("PrintSDigits","List of HMPID sdigits for event %i",iEvt);
123	if(GetLoader()->LoadSDigits()) return;
124
125	GetLoader()->TreeS()->GetEntry(0);
126	SdiLst()->Print();
127	GetLoader()->UnloadSDigits();
128	Printf("totally %i sdigits",SdiLst()->GetEntries());
129	}
130	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
131	void AliHMPID::DigPrint(Int_t iEvt)const
132	{
133	//prints a list of HMPID digits for a given event
134	if(GetLoader()->GetRunLoader()->GetEvent(iEvt)) return;
135	Printf("List of HMPID digits for event %i",iEvt);
136	if(GetLoader()->LoadDigits()) return;
137
138	GetLoader()->TreeD()->GetEntry(0);
139	DigLst()->Print();
140	GetLoader()->UnloadDigits();
141	Printf("totally %i Digits",DigLst()->GetEntries());
142	}
143	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
144	void AliHMPID::OccupancyPrint(Int_t iEvtNreq)
145	{
146	//prints occupancy for each chamber in a given event
147	Int_t iEvtNmin,iEvtNmax;
148	if(iEvtNreq==-1){
149	iEvtNmin=0;
150	iEvtNmax=gAlice->GetEventsPerRun();
151	} else {
152	iEvtNmin=iEvtNreq;iEvtNmax=iEvtNreq+1;
153	}
154
155	if(GetLoader()->GetRunLoader()->LoadHeader()) return;
156	if(GetLoader()->GetRunLoader()->LoadKinematics()) return;
157
158	// Info("Occupancy","for event %i",iEvtN);
159	if(GetLoader()->LoadHits()) return;
160	if(GetLoader()->LoadDigits()) return;
161
162
163	for(Int_t iEvtN=iEvtNmin;iEvtN<iEvtNmax;iEvtN++){
164	Int_t nDigCh[7]={0,0,0,0,0,0,0};
165	Int_t iChHits[7]={0,0,0,0,0,0,0};
166	Int_t nPrim[7]={0,0,0,0,0,0,0};
167	Int_t nSec[7]={0,0,0,0,0,0,0};
168	AliInfo(Form("events processed %i",iEvtN));
169	if(GetLoader()->GetRunLoader()->GetEvent(iEvtN)) return;
170	AliStack *pStack = GetLoader()->GetRunLoader()->Stack();
171	for(Int_t iPrimN=0;iPrimN<GetLoader()->TreeH()->GetEntries();iPrimN++){//prims loop
172	GetLoader()->TreeH()->GetEntry(iPrimN);
173	for(Int_t iHitN=0;iHitN<Hits()->GetEntries();iHitN++){
174	AliHMPIDHit pHit = (AliHMPIDHit)Hits()->At(iHitN);
175	if(pHit->E()>0){
176	iChHits[pHit->Ch()]++;
177	if(pStack->Particle(pHit->GetTrack())->Rho()<0.01) nPrim[pHit->Ch()]++;else nSec[pHit->Ch()]++;
178	}
179	}
180	}
181
182	GetLoader()->TreeD()->GetEntry(0);
183	for(Int_t iCh=0;iCh<7;iCh++){
184	for(Int_t iDig=0;iDig<DigLst(iCh)->GetEntries();iDig++){
185	AliHMPIDDigit pDig=(AliHMPIDDigit)DigLst(iCh)->At(iDig);
186	nDigCh[pDig->Ch()]++;
187	}
188	Printf("Occupancy for chamber %i = %4.2f %% and charged prim tracks %i and sec. tracks %i with total %i",
189	iCh,Float_t(nDigCh[iCh])*100/AliHMPIDDigit::kPadAll,nPrim[iCh],nSec[iCh],iChHits[iCh]);
190	}
191
192
193	}//events loop
194	GetLoader()->UnloadHits();
195	GetLoader()->UnloadDigits();
196	GetLoader()->GetRunLoader()->UnloadHeader();
197	GetLoader()->GetRunLoader()->UnloadKinematics();
198	}
199	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
200	void AliHMPID::CluPrint(Int_t iEvtN)const
201	{
202	//prints a list of HMPID clusters for a given event
203	Printf("List of HMPID clusters for event %i",iEvtN);
204	GetLoader()->GetRunLoader()->GetEvent(iEvtN);
205	if(GetLoader()->LoadRecPoints()) return;
206
207	Int_t iCluCnt=0;
208	GetLoader()->TreeR()->GetEntry(0);
209	for(Int_t iCh=0;iCh<7;iCh++){
210	TClonesArray pCluLst=(TClonesArray)fClu->At(iCh); iCluCnt+=pCluLst->GetEntries(); pCluLst->Print();
211	}
212	GetLoader()->UnloadRecPoints();
213	Printf("totally %i clusters for event %i",iCluCnt,iEvtN);
214	}
215	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
216	void AliHMPID::SummaryOfEvent(Int_t iEvtN) const
217	{
218	//prints a summary for a given event
219	AliInfo(Form("Summary of event %i",iEvtN));
220	GetLoader()->GetRunLoader()->GetEvent(iEvtN);
221	if(GetLoader()->GetRunLoader()->LoadHeader()) return;
222	if(GetLoader()->GetRunLoader()->LoadKinematics()) return;
223	AliStack *pStack=GetLoader()->GetRunLoader()->Stack();
224
225	AliGenEventHeader* pGenHeader = gAlice->GetHeader()->GenEventHeader();
226	if(pGenHeader->InheritsFrom("AliGenHijingEventHeader")) {
227	AliInfo(Form(" Hijing event with impact parameter b = %.2f (fm)",((AliGenHijingEventHeader*) pGenHeader)->ImpactParameter()));
228	}
229	Int_t nChargedPrimaries=0;
230	for(Int_t i=0;i<pStack->GetNtrack();i++) {
231	TParticle *pParticle = pStack->Particle(i);
232	if(pParticle->IsPrimary()&&pParticle->GetPDG()->Charge()!=0) nChargedPrimaries++;
233	}
234	AliInfo(Form("Total number of primaries %i",pStack->GetNprimary()));
235	AliInfo(Form("Total number of charged primaries %i",nChargedPrimaries));
236	AliInfo(Form("Total n. of tracks in stack(+sec) %i",pStack->GetNtrack()));
237	GetLoader()->GetRunLoader()->UnloadHeader();
238	GetLoader()->GetRunLoader()->UnloadKinematics();
239	}
240	//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++