Print

---

Good Morning Holly,
Assuming you have access to the org.lcsim.geometry.Detector the following will give you the position of the front of the crystal.

               HPSEcal3 ecal = (HPSEcal3) detector.getSubdetector("Ecal");
                HPSEcalDetectorElement detElement = (HPSEcalDetectorElement) ecal.getDetectorElement();
                EcalCrystal crystal = detElement.getCrystal(ix, iy);
                Hep3Vector pos = crystal.getPositionFront();

If you need the center of the crystal I would suggest that you add that method to the EcalCrystal class.
Let me know if you encounter any difficulties or have any further questions,
Norman

From: [log in to unmask] [mailto:[log in to unmask]] On Behalf Of Holly Vance
Sent: Tuesday, December 01, 2015 8:33 AM
To: hps-software
Subject: Another geometry question

Hi,

In addition to my previous e-mail, I would like to access the x,y vales (mm) of certain crystals (I would like to avoid hard-coding coordinates). I 've done this previously in clustering by using a hit (CalorimeterHit) and then getting the geometry:
IGeometryInfo geom = hit.getDetectorElement().getGeometry();
        double[] pos = geom.transformLocalToGlobal(VecOp.add(geom.transformGlobalToLocal(geom.getPosition()), (Hep3Vector) new BasicHep3Vector(0, 0, -1 * ((Trd) geom.getLogicalVolume().getSolid()).getZHalfLength()))).v();
Is there a way to do this without using a CalorimeterHit? Perhaps just using a crystal index like ix,iy?

Thanks,
Holly

________________________________

Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1

########################################################################
Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1