Skip to content
This repository has been archived by the owner on May 18, 2023. It is now read-only.

Commit

Permalink
Fix to get cosmics; bypassing IBdl calculation to reconstruct straigh…
Browse files Browse the repository at this point in the history
…t tracks.
  • Loading branch information
zieglerv committed Jan 20, 2018
1 parent ef11315 commit 43c0b00
Showing 1 changed file with 11 additions and 11 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -76,18 +76,26 @@ public List<Track> getTrackCands(CrossList crossList, DCGeant4Factory DcDetector
continue;

if(crossesInTrk.size()==3 && this.PassNSuperlayerTracking(crossesInTrk)==true) {

cand.addAll(crossesInTrk);

cand.set_Sector(crossesInTrk.get(0).get_Sector());

if(Math.abs(TORSCALE)<0.001) {
//no field --> fit straight track
this.getStraightTrack(cand);
if(cand.get_pAtOrig()!=null) {
cand.set_Id(cands.size()+1);
cands.add(cand);
}
}
//cand.set_Region3CrossPoint();
//cand.set_Region3CrossDir();

cand.set_Trajectory(traj.get_Trajectory());
cand.set_IntegralBdl(traj.get_IntegralBdl());


if(cand.size()==3) {
double theta3 = Math.atan(cand.get(2).get_Segment2().get_fittedCluster().get_clusterLineFitSlope());
double theta1 = Math.atan(cand.get(0).get_Segment2().get_fittedCluster().get_clusterLineFitSlope());
Expand All @@ -100,15 +108,7 @@ public List<Track> getTrackCands(CrossList crossList, DCGeant4Factory DcDetector

double iBdl = traj.get_IntegralBdl();

if(iBdl == 0) {
//no field --> fit straight track
//System.out.println(trking+" FitChisq "+cand.get_FitChi2());
this.getStraightTrack(cand);
if(cand.get_pAtOrig()!=null) {
cand.set_Id(cands.size()+1);
cands.add(cand);
}
}

if(iBdl != 0 || (deltaTheta != 0)) {

double pxz = Math.abs(Constants.LIGHTVEL*iBdl/deltaTheta);
Expand Down

0 comments on commit 43c0b00

Please sign in to comment.