#include <HBHETimeProfileStatusBitSetter.h>

Public Member Functions
double	EnergyThreshold ()

void	hbheSetTimeFlagsFromDigi (HBHERecHitCollection *, const std::vector< HBHEDataFrame > &, const std::vector< int > &)

	HBHETimeProfileStatusBitSetter ()

	HBHETimeProfileStatusBitSetter (double R1Min, double R1Max, double R2Min, double R2Max, double FracLeaderMin, double FracLeaderMax, double SlopeMin, double SlopeMax, double OuterMin, double OuterMax, double EnergyThreshold)

void	SetExpLimits (double R1Min, double R1Max, double R2Min, double R2Max)

void	SetFracLeaderLimits (double FracLeaderMin, double FracLeaderMax)

void	SetOuterLimits (double OuterMin, double OuterMax)

void	SetSlopeLimits (double SlopeMin, double SlopeMax)

	~HBHETimeProfileStatusBitSetter ()

Private Member Functions
double	TotalEnergyInDataFrame (const HBHEDataFrame &x)

Private Attributes
double	EnergyThreshold_

double	FracLeaderMax_

double	FracLeaderMin_

double	OuterMax_

double	OuterMin_

double	R1Max_

double	R1Min_

double	R2Max_

double	R2Min_

double	SlopeMax_

double	SlopeMin_

Detailed Description

This class sets status bit in the status words for the revised CaloRecHit objets according to informatino from the digi associated to the hit.

Author: B Jones – University of Bristol / University of Maryland

Definition at line 15 of file HBHETimeProfileStatusBitSetter.h.

Constructor & Destructor Documentation

◆ HBHETimeProfileStatusBitSetter() [1/2]

HBHETimeProfileStatusBitSetter::HBHETimeProfileStatusBitSetter ( )

Full featured constructor for HB/HE and HO (HPD-based detectors)

Definition at line 10 of file HBHETimeProfileStatusBitSetter.cc.

                                                                {
   // use simple values in default constructor
   R1Min_ = 0.1;
   R1Max_ = 0.7;
   R2Min_ = 0.2;
   R2Max_ = 0.5;
   FracLeaderMin_ = 0.4;
   FracLeaderMax_ = 0.7;
   SlopeMin_ = -1.5;
   SlopeMax_ = -0.6;
   OuterMin_ = 0.9;
   OuterMax_ = 1.0;
   EnergyThreshold_ = 30;
 }

References EnergyThreshold_, FracLeaderMax_, FracLeaderMin_, OuterMax_, OuterMin_, R1Max_, R1Min_, R2Max_, R2Min_, SlopeMax_, and SlopeMin_.

◆ HBHETimeProfileStatusBitSetter() [2/2]

HBHETimeProfileStatusBitSetter::HBHETimeProfileStatusBitSetter	(	double	R1Min,
		double	R1Max,
		double	R2Min,
		double	R2Max,
		double	FracLeaderMin,
		double	FracLeaderMax,
		double	SlopeMin,
		double	SlopeMax,
		double	OuterMin,
		double	OuterMax,
		double	EnergyThreshold
	)

Definition at line 25 of file HBHETimeProfileStatusBitSetter.cc.

                                                                                        {
   R1Min_ = R1Min;
   R1Max_ = R1Max;
   R2Min_ = R2Min;
   R2Max_ = R2Max;
   FracLeaderMin_ = FracLeaderMin;
   FracLeaderMax_ = FracLeaderMax;
   SlopeMin_ = SlopeMin;
   SlopeMax_ = SlopeMax;
   OuterMin_ = OuterMin;
   OuterMax_ = OuterMax;
   EnergyThreshold_ = EnergyThreshold;
 }

References EnergyThreshold(), EnergyThreshold_, FracLeaderMax_, FracLeaderMin_, OuterMax_, OuterMin_, R1Max_, R1Min_, R2Max_, R2Min_, SlopeMax_, and SlopeMin_.

◆ ~HBHETimeProfileStatusBitSetter()

