LISTSERV mailing list manager LISTSERV 16.5

Help for HPS-SVN Archives


HPS-SVN Archives

HPS-SVN Archives


HPS-SVN@LISTSERV.SLAC.STANFORD.EDU


View:

Message:

[

First

|

Previous

|

Next

|

Last

]

By Topic:

[

First

|

Previous

|

Next

|

Last

]

By Author:

[

First

|

Previous

|

Next

|

Last

]

Font:

Proportional Font

LISTSERV Archives

LISTSERV Archives

HPS-SVN Home

HPS-SVN Home

HPS-SVN  February 2015

HPS-SVN February 2015

Subject:

r2021 - /java/trunk/users/src/main/java/org/hps/users/baltzell/ECalRunningPedestalDriver.java

From:

[log in to unmask]

Reply-To:

Notification of commits to the hps svn repository <[log in to unmask]>

Date:

Sun, 1 Feb 2015 14:31:40 -0000

Content-Type:

text/plain

Parts/Attachments:

Parts/Attachments

text/plain (143 lines)

Author: [log in to unmask]
Date: Sun Feb  1 06:31:34 2015
New Revision: 2021

Log:
Changing initialization

Modified:
    java/trunk/users/src/main/java/org/hps/users/baltzell/ECalRunningPedestalDriver.java

Modified: java/trunk/users/src/main/java/org/hps/users/baltzell/ECalRunningPedestalDriver.java
 =============================================================================
--- java/trunk/users/src/main/java/org/hps/users/baltzell/ECalRunningPedestalDriver.java	(original)
+++ java/trunk/users/src/main/java/org/hps/users/baltzell/ECalRunningPedestalDriver.java	Sun Feb  1 06:31:34 2015
@@ -53,6 +53,7 @@
     private List<Integer>[] eventPedestals = (ArrayList<Integer>[]) new ArrayList[nChannels];
     private List<Long>[] eventTimestamps = (ArrayList<Long>[]) new ArrayList[nChannels];
 
+    private boolean debug = false;
     private EcalConditions ecalConditions = null;
 
     
@@ -61,7 +62,7 @@
         for (int ii = 0; ii < nChannels; ii++) {
             eventPedestals[ii] = new ArrayList<>();
             eventTimestamps[ii] = new ArrayList<>();
-            runningPedestals.add(-1.); // would like to initialize with DB pedestals, but they're not available yet
+            runningPedestals.add(-1.);
         }
     }
     @Override
@@ -73,27 +74,38 @@
         ecalConditions = ConditionsManager.defaultInstance()
                 .getCachedConditions(EcalConditions.class,TableConstants.ECAL_CONDITIONS)
                 .getCachedData();
+        for (int ii=0; ii<nChannels; ii++) {
+            runningPedestals.set(ii,getStaticPedestal(ii+1));
+        }
+        if (debug){
+            printPedestals();
+        }
     }
     public void setMinLookbackEvents(int nev) {
+        if (nev<1){
+            System.err.println(
+                    "ECalRunningPedestalDriver:  MinLookbackEvents too small.  Setting to 1.");
+            nev=1;
+        }
         minLookbackEvents = nev;
     }
     public void setMaxLookbackEvents(int nev) {
+        if (nev>limitLookbackEvents){
+            System.err.println(
+                    "ECalRunningPedestalDriver:  MaxLookbackEvents too big.  Setting to "
+                     +limitLookbackEvents+".");
+            nev=limitLookbackEvents;
+        }
         maxLookbackEvents = nev;
     }
     public void setMaxLookbackTime(int time) {
         maxLookbackTime = time;
     }
 
