diff --git a/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/HitReader.java b/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/HitReader.java index 6fca3cdab2..b41ec179d4 100644 --- a/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/HitReader.java +++ b/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/HitReader.java @@ -17,6 +17,9 @@ import org.jlab.detector.banks.RawDataBank; import org.jlab.detector.calib.utils.ConstantsManager; import org.jlab.detector.geant4.v2.DCGeant4Factory; +import org.jlab.io.banks.DC_tot; +import org.jlab.io.banks.HitBasedTrkg_Hits; +import org.jlab.io.banks.ai_tracks; import org.jlab.rec.dc.Constants; import org.jlab.utils.groups.IndexedTable; @@ -148,11 +151,11 @@ private int getDCRBJitter(int sector, int layer, int wire, int order) { if(dcrbjitters!=null && reverseTT!=null) { long hash = IndexedTable.DEFAULT_GENERATOR.hashCode(sector, layer, wire, order); if(!reverseTT.hasEntryByHash(hash)) - LOGGER.log(Level.SEVERE, "Missing TT entry for slco = " + sector + " " + layer + " " + wire + " " + order); + LOGGER.log(Level.SEVERE, () -> "Missing TT entry for slco = " + sector + " " + layer + " " + wire + " " + order); int crate = reverseTT.getIntValueByHash("crate", hash); int slot = reverseTT.getIntValueByHash("slot", hash); if(!dcrbjitters.hasEntry(crate, slot)) - LOGGER.log(Level.SEVERE, "Missing DC::jitter entry for crate/slot = " + crate + "/" + slot); + LOGGER.log(Level.SEVERE, () -> "Missing DC::jitter entry for crate/slot = " + crate + "/" + slot); jitter = dcrbjitters.getIntValue("jitter", crate, slot); } return jitter; @@ -185,13 +188,13 @@ public void fetch_DCHits(DataEvent event) { this.set_NumTDCBankRows(bankFiltered.rows()); for (int i = 0; i < bankFiltered.rows(); i++) { - int sector = bankFiltered.getByte("sector", i); - int layer = (bankFiltered.getByte("layer", i)-1)%6 + 1; - int superlayer = (bankFiltered.getByte("layer", i)-1)/6 + 1; - int wire = bankFiltered.getShort("component", i); + int sector = bankFiltered.getByte(DC_tot.sector, i); + int layer = (bankFiltered.getByte(DC_tot.layer, i)-1)%6 + 1; + int superlayer = (bankFiltered.getByte(DC_tot.layer, i)-1)/6 + 1; + int wire = bankFiltered.getShort(DC_tot.component, i); int order = bankFiltered.trueOrder(i); - int jitter = this.getJitter(sector, bankFiltered.getByte("layer", i), wire, order); - int tdc = bankFiltered.getInt("TDC", i) - jitter; + int jitter = this.getJitter(sector, bankFiltered.getByte(DC_tot.layer, i), wire, order); + int tdc = bankFiltered.getInt(DC_tot.TDC, i) - jitter; int index = bankFiltered.trueIndex(i); boolean passHit = true; @@ -245,15 +248,14 @@ private Map> read_Hits() { List hits = new ArrayList<>(); for (int i = 0; i < rows; i++) { - int id = bank.getShort("id", i); - int sector = bank.getByte("sector", i); - int slayer = bank.getByte("superlayer", i); - int layer = bank.getByte("layer", i); - int wire = bank.getShort("wire", i); - int tdc = bank.getInt("TDC", i); - int jitter = bank.getByte("jitter", i); - int LR = bank.getByte("LR", i); - int clusterID = bank.getShort("clusterID", i); + int id = bank.getShort(HitBasedTrkg_Hits.id, i); + int sector = bank.getByte(HitBasedTrkg_Hits.sector, i); + int slayer = bank.getByte(HitBasedTrkg_Hits.superlayer, i); + int layer = bank.getByte(HitBasedTrkg_Hits.layer, i); + int wire = bank.getShort(HitBasedTrkg_Hits.wire, i); + int tdc = bank.getInt(HitBasedTrkg_Hits.TDC, i); + int jitter = bank.getByte(HitBasedTrkg_Hits.jitter, i); + int clusterID = bank.getShort(HitBasedTrkg_Hits.clusterID, i); //use only hits that have been fit to a track if (clusterID == -1) { @@ -357,19 +359,19 @@ private void read_HBHits(TimeToDistanceEstimator tde) { double[] trkDoca = new double[rows]; for (int i = 0; i < rows; i++) { - id[i] = bank.getShort("id", i); - status[i] = bank.getShort("status", i); - sector[i] = bank.getByte("sector", i); - slayer[i] = bank.getByte("superlayer", i); - layer[i] = bank.getByte("layer", i); - wire[i] = bank.getShort("wire", i); - tdc[i] = bank.getInt("TDC", i); - jitter[i] = bank.getByte("jitter", i); - id[i] = bank.getShort("id", i); - LR[i] = bank.getByte("LR", i); + id[i] = bank.getShort(HitBasedTrkg_Hits.id, i); + status[i] = bank.getShort(HitBasedTrkg_Hits.status, i); + sector[i] = bank.getByte(HitBasedTrkg_Hits.sector, i); + slayer[i] = bank.getByte(HitBasedTrkg_Hits.superlayer, i); + layer[i] = bank.getByte(HitBasedTrkg_Hits.layer, i); + wire[i] = bank.getShort(HitBasedTrkg_Hits.wire, i); + tdc[i] = bank.getInt(HitBasedTrkg_Hits.TDC, i); + jitter[i] = bank.getByte(HitBasedTrkg_Hits.jitter, i); + id[i] = bank.getShort(HitBasedTrkg_Hits.id, i); + LR[i] = bank.getByte(HitBasedTrkg_Hits.LR, i); - trkDoca[i] = bank.getFloat("trkDoca", i); - clusterID[i] = bank.getShort("clusterID", i); + trkDoca[i] = bank.getFloat(HitBasedTrkg_Hits.trkDoca, i); + clusterID[i] = bank.getShort(HitBasedTrkg_Hits.clusterID, i); trkID[i] = -1; if(this.id2tid.containsKey(id[i]) ){ trkID[i] = this.id2tid.get(id[i]); @@ -479,17 +481,17 @@ private void read_NNHits() { for (int j = 0; j < bankAI.rows(); j++) { Ids = new int[6]; tPars = new double[4]; - Ids[0] = (int)bankAI.getShort("c1", j); // clusId in superlayer 1 - Ids[1] = (int)bankAI.getShort("c2", j); - Ids[2] = (int)bankAI.getShort("c3", j); - Ids[3] = (int)bankAI.getShort("c4", j); - Ids[4] = (int)bankAI.getShort("c5", j); - Ids[5] = (int)bankAI.getShort("c6", j); // clusId in superlayer 6 + Ids[0] = (int)bankAI.getShort(ai_tracks.c1, j); // clusId in superlayer 1 + Ids[1] = (int)bankAI.getShort(ai_tracks.c2, j); + Ids[2] = (int)bankAI.getShort(ai_tracks.c3, j); + Ids[3] = (int)bankAI.getShort(ai_tracks.c4, j); + Ids[4] = (int)bankAI.getShort(ai_tracks.c5, j); + Ids[5] = (int)bankAI.getShort(ai_tracks.c6, j); // clusId in superlayer 6 - tPars[0] = (double)bankAI.getFloat("p", j); - tPars[1] = (double)bankAI.getFloat("theta", j); - tPars[2] = (double)bankAI.getFloat("phi", j); - tPars[3] = (double)bankAI.getByte("id", j); + tPars[0] = (double)bankAI.getFloat(ai_tracks.p, j); + tPars[1] = (double)bankAI.getFloat(ai_tracks.theta, j); + tPars[2] = (double)bankAI.getFloat(ai_tracks.phi, j); + tPars[3] = (double)bankAI.getByte(ai_tracks.id, j); aimatch.clear(); for (int k = 0; k < 6; k++) { @@ -501,9 +503,14 @@ private void read_NNHits() { if(clusterID>0) { if(this.aimatch.containsKey(clusterID)) { - Hit hit = new Hit(bank.getByte("sector", i), bank.getByte("superlayer", i), - bank.getByte("layer", i), bank.getShort("wire", i), bank.getInt("TDC", i), bank.getByte("jitter", i), bank.getShort("id", i)); - hit.set_Id(bank.getShort("id", i)); + Hit hit = new Hit(bank.getByte(HitBasedTrkg_Hits.sector, i), + bank.getByte(HitBasedTrkg_Hits.superlayer, i), + bank.getByte(HitBasedTrkg_Hits.layer, i), + bank.getShort(HitBasedTrkg_Hits.wire, i), + bank.getInt(HitBasedTrkg_Hits.TDC, i), + bank.getByte(HitBasedTrkg_Hits.jitter, i), + bank.getShort(HitBasedTrkg_Hits.id, i)); + hit.set_Id(bank.getShort(HitBasedTrkg_Hits.id, i)); hit.calc_CellSize(detector); double posError = hit.get_CellSize() / Math.sqrt(12.); hit.set_DocaErr(posError); @@ -512,7 +519,7 @@ private void read_NNHits() { hit.NNTrkP = this.aimatch.get(clusterID)[0]; hit.NNTrkTheta = this.aimatch.get(clusterID)[1]; hit.NNTrkPhi = this.aimatch.get(clusterID)[2]; - LOGGER.log(Level.FINEST, "NN"+hit.printInfo()); + LOGGER.log(Level.FINEST, () -> "NN"+hit.printInfo()); this._DCHits.add(hit); } } diff --git a/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/RecoBankWriter.java b/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/RecoBankWriter.java index d018dd3c8d..e29a69b304 100644 --- a/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/RecoBankWriter.java +++ b/reconstruction/dc/src/main/java/org/jlab/rec/dc/banks/RecoBankWriter.java @@ -3,7 +3,19 @@ import java.util.ArrayList; import java.util.List; import org.jlab.detector.base.DetectorType; - +import org.jlab.io.banks.HitBasedTrkg_HBClusters; +import org.jlab.io.banks.HitBasedTrkg_HBCrosses; +import org.jlab.io.banks.HitBasedTrkg_HBHitTrkId; +import org.jlab.io.banks.HitBasedTrkg_HBHits; +import org.jlab.io.banks.HitBasedTrkg_HBSegments; +import org.jlab.io.banks.HitBasedTrkg_HBTracks; +import org.jlab.io.banks.HitBasedTrkg_Hits; +import org.jlab.io.banks.TimeBasedTrkg_TBClusters; +import org.jlab.io.banks.TimeBasedTrkg_TBCovMat; +import org.jlab.io.banks.TimeBasedTrkg_TBCrosses; +import org.jlab.io.banks.TimeBasedTrkg_TBSegments; +import org.jlab.io.banks.TimeBasedTrkg_TBTracks; +import org.jlab.io.banks.TimeBasedTrkg_Trajectory; import org.jlab.io.base.DataBank; import org.jlab.io.base.DataEvent; import org.jlab.rec.dc.cluster.FittedCluster; @@ -74,23 +86,23 @@ public DataBank fillHitsBank(DataEvent event, List hitlist) { rejCnt++; continue; } - bank.setShort("id", i-rejCnt, (short) hitlist.get(i).get_Id()); - bank.setShort("indexTDC", i-rejCnt, (short) hitlist.get(i).get_IndexTDC()); - bank.setShort("status", i-rejCnt, (short) hitlist.get(i).get_QualityFac()); - bank.setByte("superlayer", i-rejCnt, (byte) hitlist.get(i).get_Superlayer()); - bank.setByte("layer", i-rejCnt, (byte) hitlist.get(i).get_Layer()); - bank.setByte("sector", i-rejCnt, (byte) hitlist.get(i).get_Sector()); - bank.setShort("wire", i-rejCnt, (short) hitlist.get(i).get_Wire()); - bank.setFloat("docaError", i-rejCnt, (float) hitlist.get(i).get_DocaErr()); - bank.setFloat("trkDoca", i-rejCnt, (float) hitlist.get(i).get_ClusFitDoca()); - bank.setFloat("LocX", i-rejCnt, (float) hitlist.get(i).get_lX()); - bank.setFloat("LocY", i-rejCnt, (float) hitlist.get(i).get_lY()); - bank.setFloat("X", i-rejCnt, (float) hitlist.get(i).get_X()); - bank.setFloat("Z", i-rejCnt, (float) hitlist.get(i).get_Z()); - bank.setByte("LR", i-rejCnt, (byte) hitlist.get(i).get_LeftRightAmb()); - bank.setShort("clusterID", i-rejCnt, (short) hitlist.get(i).get_AssociatedClusterID()); - bank.setInt("TDC",i-rejCnt,hitlist.get(i).get_TDC()); - bank.setByte("jitter",i, (byte) hitlist.get(i).getJitter()); + bank.setShort(HitBasedTrkg_Hits.id, i-rejCnt, (short) hitlist.get(i).get_Id()); + bank.setShort(HitBasedTrkg_Hits.indexTDC, i-rejCnt, (short) hitlist.get(i).get_IndexTDC()); + bank.setShort(HitBasedTrkg_Hits.status, i-rejCnt, (short) hitlist.get(i).get_QualityFac()); + bank.setByte(HitBasedTrkg_Hits.superlayer, i-rejCnt, (byte) hitlist.get(i).get_Superlayer()); + bank.setByte(HitBasedTrkg_Hits.layer, i-rejCnt, (byte) hitlist.get(i).get_Layer()); + bank.setByte(HitBasedTrkg_Hits.sector, i-rejCnt, (byte) hitlist.get(i).get_Sector()); + bank.setShort(HitBasedTrkg_Hits.wire, i-rejCnt, (short) hitlist.get(i).get_Wire()); + bank.setFloat(HitBasedTrkg_Hits.docaError, i-rejCnt, (float) hitlist.get(i).get_DocaErr()); + bank.setFloat(HitBasedTrkg_Hits.trkDoca, i-rejCnt, (float) hitlist.get(i).get_ClusFitDoca()); + bank.setFloat(HitBasedTrkg_Hits.LocX, i-rejCnt, (float) hitlist.get(i).get_lX()); + bank.setFloat(HitBasedTrkg_Hits.LocY, i-rejCnt, (float) hitlist.get(i).get_lY()); + bank.setFloat(HitBasedTrkg_Hits.X, i-rejCnt, (float) hitlist.get(i).get_X()); + bank.setFloat(HitBasedTrkg_Hits.Z, i-rejCnt, (float) hitlist.get(i).get_Z()); + bank.setByte(HitBasedTrkg_Hits.LR, i-rejCnt, (byte) hitlist.get(i).get_LeftRightAmb()); + bank.setShort(HitBasedTrkg_Hits.clusterID, i-rejCnt, (short) hitlist.get(i).get_AssociatedClusterID()); + bank.setInt(HitBasedTrkg_Hits.TDC,i-rejCnt,hitlist.get(i).get_TDC()); + bank.setByte(HitBasedTrkg_Hits.jitter,i, (byte) hitlist.get(i).getJitter()); } @@ -117,22 +129,22 @@ public DataBank fillHBHitsBank(DataEvent event, List hitlist) { continue; } - bank.setShort("id", i-rejCnt, (short) hitlist.get(i).get_Id()); - bank.setShort("status", i-rejCnt, (short) hitlist.get(i).get_QualityFac()); - bank.setByte("superlayer", i-rejCnt, (byte) hitlist.get(i).get_Superlayer()); - bank.setByte("layer", i-rejCnt, (byte) hitlist.get(i).get_Layer()); - bank.setByte("sector", i-rejCnt, (byte) hitlist.get(i).get_Sector()); - bank.setShort("wire", i-rejCnt, (short) hitlist.get(i).get_Wire()); - bank.setFloat("docaError", i-rejCnt, (float) hitlist.get(i).get_DocaErr()); - bank.setFloat("trkDoca", i-rejCnt, (float) hitlist.get(i).get_ClusFitDoca()); - bank.setFloat("LocX", i-rejCnt, (float) hitlist.get(i).get_lX()); - bank.setFloat("LocY", i-rejCnt, (float) hitlist.get(i).get_lY()); - bank.setFloat("X", i-rejCnt, (float) hitlist.get(i).get_X()); - bank.setFloat("Z", i-rejCnt, (float) hitlist.get(i).get_Z()); - bank.setByte("LR", i-rejCnt, (byte) hitlist.get(i).get_LeftRightAmb()); - bank.setShort("clusterID", i-rejCnt, (short) hitlist.get(i).get_AssociatedClusterID()); - bank.setInt("TDC",i-rejCnt,hitlist.get(i).get_TDC()); - bank.setByte("jitter",i, (byte) hitlist.get(i).getJitter()); + bank.setShort(HitBasedTrkg_HBHits.id, i-rejCnt, (short) hitlist.get(i).get_Id()); + bank.setShort(HitBasedTrkg_HBHits.status, i-rejCnt, (short) hitlist.get(i).get_QualityFac()); + bank.setByte(HitBasedTrkg_HBHits.superlayer, i-rejCnt, (byte) hitlist.get(i).get_Superlayer()); + bank.setByte(HitBasedTrkg_HBHits.layer, i-rejCnt, (byte) hitlist.get(i).get_Layer()); + bank.setByte(HitBasedTrkg_HBHits.sector, i-rejCnt, (byte) hitlist.get(i).get_Sector()); + bank.setShort(HitBasedTrkg_HBHits.wire, i-rejCnt, (short) hitlist.get(i).get_Wire()); + bank.setFloat(HitBasedTrkg_HBHits.docaError, i-rejCnt, (float) hitlist.get(i).get_DocaErr()); + bank.setFloat(HitBasedTrkg_HBHits.trkDoca, i-rejCnt, (float) hitlist.get(i).get_ClusFitDoca()); + bank.setFloat(HitBasedTrkg_HBHits.LocX, i-rejCnt, (float) hitlist.get(i).get_lX()); + bank.setFloat(HitBasedTrkg_HBHits.LocY, i-rejCnt, (float) hitlist.get(i).get_lY()); + bank.setFloat(HitBasedTrkg_HBHits.X, i-rejCnt, (float) hitlist.get(i).get_X()); + bank.setFloat(HitBasedTrkg_HBHits.Z, i-rejCnt, (float) hitlist.get(i).get_Z()); + bank.setByte(HitBasedTrkg_HBHits.LR, i-rejCnt, (byte) hitlist.get(i).get_LeftRightAmb()); + bank.setShort(HitBasedTrkg_HBHits.clusterID, i-rejCnt, (short) hitlist.get(i).get_AssociatedClusterID()); + bank.setInt(HitBasedTrkg_HBHits.TDC,i-rejCnt,hitlist.get(i).get_TDC()); + bank.setByte(HitBasedTrkg_HBHits.jitter,i, (byte) hitlist.get(i).getJitter()); } @@ -162,11 +174,11 @@ public DataBank fillHBHitsTrkIdBank(DataEvent event, List hitlist) { rejCnt++; continue; } - bank.setShort("id", i-rejCnt, (short) hitlist.get(i).get_Id()); - bank.setShort("tid", i-rejCnt, (short) hitlist.get(i).get_AssociatedHBTrackID()); - bank.setFloat("B", i-rejCnt, (float) hitlist.get(i).getB()); - bank.setFloat("TProp", i-rejCnt, (float) hitlist.get(i).getTProp()); - bank.setFloat("TFlight", i-rejCnt, (float) hitlist.get(i).getTFlight()); + bank.setShort(HitBasedTrkg_HBHitTrkId.id, i-rejCnt, (short) hitlist.get(i).get_Id()); + bank.setShort(HitBasedTrkg_HBHitTrkId.tid, i-rejCnt, (short) hitlist.get(i).get_AssociatedHBTrackID()); + bank.setFloat(HitBasedTrkg_HBHitTrkId.B, i-rejCnt, (float) hitlist.get(i).getB()); + bank.setFloat(HitBasedTrkg_HBHitTrkId.TProp, i-rejCnt, (float) hitlist.get(i).getTProp()); + bank.setFloat(HitBasedTrkg_HBHitTrkId.TFlight, i-rejCnt, (float) hitlist.get(i).getTFlight()); // if(i>0 && hitlist.get(i).getSector()==hitlist.get(i-1).getSector() && // hitlist.get(i).get_Superlayer()==hitlist.get(i-1).get_Superlayer() && // hitlist.get(i).get_Layer()==hitlist.get(i-1).get_Layer() && @@ -215,24 +227,24 @@ public DataBank fillHBClustersBank(DataEvent event, List cluslist } double chi2 = 0; - bank.setShort("id", i, (short) cluslist.get(i).get_Id()); + bank.setShort(HitBasedTrkg_HBClusters.id, i, (short) cluslist.get(i).get_Id()); int status = 0; if(cluslist.get(i).size()<6) status = 1; - bank.setShort("status", i, (short) status); - bank.setByte("superlayer", i, (byte) cluslist.get(i).get_Superlayer()); - bank.setByte("sector", i, (byte) cluslist.get(i).get_Sector()); + bank.setShort(HitBasedTrkg_HBClusters.status, i, (short) status); + bank.setByte(HitBasedTrkg_HBClusters.superlayer, i, (byte) cluslist.get(i).get_Superlayer()); + bank.setByte(HitBasedTrkg_HBClusters.sector, i, (byte) cluslist.get(i).get_Sector()); - bank.setFloat("avgWire", i, (float) cluslist.get(i).getAvgwire()); - bank.setByte("size", i, (byte) cluslist.get(i).size()); + bank.setFloat(HitBasedTrkg_HBClusters.avgWire, i, (float) cluslist.get(i).getAvgwire()); + bank.setByte(HitBasedTrkg_HBClusters.size, i, (byte) cluslist.get(i).size()); double fitSlope = cluslist.get(i).get_clusterLineFitSlope(); double fitInterc = cluslist.get(i).get_clusterLineFitIntercept(); - bank.setFloat("fitSlope", i, (float) fitSlope); - bank.setFloat("fitSlopeErr", i, (float) cluslist.get(i).get_clusterLineFitSlopeErr()); - bank.setFloat("fitInterc", i, (float) fitInterc); - bank.setFloat("fitIntercErr", i, (float) cluslist.get(i).get_clusterLineFitInterceptErr()); + bank.setFloat(HitBasedTrkg_HBClusters.fitSlope, i, (float) fitSlope); + bank.setFloat(HitBasedTrkg_HBClusters.fitSlopeErr, i, (float) cluslist.get(i).get_clusterLineFitSlopeErr()); + bank.setFloat(HitBasedTrkg_HBClusters.fitInterc, i, (float) fitInterc); + bank.setFloat(HitBasedTrkg_HBClusters.fitIntercErr, i, (float) cluslist.get(i).get_clusterLineFitInterceptErr()); for (int j = 0; j < cluslist.get(i).size(); j++) { if (j < hitIdxArray.length) { @@ -242,7 +254,7 @@ public DataBank fillHBClustersBank(DataEvent event, List cluslist double residual = cluslist.get(i).get(j).get_ClusFitDoca() / (cluslist.get(i).get(j).get_CellSize() / Math.sqrt(12.)); chi2 += residual * residual; } - bank.setFloat("fitChisqProb", i, (float) ProbChi2perNDF.prob(chi2, cluslist.get(i).size() - 2)); + bank.setFloat(HitBasedTrkg_HBClusters.fitChisqProb, i, (float) ProbChi2perNDF.prob(chi2, cluslist.get(i).size() - 2)); for (int j = 0; j < hitIdxArray.length; j++) { String hitStrg = "Hit"; @@ -279,25 +291,25 @@ public DataBank fillHBSegmentsBank(DataEvent event, List seglist) { double chi2 = 0; - bank.setShort("id", i, (short) seglist.get(i).get_Id()); - bank.setByte("superlayer", i, (byte) seglist.get(i).get_Superlayer()); - bank.setByte("sector", i, (byte) seglist.get(i).get_Sector()); + bank.setShort(HitBasedTrkg_HBSegments.id, i, (short) seglist.get(i).get_Id()); + bank.setByte(HitBasedTrkg_HBSegments.superlayer, i, (byte) seglist.get(i).get_Superlayer()); + bank.setByte(HitBasedTrkg_HBSegments.sector, i, (byte) seglist.get(i).get_Sector()); FittedCluster cls = seglist.get(i).get_fittedCluster(); - bank.setShort("Cluster_ID", i, (short) cls.get_Id()); + bank.setShort(HitBasedTrkg_HBSegments.Cluster_ID, i, (short) cls.get_Id()); - bank.setFloat("avgWire", i, (float) cls.getAvgwire()); - bank.setByte("size", i, (byte) seglist.get(i).size()); + bank.setFloat(HitBasedTrkg_HBSegments.avgWire, i, (float) cls.getAvgwire()); + bank.setByte(HitBasedTrkg_HBSegments.size, i, (byte) seglist.get(i).size()); - bank.setFloat("fitSlope", i, (float) cls.get_clusterLineFitSlope()); - bank.setFloat("fitSlopeErr", i, (float) cls.get_clusterLineFitSlopeErr()); - bank.setFloat("fitInterc", i, (float) cls.get_clusterLineFitIntercept()); - bank.setFloat("fitIntercErr", i, (float) cls.get_clusterLineFitInterceptErr()); + bank.setFloat(HitBasedTrkg_HBSegments.fitSlope, i, (float) cls.get_clusterLineFitSlope()); + bank.setFloat(HitBasedTrkg_HBSegments.fitSlopeErr, i, (float) cls.get_clusterLineFitSlopeErr()); + bank.setFloat(HitBasedTrkg_HBSegments.fitInterc, i, (float) cls.get_clusterLineFitIntercept()); + bank.setFloat(HitBasedTrkg_HBSegments.fitIntercErr, i, (float) cls.get_clusterLineFitInterceptErr()); - bank.setFloat("SegEndPoint1X", i, (float) seglist.get(i).get_SegmentEndPoints()[0]); - bank.setFloat("SegEndPoint1Z", i, (float) seglist.get(i).get_SegmentEndPoints()[1]); - bank.setFloat("SegEndPoint2X", i, (float) seglist.get(i).get_SegmentEndPoints()[2]); - bank.setFloat("SegEndPoint2Z", i, (float) seglist.get(i).get_SegmentEndPoints()[3]); + bank.setFloat(HitBasedTrkg_HBSegments.SegEndPoint1X, i, (float) seglist.get(i).get_SegmentEndPoints()[0]); + bank.setFloat(HitBasedTrkg_HBSegments.SegEndPoint1Z, i, (float) seglist.get(i).get_SegmentEndPoints()[1]); + bank.setFloat(HitBasedTrkg_HBSegments.SegEndPoint2X, i, (float) seglist.get(i).get_SegmentEndPoints()[2]); + bank.setFloat(HitBasedTrkg_HBSegments.SegEndPoint2Z, i, (float) seglist.get(i).get_SegmentEndPoints()[3]); for (int j = 0; j < seglist.get(i).size(); j++) { if (seglist.get(i).get_Id() == -1) { @@ -310,7 +322,7 @@ public DataBank fillHBSegmentsBank(DataEvent event, List seglist) { double residual = seglist.get(i).get(j).get_ClusFitDoca() / (seglist.get(i).get(j).get_CellSize() / Math.sqrt(12.)); chi2 += residual * residual; } - bank.setFloat("fitChisqProb", i, (float) ProbChi2perNDF.prob(chi2, seglist.get(i).size() - 2)); + bank.setFloat(HitBasedTrkg_HBSegments.fitChisqProb, i, (float) ProbChi2perNDF.prob(chi2, seglist.get(i).size() - 2)); for (int j = 0; j < hitIdxArray.length; j++) { String hitStrg = "Hit"; @@ -343,24 +355,24 @@ public DataBank fillHBCrossesBank(DataEvent event, List crosslist) { int index=0; for (Cross aCrosslist : crosslist) { if (aCrosslist.get_Id() != -1) { - bank.setShort("id", index, (short) aCrosslist.get_Id()); - bank.setShort("status", index, (short) 0); - bank.setByte("sector", index, (byte) aCrosslist.get_Sector()); - bank.setByte("region", index, (byte) aCrosslist.get_Region()); - bank.setFloat("x", index, (float) aCrosslist.get_Point().x()); - bank.setFloat("y", index, (float) aCrosslist.get_Point().y()); - bank.setFloat("z", index, (float) aCrosslist.get_Point().z()); - bank.setFloat("err_x", index, (float) aCrosslist.get_PointErr().x()); - bank.setFloat("err_y", index, (float) aCrosslist.get_PointErr().y()); - bank.setFloat("err_z", index, (float) aCrosslist.get_PointErr().z()); - bank.setFloat("ux", index, (float) aCrosslist.get_Dir().x()); - bank.setFloat("uy", index, (float) aCrosslist.get_Dir().y()); - bank.setFloat("uz", index, (float) aCrosslist.get_Dir().z()); - bank.setFloat("err_ux", index, (float) aCrosslist.get_DirErr().x()); - bank.setFloat("err_uy", index, (float) aCrosslist.get_DirErr().y()); - bank.setFloat("err_uz", index, (float) aCrosslist.get_DirErr().z()); - bank.setShort("Segment1_ID", index, (short) aCrosslist.get_Segment1().get_Id()); - bank.setShort("Segment2_ID", index, (short) aCrosslist.get_Segment2().get_Id()); + bank.setShort(HitBasedTrkg_HBCrosses.id, index, (short) aCrosslist.get_Id()); + bank.setShort(HitBasedTrkg_HBCrosses.status, index, (short) 0); + bank.setByte(HitBasedTrkg_HBCrosses.sector, index, (byte) aCrosslist.get_Sector()); + bank.setByte(HitBasedTrkg_HBCrosses.region, index, (byte) aCrosslist.get_Region()); + bank.setFloat(HitBasedTrkg_HBCrosses.x, index, (float) aCrosslist.get_Point().x()); + bank.setFloat(HitBasedTrkg_HBCrosses.y, index, (float) aCrosslist.get_Point().y()); + bank.setFloat(HitBasedTrkg_HBCrosses.z, index, (float) aCrosslist.get_Point().z()); + bank.setFloat(HitBasedTrkg_HBCrosses.err_x, index, (float) aCrosslist.get_PointErr().x()); + bank.setFloat(HitBasedTrkg_HBCrosses.err_y, index, (float) aCrosslist.get_PointErr().y()); + bank.setFloat(HitBasedTrkg_HBCrosses.err_z, index, (float) aCrosslist.get_PointErr().z()); + bank.setFloat(HitBasedTrkg_HBCrosses.ux, index, (float) aCrosslist.get_Dir().x()); + bank.setFloat(HitBasedTrkg_HBCrosses.uy, index, (float) aCrosslist.get_Dir().y()); + bank.setFloat(HitBasedTrkg_HBCrosses.uz, index, (float) aCrosslist.get_Dir().z()); + bank.setFloat(HitBasedTrkg_HBCrosses.err_ux, index, (float) aCrosslist.get_DirErr().x()); + bank.setFloat(HitBasedTrkg_HBCrosses.err_uy, index, (float) aCrosslist.get_DirErr().y()); + bank.setFloat(HitBasedTrkg_HBCrosses.err_uz, index, (float) aCrosslist.get_DirErr().z()); + bank.setShort(HitBasedTrkg_HBCrosses.Segment1_ID, index, (short) aCrosslist.get_Segment1().get_Id()); + bank.setShort(HitBasedTrkg_HBCrosses.Segment2_ID, index, (short) aCrosslist.get_Segment2().get_Id()); index++; } } @@ -372,42 +384,41 @@ public DataBank fillHBTracksBank(DataEvent event, List candlist) { DataBank bank = event.createBank(name, candlist.size()); for (int i = 0; i < candlist.size(); i++) { - bank.setShort("id", i, (short) candlist.get(i).get_Id()); - bank.setByte("sector", i, (byte) candlist.get(i).getSector()); - bank.setByte("q", i, (byte) candlist.get(i).get_Q()); - //bank.setShort("status", i, (short) (100+candlist.get(i).get_Status()*10+candlist.get(i).get_MissingSuperlayer())); - bank.setShort("status", i, (short) candlist.get(i).getBitStatus()); + bank.setShort(HitBasedTrkg_HBTracks.id, i, (short) candlist.get(i).get_Id()); + bank.setByte(HitBasedTrkg_HBTracks.sector, i, (byte) candlist.get(i).getSector()); + bank.setByte(HitBasedTrkg_HBTracks.q, i, (byte) candlist.get(i).get_Q()); + bank.setShort(HitBasedTrkg_HBTracks.status, i, (short) candlist.get(i).getBitStatus()); if(candlist.get(i).get_PreRegion1CrossPoint()!=null) { - bank.setFloat("c1_x", i, (float) candlist.get(i).get_PreRegion1CrossPoint().x()); - bank.setFloat("c1_y", i, (float) candlist.get(i).get_PreRegion1CrossPoint().y()); - bank.setFloat("c1_z", i, (float) candlist.get(i).get_PreRegion1CrossPoint().z()); - bank.setFloat("c1_ux", i, (float) candlist.get(i).get_PreRegion1CrossDir().x()); - bank.setFloat("c1_uy", i, (float) candlist.get(i).get_PreRegion1CrossDir().y()); - bank.setFloat("c1_uz", i, (float) candlist.get(i).get_PreRegion1CrossDir().z()); + bank.setFloat(HitBasedTrkg_HBTracks.c1_x, i, (float) candlist.get(i).get_PreRegion1CrossPoint().x()); + bank.setFloat(HitBasedTrkg_HBTracks.c1_y, i, (float) candlist.get(i).get_PreRegion1CrossPoint().y()); + bank.setFloat(HitBasedTrkg_HBTracks.c1_z, i, (float) candlist.get(i).get_PreRegion1CrossPoint().z()); + bank.setFloat(HitBasedTrkg_HBTracks.c1_ux, i, (float) candlist.get(i).get_PreRegion1CrossDir().x()); + bank.setFloat(HitBasedTrkg_HBTracks.c1_uy, i, (float) candlist.get(i).get_PreRegion1CrossDir().y()); + bank.setFloat(HitBasedTrkg_HBTracks.c1_uz, i, (float) candlist.get(i).get_PreRegion1CrossDir().z()); } if(candlist.get(i).get_PostRegion3CrossPoint()!=null) { - bank.setFloat("c3_x", i, (float) candlist.get(i).get_PostRegion3CrossPoint().x()); - bank.setFloat("c3_y", i, (float) candlist.get(i).get_PostRegion3CrossPoint().y()); - bank.setFloat("c3_z", i, (float) candlist.get(i).get_PostRegion3CrossPoint().z()); - bank.setFloat("c3_ux", i, (float) candlist.get(i).get_PostRegion3CrossDir().x()); - bank.setFloat("c3_uy", i, (float) candlist.get(i).get_PostRegion3CrossDir().y()); - bank.setFloat("c3_uz", i, (float) candlist.get(i).get_PostRegion3CrossDir().z()); + bank.setFloat(HitBasedTrkg_HBTracks.c3_x, i, (float) candlist.get(i).get_PostRegion3CrossPoint().x()); + bank.setFloat(HitBasedTrkg_HBTracks.c3_y, i, (float) candlist.get(i).get_PostRegion3CrossPoint().y()); + bank.setFloat(HitBasedTrkg_HBTracks.c3_z, i, (float) candlist.get(i).get_PostRegion3CrossPoint().z()); + bank.setFloat(HitBasedTrkg_HBTracks.c3_ux, i, (float) candlist.get(i).get_PostRegion3CrossDir().x()); + bank.setFloat(HitBasedTrkg_HBTracks.c3_uy, i, (float) candlist.get(i).get_PostRegion3CrossDir().y()); + bank.setFloat(HitBasedTrkg_HBTracks.c3_uz, i, (float) candlist.get(i).get_PostRegion3CrossDir().z()); } if(candlist.get(i).get_Region1TrackX()!=null) { - bank.setFloat("t1_x", i, (float) candlist.get(i).get_Region1TrackX().x()); - bank.setFloat("t1_y", i, (float) candlist.get(i).get_Region1TrackX().y()); - bank.setFloat("t1_z", i, (float) candlist.get(i).get_Region1TrackX().z()); - bank.setFloat("t1_px", i, (float) candlist.get(i).get_Region1TrackP().x()); - bank.setFloat("t1_py", i, (float) candlist.get(i).get_Region1TrackP().y()); - bank.setFloat("t1_pz", i, (float) candlist.get(i).get_Region1TrackP().z()); - } - bank.setFloat("pathlength", i, (float) candlist.get(i).get_TotPathLen()); - bank.setFloat("Vtx0_x", i, (float) candlist.get(i).get_Vtx0().x()); - bank.setFloat("Vtx0_y", i, (float) candlist.get(i).get_Vtx0().y()); - bank.setFloat("Vtx0_z", i, (float) candlist.get(i).get_Vtx0().z()); - bank.setFloat("p0_x", i, (float) candlist.get(i).get_pAtOrig().x()); - bank.setFloat("p0_y", i, (float) candlist.get(i).get_pAtOrig().y()); - bank.setFloat("p0_z", i, (float) candlist.get(i).get_pAtOrig().z()); + bank.setFloat(HitBasedTrkg_HBTracks.t1_x, i, (float) candlist.get(i).get_Region1TrackX().x()); + bank.setFloat(HitBasedTrkg_HBTracks.t1_y, i, (float) candlist.get(i).get_Region1TrackX().y()); + bank.setFloat(HitBasedTrkg_HBTracks.t1_z, i, (float) candlist.get(i).get_Region1TrackX().z()); + bank.setFloat(HitBasedTrkg_HBTracks.t1_px, i, (float) candlist.get(i).get_Region1TrackP().x()); + bank.setFloat(HitBasedTrkg_HBTracks.t1_py, i, (float) candlist.get(i).get_Region1TrackP().y()); + bank.setFloat(HitBasedTrkg_HBTracks.t1_pz, i, (float) candlist.get(i).get_Region1TrackP().z()); + } + bank.setFloat(HitBasedTrkg_HBTracks.pathlength, i, (float) candlist.get(i).get_TotPathLen()); + bank.setFloat(HitBasedTrkg_HBTracks.Vtx0_x, i, (float) candlist.get(i).get_Vtx0().x()); + bank.setFloat(HitBasedTrkg_HBTracks.Vtx0_y, i, (float) candlist.get(i).get_Vtx0().y()); + bank.setFloat(HitBasedTrkg_HBTracks.Vtx0_z, i, (float) candlist.get(i).get_Vtx0().z()); + bank.setFloat(HitBasedTrkg_HBTracks.p0_x, i, (float) candlist.get(i).get_pAtOrig().x()); + bank.setFloat(HitBasedTrkg_HBTracks.p0_y, i, (float) candlist.get(i).get_pAtOrig().y()); + bank.setFloat(HitBasedTrkg_HBTracks.p0_z, i, (float) candlist.get(i).get_pAtOrig().z()); //fill associated IDs for(int r = 0; r < 3; r++) { bank.setShort("Cross"+String.valueOf(r+1)+"_ID", @@ -430,13 +441,13 @@ public DataBank fillHBTracksBank(DataEvent event, List candlist) { i, (short) candlist.get(i).getSingleSuperlayer().get_fittedCluster().get_Id()); } - bank.setFloat("chi2", i, (float) candlist.get(i).get_FitChi2()); - bank.setShort("ndf", i, (short) candlist.get(i).get_FitNDF()); - bank.setFloat("x", i, (float) candlist.get(i).getFinalStateVec().x()); - bank.setFloat("y", i, (float) candlist.get(i).getFinalStateVec().y()); - bank.setFloat("z", i, (float) candlist.get(i).getFinalStateVec().getZ()); - bank.setFloat("tx", i, (float) candlist.get(i).getFinalStateVec().tanThetaX()); - bank.setFloat("ty", i, (float) candlist.get(i).getFinalStateVec().tanThetaY()); + bank.setFloat(HitBasedTrkg_HBTracks.chi2, i, (float) candlist.get(i).get_FitChi2()); + bank.setShort(HitBasedTrkg_HBTracks.ndf, i, (short) candlist.get(i).get_FitNDF()); + bank.setFloat(HitBasedTrkg_HBTracks.x, i, (float) candlist.get(i).getFinalStateVec().x()); + bank.setFloat(HitBasedTrkg_HBTracks.y, i, (float) candlist.get(i).getFinalStateVec().y()); + bank.setFloat(HitBasedTrkg_HBTracks.z, i, (float) candlist.get(i).getFinalStateVec().getZ()); + bank.setFloat(HitBasedTrkg_HBTracks.tx, i, (float) candlist.get(i).getFinalStateVec().tanThetaX()); + bank.setFloat(HitBasedTrkg_HBTracks.ty, i, (float) candlist.get(i).getFinalStateVec().tanThetaY()); } //bank.show(); @@ -458,50 +469,50 @@ private DataBank fillTrackCovMatLabBank(DataEvent event, List candlist) { DataBank bank = event.createBank(bankNames.getCovmatBank(), candlist.size()); for (int i = 0; i < candlist.size(); i++) { - bank.setShort("id", i, (short) candlist.get(i).get_Id()); + bank.setShort(TimeBasedTrkg_TBCovMat.id, i, (short) candlist.get(i).get_Id()); if(candlist.get(i).get_CMInLab()!=null) { double[][] CM = candlist.get(i).get_CMInLab(); - bank.setFloat("C11", i, (float) CM[0][0]); - bank.setFloat("C12", i, (float) CM[0][1]); - bank.setFloat("C13", i, (float) CM[0][2]); - bank.setFloat("C14", i, (float) CM[0][3]); - bank.setFloat("C15", i, (float) CM[0][4]); - bank.setFloat("C16", i, (float) CM[0][5]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C11, i, (float) CM[0][0]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C12, i, (float) CM[0][1]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C13, i, (float) CM[0][2]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C14, i, (float) CM[0][3]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C15, i, (float) CM[0][4]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C16, i, (float) CM[0][5]); - bank.setFloat("C21", i, (float) CM[1][0]); - bank.setFloat("C22", i, (float) CM[1][1]); - bank.setFloat("C23", i, (float) CM[1][2]); - bank.setFloat("C24", i, (float) CM[1][3]); - bank.setFloat("C25", i, (float) CM[1][4]); - bank.setFloat("C26", i, (float) CM[1][5]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C21, i, (float) CM[1][0]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C22, i, (float) CM[1][1]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C23, i, (float) CM[1][2]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C24, i, (float) CM[1][3]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C25, i, (float) CM[1][4]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C26, i, (float) CM[1][5]); - bank.setFloat("C31", i, (float) CM[2][0]); - bank.setFloat("C32", i, (float) CM[2][1]); - bank.setFloat("C33", i, (float) CM[2][2]); - bank.setFloat("C34", i, (float) CM[2][3]); - bank.setFloat("C35", i, (float) CM[2][4]); - bank.setFloat("C36", i, (float) CM[2][5]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C31, i, (float) CM[2][0]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C32, i, (float) CM[2][1]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C33, i, (float) CM[2][2]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C34, i, (float) CM[2][3]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C35, i, (float) CM[2][4]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C36, i, (float) CM[2][5]); - bank.setFloat("C41", i, (float) CM[3][0]); - bank.setFloat("C42", i, (float) CM[3][1]); - bank.setFloat("C43", i, (float) CM[3][2]); - bank.setFloat("C44", i, (float) CM[3][3]); - bank.setFloat("C45", i, (float) CM[3][4]); - bank.setFloat("C46", i, (float) CM[3][5]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C41, i, (float) CM[3][0]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C42, i, (float) CM[3][1]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C43, i, (float) CM[3][2]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C44, i, (float) CM[3][3]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C45, i, (float) CM[3][4]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C46, i, (float) CM[3][5]); - bank.setFloat("C51", i, (float) CM[4][0]); - bank.setFloat("C52", i, (float) CM[4][1]); - bank.setFloat("C53", i, (float) CM[4][2]); - bank.setFloat("C54", i, (float) CM[4][3]); - bank.setFloat("C55", i, (float) CM[4][4]); - bank.setFloat("C56", i, (float) CM[4][5]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C51, i, (float) CM[4][0]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C52, i, (float) CM[4][1]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C53, i, (float) CM[4][2]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C54, i, (float) CM[4][3]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C55, i, (float) CM[4][4]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C56, i, (float) CM[4][5]); - bank.setFloat("C61", i, (float) CM[5][0]); - bank.setFloat("C62", i, (float) CM[5][1]); - bank.setFloat("C63", i, (float) CM[5][2]); - bank.setFloat("C64", i, (float) CM[5][3]); - bank.setFloat("C65", i, (float) CM[5][4]); - bank.setFloat("C66", i, (float) CM[5][5]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C61, i, (float) CM[5][0]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C62, i, (float) CM[5][1]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C63, i, (float) CM[5][2]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C64, i, (float) CM[5][3]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C65, i, (float) CM[5][4]); + bank.setFloat(TimeBasedTrkg_TBCovMat.C66, i, (float) CM[5][5]); } } @@ -608,24 +619,24 @@ private DataBank fillTBClustersBank(DataEvent event, List cluslis } double chi2 = 0; - bank.setShort("id", i, (short) cluslist.get(i).get_Id()); + bank.setShort(TimeBasedTrkg_TBClusters.id, i, (short) cluslist.get(i).get_Id()); // int status =0; // if(cluslist.get(i).size()<6) // status = 1; - bank.setShort("status", i, (short) 0); - bank.setByte("superlayer", i, (byte) cluslist.get(i).get_Superlayer()); - bank.setByte("sector", i, (byte) cluslist.get(i).get_Sector()); + bank.setShort(TimeBasedTrkg_TBClusters.status, i, (short) 0); + bank.setByte(TimeBasedTrkg_TBClusters.superlayer, i, (byte) cluslist.get(i).get_Superlayer()); + bank.setByte(TimeBasedTrkg_TBClusters.sector, i, (byte) cluslist.get(i).get_Sector()); - bank.setFloat("avgWire", i, (float) cluslist.get(i).getAvgwire()); - bank.setByte("size", i, (byte) cluslist.get(i).size()); + bank.setFloat(TimeBasedTrkg_TBClusters.avgWire, i, (float) cluslist.get(i).getAvgwire()); + bank.setByte(TimeBasedTrkg_TBClusters.size, i, (byte) cluslist.get(i).size()); double fitSlope = cluslist.get(i).get_clusterLineFitSlope(); double fitInterc = cluslist.get(i).get_clusterLineFitIntercept(); - bank.setFloat("fitSlope", i, (float) fitSlope); - bank.setFloat("fitSlopeErr", i, (float) cluslist.get(i).get_clusterLineFitSlopeErr()); - bank.setFloat("fitInterc", i, (float) fitInterc); - bank.setFloat("fitIntercErr", i, (float) cluslist.get(i).get_clusterLineFitInterceptErr()); + bank.setFloat(TimeBasedTrkg_TBClusters.fitSlope, i, (float) fitSlope); + bank.setFloat(TimeBasedTrkg_TBClusters.fitSlopeErr, i, (float) cluslist.get(i).get_clusterLineFitSlopeErr()); + bank.setFloat(TimeBasedTrkg_TBClusters.fitInterc, i, (float) fitInterc); + bank.setFloat(TimeBasedTrkg_TBClusters.fitIntercErr, i, (float) cluslist.get(i).get_clusterLineFitInterceptErr()); for (int j = 0; j < cluslist.get(i).size(); j++) { if (j < hitIdxArray.length) { @@ -635,7 +646,7 @@ private DataBank fillTBClustersBank(DataEvent event, List cluslis double residual = cluslist.get(i).get(j).get_ClusFitDoca() / (cluslist.get(i).get(j).get_CellSize() / Math.sqrt(12.)); chi2 += residual * residual; } - bank.setFloat("fitChisqProb", i, (float) ProbChi2perNDF.prob(chi2, cluslist.get(i).size() - 2)); + bank.setFloat(TimeBasedTrkg_TBClusters.fitChisqProb, i, (float) ProbChi2perNDF.prob(chi2, cluslist.get(i).size() - 2)); for (int j = 0; j < hitIdxArray.length; j++) { String hitStrg = "Hit"; @@ -671,25 +682,25 @@ private DataBank fillTBSegmentsBank(DataEvent event, List seglist) { double chi2 = 0; - bank.setShort("id", i, (short) seglist.get(i).get_Id()); - bank.setShort("status", i, (short) seglist.get(i).get_Status()); - bank.setByte("superlayer", i, (byte) seglist.get(i).get_Superlayer()); - bank.setByte("sector", i, (byte) seglist.get(i).get_Sector()); + bank.setShort(TimeBasedTrkg_TBSegments.id, i, (short) seglist.get(i).get_Id()); + bank.setShort(TimeBasedTrkg_TBSegments.status, i, (short) seglist.get(i).get_Status()); + bank.setByte(TimeBasedTrkg_TBSegments.superlayer, i, (byte) seglist.get(i).get_Superlayer()); + bank.setByte(TimeBasedTrkg_TBSegments.sector, i, (byte) seglist.get(i).get_Sector()); FittedCluster cls = seglist.get(i).get_fittedCluster(); - bank.setShort("Cluster_ID", i, (short) cls.get_Id()); - - bank.setFloat("avgWire", i, (float) cls.getAvgwire()); - bank.setByte("size", i, (byte) seglist.get(i).size()); - bank.setFloat("fitSlope", i, (float) cls.get_clusterLineFitSlope()); - bank.setFloat("fitSlopeErr", i, (float) cls.get_clusterLineFitSlopeErr()); - bank.setFloat("fitInterc", i, (float) cls.get_clusterLineFitIntercept()); - bank.setFloat("fitIntercErr", i, (float) cls.get_clusterLineFitInterceptErr()); - bank.setFloat("resiSum", i, (float) seglist.get(i).get_ResiSum()); - bank.setFloat("timeSum", i, (float) seglist.get(i).get_TimeSum()); - bank.setFloat("SegEndPoint1X", i, (float) seglist.get(i).get_SegmentEndPoints()[0]); - bank.setFloat("SegEndPoint1Z", i, (float) seglist.get(i).get_SegmentEndPoints()[1]); - bank.setFloat("SegEndPoint2X", i, (float) seglist.get(i).get_SegmentEndPoints()[2]); - bank.setFloat("SegEndPoint2Z", i, (float) seglist.get(i).get_SegmentEndPoints()[3]); + bank.setShort(TimeBasedTrkg_TBSegments.Cluster_ID, i, (short) cls.get_Id()); + + bank.setFloat(TimeBasedTrkg_TBSegments.avgWire, i, (float) cls.getAvgwire()); + bank.setByte(TimeBasedTrkg_TBSegments.size, i, (byte) seglist.get(i).size()); + bank.setFloat(TimeBasedTrkg_TBSegments.fitSlope, i, (float) cls.get_clusterLineFitSlope()); + bank.setFloat(TimeBasedTrkg_TBSegments.fitSlopeErr, i, (float) cls.get_clusterLineFitSlopeErr()); + bank.setFloat(TimeBasedTrkg_TBSegments.fitInterc, i, (float) cls.get_clusterLineFitIntercept()); + bank.setFloat(TimeBasedTrkg_TBSegments.fitIntercErr, i, (float) cls.get_clusterLineFitInterceptErr()); + bank.setFloat(TimeBasedTrkg_TBSegments.resiSum, i, (float) seglist.get(i).get_ResiSum()); + bank.setFloat(TimeBasedTrkg_TBSegments.timeSum, i, (float) seglist.get(i).get_TimeSum()); + bank.setFloat(TimeBasedTrkg_TBSegments.SegEndPoint1X, i, (float) seglist.get(i).get_SegmentEndPoints()[0]); + bank.setFloat(TimeBasedTrkg_TBSegments.SegEndPoint1Z, i, (float) seglist.get(i).get_SegmentEndPoints()[1]); + bank.setFloat(TimeBasedTrkg_TBSegments.SegEndPoint2X, i, (float) seglist.get(i).get_SegmentEndPoints()[2]); + bank.setFloat(TimeBasedTrkg_TBSegments.SegEndPoint2Z, i, (float) seglist.get(i).get_SegmentEndPoints()[3]); for (int j = 0; j < seglist.get(i).size(); j++) { if (j < hitIdxArray.length) { @@ -699,7 +710,7 @@ private DataBank fillTBSegmentsBank(DataEvent event, List seglist) { double residual = seglist.get(i).get(j).get_ClusFitDoca() / (seglist.get(i).get(j).get_CellSize() / Math.sqrt(12.)); chi2 += residual * residual; } - bank.setFloat("fitChisqProb", i, (float) ProbChi2perNDF.prob(chi2, seglist.get(i).size() - 2)); + bank.setFloat(TimeBasedTrkg_TBSegments.fitChisqProb, i, (float) ProbChi2perNDF.prob(chi2, seglist.get(i).size() - 2)); for (int j = 0; j < hitIdxArray.length; j++) { String hitStrg = "Hit"; @@ -731,24 +742,24 @@ private DataBank fillTBCrossesBank(DataEvent event, List crosslist) { int index=0; for (Cross aCrosslist : crosslist) { if (aCrosslist.get_Id() != -1) { - bank.setShort("id", index, (short) aCrosslist.get_Id()); - bank.setShort("status", index, (short) (aCrosslist.get_Segment1().get_Status() + aCrosslist.get_Segment2().get_Status())); - bank.setByte("sector", index, (byte) aCrosslist.get_Sector()); - bank.setByte("region", index, (byte) aCrosslist.get_Region()); - bank.setFloat("x", index, (float) aCrosslist.get_Point().x()); - bank.setFloat("y", index, (float) aCrosslist.get_Point().y()); - bank.setFloat("z", index, (float) aCrosslist.get_Point().z()); - bank.setFloat("err_x", index, (float) aCrosslist.get_PointErr().x()); - bank.setFloat("err_y", index, (float) aCrosslist.get_PointErr().y()); - bank.setFloat("err_z", index, (float) aCrosslist.get_PointErr().z()); - bank.setFloat("ux", index, (float) aCrosslist.get_Dir().x()); - bank.setFloat("uy", index, (float) aCrosslist.get_Dir().y()); - bank.setFloat("uz", index, (float) aCrosslist.get_Dir().z()); - bank.setFloat("err_ux", index, (float) aCrosslist.get_DirErr().x()); - bank.setFloat("err_uy", index, (float) aCrosslist.get_DirErr().y()); - bank.setFloat("err_uz", index, (float) aCrosslist.get_DirErr().z()); - bank.setShort("Segment1_ID", index, (short) aCrosslist.get_Segment1().get_Id()); - bank.setShort("Segment2_ID", index, (short) aCrosslist.get_Segment2().get_Id()); + bank.setShort(TimeBasedTrkg_TBCrosses.id, index, (short) aCrosslist.get_Id()); + bank.setShort(TimeBasedTrkg_TBCrosses.status, index, (short) (aCrosslist.get_Segment1().get_Status() + aCrosslist.get_Segment2().get_Status())); + bank.setByte(TimeBasedTrkg_TBCrosses.sector, index, (byte) aCrosslist.get_Sector()); + bank.setByte(TimeBasedTrkg_TBCrosses.region, index, (byte) aCrosslist.get_Region()); + bank.setFloat(TimeBasedTrkg_TBCrosses.x, index, (float) aCrosslist.get_Point().x()); + bank.setFloat(TimeBasedTrkg_TBCrosses.y, index, (float) aCrosslist.get_Point().y()); + bank.setFloat(TimeBasedTrkg_TBCrosses.z, index, (float) aCrosslist.get_Point().z()); + bank.setFloat(TimeBasedTrkg_TBCrosses.err_x, index, (float) aCrosslist.get_PointErr().x()); + bank.setFloat(TimeBasedTrkg_TBCrosses.err_y, index, (float) aCrosslist.get_PointErr().y()); + bank.setFloat(TimeBasedTrkg_TBCrosses.err_z, index, (float) aCrosslist.get_PointErr().z()); + bank.setFloat(TimeBasedTrkg_TBCrosses.ux, index, (float) aCrosslist.get_Dir().x()); + bank.setFloat(TimeBasedTrkg_TBCrosses.uy, index, (float) aCrosslist.get_Dir().y()); + bank.setFloat(TimeBasedTrkg_TBCrosses.uz, index, (float) aCrosslist.get_Dir().z()); + bank.setFloat(TimeBasedTrkg_TBCrosses.err_ux, index, (float) aCrosslist.get_DirErr().x()); + bank.setFloat(TimeBasedTrkg_TBCrosses.err_uy, index, (float) aCrosslist.get_DirErr().y()); + bank.setFloat(TimeBasedTrkg_TBCrosses.err_uz, index, (float) aCrosslist.get_DirErr().z()); + bank.setShort(TimeBasedTrkg_TBCrosses.Segment1_ID, index, (short) aCrosslist.get_Segment1().get_Id()); + bank.setShort(TimeBasedTrkg_TBCrosses.Segment2_ID, index, (short) aCrosslist.get_Segment2().get_Id()); index++; } } @@ -765,45 +776,42 @@ private DataBank fillTBTracksBank(DataEvent event, List candlist) { String name = bankNames.getTracksBank(); DataBank bank = event.createBank(name, candlist.size()); for (int i = 0; i < candlist.size(); i++) { - bank.setShort("id", i, (short) candlist.get(i).get_Id()); - - //bank.setShort("status", i, (short) (100+candlist.get(i).get_Status()*10+candlist.get(i).get_MissingSuperlayer())); - bank.setShort("status", i, (short) candlist.get(i).getBitStatus()); - bank.setByte("sector", i, (byte) candlist.get(i).getSector()); - bank.setByte("q", i, (byte) candlist.get(i).get_Q()); - //bank.setFloat("p", i, (float) candlist.get(i).get_P()); + bank.setShort(TimeBasedTrkg_TBTracks.id, i, (short) candlist.get(i).get_Id()); + bank.setShort(TimeBasedTrkg_TBTracks.status, i, (short) candlist.get(i).getBitStatus()); + bank.setByte(TimeBasedTrkg_TBTracks.sector, i, (byte) candlist.get(i).getSector()); + bank.setByte(TimeBasedTrkg_TBTracks.q, i, (byte) candlist.get(i).get_Q()); if(candlist.get(i).get_PreRegion1CrossPoint()!=null) { - bank.setFloat("c1_x", i, (float) candlist.get(i).get_PreRegion1CrossPoint().x()); - bank.setFloat("c1_y", i, (float) candlist.get(i).get_PreRegion1CrossPoint().y()); - bank.setFloat("c1_z", i, (float) candlist.get(i).get_PreRegion1CrossPoint().z()); - bank.setFloat("c1_ux", i, (float) candlist.get(i).get_PreRegion1CrossDir().x()); - bank.setFloat("c1_uy", i, (float) candlist.get(i).get_PreRegion1CrossDir().y()); - bank.setFloat("c1_uz", i, (float) candlist.get(i).get_PreRegion1CrossDir().z()); + bank.setFloat(TimeBasedTrkg_TBTracks.c1_x, i, (float) candlist.get(i).get_PreRegion1CrossPoint().x()); + bank.setFloat(TimeBasedTrkg_TBTracks.c1_y, i, (float) candlist.get(i).get_PreRegion1CrossPoint().y()); + bank.setFloat(TimeBasedTrkg_TBTracks.c1_z, i, (float) candlist.get(i).get_PreRegion1CrossPoint().z()); + bank.setFloat(TimeBasedTrkg_TBTracks.c1_ux, i, (float) candlist.get(i).get_PreRegion1CrossDir().x()); + bank.setFloat(TimeBasedTrkg_TBTracks.c1_uy, i, (float) candlist.get(i).get_PreRegion1CrossDir().y()); + bank.setFloat(TimeBasedTrkg_TBTracks.c1_uz, i, (float) candlist.get(i).get_PreRegion1CrossDir().z()); } if(candlist.get(i).get_PostRegion3CrossPoint()!=null) { - bank.setFloat("c3_x", i, (float) candlist.get(i).get_PostRegion3CrossPoint().x()); - bank.setFloat("c3_y", i, (float) candlist.get(i).get_PostRegion3CrossPoint().y()); - bank.setFloat("c3_z", i, (float) candlist.get(i).get_PostRegion3CrossPoint().z()); - bank.setFloat("c3_ux", i, (float) candlist.get(i).get_PostRegion3CrossDir().x()); - bank.setFloat("c3_uy", i, (float) candlist.get(i).get_PostRegion3CrossDir().y()); - bank.setFloat("c3_uz", i, (float) candlist.get(i).get_PostRegion3CrossDir().z()); + bank.setFloat(TimeBasedTrkg_TBTracks.c3_x, i, (float) candlist.get(i).get_PostRegion3CrossPoint().x()); + bank.setFloat(TimeBasedTrkg_TBTracks.c3_y, i, (float) candlist.get(i).get_PostRegion3CrossPoint().y()); + bank.setFloat(TimeBasedTrkg_TBTracks.c3_z, i, (float) candlist.get(i).get_PostRegion3CrossPoint().z()); + bank.setFloat(TimeBasedTrkg_TBTracks.c3_ux, i, (float) candlist.get(i).get_PostRegion3CrossDir().x()); + bank.setFloat(TimeBasedTrkg_TBTracks.c3_uy, i, (float) candlist.get(i).get_PostRegion3CrossDir().y()); + bank.setFloat(TimeBasedTrkg_TBTracks.c3_uz, i, (float) candlist.get(i).get_PostRegion3CrossDir().z()); } if(candlist.get(i).get_Region1TrackX()!=null) { - bank.setFloat("t1_x", i, (float) candlist.get(i).get_Region1TrackX().x()); - bank.setFloat("t1_y", i, (float) candlist.get(i).get_Region1TrackX().y()); - bank.setFloat("t1_z", i, (float) candlist.get(i).get_Region1TrackX().z()); - bank.setFloat("t1_px", i, (float) candlist.get(i).get_Region1TrackP().x()); - bank.setFloat("t1_py", i, (float) candlist.get(i).get_Region1TrackP().y()); - bank.setFloat("t1_pz", i, (float) candlist.get(i).get_Region1TrackP().z()); + bank.setFloat(TimeBasedTrkg_TBTracks.t1_x, i, (float) candlist.get(i).get_Region1TrackX().x()); + bank.setFloat(TimeBasedTrkg_TBTracks.t1_y, i, (float) candlist.get(i).get_Region1TrackX().y()); + bank.setFloat(TimeBasedTrkg_TBTracks.t1_z, i, (float) candlist.get(i).get_Region1TrackX().z()); + bank.setFloat(TimeBasedTrkg_TBTracks.t1_px, i, (float) candlist.get(i).get_Region1TrackP().x()); + bank.setFloat(TimeBasedTrkg_TBTracks.t1_py, i, (float) candlist.get(i).get_Region1TrackP().y()); + bank.setFloat(TimeBasedTrkg_TBTracks.t1_pz, i, (float) candlist.get(i).get_Region1TrackP().z()); } - bank.setFloat("pathlength", i, (float) candlist.get(i).get_TotPathLen()); - bank.setFloat("Vtx0_x", i, (float) candlist.get(i).get_Vtx0().x()); - bank.setFloat("Vtx0_y", i, (float) candlist.get(i).get_Vtx0().y()); - bank.setFloat("Vtx0_z", i, (float) candlist.get(i).get_Vtx0().z()); - bank.setFloat("p0_x", i, (float) candlist.get(i).get_pAtOrig().x()); - bank.setFloat("p0_y", i, (float) candlist.get(i).get_pAtOrig().y()); - bank.setFloat("p0_z", i, (float) candlist.get(i).get_pAtOrig().z()); + bank.setFloat(TimeBasedTrkg_TBTracks.pathlength, i, (float) candlist.get(i).get_TotPathLen()); + bank.setFloat(TimeBasedTrkg_TBTracks.Vtx0_x, i, (float) candlist.get(i).get_Vtx0().x()); + bank.setFloat(TimeBasedTrkg_TBTracks.Vtx0_y, i, (float) candlist.get(i).get_Vtx0().y()); + bank.setFloat(TimeBasedTrkg_TBTracks.Vtx0_z, i, (float) candlist.get(i).get_Vtx0().z()); + bank.setFloat(TimeBasedTrkg_TBTracks.p0_x, i, (float) candlist.get(i).get_pAtOrig().x()); + bank.setFloat(TimeBasedTrkg_TBTracks.p0_y, i, (float) candlist.get(i).get_pAtOrig().y()); + bank.setFloat(TimeBasedTrkg_TBTracks.p0_z, i, (float) candlist.get(i).get_pAtOrig().z()); //fill associated IDs for(int r = 0; r < 3; r++) { bank.setShort("Cross"+String.valueOf(r+1)+"_ID", @@ -825,7 +833,7 @@ private DataBank fillTBTracksBank(DataEvent event, List candlist) { bank.setShort("Cluster"+String.valueOf(candlist.get(i).getSingleSuperlayer().get_Superlayer())+"_ID", i, (short) candlist.get(i).getSingleSuperlayer().get_fittedCluster().get_Id()); } - bank.setFloat("chi2", i, (float) candlist.get(i).get_FitChi2()); + bank.setFloat(TimeBasedTrkg_TBTracks.chi2, i, (float) candlist.get(i).get_FitChi2()); // To not interrupt current type of ndf, ndf weighted by DAF is converted from float to interger int ndfDAF = 999; if(candlist.get(i).get_NDFDAF() > 0){ @@ -834,9 +842,9 @@ private DataBank fillTBTracksBank(DataEvent event, List candlist) { else if (candlist.get(i).get_NDFDAF() < 0){ ndfDAF = (int) Math.floor(candlist.get(i).get_NDFDAF()); } - bank.setShort("ndf", i, (short) ndfDAF); + bank.setShort(TimeBasedTrkg_TBTracks.ndf, i, (short) ndfDAF); // ndf0 is for traditional ndf for the track; # of hits can be obtained through it - bank.setShort("ndf0", i, (short) candlist.get(i).get_FitNDF()); + bank.setShort(TimeBasedTrkg_TBTracks.ndf0, i, (short) candlist.get(i).get_FitNDF()); } return bank; @@ -863,20 +871,20 @@ private DataBank fillTrajectoryBank(DataEvent event, List tracks) { if (track.getTrajectory().get(j).getDetector() == DetectorType.DC.getDetectorId() && (track.getTrajectory().get(j).getLayer() - 6) % 6 != 0) continue; // save the last layer in a superlayer - bank.setShort("id", i1, (short) track.get_Id()); - bank.setByte("detector", i1, (byte) track.getTrajectory().get(j).getDetector()); - bank.setByte("sector", i1, (byte) track.getSector()); - bank.setByte("layer", i1, (byte) track.getTrajectory().get(j).getLayer()); - bank.setFloat("x", i1, (float) track.getTrajectory().get(j).getPoint().x()); - bank.setFloat("y", i1, (float) track.getTrajectory().get(j).getPoint().y()); - bank.setFloat("z", i1, (float) track.getTrajectory().get(j).getPoint().z()); - bank.setFloat("tx", i1, (float) track.getTrajectory().get(j).getDirection().x()); - bank.setFloat("ty", i1, (float) track.getTrajectory().get(j).getDirection().y()); - bank.setFloat("tz", i1, (float) track.getTrajectory().get(j).getDirection().z()); - bank.setFloat("B", i1, (float) track.getTrajectory().get(j).getiBdl()); - bank.setFloat("path", i1, (float) track.getTrajectory().get(j).getPath()); - bank.setFloat("dx", i1, (float) track.getTrajectory().get(j).getDx()); - bank.setFloat("edge", i1, (float) track.getTrajectory().get(j).getEdge()); + bank.setShort(TimeBasedTrkg_Trajectory.id, i1, (short) track.get_Id()); + bank.setByte(TimeBasedTrkg_Trajectory.detector, i1, (byte) track.getTrajectory().get(j).getDetector()); + bank.setByte(TimeBasedTrkg_Trajectory.sector, i1, (byte) track.getSector()); + bank.setByte(TimeBasedTrkg_Trajectory.layer, i1, (byte) track.getTrajectory().get(j).getLayer()); + bank.setFloat(TimeBasedTrkg_Trajectory.x, i1, (float) track.getTrajectory().get(j).getPoint().x()); + bank.setFloat(TimeBasedTrkg_Trajectory.y, i1, (float) track.getTrajectory().get(j).getPoint().y()); + bank.setFloat(TimeBasedTrkg_Trajectory.z, i1, (float) track.getTrajectory().get(j).getPoint().z()); + bank.setFloat(TimeBasedTrkg_Trajectory.tx, i1, (float) track.getTrajectory().get(j).getDirection().x()); + bank.setFloat(TimeBasedTrkg_Trajectory.ty, i1, (float) track.getTrajectory().get(j).getDirection().y()); + bank.setFloat(TimeBasedTrkg_Trajectory.tz, i1, (float) track.getTrajectory().get(j).getDirection().z()); + bank.setFloat(TimeBasedTrkg_Trajectory.B, i1, (float) track.getTrajectory().get(j).getiBdl()); + bank.setFloat(TimeBasedTrkg_Trajectory.path, i1, (float) track.getTrajectory().get(j).getPath()); + bank.setFloat(TimeBasedTrkg_Trajectory.dx, i1, (float) track.getTrajectory().get(j).getDx()); + bank.setFloat(TimeBasedTrkg_Trajectory.edge, i1, (float) track.getTrajectory().get(j).getEdge()); i1++; } }