HBHETimeProfileStatusBitSetter::~HBHETimeProfileStatusBitSetter ( )

Definition at line 49 of file HBHETimeProfileStatusBitSetter.cc.

49 {}

Member Function Documentation

◆ EnergyThreshold()

double HBHETimeProfileStatusBitSetter::EnergyThreshold ( )

inline

Definition at line 56 of file HBHETimeProfileStatusBitSetter.h.

56 { return EnergyThreshold_; }

References EnergyThreshold_.

Referenced by HBHETimeProfileStatusBitSetter().

◆ hbheSetTimeFlagsFromDigi()

void HBHETimeProfileStatusBitSetter::hbheSetTimeFlagsFromDigi	(	HBHERecHitCollection *	hbhe,
		const std::vector< HBHEDataFrame > &	udigi,
		const std::vector< int > &	RecHitIndex
	)

Definition at line 64 of file HBHETimeProfileStatusBitSetter.cc.

                                                                                                  {
   bool Bits[4] = {false, false, false, false};
   std::vector<HBHEDataFrame> digi(udigi);
   std::sort(digi.begin(), digi.end(), compareDigiEnergy);  // sort digis according to energies
   std::vector<double> PulseShape;                          // store fC values for each time slice
   int DigiSize = 0;
   //  int LeadingEta=0;
   int LeadingPhi = 0;
   bool FoundLeadingChannel = false;
   for (std::vector<HBHEDataFrame>::const_iterator itDigi = digi.begin(); itDigi != digi.end(); itDigi++) {
     if (!FoundLeadingChannel) {
       //      LeadingEta = itDigi->id().ieta();
       LeadingPhi = itDigi->id().iphi();
       DigiSize = (*itDigi).size();
       PulseShape.clear();
       PulseShape.resize(DigiSize, 0);
       FoundLeadingChannel = true;
     }
     if (abs(LeadingPhi - itDigi->id().iphi()) < 2)
       for (int i = 0; i != DigiSize; i++)
         PulseShape[i] += itDigi->sample(i).nominal_fC();
   }
  
   if (!RecHitIndex.empty()) {
     double FracInLeader = -1;
     //double Slope=0; // not currently used
     double R1 = -1;
     double R2 = -1;
     double OuterEnergy = -1;
     double TotalEnergy = 0;
     int PeakPosition = 0;
  
     for (int i = 0; i != DigiSize; i++) {
       if (PulseShape[i] > PulseShape[PeakPosition])
         PeakPosition = i;
       TotalEnergy += PulseShape[i];
     }
  
     if (PeakPosition < (DigiSize - 2)) {
       R1 = PulseShape[PeakPosition + 1] / PulseShape[PeakPosition];
       R2 = PulseShape[PeakPosition + 2] / PulseShape[PeakPosition + 1];
     }
  
     FracInLeader = PulseShape[PeakPosition] / TotalEnergy;
  
     if ((PeakPosition > 0) && (PeakPosition < (DigiSize - 2))) {
       OuterEnergy = 1. - ((PulseShape[PeakPosition - 1] + PulseShape[PeakPosition] + PulseShape[PeakPosition + 1] +
                            PulseShape[PeakPosition + 2]) /
                           TotalEnergy);
     }
  
     /*      TH1D * HistForFit = new TH1D("HistForFit","HistForFit",DigiSize,0,DigiSize);
       for(int i=0; i!=DigiSize; i++)
     {
       HistForFit->Fill(i,PulseShape[i]);
       HistForFit->Fit("expo","WWQ","",PeakPosition, DigiSize-1);
       TF1 * Fit = HistForFit->GetFunction("expo");
       Slope = Fit->GetParameter("Slope");
     }
       delete HistForFit;
       */
     if (R1 != -1 && R2 != -1)
       Bits[0] = (R1Min_ > R1) || (R1Max_ < R1) || (R2Min_ > R2) || (R2Max_ < R2);
     if (FracInLeader != -1)
       Bits[1] = (FracInLeader < FracLeaderMin_) || (FracInLeader > FracLeaderMax_);
     if (OuterEnergy != -1)
       Bits[2] = (OuterEnergy < OuterMin_) || (OuterEnergy > OuterMax_);
     //  Bits[3] = (SlopeMin_ > Slope) || (SlopeMax_ < Slope);
     Bits[3] = false;
  
   }  // if (RecHitIndex.size()>0)
   else {
     Bits[0] = false;
     Bits[1] = false;
     Bits[2] = false;
     Bits[3] = true;
  
   }  // (RecHitIndex.size()==0; no need to set Bit3 true?)
  
   for (unsigned int i = 0; i != RecHitIndex.size(); i++) {
     // Write calculated bit values starting from position FirstBit
     (*hbhe)[RecHitIndex.at(i)].setFlagField(Bits[0], HcalCaloFlagLabels::HSCP_R1R2);
     (*hbhe)[RecHitIndex.at(i)].setFlagField(Bits[1], HcalCaloFlagLabels::HSCP_FracLeader);
     (*hbhe)[RecHitIndex.at(i)].setFlagField(Bits[2], HcalCaloFlagLabels::HSCP_OuterEnergy);
     (*hbhe)[RecHitIndex.at(i)].setFlagField(Bits[3], HcalCaloFlagLabels::HSCP_ExpFit);
   }
 }

