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

LP new signal simulation #20

Merged
merged 16 commits into from
Apr 3, 2024
Merged

LP new signal simulation #20

merged 16 commits into from
Apr 3, 2024

Conversation

lopezzot
Copy link
Member

@lopezzot lopezzot commented Apr 3, 2024

This PR includes several improvements in signal simulation and event tagging.

  • Add geometry constants for small silicon cells (e4e6ff6)
  • Include signal in silicon half cells (34e76a9)
  • Include class to compute si cell position and tag events (c4ed633 53b3761 9641134 9f6dc6d d7f97ee 7fdd87d 1de6dff c954cc4)
  • Save primary particle info in ROOT file (ba79b0e)
  • Add version control for method checking geometry overlaps (db9f6b6)
  • Add method (not used) for simplified Birk treatment (c08653f)
  • Fix bug in CMS Birk treatment by using correct unit in Birk parameters (5c7dd13)
  • Update MIP calibration constants for silicon and tiles (6d5f944)
  • Update ROOT analysis for pion data. Use new variables to tag CEE and CHE events, and new gaussian fitting technique (54f927d)

Add apothem and radius of small silicon cells.
Take into account also half cells in silicon wafers by using the same
sensitive detector of full cells. Extend number of cells to 132 for both
CEE and CEH.
Include auxiliary class to calculate x and y position of small si cells.
Add functions to calculate MIP signal per layer and tag pion interacting
layer in CEE according to CMS publication.
Access PDGID and energy of primary particle using the PrimaryGenAction
and save info in root format.
Use SignalHelper class methods to tag interacting pions in CEE and save
interaction layer in ROOT output file.
Add version control for overlap check method. Define it only for G4-11.0
and up.
In order to find layer of corresponding nuclear interaction in CEE, pass
vector of signals already calibrated at the MIP scale. Use an auxiliary
lambda function for in EventAction.
Add printout information for SignalHelper IsInteraction() method.
To identify CEE interaction layer check also if the layer has at least 3
Si cells above threshold.
Apply the same stragy to tag interacting pions in CHE as done in CEE.
Simply consider the central Si pad for CHE layers. Save corresponding
layer number in ROOT file.
Fix bug when accessing signals per CHE layer in EventAction EndOfEvent
method.
Add simplified method for Birk correction. To be used only for
comparison with the CMS one.
Add units for Brik1 coefficient as done in CMSSW. It restores
a-dimensionality of Birks Law. Remove density convertion to g/cm3 for
consistency.
Update value of MIP calibration constants for silicon and tiles.
Discard events interacting in the first two layers of CEE, or not
interacing both in CEE and CHE or with interacting layer not found. Use
CEEInteracted variable to select CEE and CHE events. Perform gaussian fit to
extract response as the pick value of distributions (as written in TB
paper) by fitting around maximum of histograms.
@lopezzot
Copy link
Member Author

lopezzot commented Apr 3, 2024

Accepting PR #20.

@lopezzot lopezzot merged commit d1d1dd1 into main Apr 3, 2024
@lopezzot lopezzot deleted the lp-mapcells branch April 3, 2024 15:19
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.

1 participant