@@ -414,18 +414,18 @@ void sbn::TrackCaloSkimmer::analyze(art::Event const& e)
414414 else if (hasPFPT0) whicht0 = 0 ;
415415 else if (hasCRTTrackT0)
416416 {
417- whicht0 = 1 ;
418- crtMatchingScore = t0CRTTrackScore;
417+ whicht0 = 1 ;
418+ crtMatchingScore = t0CRTTrackScore;
419419 }
420420 else if (hasCRTHitT0)
421421 {
422- whicht0 = 2 ;
423- crtMatchingScore = t0CRTHitScore;
422+ whicht0 = 2 ;
423+ crtMatchingScore = t0CRTHitScore;
424424 }
425425 else if (hasCRTSpacePointT0)
426426 {
427- whicht0 = 3 ;
428- crtMatchingScore = t0CRTSpacePointScore;
427+ whicht0 = 3 ;
428+ crtMatchingScore = t0CRTSpacePointScore;
429429 }
430430
431431 T0TimingInfo thisTrackTimingInfo = {t0PFP, t0CRTTrack, t0CRTHit, t0CRTSpacePoint, hasPFPT0, hasCRTTrackT0, hasCRTHitT0, hasCRTSpacePointT0, crtMatchingScore};
@@ -1186,26 +1186,26 @@ void sbn::TrackCaloSkimmer::FillTrack(const recob::Track &track,
11861186 const double driftv (dprop.DriftVelocity (dprop.Efield (), dprop.Temperature ()));
11871187
11881188 if (t0Info.hasT0CRTTrack )
1189- {
1190- const double xshift = driftDir*driftv*t0Info.t0CRTTrack *1e-3 ;
1191- fTrack ->start .x = track.Start ().X () + xshift;
1192- fTrack ->end .x = track.End ().X () + xshift;
1193- fTrack ->xShiftCRT = xshift;
1194- }
1189+ {
1190+ const double xshift = driftDir*driftv*t0Info.t0CRTTrack *1e-3 ;
1191+ fTrack ->start .x = track.Start ().X () + xshift;
1192+ fTrack ->end .x = track.End ().X () + xshift;
1193+ fTrack ->xShiftCRT = xshift;
1194+ }
11951195 else if (t0Info.hasT0CRTHit )
1196- {
1197- const double xshift = driftDir*driftv*t0Info.t0CRTHit *1e-3 ;
1198- fTrack ->start .x = track.Start ().X () + xshift;
1199- fTrack ->end .x = track.End ().X () + xshift;
1200- fTrack ->xShiftCRT = xshift;
1201- }
1196+ {
1197+ const double xshift = driftDir*driftv*t0Info.t0CRTHit *1e-3 ;
1198+ fTrack ->start .x = track.Start ().X () + xshift;
1199+ fTrack ->end .x = track.End ().X () + xshift;
1200+ fTrack ->xShiftCRT = xshift;
1201+ }
12021202 else if (t0Info.hasT0CRTSpacePoint )
1203- {
1204- const double xshift = driftDir*driftv*t0Info.t0CRTSpacePoint *1e-3 ;
1205- fTrack ->start .x = track.Start ().X () + xshift;
1206- fTrack ->end .x = track.End ().X () + xshift;
1207- fTrack ->xShiftCRT = xshift;
1208- }
1203+ {
1204+ const double xshift = driftDir*driftv*t0Info.t0CRTSpacePoint *1e-3 ;
1205+ fTrack ->start .x = track.Start ().X () + xshift;
1206+ fTrack ->end .x = track.End ().X () + xshift;
1207+ fTrack ->xShiftCRT = xshift;
1208+ }
12091209 }
12101210
12111211 if (hits.size () > 0 ) {
@@ -1468,11 +1468,11 @@ sbn::TrackHitInfo sbn::TrackCaloSkimmer::MakeHit(const recob::Hit &hit,
14681468 double time = std::numeric_limits<float >::signaling_NaN ();
14691469
14701470 if (t0Info.hasT0CRTTrack )
1471- time = t0Info.t0CRTTrack *1e-3 ;
1471+ time = t0Info.t0CRTTrack *1e-3 ;
14721472 else if (t0Info.hasT0CRTHit )
1473- time = t0Info.t0CRTHit *1e-3 ;
1473+ time = t0Info.t0CRTHit *1e-3 ;
14741474 else if (t0Info.hasT0CRTSpacePoint )
1475- time = t0Info.t0CRTSpacePoint *1e-3 ;
1475+ time = t0Info.t0CRTSpacePoint *1e-3 ;
14761476
14771477 double anodeDistance = (hit.PeakTime ()-dclock.Time2Tick (dclock.TriggerTime ())-time/dclock.TPCClock ().TickPeriod ())*dclock.TPCClock ().TickPeriod ()*driftv;
14781478 double wirePlaneX = wireReadout->Plane (geo::PlaneID (hit.WireID ().Cryostat , hit.WireID ().TPC , hit.WireID ().Plane )).GetCenter ().X ();
0 commit comments