References funct::abs(), FracLeaderMax_, FracLeaderMin_, HcalCaloFlagLabels::HSCP_ExpFit, HcalCaloFlagLabels::HSCP_FracLeader, HcalCaloFlagLabels::HSCP_OuterEnergy, HcalCaloFlagLabels::HSCP_R1R2, mps_fire::i, OuterMax_, OuterMin_, R1Max_, R1Min_, R2Max_, R2Min_, and jetUpdater_cfi::sort.

◆ SetExpLimits()

void HBHETimeProfileStatusBitSetter::SetExpLimits	(	double	R1Min,
		double	R1Max,
		double	R2Min,
		double	R2Max
	)

inline

Definition at line 38 of file HBHETimeProfileStatusBitSetter.h.

                                                                             {
     R1Min_ = R1Min;
     R1Max_ = R1Max;
     R2Min_ = R2Max;
     R2Max_ = R2Max;
   }

References R1Max_, R1Min_, R2Max_, and R2Min_.

◆ SetFracLeaderLimits()

void HBHETimeProfileStatusBitSetter::SetFracLeaderLimits	(	double	FracLeaderMin,
		double	FracLeaderMax
	)

inline

Definition at line 44 of file HBHETimeProfileStatusBitSetter.h.

                                                                        {
     FracLeaderMin_ = FracLeaderMin;
     FracLeaderMax_ = FracLeaderMax;
   }

References FracLeaderMax_, and FracLeaderMin_.

◆ SetOuterLimits()

void HBHETimeProfileStatusBitSetter::SetOuterLimits	(	double	OuterMin,
		double	OuterMax
	)

inline

Definition at line 52 of file HBHETimeProfileStatusBitSetter.h.

                                                         {
     OuterMin_ = OuterMin;
     OuterMax_ = OuterMax;
   }

References OuterMax_, and OuterMin_.

◆ SetSlopeLimits()

void HBHETimeProfileStatusBitSetter::SetSlopeLimits	(	double	SlopeMin,
		double	SlopeMax
	)

inline

Definition at line 48 of file HBHETimeProfileStatusBitSetter.h.

                                                         {
     SlopeMin_ = SlopeMin;
     SlopeMax_ = SlopeMax;
   }

References SlopeMax_, and SlopeMin_.

◆ TotalEnergyInDataFrame()

double HBHETimeProfileStatusBitSetter::TotalEnergyInDataFrame ( const HBHEDataFrame & x )

inlineprivate

Definition at line 66 of file HBHETimeProfileStatusBitSetter.h.

                                                         {
     double Total = 0;
     for (int i = 0; i != x.size(); Total += x.sample(i++).nominal_fC())
       ;
     return Total;
   }