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.

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

Add a steering file for monitoring using KF tracks and add a tool to remove hits from pattern reco in LCIO #935

Merged
merged 3 commits into from
Nov 28, 2022

Conversation

pbutti
Copy link
Contributor

@pbutti pbutti commented Nov 22, 2022

Title says all.

Copy link
Collaborator

@cbravo135 cbravo135 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Code added to SensorSetup should be safe since existing steering files won't have "SkipSensors" parameter list.

SiSensor sensor = (SiSensor) rawHit.getDetectorElement();
String name = sensor.getName();

if (sensorNames_.size() > 0) {
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This check is unnecessary. If sensorNames_ is empty, the whole for loop will be skipped.

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The loop is over all the hits, this check is on the list of sensors we want to skip. I do think sensorNames_ is potentially a confusing variable name since it is really the sensors to be skipped.

String name = sensor.getName();

if (sensorNames_.size() > 0) {

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think the following works as well ...

for (LCRelation fittedHit : fittedHits) { 

  RawTrackerHit rawHit = FittedRawTrackerHit.getRawTrackerHit(fittedHit);
  SiSensor sensor = (SiSensor) rawHit.getDetectorElement();
  if (sensorNames_.contains(sensor.name()) continue; 

  sensor.getReadout().addHit(fittedHit); 

}

Note the use of the contains method. When using a HashSet, this gives you O(1) lookup time vs O(n) when using an ArrayList or just iterating over the elements. HashSet is appropriate here since all of the sensor names will be unique. It also might be worth changing the name of the set sensorNames_ to maskedSensorNames_ so it's clearer what it does.

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The behavior of this will be different. contains behavior is different when called on an array vs a string. I think we do indeed want contains called on the string and not the array, which gives us easier control over which sensors we drop. This is only called once, so the complexity of the code doesn't really matter, and either way will effectively run in the same time. We would like to be able to remove multiple sensors with a single string if they all contain the same sub-string.

@@ -16,6 +18,8 @@ public class SensorSetup extends RawTrackerHitSensorSetup {

/// Name of the collection of fitted hits
private String fittedHitColName_ = "";

private List<String> sensorNames_ = new ArrayList<String>();
Copy link
Collaborator

@omar-moreno omar-moreno Nov 23, 2022

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Change this to a HashSet. You don't need to initialize the array since you do so in the setSkipSensors method below.

@@ -29,6 +33,10 @@ public SensorSetup() { }
*/
public void setFittedHitCollection(String fittedHitColName) { fittedHitColName_ = fittedHitColName; }

public void setSkipSensors(String[] sensorNames) {
this.sensorNames_ = new ArrayList<String>(Arrays.asList(sensorNames));
Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Change this to Hash Map.

On a separate note, I wonder if this method should check if sensorNames was already initializing before doing it again. I can see a situation where this gets set twice on accident in a steering file and could lead to unexpected behavior. Maybe if sensorNames has already been initialized, throw an exception?

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I think the framework already checks you only set a specific variable in the steering file once and throws an error if it is there more than once.

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Given the behavior we would like from this code, not clear to me HashMap really buys us anything.

Copy link
Collaborator

@omar-moreno omar-moreno left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I don't agree with these changes but will approve it anyway. The performance issues can be taken care of in the future.

@omar-moreno omar-moreno merged commit d20c28d into master Nov 28, 2022
This pull request was closed.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

3 participants