-
-    public double getPedestal(int channel_id) {
-        final int nped = eventPedestals[channel_id - 1].size();
-        if (nped < minLookbackEvents) return getStaticPedestal(channel_id);
-        else                          return runningPedestals.get(channel_id - 1);
-    }
-
     public void printPedestals() {
         for (int ii = 0; ii < nChannels; ii++) {
-            System.out.printf("(%d,%.2f,%.2f) ",ii,runningPedestals.get(ii),getStaticPedestal(ii + 1));
+            System.out.printf("(%d,%.2f,%.2f) ",ii,runningPedestals.get(ii),
+                    getStaticPedestal(ii + 1));
         }
         System.out.printf("\n");
     }
@@ -109,25 +121,19 @@
             GenericObject extraData = (GenericObject) rel.getTo();
             updatePedestal(event,hit,extraData);
         }
-
-        // quick fix until I know how to read from DB before 'process':
-        List<Double> peds = new ArrayList<Double>(nChannels);
-        for (int ii=0; ii<nChannels; ii++){
-            peds.add(ii,getPedestal(ii+1));
-        }
         
         //
-        // don't care right now whether this persists in output slcio,
+        // Don't care right now whether this persists in output slcio,
         // just that it is accessible during reconstruction (and it is)
         //
         // Another option would be to put hits' running pedestals into HitExtraData.Mode7Data
-        // Or create another LCRelation
-        // Either would also remove the need for indexing later.
+        // Or create another LCRelation.  Either would also remove the need for indexing later.
         //
-        event.put(runningPedestalsName,peds,Double.class,1,"dog");
-//        event.put(runningPedestalsName,runningPedestals,Double.class,1,"dog");
-
-//        printPedestals();
+        event.put(runningPedestalsName,runningPedestals,Double.class,1,"dog");
+
+        if (debug) {
+            printPedestals();
+        }
     }
 
     public void updatePedestal(EventHeader event, RawCalorimeterHit hit, GenericObject mode7data) {
@@ -160,7 +166,7 @@
 
         if (eventPedestals[ii].size() > 1) {
 
-            // remove oldest pedestal if we surpassed limit on #events:
+            // remove oldest pedestal if surpassed limit on #events:
             if (eventPedestals[ii].size() > limitLookbackEvents
                     || (maxLookbackEvents > 0 && eventPedestals[ii].size() > maxLookbackEvents)) {
                 eventPedestals[ii].remove(0);
@@ -181,7 +187,7 @@
         }
 
         // Update running pedestal average:
-        if (eventPedestals[ii].size() > 0) {
+        if (eventPedestals[ii].size() >= minLookbackEvents) {
             double avg = 0;
             for (int jj = 0; jj < eventPedestals[ii].size(); jj++) {
                 avg += eventPedestals[ii].get(jj);
@@ -192,6 +198,12 @@
         }
     }
 
+    public double getPedestal(int channel_id) {
+        return runningPedestals.get(channel_id-1);
+//        final int nped = eventPedestals[channel_id - 1].size();
+//        if (nped < minLookbackEvents) return getStaticPedestal(channel_id);
+//        else                          return runningPedestals.get(channel_id - 1);
+    }
     public double getPedestal(RawCalorimeterHit hit) {
         return getPedestal(getChannelID(hit));
     }

Top of Message | Previous Page | Permalink

Advanced Options


Options

Log In

Log In

Get Password

Get Password


Search Archives

Search Archives


Subscribe or Unsubscribe

Subscribe or Unsubscribe


Archives

November 2017
August 2017
July 2017
January 2017
December 2016
November 2016
October 2016
September 2016
August 2016
July 2016
June 2016
May 2016
April 2016
March 2016
February 2016
January 2016
December 2015
November 2015
October 2015
September 2015
August 2015
July 2015
June 2015
May 2015
April 2015
March 2015
February 2015
January 2015
December 2014
November 2014
October 2014
September 2014
August 2014
July 2014
June 2014
May 2014
April 2014
March 2014
February 2014
January 2014
December 2013
November 2013

ATOM RSS1 RSS2



LISTSERV.SLAC.STANFORD.EDU

Secured by F-Secure Anti-Virus CataList Email List Search Powered by the LISTSERV Email List Manager

Privacy Notice, Security Notice and Terms of Use