Skip to content
Snippets Groups Projects
Commit fc487da6 authored by scott snyder's avatar scott snyder
Browse files

FTKStandaloneMonitoring: Comply with ATLAS naming conventions. 

Private data members should start with m_.
(Constants can be all caps.)
parent 7c6214b4
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Trigger/TrigFTK/FTKStandaloneMonitoring/FTKStandaloneMonitoring/CompareFTKEvents.h
+ 16
− 16
View file @ fc487da6
/* /*
Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/ */
#ifndef __CompareFTKEvents__ #ifndef FTKSTANDALONEMONITORING_COMPAREFTKEVENTS_H
#define __CompareFTKEvents__ #define FTKSTANDALONEMONITORING_COMPAREFTKEVENTS_H
#include <iostream> #include <iostream>
...@@ -45,18 +45,18 @@ class CompareFTKEvents{ ...@@ -45,18 +45,18 @@ class CompareFTKEvents{
void unpackPixCluster(OFFLINE_FRAGMENTS_NAMESPACE::PointerType data, FTK_RawPixelCluster& cluster); void unpackPixCluster(OFFLINE_FRAGMENTS_NAMESPACE::PointerType data, FTK_RawPixelCluster& cluster);
void unpackSCTCluster(OFFLINE_FRAGMENTS_NAMESPACE::PointerType data, FTK_RawSCT_Cluster& cluster) { cluster.setWord(*data); } void unpackSCTCluster(OFFLINE_FRAGMENTS_NAMESPACE::PointerType data, FTK_RawSCT_Cluster& cluster) { cluster.setWord(*data); }
int m_nevtmax=0; int m_nevtmax=0;
const size_t TrackBlobSize = 22; // magic number from BS specification static const size_t TRACKBLOBSIZE = 22; // magic number from BS specification
const size_t TrackParamsBlobSize = 6; // --||-- static const size_t TRACKPARAMSBLOBSIZE = 6; // --||--
const size_t PixHitParamsBlobSize = 2; // --||-- static const size_t PIXHITPARAMSBLOBSIZE = 2; // --||--
const size_t SCTHitParamsBlobSize = 1; // --||-- static const size_t SCTHITPARAMSBLOBSIZE = 1; // --||--
const size_t NPixLayers=4; static const size_t NPIXLAYERS=4;
const size_t NSCTLayers=8; static const size_t NSCTLAYERS=8;
std::vector<uint32_t> m_data; std::vector<uint32_t> m_data;
std::string m_BSfile; std::string m_BSfile;
std::string m_NTUP_FTK_file; std::string m_NTUP_FTK_file;
uint32_t m_checkword=0xaa1234aa; uint32_t m_checkword=0xaa1234aa;
uint32_t trackBlockOffsetStart=0; // Nr of words in ROB data block before track blocks start uint32_t m_trackBlockOffsetStart=0; // Nr of words in ROB data block before track blocks start
uint32_t trackBlockOffsetEnd=0; // Nr of words in ROB data block after track blocks end uint32_t m_trackBlockOffsetEnd=0; // Nr of words in ROB data block after track blocks end
int m_Nevents_BS; int m_Nevents_BS;
int m_Nevents_NTUPFTK; int m_Nevents_NTUPFTK;
std::ifstream m_myBSfile; std::ifstream m_myBSfile;
...@@ -64,15 +64,15 @@ class CompareFTKEvents{ ...@@ -64,15 +64,15 @@ class CompareFTKEvents{
TTree *m_tevtinfo; TTree *m_tevtinfo;
TTree *m_theTree; TTree *m_theTree;
TFile *m_fntupftk; TFile *m_fntupftk;
int RN, EN, LB, BCID, EL1ID, l1TT; int m_RN, m_EN, m_LB, m_BCID, m_EL1ID, m_l1TT;
FTKTrackStream *m_ft; FTKTrackStream *m_ft;
bool m_allmatched=true; bool m_allmatched=true;
bool m_verbose=false; bool m_verbose=false;
bool m_setup_partition=false; bool m_setup_partition=false;
CompareFTKTracks * m_compTrk; CompareFTKTracks * m_compTrk;
std::vector<std::string> m_histo_list; std::vector<std::string> m_histo_list;
std::map<std::string , TH1D * > map_histo; std::map<std::string , TH1D * > m_map_histo;
std::map<std::string , TH2D * > map_histo_2D; std::map<std::string , TH2D * > m_map_histo_2D;
TFile * m_fout; TFile * m_fout;
void CreateHistos(); void CreateHistos();
void WriteHistos(); void WriteHistos();
...@@ -80,13 +80,13 @@ class CompareFTKEvents{ ...@@ -80,13 +80,13 @@ class CompareFTKEvents{
std::string m_str_tree_ftkdata="ftkdata"; std::string m_str_tree_ftkdata="ftkdata";
std::string m_str_tree_ftkstream="FTKMergedTracksStream"; std::string m_str_tree_ftkstream="FTKMergedTracksStream";
std::string m_partition_name=""; std::string m_partition_name="";
std::map<std::string, std::vector<double> > histo_param={{"pt",{1000.,0.,100000.}}, std::map<std::string, std::vector<double> > m_histo_param={{"pt",{1000.,0.,100000.}},
{"eta",{100.,-2.5,2.5}}, {"eta",{100.,-2.5,2.5}},
{"phi",{100.,-3.2,3.2}}, {"phi",{100.,-3.2,3.2}},
{"d0",{100.,-5.,5.}}, {"d0",{100.,-5.,5.}},
{"z0",{100,-50,50}}, {"z0",{100,-50,50}},
{"chi2",{100,0,50}}, {"chi2",{100,0,50}},
{"ETA_PHI",{100,-2.5,2.5,100,-3.2,3.2}}}; {"ETA_PHI",{100,-2.5,2.5,100,-3.2,3.2}}};
std::vector<std::string> variable_list={"pt","eta","phi","d0","z0","chi2","ETA_PHI"}; std::vector<std::string> m_variable_list={"pt","eta","phi","d0","z0","chi2","ETA_PHI"};
}; };
#endif //__CompareFTKEvents__ #endif //FTKSTANDALONEMONITORING_COMPAREFTKEVENTS_H
Trigger/TrigFTK/FTKStandaloneMonitoring/FTKStandaloneMonitoring/CompareFTKTracks.h
+ 4
− 4
View file @ fc487da6
/* /*
Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/ */
#ifndef __CompareFTKTracks__ #ifndef FTKSTANDALONEMONITORING_COMPAREFTKTRACKS_H
#define __CompareFTKTracks__ #define FTKSTANDALONEMONITORING_COMPAREFTKTRACKS_H
#include <iostream> #include <iostream>
...@@ -34,6 +34,6 @@ class CompareFTKTracks{ ...@@ -34,6 +34,6 @@ class CompareFTKTracks{
FTKTrkAssoc* m_associator; FTKTrkAssoc* m_associator;
double m_dmax=0.2; double m_dmax=0.2;
bool m_allmatched=true; bool m_allmatched=true;
std::vector<std::string> variable_list={"pt","eta","phi","d0","z0","chi2","ETA_PHI"}; std::vector<std::string> m_variable_list={"pt","eta","phi","d0","z0","chi2","ETA_PHI"};
}; };
#endif //__CompareFTKTracks__ #endif //FTKSTANDALONEMONITORING_COMPAREFTKTRACKS_H
Trigger/TrigFTK/FTKStandaloneMonitoring/src/CompareFTKEvents.cxx
+ 42
− 42
View file @ fc487da6
...@@ -38,36 +38,36 @@ void CompareFTKEvents::SetHistos(std::vector<std::string> histo_list){ ...@@ -38,36 +38,36 @@ void CompareFTKEvents::SetHistos(std::vector<std::string> histo_list){
std::string title = ss.str(); std::string title = ss.str();
//if histo is not eta vs phi 2D histo //if histo is not eta vs phi 2D histo
if (istr.find("ETA_PHI")==std::string::npos){ if (istr.find("ETA_PHI")==std::string::npos){
if (istr.find("nTrk")!=std::string::npos) map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),title.c_str(),100,0,100))); if (istr.find("nTrk")!=std::string::npos) m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),title.c_str(),100,0,100)));
else if (istr.find("HWSW")!=std::string::npos|| istr.find("HWonly")!=std::string::npos||istr.find("SWonly")!=std::string::npos||(istr.find("HWvsSW")!=std::string::npos&&istr.find("diff")==std::string::npos)) else if (istr.find("HWSW")!=std::string::npos|| istr.find("HWonly")!=std::string::npos||istr.find("SWonly")!=std::string::npos||(istr.find("HWvsSW")!=std::string::npos&&istr.find("diff")==std::string::npos))
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),title.c_str(),histo_param[str3].at(0),histo_param[str3].at(1),histo_param[str3].at(2)))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),title.c_str(),m_histo_param[str3].at(0),m_histo_param[str3].at(1),m_histo_param[str3].at(2))));
else if (istr.find("HWvsSWdiff_pt")!=std::string::npos) else if (istr.find("HWvsSWdiff_pt")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta p_{T} (test-ref); p_{T}^{test} - p_{T}^{ref} ;FTK Tracks",1000,-5000,5000))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta p_{T} (test-ref); p_{T}^{test} - p_{T}^{ref} ;FTK Tracks",1000,-5000,5000)));
else if (istr.find("HWvsSWdiff_eta")!=std::string::npos) else if (istr.find("HWvsSWdiff_eta")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta #eta (test-ref); #eta^{test} - #eta^{ref} ;FTK Tracks",2000,-0.1,0.1))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta #eta (test-ref); #eta^{test} - #eta^{ref} ;FTK Tracks",2000,-0.1,0.1)));
else if (istr.find("HWvsSWdiff_phi")!=std::string::npos) else if (istr.find("HWvsSWdiff_phi")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta #phi (test-ref); #phi^{test} - #phi^{ref} ;FTK Tracks",2000,-0.1,0.1))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta #phi (test-ref); #phi^{test} - #phi^{ref} ;FTK Tracks",2000,-0.1,0.1)));
else if (istr.find("HWvsSWdiff_d0")!=std::string::npos) else if (istr.find("HWvsSWdiff_d0")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta d_{0} (test-ref); d_{0}^{test} - d_{0}^{ref} ;FTK Tracks",4000,-5.,5.))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta d_{0} (test-ref); d_{0}^{test} - d_{0}^{ref} ;FTK Tracks",4000,-5.,5.)));
else if (istr.find("HWvsSWdiff_z0")!=std::string::npos) else if (istr.find("HWvsSWdiff_z0")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta z_{0} (test-ref); z_{0}^{test} - z_{0}^{ref} ;FTK Tracks",4000,-2.,2.))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta z_{0} (test-ref); z_{0}^{test} - z_{0}^{ref} ;FTK Tracks",4000,-2.,2.)));
else if (istr.find("HWvsSWdiff_chi2")!=std::string::npos) else if (istr.find("HWvsSWdiff_chi2")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta chi2 (test-ref); chi2^{test} - chi2^{ref} ;FTK Tracks",4000,-5.,5.))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta chi2 (test-ref); chi2^{test} - chi2^{ref} ;FTK Tracks",4000,-5.,5.)));
else if (istr.find("res_pt")!=std::string::npos) else if (istr.find("res_pt")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta p_{T} (test-ref)/ p_{T}^{ref};( p_{T}^{test} - p_{T}^{ref} )/ p_{T}^{ref};FTK Tracks",2000,-1.,1.))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta p_{T} (test-ref)/ p_{T}^{ref};( p_{T}^{test} - p_{T}^{ref} )/ p_{T}^{ref};FTK Tracks",2000,-1.,1.)));
else if (istr.find("res_eta")!=std::string::npos) else if (istr.find("res_eta")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta #eta (test-ref)/ #eta^{ref};( #eta^{test} - #eta^{ref} )/ #eta^{ref};FTK Tracks",2000,-1.,1.))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta #eta (test-ref)/ #eta^{ref};( #eta^{test} - #eta^{ref} )/ #eta^{ref};FTK Tracks",2000,-1.,1.)));
else if (istr.find("res_phi")!=std::string::npos) else if (istr.find("res_phi")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta #phi (test-ref)/ #phi^{ref};( #phi^{test} - #phi^{ref} )/ #phi^{ref};FTK Tracks",2000,-1.,1.))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta #phi (test-ref)/ #phi^{ref};( #phi^{test} - #phi^{ref} )/ #phi^{ref};FTK Tracks",2000,-1.,1.)));
else if (istr.find("res_d0")!=std::string::npos) else if (istr.find("res_d0")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta d_{0} (test-ref)/ d_{0}^{ref};( d_{0}^{test} - d_{0}^{ref} )/ d_{0}^{ref};FTK Tracks",4000,-2.,2.))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta d_{0} (test-ref)/ d_{0}^{ref};( d_{0}^{test} - d_{0}^{ref} )/ d_{0}^{ref};FTK Tracks",4000,-2.,2.)));
else if (istr.find("res_z0")!=std::string::npos) else if (istr.find("res_z0")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta z_{0} (test-ref)/ z_{0}^{ref};( z_{0}^{test} - z_{0}^{ref} )/ z_{0}^{ref};FTK Tracks",4000,-2.,2.))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta z_{0} (test-ref)/ z_{0}^{ref};( z_{0}^{test} - z_{0}^{ref} )/ z_{0}^{ref};FTK Tracks",4000,-2.,2.)));
else if (istr.find("res_chi2")!=std::string::npos) else if (istr.find("res_chi2")!=std::string::npos)
map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta chi2 (test-ref)/ chi2^{ref};( chi2^{test} - chi2^{ref} )/ chi2^{ref};FTK Tracks",4000,-2.,2.))); m_map_histo.insert(std::map<std::string, TH1D *>::value_type(istr, new TH1D(istr.c_str(),"#Delta chi2 (test-ref)/ chi2^{ref};( chi2^{test} - chi2^{ref} )/ chi2^{ref};FTK Tracks",4000,-2.,2.)));
} }
else { else {
map_histo_2D.insert(std::map<std::string, TH2D *>::value_type(istr, new TH2D(istr.c_str(),title.c_str(),100,-2.5,2.5,100,-3.2,3.2))); m_map_histo_2D.insert(std::map<std::string, TH2D *>::value_type(istr, new TH2D(istr.c_str(),title.c_str(),100,-2.5,2.5,100,-3.2,3.2)));
} }
} }
} }
...@@ -92,12 +92,12 @@ void CompareFTKEvents::readNTUP_FTKfile() ...@@ -92,12 +92,12 @@ void CompareFTKEvents::readNTUP_FTKfile()
if (!m_fntupftk) std::cout<<"ERROR!! NTUP_FTK file could not be opened: "<<m_NTUP_FTK_file.c_str()<<std::endl; if (!m_fntupftk) std::cout<<"ERROR!! NTUP_FTK file could not be opened: "<<m_NTUP_FTK_file.c_str()<<std::endl;
if (!m_fntupftk->GetListOfKeys()->Contains(m_str_tree_evtinfo.c_str())) std::cout<<"ERROR!! NTUP_FTK file does not contain tree: "<<m_str_tree_evtinfo.c_str()<<std::endl; if (!m_fntupftk->GetListOfKeys()->Contains(m_str_tree_evtinfo.c_str())) std::cout<<"ERROR!! NTUP_FTK file does not contain tree: "<<m_str_tree_evtinfo.c_str()<<std::endl;
m_tevtinfo = (TTree*)m_fntupftk->Get(m_str_tree_evtinfo.c_str()); m_tevtinfo = (TTree*)m_fntupftk->Get(m_str_tree_evtinfo.c_str());
m_tevtinfo->SetBranchAddress("RunNumber",&RN); m_tevtinfo->SetBranchAddress("RunNumber",&m_RN);
m_tevtinfo->SetBranchAddress("EventNumber",&EN); m_tevtinfo->SetBranchAddress("EventNumber",&m_EN);
m_tevtinfo->SetBranchAddress("LB",&LB); m_tevtinfo->SetBranchAddress("LB",&m_LB);
m_tevtinfo->SetBranchAddress("BCID",&BCID); m_tevtinfo->SetBranchAddress("BCID",&m_BCID);
m_tevtinfo->SetBranchAddress("ExtendedLevel1ID",&EL1ID); m_tevtinfo->SetBranchAddress("ExtendedLevel1ID",&m_EL1ID);
m_tevtinfo->SetBranchAddress("Level1TriggerType",&l1TT); m_tevtinfo->SetBranchAddress("Level1TriggerType",&m_l1TT);
if (!m_fntupftk->GetListOfKeys()->Contains(m_str_tree_ftkdata.c_str())) std::cout<<"ERROR!! NTUP_FTK file does not contain tree: "<<m_str_tree_ftkdata.c_str()<<std::endl; if (!m_fntupftk->GetListOfKeys()->Contains(m_str_tree_ftkdata.c_str())) std::cout<<"ERROR!! NTUP_FTK file does not contain tree: "<<m_str_tree_ftkdata.c_str()<<std::endl;
m_theTree=(TTree *)m_fntupftk->Get(m_str_tree_ftkdata.c_str()); m_theTree=(TTree *)m_fntupftk->Get(m_str_tree_ftkdata.c_str());
m_ft=new FTKTrackStream(); m_ft=new FTKTrackStream();
...@@ -112,27 +112,27 @@ StatusCode CompareFTKEvents::decode(uint32_t nTracks, OFFLINE_FRAGMENTS_NAMESPAC ...@@ -112,27 +112,27 @@ StatusCode CompareFTKEvents::decode(uint32_t nTracks, OFFLINE_FRAGMENTS_NAMESPAC
result->reserve(result->size() + nTracks); result->reserve(result->size() + nTracks);
for ( size_t i = 0; i < nTracks; ++i ) { for ( size_t i = 0; i < nTracks; ++i ) {
FTK_RawTrack* track = unpackFTTrack( rodData ); FTK_RawTrack* track = unpackFTTrack( rodData );
rodData += TrackBlobSize; rodData += TRACKBLOBSIZE;
result->push_back(track); result->push_back(track);
} }
return StatusCode::SUCCESS; return StatusCode::SUCCESS;
} }
FTK_RawTrack* CompareFTKEvents::unpackFTTrack( OFFLINE_FRAGMENTS_NAMESPACE::PointerType data) { FTK_RawTrack* CompareFTKEvents::unpackFTTrack( OFFLINE_FRAGMENTS_NAMESPACE::PointerType data) {
FTK_RawTrack* track = new FTK_RawTrack(data[0], data[1], data[2], data[3], data[4], data[5]); // first six words are track params FTK_RawTrack* track = new FTK_RawTrack(data[0], data[1], data[2], data[3], data[4], data[5]); // first six words are track params
data += TrackParamsBlobSize; data += TRACKPARAMSBLOBSIZE;
// get pixel hits // get pixel hits
track->getPixelClusters().resize(NPixLayers); track->getPixelClusters().resize(NPIXLAYERS);
for ( size_t i = 0; i < size_t(NPixLayers); ++i) { for ( size_t i = 0; i < size_t(NPIXLAYERS); ++i) {
size_t offset = PixHitParamsBlobSize*i; size_t offset = PIXHITPARAMSBLOBSIZE*i;
unpackPixCluster(data+offset, track->getPixelCluster(i) ); unpackPixCluster(data+offset, track->getPixelCluster(i) );
} }
data += PixHitParamsBlobSize*NPixLayers; data += PIXHITPARAMSBLOBSIZE*NPIXLAYERS;
// gets SCT hits // gets SCT hits
track->getSCTClusters().resize(NSCTLayers); track->getSCTClusters().resize(NSCTLAYERS);
for ( size_t i = 0; i < size_t(NSCTLayers); ++i) { for ( size_t i = 0; i < size_t(NSCTLAYERS); ++i) {
size_t offset = SCTHitParamsBlobSize*i; size_t offset = SCTHITPARAMSBLOBSIZE*i;
unpackSCTCluster(data+offset, track->getSCTCluster(i) ); unpackSCTCluster(data+offset, track->getSCTCluster(i) );
} }
// no more shifts needed // no more shifts needed
...@@ -174,13 +174,13 @@ void CompareFTKEvents::EventLoop() ...@@ -174,13 +174,13 @@ void CompareFTKEvents::EventLoop()
} }
const uint32_t *rod_data = ROBFragment.rod_data(); const uint32_t *rod_data = ROBFragment.rod_data();
uint32_t rod_ndata = ROBFragment.rod_ndata(); uint32_t rod_ndata = ROBFragment.rod_ndata();
if( ( rod_ndata - trackBlockOffsetStart - trackBlockOffsetEnd ) % TrackBlobSize != 0 ){ if( ( rod_ndata - m_trackBlockOffsetStart - m_trackBlockOffsetEnd ) % TRACKBLOBSIZE != 0 ){
std::cout<<"Error: wrong size of rod"<<rod_ndata<<std::endl; std::cout<<"Error: wrong size of rod"<<rod_ndata<<std::endl;
continue; continue;
} }
niTrks = ( rod_ndata - trackBlockOffsetStart - trackBlockOffsetEnd ) / TrackBlobSize; niTrks = ( rod_ndata - m_trackBlockOffsetStart - m_trackBlockOffsetEnd ) / TRACKBLOBSIZE;
nTrks+=niTrks; nTrks+=niTrks;
rod_data += trackBlockOffsetStart; // moving data pointer to first track block rod_data += m_trackBlockOffsetStart; // moving data pointer to first track block
if (m_verbose){ if (m_verbose){
std::cout<<"test of working rod l1 "<<ROBFragment.rod_lvl1_id()<<" ndata "<<rod_ndata<<std::endl; std::cout<<"test of working rod l1 "<<ROBFragment.rod_lvl1_id()<<" ndata "<<rod_ndata<<std::endl;
} }
...@@ -201,8 +201,8 @@ void CompareFTKEvents::EventLoop() ...@@ -201,8 +201,8 @@ void CompareFTKEvents::EventLoop()
} }
nrob+=1; nrob+=1;
} }
auto search = map_histo.find("nTrk_HW"); auto search = m_map_histo.find("nTrk_HW");
if(search != map_histo.end()) search->second->Fill(nTrks); if(search != m_map_histo.end()) search->second->Fill(nTrks);
else std::cout << "Not found "<<search->first<<std::endl; else std::cout << "Not found "<<search->first<<std::endl;
if (nrob>1){std::cout<<"!!!!!!!!!!SOMETHING WRONG number of robs >1:"<<nrob<<std::endl;} if (nrob>1){std::cout<<"!!!!!!!!!!SOMETHING WRONG number of robs >1:"<<nrob<<std::endl;}
std::cout<<"Event "<<ievent<<" N tracks BS file "<<trkcontainer->size()<<std::endl; std::cout<<"Event "<<ievent<<" N tracks BS file "<<trkcontainer->size()<<std::endl;
...@@ -214,8 +214,8 @@ void CompareFTKEvents::EventLoop() ...@@ -214,8 +214,8 @@ void CompareFTKEvents::EventLoop()
std::vector<const FTK_RawTrack *> ftkNTUPtest; std::vector<const FTK_RawTrack *> ftkNTUPtest;
int NTracksNTUP; int NTracksNTUP;
NTracksNTUP=m_ft->getNTracks(); NTracksNTUP=m_ft->getNTracks();
auto search_SW = map_histo.find("nTrk_SW"); auto search_SW = m_map_histo.find("nTrk_SW");
if(search_SW != map_histo.end()) search_SW->second->Fill(NTracksNTUP); if(search_SW != m_map_histo.end()) search_SW->second->Fill(NTracksNTUP);
else std::cout << "Not found "<<search_SW->first<<std::endl; else std::cout << "Not found "<<search_SW->first<<std::endl;
for(int it=0; it<NTracksNTUP;it++){ for(int it=0; it<NTracksNTUP;it++){
FTKTrack* ftktrk=m_ft->getTrack(it); FTKTrack* ftktrk=m_ft->getTrack(it);
...@@ -236,8 +236,8 @@ void CompareFTKEvents::EventLoop() ...@@ -236,8 +236,8 @@ void CompareFTKEvents::EventLoop()
// start comparing infos // start comparing infos
m_tevtinfo->GetEntry(ievent-1); m_tevtinfo->GetEntry(ievent-1);
if (NTracksNTUP!=nTrks){std::cout<<"different N tracks: "<<NTracksNTUP<<" "<<nTrks<<std::endl;m_allmatched=false;} if (NTracksNTUP!=nTrks){std::cout<<"different N tracks: "<<NTracksNTUP<<" "<<nTrks<<std::endl;m_allmatched=false;}
if (fe.run_no()!=(unsigned int)RN){std::cout<<"different RunNUmber: "<<fe.run_no()<<" "<<RN<<std::endl; m_allmatched=false;} if (fe.run_no()!=(unsigned int)m_RN){std::cout<<"different RunNUmber: "<<fe.run_no()<<" "<<m_RN<<std::endl; m_allmatched=false;}
if (fe.bc_id()!=(unsigned int)BCID){std::cout<<"different BCID: "<<fe.bc_id()<<" "<<BCID<<std::endl;m_allmatched=false;} if (fe.bc_id()!=(unsigned int)m_BCID){std::cout<<"different BCID: "<<fe.bc_id()<<" "<<m_BCID<<std::endl;m_allmatched=false;}
// end comparing infos // end comparing infos
//---------------------------------------------------------------------// //---------------------------------------------------------------------//
// start matching // start matching
...@@ -264,7 +264,7 @@ void CompareFTKEvents::EventLoop() ...@@ -264,7 +264,7 @@ void CompareFTKEvents::EventLoop()
} }
else {if (i%std::min(nTrks,NTracksNTUP)==0 && m_allmatched){std::cout<<"Track Parameters btw BS_FTK and NTUP_FTK are the same: good!"<<std::endl;}} else {if (i%std::min(nTrks,NTracksNTUP)==0 && m_allmatched){std::cout<<"Track Parameters btw BS_FTK and NTUP_FTK are the same: good!"<<std::endl;}}
} }
m_compTrk= new CompareFTKTracks(ftkBSref,ftkNTUPtest, map_histo, map_histo_2D); m_compTrk= new CompareFTKTracks(ftkBSref,ftkNTUPtest, m_map_histo, m_map_histo_2D);
m_compTrk->AssociateTracks(); m_compTrk->AssociateTracks();
m_compTrk->FillHistos(); m_compTrk->FillHistos();
delete data32; delete data32;
...@@ -273,10 +273,10 @@ void CompareFTKEvents::EventLoop() ...@@ -273,10 +273,10 @@ void CompareFTKEvents::EventLoop()
void CompareFTKEvents::WriteHistos(){ void CompareFTKEvents::WriteHistos(){
m_fout= new TFile("./out.histo.root","RECREATE"); m_fout= new TFile("./out.histo.root","RECREATE");
m_fout->cd(); m_fout->cd();
for(auto & imap : map_histo){ for(auto & imap : m_map_histo){
imap.second->Write(); imap.second->Write();
} }
for(auto & imap : map_histo_2D){ for(auto & imap : m_map_histo_2D){
imap.second->Write(); imap.second->Write();
} }
m_fout->Write(); m_fout->Write();
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment