sPHENIX discussion of software

["Huang, Jin" <jhuang AT bnl.gov>]

Also available on agenda: https://indico.bnl.gov/conferenceDisplay.py?confId=1663

 

Software update

 

·       sPHENIX default switched to new compiler of GCC 4.8. Several commit to coresoftware and analysis repository to fix compatiblity issues. Please report on any problem.

·       Pull request #101: Mike introduced shower truth compression last week via this pull request. Jin is assigned to cross check. Check not yet finished.
Action Item [Jin]: finish checking and merge as soon as possible

·       Pull request #103: Nils is fixing a style problem in the forward crystal calorimeter. Consensus is to merge once Nils submit final fix following the discussion on the GitHub thread.   

 

 

Standardized calorimeter QA: Jin Huang

 

Page 4/6:

Chris: for higher level check (E/p, jets), suggest separating from the Geant4 level hit checks for two purposes 1) separate G4-level and reco-level debugging 2) allow access larger statistics via Geant4 level production. For example. we can buffer latest production DST of PHG4hits, and run jet reco code on top of it to check jet reconstruction.

Jin: Good point. Adopt this approach.

 

 

Checking single particle simulations using Jet finder: Sarah Campbell

 

Sarah: this is still based on Dec 2015 Hcal workfest production sample. Not newer ones via lastest Geant4 version.

 

Page 3: analyze single particle events using jet finders

Sarah:

1) This is based on jet reco using calorimeter towers as input

2) dPhi and dEta cuts are difference between reconstructed “jet” and truth particle’s angle

 

Chris: Why 10^6 miss at low energy on black curve?

Sarah/Jin: “noisy” jets from towers with pedestal hits

 

Page 4:

Chris: there is a known bug that Hcal towers were shifted by a fraction of tower size in phi.

Action Item [Chris]: fix this shift

 

Page 5:

Jin: besides the charged bending, the shift in azimuthal direction also depends on the depth of the shower via tilted tile tower geometry. The shower depth in the calorimeter can be accessed via PHG4hits, such as in HCALAnalysis::process_event(PHCompositeNode *topNode):

PHG4HitContainer::ConstRange hcalout_hit_range = _hcalout_hit_container->getHits();

   70  for (PHG4HitContainer::ConstIterator hit_iter = hcalout_hit_range.first;

   71  hit_iter != hcalout_hit_range.second; hit_iter++){

   72 

   73  PHG4Hit *this_hit = hit_iter->second;

   74  assert(this_hit);

   75 ....

   78 

   79  float rin = sqrt(pow(this_hit->get_x(0),2)+pow(this_hit->get_y(0),2));

   80  float rout = sqrt(pow(this_hit->get_x(1),2)+pow(this_hit->get_y(1),2));

   81 

   82  lcg += this_hit->get_edep()*(rin + ((rout-rin)/2.0));

   83  hcalout_lcg += this_hit->get_edep()*(rin + ((rout-rin)/2.0));    

   86  }

   87 

Page 9:

Sarah: this is based on Outter HCal only

 

Additional discussions:

We should fix geometry problem on Hcal first to avoid confusions

Expect shower energy dependence via depth-phi correlation.

Good to quantify a charge dependence of dPhi width, via coincidence between bending of the track and tilting of the tower.

 

 

Next meeting:

Murad (via email): plan to present prototype simulation

Jin: plan to generate pull request for calorimeter QA and check newer Geant4 production sample.

 

 

______________________________

 

Jin HUANG

 

Brookhaven National Laboratory

Physics Department, Bldg 510 C

Upton, NY 11973-5000

 

Office: 631-344-5898

Cell:   757-604-9946

______________________________

 

 

> -----Original Message-----

> From: sphenix-l-bounces AT lists.bnl.gov [mailto:sphenix-l-

> bounces AT lists.bnl.gov] On Behalf Of Huang, Jin

> Sent: Monday, January 11, 2016 1:57 PM

> To: sphenix-l AT lists.bnl.gov; sphenix-software-l AT lists.bnl.gov

> Subject: [Sphenix-l] [One-time room 3-192 ] sPHENIX simulation meetings,

> Tue Jan 12th, 1PM EST

>

> Dear All,

>

> Although this is the week of Jet workshop, we will have a brief sPHENIX

> simulation meetings, Tue Jan 12th, 1PM EST.

>

> Currently on the agenda: following software update, I will talk about a QA

> service for calorimeters, which was an action item of the last meeting. More

> contributions will be welcomed too.

>

> Please note a one-time new room (BNL 510C 3-192) will be used for this

> meeting.

>

> The agenda page is at

> https://indico.bnl.gov/conferenceDisplay.py?confId=1645.

> *            For residents at BNL, we will meet at the 3-192 conference room

> *            To join the Meeting via phone bridge:

> https://bluejeans.com/619960505

> *            To join via Phone:

> 1) Dial:

>         +1.408.740.7256

>         +1.888.240.2560(US Toll Free)

>         +1.408.317.9253(Alternate number)

>         (see all numbers - http://bluejeans.com/numbers)

> 2) Enter Conference ID: 619960505

>

>

> Cheers,

>

> Jin

>

> ______________________________

>

> Jin HUANG

>

> Brookhaven National Laboratory

> Physics Department, Bldg 510 C

> Upton, NY 11973-5000

>

> Office: 631-344-5898

> Cell:   757-604-9946

> ______________________________

> _______________________________________________

> Sphenix-l mailing list

> Sphenix-l AT lists.bnl.gov

> https://lists.bnl.gov/mailman/listinfo/sphenix-l