Printf("List of HMPID ESD summary for event %i",fEvt);
for(Int_t iTrk=0;iTrk<fEsd->GetNumberOfTracks();iTrk++){
AliESDtrack *pTrk = fEsd->GetTrack(iTrk);
+
+ Double_t xout[3],pout[3];
+
+ pTrk->GetOuterPxPyPz(pout);
+ Double_t pMomOut = TMath::Sqrt(pout[0]*pout[0]+pout[1]*pout[1]+pout[2]*pout[2]);
+
Float_t x,y;Int_t q,nacc; pTrk->GetHMPIDmip(x,y,q,nacc);
- Printf("Track %02i with phi %7.2f theta %7.2f p %7.2f with ThetaCer %5.3f with %i photons",iTrk,pTrk->Phi()*TMath::RadToDeg(),
- pTrk->Theta()*TMath::RadToDeg(),
- pTrk->GetP(),pTrk->GetHMPIDsignal(),nacc);
+ Float_t xra,yra,th,ph; pTrk->GetHMPIDtrk(xra,yra,th,ph);
+// Printf("xra %f yra %f th %f phi %f",xra,yra,th,ph);
+ Int_t ch,idx,size;
+ Int_t word = pTrk->GetHMPIDcluIdx();
+ ch = word/1000000;
+ word = word%1000000;
+ size = word/1000;
+ idx = word%1000;
+ Double_t rout[3]; pTrk->GetOuterXYZ(rout);
+ vol = gGeoManager->FindNode(rout[0],rout[1],rout[2]);
+ Printf("Track %02i Ch. %2i with pOuter %7.2f with ThetaCer %7.3f with %3i photons with MIP Q %4i size %2i (idx %3i) in vol. %s",iTrk,ch,pMomOut,
+ pTrk->GetHMPIDsignal(),nacc,q,size,idx,vol->GetName());
}
}//PrintEsd()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
void RenderDig(TObjArray *pDigLst)
{//used by ReadEvent() to render digs to Tbox structures, one per chamber
- Printf("RENDERDIG : event n. %i",fEvt);
-
for(Int_t iCh=0;iCh<=6;iCh++){ //chambers loop
TClonesArray *pDigCham=(TClonesArray*)pDigLst->At(iCh); //get digs list for this chamber
nDigs[iCh] = 0;
void DoZoom(Int_t evt, Int_t px, Int_t py, TObject *)
{
// Printf(" px %i py%i event %i",px,py,evt);
+
+ if(evt== 1) evt = 5; //for laptop with touchpad: zoom in with left button
+ if(evt== 2) evt = 6; //for laptop with touchpad: zoom out with both buttons
+
if(evt==kMouseMotion||evt==kButton1Down) DisplayInfo(evt,px,py);
- if(evt==11) CloseInfo();
+// if(evt==11) CloseInfo();
if(evt!=5 && evt!=6) return; //5- zoom in 6-zoom out
const Int_t minZoom=64;
const Int_t maxZoom=2;
{
ClearRenders();
- Printf("get Event...... event n. %i",fEvt);
ReadEvent();
while(!Trigger()) {
if(fNevt && fEvt>=fNevt) fEvt=0; //loop over max event
if(fNevt && fEvt<0) fEvt=fNevt-1; //loop over max event
- Printf("getting event %i out of %i",fEvt,fNevt);
+// Printf("getting event %i out of %i",fEvt,fNevt);
if(gDL) {