Changeset 129 in Idarraga for mafalda

Timestamp:

Apr 14, 2011, 6:22:07 PM (13 years ago)

Author:

idarraga

Message:

 

File:

• mafalda/FEI3/FEI3Mips/FEI3Mips.cpp (modified) (10 diffs)

mafalda/FEI3/FEI3Mips/FEI3Mips.cpp

@@ -12 +12 @@
 #include "MPXAlgo/Highlighter.h"
 
-#include "TColor.h"
+#include "TH1F.h"
 
 using namespace MSG;
@@ -27 +27 @@
 {
 
+
 }
 
@@ -38 +39 @@
         getMyTree()->Branch("chargeWeights", &m_chargeWeights);
         getMyTree()->Branch("lvl1Weights", &m_lvl1Weights);
+        getMyTree()->Branch("chargePerSection", &m_chargePerSection);
+
         getMyTree()->Branch("clusterSize", &m_clusterSize, "clusterSize/I");
         getMyTree()->Branch("alphaAngle", &m_alphaAngle, "alphaAngle/D");
@@ -89 +92 @@
                                 m_chargeWeights.push_back(0.0);
                                 m_lvl1Weights.push_back(0.0);
+                                m_chargePerSection.push_back(vector<Int_t>()); // push back empty vectors
                         }
 
@@ -225 +229 @@
                                         m_lvl1Weights[__HEAD_INDX] += GetLVL1(*contItr);
 
+                                        // store charge for landau distribution per section
+                                        m_chargePerSection[__HEAD_INDX].push_back( GetMatrixElement(*contItr) );
+
                                 }else if( firstThreeBits == __RIGHT_CORNER_MASK_1
                                                 || firstThreeBits == __RIGHT_CORNER_MASK_2){
@@ -231 +238 @@
                                         m_chargeWeights[__TAIL_INDX] += GetMatrixElement(*contItr);
                                         m_lvl1Weights[__TAIL_INDX] += GetLVL1(*contItr);
+
+                                        // store charge for landau distribution per section
+                                        m_chargePerSection[__TAIL_INDX].push_back( GetMatrixElement(*contItr) );
 
                                 }else{
@@ -241 +251 @@
 
                                         Int_t nB = NumberOfBitsOnBeforeDoubleZero(locationBitWord);
-                                        Log << MSG::DEBUG << "... column " << nB << " ... " << endreq;
-
+
+                                        // Handle special case when the pixel is equidistant to two lines.
+                                        // We end up with a word  001111  for m_nDivisions = 4.  This pixel belongs
+                                        //  to the 3rd slice.
+                                        if(nB == m_nDivisions){
+                                                nB--;
+                                        }
+                                        Log << MSG::DEBUG << "... column " << nB << " ... " << GetFrameId() << endreq;
                                         m_chargeWeights[nB] += GetMatrixElement(*contItr);
                                         m_lvl1Weights[nB] += GetLVL1(*contItr);
+                                        // Store charge for landau distribution per section
+                                        m_chargePerSection[nB].push_back( GetMatrixElement(*contItr) );
+
                                 }
 
@@ -260 +279 @@
                         }
 
+
                         // WARNING, filling once per mip !!!
                         m_frameId = GetFrameId();
@@ -266 +286 @@
                         m_chargeWeights.clear();
                         m_lvl1Weights.clear();
+                        m_chargePerSection.clear();
                         m_alphaAngle = 0.;
                         m_clusterSize = 0;
@@ -281 +302 @@
         m_chargeWeights.clear();
         m_lvl1Weights.clear();
+        m_chargePerSection.clear();
 
 }