29{
30
33
34 std::vector<int> isampmap(CaloSampling::Unknown,-1);
35 unsigned int nAreas(0);
36 std::vector<CaloLocalHadCoeff::LocalHadArea> theAreas;
37 for (unsigned int iArea=0;iArea<CaloSampling::Unknown;iArea++) {
39 theAreas.push_back(std::move(theArea));
40 }
41
42 for (int iread=0;iread<2;iread++) {
43
44 TFile* theLCWeightFile =
new TFile(
file.c_str());
45 if ( !theLCWeightFile ) {
46 return StatusCode::FAILURE;
47 }
48 TList * theKeyList = theLCWeightFile->GetListOfKeys();
49 while ( theKeyList->GetEntries() ) {
50 TProfile2D * prof = (TProfile2D *)theLCWeightFile->Get(theKeyList->First()->GetTitle());
51
52 std::string sTitle(prof->GetTitle());
53
54 std::string
key(
"_isamp_");
55 size_t found = sTitle.find(key);
56 if ( found == std::string::npos ) {
57 ATH_MSG_ERROR(
"Could not find key " << key <<
" in current histogram." );
58 }
59 else {
60 std::istringstream tstr(sTitle.substr(found+
key.length()));
61 int isamp(-1);
62 tstr >> isamp;
63 if ( isamp <= 0 || isamp >= CaloSampling::Unknown ) {
64 ATH_MSG_ERROR(
"Found invalid sampling " << isamp <<
" in current histogram." );
65 }
66 else {
67
68 bool allValid(true);
69
70 std::vector<std::string>
keys;
71 keys.emplace_back(
"_iside_");
72 keys.emplace_back(
"_ieta_");
73 keys.emplace_back(
"_iphi_");
74
75 std::vector<std::string>
names;
76 names.emplace_back(
"side");
77 names.emplace_back(
"|eta|");
78 names.emplace_back(
"phi");
79 names.emplace_back(prof->GetXaxis()->GetTitle());
80 names.emplace_back(prof->GetYaxis()->GetTitle());
81
82 std::vector<int>
types;
88
89 std::vector<int> ibin(
names.size(),-1);
90 std::vector<double> rmin(
names.size(),-1);
91 std::vector<double> rmax(
names.size(),-1);
92 std::vector<int> nbin(
names.size(),-1);
93 unsigned int idim;
94
95 for (idim=0;idim<
keys.size();idim++) {
96 found = sTitle.find(keys[idim]);
97 if ( found == std::string::npos ) {
98 ATH_MSG_ERROR(
"Could not find key " << keys[idim] <<
" in current histogram." );
99 allValid = false;
100 }
101 else {
103 std::istringstream tstr(sTitle.substr(found+keys[idim].length()));
104 tstr >> ibin[idim] >>
c >>
c >> rmin[idim] >>
c >> rmax[idim] >>
c >> nbin[idim];
105 if ( ibin[idim] < 0 || ibin[idim] >= nbin[idim] ) {
106 ATH_MSG_ERROR(
"Found invalid bin number " << ibin[idim] <<
" not in valid range [0," << nbin[idim] <<
" in current histogram." );
107 allValid = false;
108 }
109 }
110 }
111 if ( allValid ) {
112
113 nbin[
names.size()-2] = prof->GetNbinsX();
114 rmin[
names.size()-2] = prof->GetXaxis()->GetXmin();
115 rmax[
names.size()-2] = prof->GetXaxis()->GetXmax();
116 nbin[
names.size()-1] = prof->GetNbinsY();
117 rmin[
names.size()-1] = prof->GetYaxis()->GetXmin();
118 rmax[
names.size()-1] = prof->GetYaxis()->GetXmax();
119
120 if ( isampmap[isamp] < 0 ) {
121 isampmap[isamp] = nAreas;
122 nAreas++;
123 for (idim = 0;idim<
names.size();idim++) {
124 CaloLocalHadCoeff::LocalHadDimension theDim(names[idim].c_str(),
types[idim],nbin[idim],rmin[idim],rmax[idim]);
125 theAreas[isamp].addDimension(theDim);
126 }
127 }
128 else if ( iread == 1 ) {
129
130 TAxis * xax = prof->GetXaxis();
131 TAxis * yax = prof->GetYaxis();
132 for (ibin[
names.size()-2]=0;ibin[
names.size()-2]<prof->GetNbinsX();ibin[
names.size()-2]++) {
133 for (ibin[
names.size()-1]=0;ibin[
names.size()-1]<prof->GetNbinsY();ibin[
names.size()-1]++) {
134 float logEClus = xax->GetBinCenter(ibin[
names.size()-2]+1);
135 float logEDensCell = yax->GetBinCenter(ibin[
names.size()-1]+1);
136 int iBin = prof->FindBin(logEClus,logEDensCell);
137
142
143 msg() << MSG::INFO <<
"Now set data for isamp=" << isamp
144 << " and bins: ";
145 for(unsigned int ii=0;ii<ibin.size();ii++)
146 msg() << ibin[ii] <<
" ";
148 int bin =
data.getBin(isampmap[isamp],ibin);
149 if (bin < 0) {
151 return StatusCode::FAILURE;
152 }
153 data.setCoeff(bin,theData);
154 }
155 }
156 }
157 }
158 }
159 }
160 theKeyList->RemoveFirst();
161 }
162 if ( iread == 0 ) {
163
164 for(unsigned int jsamp=0;jsamp<theAreas.size();jsamp++) {
165 if ( theAreas[jsamp].getNdim() > 0 ) {
166 ATH_MSG_INFO(
"adding Area for jsamp=" << jsamp <<
" with nDim = " << theAreas[jsamp].getNdim() );
167 data.addArea(theAreas[jsamp]);
168 }
169 }
170 }
171 theLCWeightFile->Close();
172 delete theLCWeightFile;
173 }
174 return StatusCode::SUCCESS;
175}
char data[hepevt_bytes_allocation_ATLAS]
static const std::vector< std::string > types
std::vector< float > LocalHadCoeff
Correction parameters for one general bin.
static const std::string & getSamplingName(const CaloSampling::CaloSample theSample)
Returns a string (name) for each CaloSampling.
static std::string find_file(const std::string &logical_file_name, const std::string &search_path)
1.14.0