166 art::ServiceHandle<art::TFileService>
tfs;
167 for(
int i = 30; i < 50 + 1; i++){
168 std::string
tagger =
"All";
169 if (i > 34 && i < 40)
continue;
170 if (i==48 || i==49)
continue;
173 std::cout <<
"tagger: " << tagger.c_str() << std::endl;
174 hDCA[tagger] = tfs->make<TH1D>(Form(
"DCA_%s", tagger.c_str()),
"", 50, 0, 100);
175 hMatchDCA[tagger] = tfs->make<TH1D>(Form(
"MatchDCA_%s", tagger.c_str()),
"", 50, 0, 100);
176 hNoMatchDCA[tagger] = tfs->make<TH1D>(Form(
"NoMatchDCA_%s", tagger.c_str()),
"", 50, 0, 100);
178 hDoL[tagger] = tfs->make<TH1D>(Form(
"DoL_%s", tagger.c_str()),
"", 100, 0, 0.25);
179 hMatchDoL[tagger] = tfs->make<TH1D>(Form(
"MatchDoL_%s", tagger.c_str()),
"", 100, 0, 0.25);
180 hNoMatchDoL[tagger] = tfs->make<TH1D>(Form(
"NoMatchDoL_%s", tagger.c_str()),
"", 100, 0, 0.25);
182 hT0[tagger] = tfs->make<TH1D>(Form(
"T0_%s", tagger.c_str()),
"", 600, -3000, 3000);
183 hMatchT0[tagger] = tfs->make<TH1D>(Form(
"MatchT0_%s", tagger.c_str()),
"", 600, -3000, 3000);
184 hNoMatchT0[tagger] = tfs->make<TH1D>(Form(
"NoMatchT0_%s", tagger.c_str()),
"", 600, -3000, 3000);
186 hEffDCATotal[tagger] = tfs->make<TH1D>(Form(
"EffDCATotal_%s", tagger.c_str()),
"", 50, 0, 100);
187 hEffDCAReco[tagger] = tfs->make<TH1D>(Form(
"EffDCAReco_%s", tagger.c_str()),
"", 50, 0, 100);
188 hEffDoLTotal[tagger] = tfs->make<TH1D>(Form(
"EffDoLTotal_%s", tagger.c_str()),
"", 100, 0, 0.25);
189 hEffDoLReco[tagger] = tfs->make<TH1D>(Form(
"EffDoLReco_%s", tagger.c_str()),
"", 100, 0, 0.25);
190 hEffLengthTotal[tagger] = tfs->make<TH1D>(Form(
"EffLengthTotal_%s", tagger.c_str()),
"", 20, 0, 600);
191 hEffLengthReco[tagger] = tfs->make<TH1D>(Form(
"EffLengthReco_%s", tagger.c_str()),
"", 20, 0, 600);
193 hPurityDCATotal[tagger] = tfs->make<TH1D>(Form(
"PurityDCATotal_%s", tagger.c_str()),
"", 50, 0, 100);
194 hPurityDCAReco[tagger] = tfs->make<TH1D>(Form(
"PurityDCAReco_%s", tagger.c_str()),
"", 50, 0, 100);
195 hPurityDoLTotal[tagger] = tfs->make<TH1D>(Form(
"PurityDoLTotal_%s", tagger.c_str()),
"", 100, 0, 0.25);
196 hPurityDoLReco[tagger] = tfs->make<TH1D>(Form(
"PurityDoLReco_%s", tagger.c_str()),
"", 100, 0, 0.25);
197 hPurityLengthTotal[tagger] = tfs->make<TH1D>(Form(
"PurityLengthTotal_%s", tagger.c_str()),
"", 20, 0, 600);
198 hPurityLengthReco[tagger] = tfs->make<TH1D>(Form(
"PurityLengthReco_%s", tagger.c_str()),
"", 20, 0, 600);
202 if(
fVerbose)
std::cout<<
"----------------- CRT T0 Matching Ana Module -------------------"<<std::endl;
bool fVerbose
print information about what's going on
std::map< std::string, TH1D * > hNoMatchDCA
std::map< std::string, TH1D * > hDoL
std::map< std::string, TH1D * > hPurityDoLReco
std::map< std::string, TH1D * > hNoMatchT0
std::map< std::string, TH1D * > hEffDoLReco
std::map< std::string, TH1D * > hEffDCAReco
std::map< std::string, TH1D * > hPurityLengthReco
std::map< std::string, TH1D * > hPurityDCATotal
std::map< std::string, TH1D * > hMatchDCA
std::map< std::string, TH1D * > hPurityDCAReco
std::map< std::string, TH1D * > hNoMatchDoL
std::map< std::string, TH1D * > hT0
std::map< std::string, TH1D * > hEffDoLTotal
string GetRegionNameFromNum(int num)
std::map< std::string, TH1D * > hEffLengthReco
std::map< std::string, TH1D * > hPurityLengthTotal
std::map< std::string, TH1D * > hEffDCATotal
std::map< std::string, TH1D * > hMatchDoL
std::map< std::string, TH1D * > hDCA
std::map< std::string, TH1D * > hPurityDoLTotal
std::map< std::string, TH1D * > hEffLengthTotal
art::ServiceHandle< art::TFileService > tfs
icarus::crt::CRTCommonUtils * fCrtutils
std::map< std::string, TH1D * > hMatchT0
BEGIN_PROLOG could also be cout