Dear Colleagues,

 

One of the few outstanding issues holding up the generation of a high-statistics MC sample to accompany the 2016 data analysis is the determination of the ECal position relative to the SVT. I have repeated the analysis pioneered by Nathan to plot the track-cluster match nSigma vs y for the 2016 data. Specifically, I used the FEE runs 7479 and 7483. I get essentially the same values (~42.5) as were found for the 2015 data. I have posted plots and some description to the slack 2pt3gev channel. I would appreciate it if someone were to review the methodology to ensure that this technique is valid, i.e. that there are no hidden corrections in the calculation of nSigma which would skew the measurement.

 

Sincerely,

Norman

 

If anyone wants to replicate this, I used the org.hps.analysis.alignment.SvtCalorimeterAlignmentDriver class with the following modification to pick up the new FEE collection,

 

    protected void process(EventHeader event) {

        List<ReconstructedParticle> rpList = event.get(ReconstructedParticle.class, "FinalStateParticles");

        if (event.hasCollection(ReconstructedParticle.class, "OtherElectrons")) {

            rpList.addAll(event.get(ReconstructedParticle.class, "OtherElectrons"));

        }

        for (ReconstructedParticle rp : rpList) {

 

            if (!TrackType.isGBL(rp.getType())) {

                continue;

            }

 

            // require both track and cluster

            if (rp.getClusters().size() != 1) {

                continue;

            }

 

            if (rp.getTracks().size() != 1) {

                continue;

            }

 

            double nSigma = rp.getGoodnessOfPID();

            Track t = rp.getTracks().get(0);

            TrackState trackAtEcal = TrackStateUtils.getTrackStateAtECal(t);

            double[] tposAtEcal = trackAtEcal.getReferencePoint();

 

            // look for calorimeter edge wrt SVT

            if (tposAtEcal[2] > 0) {

                trkAtEcalXvsNSigmaTop.fill(nSigma, tposAtEcal[2]);

            } else {

                trkAtEcalXvsNSigmaBottom.fill(nSigma, -tposAtEcal[2]);

            }

        }

    }