Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Add baseline time calculation for HGCal clusters #15940

Merged
merged 1 commit into from Sep 23, 2016
Merged

Add baseline time calculation for HGCal clusters #15940

Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
19 changes: 18 additions & 1 deletion RecoParticleFlow/PFClusterProducer/src/Cluster3DPCACalculator.cc
View file
Open in desktop
Expand Up @@ -60,15 +60,25 @@ calculateAndSetPositionActual(reco::PFCluster& cluster) {
<< " Found a cluster with no seed: " << cluster;
}
double cl_energy = 0;
double max_e = 0.0;
double max_e = 0.0;
double avg_time = 0.0;
double time_norm = 0.0;
PFLayer::Layer max_e_layer = PFLayer::NONE;
reco::PFRecHitRef refseed;
double pcavars[3];

for( const reco::PFRecHitFraction& rhf : cluster.recHitFractions() ) {
const reco::PFRecHitRef& refhit = rhf.recHitRef();
double rh_energy = refhit->energy();
double rh_time = refhit->time();
cl_energy += rh_energy * rhf.fraction();
if( rh_time > 0.0 ) { // time == -1 means no measurement
// all times are offset by one nanosecond in digitizer
// remove that here so all times of flight
// are with respect to (0,0,0)
avg_time += (rh_time - 1.0);
time_norm += 1.0;
}
if( rh_energy > max_e ) {
max_e = rh_energy;
max_e_layer = rhf.recHitRef()->layer();
Expand Down Expand Up @@ -102,10 +112,17 @@ calculateAndSetPositionActual(reco::PFCluster& cluster) {
math::XYZPoint barycenter(means[0],means[1],means[2]);
math::XYZVector axis(eigens(0,0),eigens(1,0),eigens(2,0));

if( time_norm > 0.0 ) {
avg_time = avg_time/time_norm;
} else {
avg_time = std::numeric_limits<double>::min();
}

if( axis.z()*barycenter.z() < 0.0 ) {
axis = math::XYZVector(-eigens(0,0),-eigens(1,0),-eigens(2,0));
}

cluster.setTime(avg_time);
cluster.setPosition(barycenter);
cluster.calculatePositionREP();

Expand Down