|
|
|
#include <HcalHitReconstructor.h>
Public Member Functions | |
| HcalHitReconstructor (const edm::ParameterSet &ps) | |
| virtual void | produce (edm::Event &e, const edm::EventSetup &c) |
| virtual | ~HcalHitReconstructor () |
Private Attributes | |
| bool | correctTiming_ |
| DetId::Detector | det_ |
| bool | dropZSmarkedPassed_ |
| int | firstauxTS_ |
| HBHEStatusBitSetter * | hbheFlagSetter_ |
| HBHETimeProfileStatusBitSetter * | hbheHSCPFlagSetter_ |
| HBHEPulseShapeFlagSetter * | hbhePulseShapeFlagSetter_ |
| HBHETimingShapedFlagSetter * | hbheTimingShapedFlagSetter_ |
| HcalHFStatusBitFromDigis * | hfdigibit_ |
| HcalHF_PETalgorithm * | hfPET_ |
| HcalHF_S9S1algorithm * | hfS8S1_ |
| HcalHF_S9S1algorithm * | hfS9S1_ |
| HFTimingTrustFlag * | HFTimingTrustFlagSetter_ |
| edm::InputTag | inputLabel_ |
| HcalSimpleRecAlgo | reco_ |
| HcalADCSaturationFlag * | saturationFlagSetter_ |
| bool | setHSCPFlags_ |
| bool | setNoiseFlags_ |
| bool | setPulseShapeFlags_ |
| bool | setSaturationFlags_ |
| bool | setTimingTrustFlags_ |
| int | subdet_ |
| HcalOtherSubdetector | subdetOther_ |
| HcalHitReconstructor::HcalHitReconstructor | ( | const edm::ParameterSet & | ps | ) | [explicit] |
Definition at line 21 of file HcalHitReconstructor.cc.
References DetId::Calo, gather_cfg::cout, det_, edm::ParameterSet::getParameter(), hbheFlagSetter_, hbheHSCPFlagSetter_, hbhePulseShapeFlagSetter_, hbheTimingShapedFlagSetter_, HcalBarrel, HcalCalibration, HcalForward, HcalOther, HcalOuter, hfdigibit_, hfPET_, hfS8S1_, hfS9S1_, HFTimingTrustFlagSetter_, saturationFlagSetter_, setHSCPFlags_, setNoiseFlags_, setPulseShapeFlags_, setSaturationFlags_, setTimingTrustFlags_, subdet_, HcalZDCDetId::SubdetectorId, and subdetOther_.
: reco_(conf.getParameter<int>("firstSample"), conf.getParameter<int>("samplesToAdd"), conf.getParameter<bool>("correctForTimeslew"), conf.getParameter<bool>("correctForPhaseContainment"), conf.getParameter<double>("correctionPhaseNS")), det_(DetId::Hcal), inputLabel_(conf.getParameter<edm::InputTag>("digiLabel")), correctTiming_(conf.getParameter<bool>("correctTiming")), setNoiseFlags_(conf.getParameter<bool>("setNoiseFlags")), setHSCPFlags_(conf.getParameter<bool>("setHSCPFlags")), setSaturationFlags_(conf.getParameter<bool>("setSaturationFlags")), setTimingTrustFlags_(conf.getParameter<bool>("setTimingTrustFlags")), setPulseShapeFlags_(conf.getParameter<bool>("setPulseShapeFlags")), dropZSmarkedPassed_(conf.getParameter<bool>("dropZSmarkedPassed")), firstauxTS_(conf.getParameter<int>("firstSample")+conf.getParameter<int>("firstAuxOffset")) { std::string subd=conf.getParameter<std::string>("Subdetector"); //Set all FlagSetters to 0 /* Important to do this! Otherwise, if the setters are turned off, the "if (XSetter_) delete XSetter_;" commands can crash */ hbheFlagSetter_ = 0; hbheHSCPFlagSetter_ = 0; hbhePulseShapeFlagSetter_ = 0; hbheTimingShapedFlagSetter_ = 0; hfdigibit_ = 0; hfS9S1_ = 0; hfS8S1_ = 0; hfPET_ = 0; saturationFlagSetter_ = 0; HFTimingTrustFlagSetter_ = 0; if (setSaturationFlags_) { const edm::ParameterSet& pssat = conf.getParameter<edm::ParameterSet>("saturationParameters"); saturationFlagSetter_ = new HcalADCSaturationFlag(pssat.getParameter<int>("maxADCvalue")); } if (!strcasecmp(subd.c_str(),"HBHE")) { subdet_=HcalBarrel; bool timingShapedCutsFlags = conf.getParameter<bool>("setTimingShapedCutsFlags"); if (timingShapedCutsFlags) { const edm::ParameterSet& psTshaped = conf.getParameter<edm::ParameterSet>("timingshapedcutsParameters"); hbheTimingShapedFlagSetter_ = new HBHETimingShapedFlagSetter(psTshaped.getParameter<std::vector<double> >("tfilterEnvelope"), psTshaped.getParameter<bool>("ignorelowest"), psTshaped.getParameter<bool>("ignorehighest"), psTshaped.getParameter<double>("win_offset"), psTshaped.getParameter<double>("win_gain")); } if (setNoiseFlags_) { const edm::ParameterSet& psdigi =conf.getParameter<edm::ParameterSet>("flagParameters"); hbheFlagSetter_=new HBHEStatusBitSetter(psdigi.getParameter<double>("nominalPedestal"), psdigi.getParameter<double>("hitEnergyMinimum"), psdigi.getParameter<int>("hitMultiplicityThreshold"), psdigi.getParameter<std::vector<edm::ParameterSet> >("pulseShapeParameterSets"), conf.getParameter<int>("firstSample"), conf.getParameter<int>("samplesToAdd")); } // if (setNoiseFlags_) if (setHSCPFlags_) { const edm::ParameterSet& psHSCP = conf.getParameter<edm::ParameterSet>("hscpParameters"); hbheHSCPFlagSetter_ = new HBHETimeProfileStatusBitSetter(psHSCP.getParameter<double>("r1Min"), psHSCP.getParameter<double>("r1Max"), psHSCP.getParameter<double>("r2Min"), psHSCP.getParameter<double>("r2Max"), psHSCP.getParameter<double>("fracLeaderMin"), psHSCP.getParameter<double>("fracLeaderMax"), psHSCP.getParameter<double>("slopeMin"), psHSCP.getParameter<double>("slopeMax"), psHSCP.getParameter<double>("outerMin"), psHSCP.getParameter<double>("outerMax"), psHSCP.getParameter<double>("TimingEnergyThreshold")); } // if (setHSCPFlags_) if (setPulseShapeFlags_) { const edm::ParameterSet &psPulseShape = conf.getParameter<edm::ParameterSet>("pulseShapeParameters"); hbhePulseShapeFlagSetter_ = new HBHEPulseShapeFlagSetter( psPulseShape.getParameter<double>("MinimumChargeThreshold"), psPulseShape.getParameter<unsigned int>("TrianglePeakTS"), psPulseShape.getParameter<std::vector<double> >("LinearThreshold"), psPulseShape.getParameter<std::vector<double> >("LinearCut"), psPulseShape.getParameter<std::vector<double> >("RMS8MaxThreshold"), psPulseShape.getParameter<std::vector<double> >("RMS8MaxCut"), psPulseShape.getParameter<std::vector<double> >("LeftSlopeThreshold"), psPulseShape.getParameter<std::vector<double> >("LeftSlopeCut"), psPulseShape.getParameter<std::vector<double> >("RightSlopeThreshold"), psPulseShape.getParameter<std::vector<double> >("RightSlopeCut"), psPulseShape.getParameter<std::vector<double> >("RightSlopeSmallThreshold"), psPulseShape.getParameter<std::vector<double> >("RightSlopeSmallCut"), psPulseShape.getParameter<bool>("UseDualFit"), psPulseShape.getParameter<bool>("TriangleIgnoreSlow")); } // if (setPulseShapeFlags_) produces<HBHERecHitCollection>(); } else if (!strcasecmp(subd.c_str(),"HO")) { subdet_=HcalOuter; produces<HORecHitCollection>(); } else if (!strcasecmp(subd.c_str(),"HF")) { subdet_=HcalForward; if (setTimingTrustFlags_) { const edm::ParameterSet& pstrust = conf.getParameter<edm::ParameterSet>("hfTimingTrustParameters"); HFTimingTrustFlagSetter_=new HFTimingTrustFlag(pstrust.getParameter<int>("hfTimingTrustLevel1"), pstrust.getParameter<int>("hfTimingTrustLevel2")); } if (setNoiseFlags_) { const edm::ParameterSet& psdigi =conf.getParameter<edm::ParameterSet>("digistat"); const edm::ParameterSet& psTimeWin =conf.getParameter<edm::ParameterSet>("HFInWindowStat"); hfdigibit_=new HcalHFStatusBitFromDigis(conf.getParameter<int>("firstSample"), conf.getParameter<int>("samplesToAdd"), psdigi, psTimeWin); const edm::ParameterSet& psS9S1 = conf.getParameter<edm::ParameterSet>("S9S1stat"); hfS9S1_ = new HcalHF_S9S1algorithm(psS9S1.getParameter<std::vector<double> >("short_optimumSlope"), psS9S1.getParameter<std::vector<double> >("shortEnergyParams"), psS9S1.getParameter<std::vector<double> >("shortETParams"), psS9S1.getParameter<std::vector<double> >("long_optimumSlope"), psS9S1.getParameter<std::vector<double> >("longEnergyParams"), psS9S1.getParameter<std::vector<double> >("longETParams"), psS9S1.getParameter<int>("flagsToSkip"), psS9S1.getParameter<bool>("isS8S1") ); const edm::ParameterSet& psS8S1 = conf.getParameter<edm::ParameterSet>("S8S1stat"); hfS8S1_ = new HcalHF_S9S1algorithm(psS8S1.getParameter<std::vector<double> >("short_optimumSlope"), psS8S1.getParameter<std::vector<double> >("shortEnergyParams"), psS8S1.getParameter<std::vector<double> >("shortETParams"), psS8S1.getParameter<std::vector<double> >("long_optimumSlope"), psS8S1.getParameter<std::vector<double> >("longEnergyParams"), psS8S1.getParameter<std::vector<double> >("longETParams"), psS8S1.getParameter<int>("flagsToSkip"), psS8S1.getParameter<bool>("isS8S1") ); const edm::ParameterSet& psPET = conf.getParameter<edm::ParameterSet>("PETstat"); hfPET_ = new HcalHF_PETalgorithm(psPET.getParameter<std::vector<double> >("short_R"), psPET.getParameter<std::vector<double> >("shortEnergyParams"), psPET.getParameter<std::vector<double> >("shortETParams"), psPET.getParameter<std::vector<double> >("long_R"), psPET.getParameter<std::vector<double> >("longEnergyParams"), psPET.getParameter<std::vector<double> >("longETParams"), psPET.getParameter<int>("flagsToSkip"), psPET.getParameter<std::vector<double> >("short_R_29"), psPET.getParameter<std::vector<double> >("long_R_29") ); } produces<HFRecHitCollection>(); } else if (!strcasecmp(subd.c_str(),"ZDC")) { det_=DetId::Calo; subdet_=HcalZDCDetId::SubdetectorId; produces<ZDCRecHitCollection>(); } else if (!strcasecmp(subd.c_str(),"CALIB")) { subdet_=HcalOther; subdetOther_=HcalCalibration; produces<HcalCalibRecHitCollection>(); } else { std::cout << "HcalHitReconstructor is not associated with a specific subdetector!" << std::endl; } }
| HcalHitReconstructor::~HcalHitReconstructor | ( | ) | [virtual] |
Definition at line 190 of file HcalHitReconstructor.cc.
References hbheFlagSetter_, hbheHSCPFlagSetter_, hbhePulseShapeFlagSetter_, hfdigibit_, hfPET_, and hfS9S1_.
{
if (hbheFlagSetter_) delete hbheFlagSetter_;
if (hfdigibit_) delete hfdigibit_;
if (hbheHSCPFlagSetter_) delete hbheHSCPFlagSetter_;
if (hbhePulseShapeFlagSetter_) delete hbhePulseShapeFlagSetter_;
if (hfS9S1_) delete hfS9S1_;
if (hfPET_) delete hfPET_;
}
| void HcalHitReconstructor::produce | ( | edm::Event & | e, |
| const edm::EventSetup & | c | ||
| ) | [virtual] |
Implements edm::EDProducer.
Definition at line 199 of file HcalHitReconstructor.cc.
References abs, HBHEStatusBitSetter::Clear(), HcalTimingCorrector::Correct(), correctTiming_, det_, HcalSeverityLevelComputer::dropChannel(), dropZSmarkedPassed_, HBHETimeProfileStatusBitSetter::EnergyThreshold(), firstauxTS_, edm::EventSetup::get(), edm::Event::getByLabel(), HcalChannelStatus::getValue(), HcalCondObjectContainer< Item >::getValues(), hbheFlagSetter_, hbheHSCPFlagSetter_, hbhePulseShapeFlagSetter_, HBHETimeProfileStatusBitSetter::hbheSetTimeFlagsFromDigi(), hbheTimingShapedFlagSetter_, DetId::Hcal, HcalBarrel, HcalCalibration, HcalEndcap, HcalForward, HcalOther, HcalOuter, hfdigibit_, hfPET_, hfS8S1_, hfS9S1_, HcalHFStatusBitFromDigis::hfSetFlagFromDigi(), HcalHF_PETalgorithm::HFSetFlagFromPET(), HcalHF_S9S1algorithm::HFSetFlagFromS9S1(), HFTimingTrustFlagSetter_, i, inputLabel_, edm::EventBase::isRealData(), j, gen::k, L1TEmulatorMonitor_cff::p, edm::ESHandle< T >::product(), edm::Event::put(), DetId::rawId(), reco_, HcalSimpleRecAlgo::reconstruct(), HcalHFStatusBitFromDigis::resetTimeSamples(), HcalSimpleRecAlgo::resetTimeSamples(), edm::Event::run(), saturationFlagSetter_, HBHEStatusBitSetter::SetFlagsFromDigi(), HBHEStatusBitSetter::SetFlagsFromRecHits(), HFTimingTrustFlag::setHFTimingTrustFlag(), setHSCPFlags_, setNoiseFlags_, HBHEPulseShapeFlagSetter::SetPulseShapeFlags(), setPulseShapeFlags_, HcalADCSaturationFlag::setSaturationFlag(), setSaturationFlags_, HBHETimingShapedFlagSetter::SetTimingShapedFlags(), setTimingTrustFlags_, subdet_, and subdetOther_.
{
//bool isData=e.isRealData(); // some flags should only be applied to real data
// get conditions
edm::ESHandle<HcalDbService> conditions;
eventSetup.get<HcalDbRecord>().get(conditions);
const HcalQIEShape* shape = conditions->getHcalShape (); // this one is generic
edm::ESHandle<HcalChannelQuality> p;
eventSetup.get<HcalChannelQualityRcd>().get(p);
HcalChannelQuality* myqual = new HcalChannelQuality(*p.product());
edm::ESHandle<HcalSeverityLevelComputer> mycomputer;
eventSetup.get<HcalSeverityLevelComputerRcd>().get(mycomputer);
const HcalSeverityLevelComputer* mySeverity = mycomputer.product();
if (det_==DetId::Hcal) {
if (subdet_==HcalBarrel || subdet_==HcalEndcap) {
edm::Handle<HBHEDigiCollection> digi;
e.getByLabel(inputLabel_,digi);
// create empty output
std::auto_ptr<HBHERecHitCollection> rec(new HBHERecHitCollection);
rec->reserve(digi->size());
// run the algorithm
if (setNoiseFlags_) hbheFlagSetter_->Clear();
HBHEDigiCollection::const_iterator i;
std::vector<HBHEDataFrame> HBDigis;
std::vector<int> RecHitIndex;
// Vote on majority TS0 CapId
int favorite_capid = 0;
if (correctTiming_) {
long capid_votes[4] = {0,0,0,0};
for (i=digi->begin(); i!=digi->end(); i++) {
capid_votes[(*i)[0].capid()]++;
}
for (int k = 0; k < 4; k++)
if (capid_votes[k] > capid_votes[favorite_capid])
favorite_capid = k;
}
for (i=digi->begin(); i!=digi->end(); i++) {
HcalDetId cell = i->id();
DetId detcell=(DetId)cell;
// check on cells to be ignored and dropped: (rof,20.Feb.09)
const HcalChannelStatus* mydigistatus=myqual->getValues(detcell.rawId());
if (mySeverity->dropChannel(mydigistatus->getValue() ) ) continue;
if (dropZSmarkedPassed_)
if (i->zsMarkAndPass()) continue;
const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
HcalCoderDb coder (*channelCoder, *shape);
rec->push_back(reco_.reconstruct(*i,coder,calibrations));
// Set auxiliary flag
int auxflag=0;
for (int xx=firstauxTS_;xx<firstauxTS_+4 && xx<i->size();++xx)
auxflag+=(i->sample(xx).adc())<<(7*(xx-firstauxTS_)); // store the time slices in the first 28 bits of aux, a set of 4 7-bit adc values
// bits 28 and 29 are reserved for capid of the first time slice saved in aux
auxflag+=((i->sample(firstauxTS_).capid())<<28);
(rec->back()).setAux(auxflag);
(rec->back()).setFlags(0); // this sets all flag bits to 0
if (hbheTimingShapedFlagSetter_!=0)
hbheTimingShapedFlagSetter_->SetTimingShapedFlags(rec->back());
if (setNoiseFlags_)
hbheFlagSetter_->SetFlagsFromDigi(rec->back(),*i,coder,calibrations);
if (setPulseShapeFlags_ == true)
hbhePulseShapeFlagSetter_->SetPulseShapeFlags(rec->back(), *i, coder, calibrations);
if (setSaturationFlags_)
saturationFlagSetter_->setSaturationFlag(rec->back(),*i);
if (correctTiming_)
HcalTimingCorrector::Correct(rec->back(), *i, favorite_capid);
if (setHSCPFlags_ && i->id().ietaAbs()<16)
{
double DigiEnergy=0;
for(int j=0; j!=i->size(); DigiEnergy += i->sample(j++).nominal_fC());
if(DigiEnergy > hbheHSCPFlagSetter_->EnergyThreshold())
{
HBDigis.push_back(*i);
RecHitIndex.push_back(rec->size()-1);
}
} // if (set HSCPFlags_ && |ieta|<16)
} // loop over HBHE digis
if (setNoiseFlags_) hbheFlagSetter_->SetFlagsFromRecHits(*rec);
if (setHSCPFlags_) hbheHSCPFlagSetter_->hbheSetTimeFlagsFromDigi(rec.get(), HBDigis, RecHitIndex);
// return result
e.put(rec);
} else if (subdet_==HcalOuter) {
edm::Handle<HODigiCollection> digi;
e.getByLabel(inputLabel_,digi);
// create empty output
std::auto_ptr<HORecHitCollection> rec(new HORecHitCollection);
rec->reserve(digi->size());
// run the algorithm
HODigiCollection::const_iterator i;
// Vote on majority TS0 CapId
int favorite_capid = 0;
if (correctTiming_) {
long capid_votes[4] = {0,0,0,0};
for (i=digi->begin(); i!=digi->end(); i++) {
capid_votes[(*i)[0].capid()]++;
}
for (int k = 0; k < 4; k++)
if (capid_votes[k] > capid_votes[favorite_capid])
favorite_capid = k;
}
for (i=digi->begin(); i!=digi->end(); i++) {
HcalDetId cell = i->id();
DetId detcell=(DetId)cell;
// check on cells to be ignored and dropped: (rof,20.Feb.09)
const HcalChannelStatus* mydigistatus=myqual->getValues(detcell.rawId());
if (mySeverity->dropChannel(mydigistatus->getValue() ) ) continue;
if (dropZSmarkedPassed_)
if (i->zsMarkAndPass()) continue;
const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
HcalCoderDb coder (*channelCoder, *shape);
rec->push_back(reco_.reconstruct(*i,coder,calibrations));
// Set auxiliary flag
int auxflag=0;
for (int xx=firstauxTS_;xx<firstauxTS_+4 && xx<i->size();++xx)
auxflag+=(i->sample(xx).adc())<<(7*(xx-firstauxTS_)); // store the time slices in the first 28 bits of aux, a set of 4 7-bit adc values
// bits 28 and 29 are reserved for capid of the first time slice saved in aux
auxflag+=((i->sample(firstauxTS_).capid())<<28);
(rec->back()).setAux(auxflag);
(rec->back()).setFlags(0);
if (setSaturationFlags_)
saturationFlagSetter_->setSaturationFlag(rec->back(),*i);
if (correctTiming_)
HcalTimingCorrector::Correct(rec->back(), *i, favorite_capid);
}
// return result
e.put(rec);
} else if (subdet_==HcalForward) {
edm::Handle<HFDigiCollection> digi;
e.getByLabel(inputLabel_,digi);
// ugly hack only for purposes of 3.11 HF treatment
if (e.isRealData() && e.run() <= 153943)
{
reco_.resetTimeSamples(3,4);
if (hfdigibit_) hfdigibit_->resetTimeSamples(3,4);
firstauxTS_=3; // hard-code starting position of aux word
}
else
{
reco_.resetTimeSamples(4,2);
if (hfdigibit_) hfdigibit_->resetTimeSamples(3,3); // flag uses 3 TS, even if reco uses 2 TS
firstauxTS_=3; // hard-code
}
// create empty output
std::auto_ptr<HFRecHitCollection> rec(new HFRecHitCollection);
rec->reserve(digi->size());
// run the algorithm
HFDigiCollection::const_iterator i;
// Vote on majority TS0 CapId
int favorite_capid = 0;
if (correctTiming_) {
long capid_votes[4] = {0,0,0,0};
for (i=digi->begin(); i!=digi->end(); i++) {
capid_votes[(*i)[0].capid()]++;
}
for (int k = 0; k < 4; k++)
if (capid_votes[k] > capid_votes[favorite_capid])
favorite_capid = k;
}
for (i=digi->begin(); i!=digi->end(); i++) {
HcalDetId cell = i->id();
DetId detcell=(DetId)cell;
// check on cells to be ignored and dropped: (rof,20.Feb.09)
const HcalChannelStatus* mydigistatus=myqual->getValues(detcell.rawId());
if (mySeverity->dropChannel(mydigistatus->getValue() ) ) continue;
if (dropZSmarkedPassed_)
if (i->zsMarkAndPass()) continue;
const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
HcalCoderDb coder (*channelCoder, *shape);
rec->push_back(reco_.reconstruct(*i,coder,calibrations));
// Set auxiliary flag
int auxflag=0;
for (int xx=firstauxTS_;xx<firstauxTS_+4 && xx<i->size();++xx)
auxflag+=(i->sample(xx).adc())<<(7*(xx-firstauxTS_)); // store the time slices in the first 28 bits of aux, a set of 4 7-bit adc values
// bits 28 and 29 are reserved for capid of the first time slice saved in aux
auxflag+=((i->sample(firstauxTS_).capid())<<28);
(rec->back()).setAux(auxflag);
// Clear flags
(rec->back()).setFlags(0);
// This calls the code for setting the HF noise bit determined from digi shape
if (setNoiseFlags_) hfdigibit_->hfSetFlagFromDigi(rec->back(),*i,coder,calibrations);
if (setSaturationFlags_)
saturationFlagSetter_->setSaturationFlag(rec->back(),*i);
if (setTimingTrustFlags_)
HFTimingTrustFlagSetter_->setHFTimingTrustFlag(rec->back(),*i);
if (correctTiming_)
HcalTimingCorrector::Correct(rec->back(), *i, favorite_capid);
} // for (i=digi->begin(); i!=digi->end(); i++) -- loop on all HF digis
// The following flags require the full set of rechits
// These need to be set consecutively, so an energy check should be the first
// test performed on these hits (to minimize the loop time)
if (setNoiseFlags_)
{
// Step 1: Set PET flag (short fibers of |ieta|==29)
// Neighbor/partner channels that are flagged by Pulse Shape algorithm (HFDigiTime)
// won't be considered in these calculations
for (HFRecHitCollection::iterator i = rec->begin();i!=rec->end();++i)
{
int depth=i->id().depth();
int ieta=i->id().ieta();
// Short fibers and all channels at |ieta|=29 use PET settings in Algo 3
if (depth==2 || abs(ieta)==29 )
hfPET_->HFSetFlagFromPET(*i,*rec,myqual,mySeverity);
}
// Step 2: Set S8S1 flag (short fibers or |ieta|==29)
for (HFRecHitCollection::iterator i = rec->begin();i!=rec->end();++i)
{
int depth=i->id().depth();
int ieta=i->id().ieta();
// Short fibers and all channels at |ieta|=29 use PET settings in Algo 3
if (depth==2 || abs(ieta)==29 )
hfS8S1_->HFSetFlagFromS9S1(*i,*rec,myqual,mySeverity);
}
// Set 3: Set S9S1 flag (long fibers)
for (HFRecHitCollection::iterator i = rec->begin();i!=rec->end();++i)
{
int depth=i->id().depth();
int ieta=i->id().ieta();
// Short fibers and all channels at |ieta|=29 use PET settings in Algo 3
if (depth==1 && abs(ieta)!=29 )
hfS9S1_->HFSetFlagFromS9S1(*i,*rec,myqual, mySeverity);
}
}
// return result
e.put(rec);
} else if (subdet_==HcalOther && subdetOther_==HcalCalibration) {
edm::Handle<HcalCalibDigiCollection> digi;
e.getByLabel(inputLabel_,digi);
// create empty output
std::auto_ptr<HcalCalibRecHitCollection> rec(new HcalCalibRecHitCollection);
rec->reserve(digi->size());
// run the algorithm
HcalCalibDigiCollection::const_iterator i;
for (i=digi->begin(); i!=digi->end(); i++) {
HcalCalibDetId cell = i->id();
DetId detcell=(DetId)cell;
// check on cells to be ignored and dropped: (rof,20.Feb.09)
const HcalChannelStatus* mydigistatus=myqual->getValues(detcell.rawId());
if (mySeverity->dropChannel(mydigistatus->getValue() ) ) continue;
if (dropZSmarkedPassed_)
if (i->zsMarkAndPass()) continue;
const HcalCalibrations& calibrations=conditions->getHcalCalibrations(cell);
const HcalQIECoder* channelCoder = conditions->getHcalCoder (cell);
HcalCoderDb coder (*channelCoder, *shape);
rec->push_back(reco_.reconstruct(*i,coder,calibrations));
/*
// Flag setting not available for calibration rechits
// Set auxiliary flag
int auxflag=0;
for (int xx=firstauxTS_;xx<firstauxTS_+4 && xx<i->size();++xx)
auxflag+=(i->sample(xx).adc())<<(7*(xx-firstauxTS_)); // store the time slices in the first 28 bits of aux, a set of 4 7-bit adc values
// bits 28 and 29 are reserved for capid of the first time slice saved in aux
auxflag+=((i->sample(firstauxTS_).capid())<<28);
(rec->back()).setAux(auxflag);
(rec->back()).setFlags(0); // Not yet implemented for HcalCalibRecHit
*/
}
// return result
e.put(rec);
}
}
delete myqual;
} // void HcalHitReconstructor::produce(...)
bool HcalHitReconstructor::correctTiming_ [private] |
Definition at line 58 of file HcalHitReconstructor.h.
Referenced by produce().
DetId::Detector HcalHitReconstructor::det_ [private] |
Definition at line 53 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
bool HcalHitReconstructor::dropZSmarkedPassed_ [private] |
Definition at line 64 of file HcalHitReconstructor.h.
Referenced by produce().
int HcalHitReconstructor::firstauxTS_ [private] |
Definition at line 66 of file HcalHitReconstructor.h.
Referenced by produce().
Definition at line 44 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), produce(), and ~HcalHitReconstructor().
Definition at line 45 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), produce(), and ~HcalHitReconstructor().
Definition at line 47 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), produce(), and ~HcalHitReconstructor().
Definition at line 46 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
Definition at line 48 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), produce(), and ~HcalHitReconstructor().
HcalHF_PETalgorithm* HcalHitReconstructor::hfPET_ [private] |
Definition at line 51 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), produce(), and ~HcalHitReconstructor().
Definition at line 50 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
Definition at line 49 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), produce(), and ~HcalHitReconstructor().
Definition at line 43 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
Definition at line 56 of file HcalHitReconstructor.h.
Referenced by produce().
HcalSimpleRecAlgo HcalHitReconstructor::reco_ [private] |
Definition at line 41 of file HcalHitReconstructor.h.
Referenced by produce().
Definition at line 42 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
bool HcalHitReconstructor::setHSCPFlags_ [private] |
Definition at line 60 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
bool HcalHitReconstructor::setNoiseFlags_ [private] |
Definition at line 59 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
bool HcalHitReconstructor::setPulseShapeFlags_ [private] |
Definition at line 63 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
bool HcalHitReconstructor::setSaturationFlags_ [private] |
Definition at line 61 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
bool HcalHitReconstructor::setTimingTrustFlags_ [private] |
Definition at line 62 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
int HcalHitReconstructor::subdet_ [private] |
Definition at line 54 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
Definition at line 55 of file HcalHitReconstructor.h.
Referenced by HcalHitReconstructor(), and produce().
1.7.3