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);
+ 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;
- Printf("Track %02i Ch. %2i with pOuter %7.2f with ThetaCer %7.3f with %3i photons with MIP Q %4i size %2i (idx %3i)",iTrk,ch,pMomOut,
- pTrk->GetHMPIDsignal(),nacc,q,size,idx);
+ 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()
//++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++