Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Ecalpf cluster corrections93 x - taking care of events with flag==5 #19827

Merged
merged 6 commits into from Aug 1, 2017
Merged

Ecalpf cluster corrections93 x - taking care of events with flag==5 #19827

Changes from all commits
Commits
Show all changes
6 commits
Select commit Hold shift + click to select a range
bff41de
Taken care of flag==5
jainshilpi Jul 19, 2017
6f95120
Taken care of flag==5
jainshilpi Jul 19, 2017
f9ce512
Taken care of flags 2,6 and > 7
jainshilpi Jul 26, 2017
6472198
Changed the line not ZS_bit to ZS_bit ne 0
jainshilpi Jul 28, 2017
bbebe81
Changed the line not FR_bit to FR_bit ne 0
jainshilpi Jul 28, 2017
d067060
corrected the ZS and FR conditions - bug which was introduced earlier
jainshilpi Jul 31, 2017
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Jump to
Jump to file
Failed to load files.
Diff view
Diff view
25 changes: 20 additions & 5 deletions RecoParticleFlow/PFClusterProducer/src/PFClusterEMEnergyCorrector.cc
View file
Open in desktop
Expand Up @@ -335,18 +335,33 @@ void PFClusterEMEnergyCorrector::correctEnergies(const edm::Event &evt,
int coridx = 0;
int regind = 0;
if (!iseb) regind = 4;
if (clusFlag==1)

///////////////////////////////////////////////////////////////////////////////////
///a hit can be ZS or forced ZS. A hit can be in Full readout or Forced to be FULL readout
///if it is ZS, then clusFlag (in binary) = 0001
///if it is forced ZS, then clusFlag (in binary) = 0101
///if it is FR, then clusFlag (in binary) = 0011
///if it is forced FR, then clusFlag (in binary) = 0111
///i.e 3rd bit is set.
///Even if it is forced, we should mark it is as ZS or FR. To take care of it, just check the LSB and second LSB(SLSB)
///////////////////////////////////////////////////////////////////////////////////
int ZS_bit = clusFlag>>0&1;
int FR_bit = clusFlag>>1&1;



if (ZS_bit!=0 && FR_bit==0) ///it is clusFlag==1, 5
coridx = 0 + regind;
else {
else{
if (pt<2.5) coridx = 1 + regind;
else if (pt>=2.5 && pt<6.) coridx = 2 + regind;
else if (pt>=6.) coridx = 3 + regind;
}
if (clusFlag!=1 && clusFlag!=3) {
edm::LogWarning("PFClusterEMEnergyCorrector") << "We can only correct regions readout in ZS (flag 1) or FULL readout (flag 3). Flag " << clusFlag << " is not recognized."
if (ZS_bit==0 || clusFlag > 7) {
edm::LogWarning("PFClusterEMEnergyCorrector") << "We can only correct regions readout in ZS (flag 1,5) or FULL readout (flag 3,7). Flag " << clusFlag << " is not recognized."
<< "\n" << "Assuming FULL readout and continuing";
}

const GBRForestD &meanforest = *forestH_mean[coridx].product();
const GBRForestD &sigmaforest = *forestH_sigma[coridx].product();

Expand Down