LISTSERV 16.5 - HPS-SOFTWARE Archives

Hi, Luca.

I think I know the problem....

The first ".evio.0" file has a pre start event in it so the run number is read from this and written correctly. All subsequent files in the run do not have the pre start event, and so the run number was (incorrectly) set to zero. These files should have been processed in a way that avoided this problem by using "-R 3261" to force a run number into the job externally. This is a problem with all pass0 LCIO files probably that are not the first file in the run.

But I think you can get around this problem by using ConditionsDriver, which I have covered in previous emails.

Basically, you put this Driver in your steering list as the first one....

<detectorName>HPS-ECalCommissioning</detectorName>

Then when running your job you should see that the conditions system is forced to use this run number for the whole job.

Let me know if that works....

--Jeremy

From: Luca Colaneri <[log in to unmask]>
Sent: Wednesday, January 14, 2015 2:48 AM
To: McCormick, Jeremy I.
Subject: Re: channel not found

Hi Jeremy, I don't know if you recived my last mail, so here it is again.

Sorry to push but it is kind of important.

In the attachment you have the log of the output and the steering I used.
The problem is that when passing from file 0 to file 1 it reads file 1 as "run#0" (should be 3261) and crashes.

I need to be able to read run files all together and obtain from each channel the gain values.

Luca

Il 12/01/2015 18:19, McCormick, Jeremy I. ha scritto:

Can you send me the full run log as an attachment?
On Jan 12, 2015, at 9:15 AM, Luca Colaneri <[log in to unmask]> wrote:
I'm not running any converter, I'm just reading the recon data available on tape, run 3261 to be precise..

Il 12/01/15 18:14, McCormick, Jeremy I. ha scritto:

Hi, Luca.

If you are running EvioToLcio then you may need to specify the run number manually if there is not an EVIO header record first.

java -cp ~/.m2/repository/org/hps/hps-distribution/3.1-SNAPSHOT/hps-distribution-3.1-SNAPSHOT-bin.jar org.hps.evio.EvioToLcio -R 3261 [args]

This should not be necessary if you are running off of the LCIO files as the run number is extracted from each event in that format, and it is written there by the EVIO converter.

—Jeremy

On Jan 12, 2015, at 6:44 AM, Luca Colaneri <[log in to unmask]> wrote:

Hello all,

I may have found the issue. I was running the code over the full batch of files of run 3261. When starting, the condition manager recognize the run and everything work, but when passing to the second file, it doesn't and says "run#0". I run the code separetly on the three files and this happens for the second and the third, the output message is in the attachment.

Am I doing something wrong or the files just don't have the run info and conditions?

Luca

Il 12/01/2015 12:32, Luca Colaneri ha scritto:

Il 10/01/2015 18:57, Nathan Baltzell ha scritto:

Hi Luca and All,

I have a few comments:

Why would it say it is using run #0?

I have no idea. Jeremy can you comment that?

Can you print out the EcalChannel you are passing to getChannelConstants before this error?
(Looks like there are EcalChannel.getChannelId/Slot/Crate/X/Y methods.)

This is the piece of code:

EcalConditions ecalConditions = ConditionsManager.defaultInstance().getCachedConditions(EcalConditions.class, TableConstants.ECAL_CONDITIONS).getCachedData();
EcalChannelCollection channels = ecalConditions.getChannelCollection();

then i collect the cluster collection and in a for loop over the clusters:

EcalChannel channel = channels.findGeometric(cluster.getCalorimeterHits().get(0).getCellID());
EcalChannelConstants channelConstants = ecalConditions.getChannelConstants(channel);
channelConstants.getGain().getGain()

Do you need to divide by the gain to get the ADC? I thought that raw energy was directly available?
Alternatively, you could just multiply your new scale factors by the cosmic gains in the end.

Raw energy is available, but then I would have to rerun the clustering and i think it's easier to just divide out the cosmic gains.

Luca

-Nathan

On Jan 10, 2015, at 12:17 PM, Luca Colaneri <[log in to unmask]> wrote:

Hello all,
In order to make my calibration "cosmics" independent I was trying to get out from the data the gain values to divide out from the corrected energies and get the raw data.
After around 31M events I got the error you can read in the attached screenshot.

Suggestions?

L.

--
PhD Student,
Università degli Studi di Roma Tor Vergata
INFN sez. Roma2
Via della Ricerca Scientifica 1, 00133 Rome, Italy
+39 0672594562

########################################################################
Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1
<channelnotfound.png>

########################################################################
Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1

########################################################################
Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1

Use REPLY-ALL to reply to list

To unsubscribe from the HPS-SOFTWARE list, click the following link:
https://listserv.slac.stanford.edu/cgi-bin/wa?SUBED1=HPS-SOFTWARE&A=1
-- 
Luca Colaneri
Dipartimento di Fisica 
Universita' degli Studi di Roma Tor Vergata
Via della Ricerca Scientifica s.n.c. 00133 Roma
+39 06 72594560

-- 
PhD Student,
Università degli Studi di Roma Tor Vergata 
INFN sez. Roma2
Via della Ricerca Scientifica 1, 00133 Rome, Italy
+39 0672594562