sPHENIX EMCal discussion

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

Hi, Guys

 

Following discussion on the EMCal meeting today, I am writing to send an example macro to read Chris’ test production data with sPHENIX software:

 

Location of the production data:

The last test production data with single particle per file is located here

·       rcf:/gpfs02/phenix/prod/sPHENIX/preCDR/pro.1-beta.2/spacal1d

·       rcf:/gpfs02/phenix/prod/sPHENIX/preCDR/pro.1-beta.2/spacal2d

Chris is running another new production (pro.1-beta.3) with the SPACAL tower code as mentioned in today’s meeting.

 

These files contain Geant4 hits and truth info. The following steps would run the default analysis software on them:

 

1.       Checkout the macro repo https://github.com/sPHENIX-Collaboration/macros

2.       Edit the Fun4All_G4_sPHENIX.C macro in order to:

a.       Read in the production file

b.      Turn on the evaluator for EMCal.

one example edit can would look like this, with changes highlighted:

 

int Fun4All_G4_sPHENIX(

                       const int nEvents = 10,

                       const char * inputFile = "/gpfs02/phenix/prod/sPHENIX/preCDR/pro.1-beta.2/spacal2d/fieldmap/G4Hits_sPHENIX_gamma_eta0_24GeV.root", // <<<<<<< One output with gamma of 24 GeV

                       const char * outputFile = "G4sPHENIXCells.root"

                       )

{

  //===============                                                                                                                    

  // Input options                                                                                                                      

  //===============                                                                                                                    

 

  // Either:                                                                                                                            

  // read previously generated g4-hits files, in this case it opens a DST and skips                                                    

  // the simulations step completely. The G4Setup macro is only loaded to get information                                              

  // about the number of layers used for the cell reco code                                                                            

  const bool readhits = true; //<<<<<<<<<<<<<<< Read production output

  // Or:                                                                                                                                

  // read files in HepMC format (typically output from event generators like hijing or pythia)                                         

  const bool readhepmc = false; // read HepMC files                                                                                    

  // Or:                                                                                                                               

  // Use particle generator                                                                                                            

 

  //======================                                                                                                              

  // What to run                                                                                                                       

  //======================                                                                                                              

 

  bool do_pipe = true;

 

  bool do_svtx = true;

  bool do_svtx_cell = true;

  bool do_svtx_track = true;

  bool do_svtx_eval = true;

 

  bool do_preshower = false;

 

  bool do_cemc = true;

  bool do_cemc_cell = true;

  bool do_cemc_twr = true;

  bool do_cemc_cluster = true;

  bool do_cemc_eval = true; ; //<<<<<<<<<<<<<<< Run evaluator.

......

 

3.       Now we run the code :

a.       Only if have to: use the production build, to avoid compatibility issue from future newer software. Right not I do not need this line to read the production data with the current build.
source /opt/sphenix/core/bin/sphenix_setup.csh -n pro.1  

b.      root  -b -q Fun4All_G4_sPHENIX.C

4.       Then 10 events should be processed in a few seconds, and an output evaluation file with clusters in them would become available:

a.       root g4cemc_eval.root

b.      root [6] ntp_cluster->Draw("eta","e") // cluster distribution in eta weighted by energy. Note energy here is about 2% of the incoming photon from the sampling fraction of EMCal.

 

Please let me know if there is any question or suggestion in running these data. You are welcome to bring the discussion to the simulation meeting too (Tomorrow 1PM EST https://indico.bnl.gov/conferenceDisplay.py?confId=1337 ).

 

Cheers

 

Jin

 

 

 

 

 

______________________________

 

Jin HUANG

 

Brookhaven National Laboratory

Physics Department, Bldg 510 C

Upton, NY 11973-5000

 

Office: 631-344-5898

Cell:   757-604-9946

______________________________

 

From: sphenix-emcal-l-bounces AT lists.bnl.gov [mailto:sphenix-emcal-l-bounces AT lists.bnl.gov] On Behalf Of Sickles, Anne M
Sent: Friday, August 28, 2015 12:46 PM
To: sphenix-emcal-l AT lists.bnl.gov
Subject: [Sphenix-emcal-l] bi-weekly sPHENIX EMCal meeting Monday 8/31

 

Hi,

 

We’ll meet again on Monday, August 31, at 4 PM ET for the sPHENIX EMCal meeting.  The indico page is here:

 

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

 

In addition to the usual topics, I’d like to discuss progress on the pCDR.

 

Best,

Anne

 

--------------------------------------------------------------

Anne Sickles

Assistant Professor, Department of Physics

University of Illinois at Urbana-Champaign

--------------------------------------------------------------