diff --git a/Trigger/TrigAnalysis/TrigInDetAnalysisExample/TrigInDetAnalysisExample/Analysis_Tier0.h b/Trigger/TrigAnalysis/TrigInDetAnalysisExample/TrigInDetAnalysisExample/Analysis_Tier0.h index 44800f5d9661aba88afca3101c9006096c179d2c..f928082d62921af5a85dd7817e7a479735974499 100644 --- a/Trigger/TrigAnalysis/TrigInDetAnalysisExample/TrigInDetAnalysisExample/Analysis_Tier0.h +++ b/Trigger/TrigAnalysis/TrigInDetAnalysisExample/TrigInDetAnalysisExample/Analysis_Tier0.h @@ -151,6 +151,8 @@ class Analysis_Tier0 : public TrackAnalysis { TH1F* h_ntrk_rec; + TH1F* m_h_layer; + TH1F* h_trkpT_rec; TH1F* h_trketa_rec; @@ -169,6 +171,9 @@ class Analysis_Tier0 : public TrackAnalysis { TH1F* h_ntrt_rec; + TH1F* m_h_layer_rec; + + TH1F* h_trkpT_residual; TH1F* h_trkipT_residual; TH1F* h_trketa_residual; diff --git a/Trigger/TrigAnalysis/TrigInDetAnalysisExample/src/Analysis_Tier0.cxx b/Trigger/TrigAnalysis/TrigInDetAnalysisExample/src/Analysis_Tier0.cxx index 3c02ec9a8a837b3269a9bf7f080b3815d7af3037..77513beaa64ee58edcf4112a388409d695cedddf 100644 --- a/Trigger/TrigAnalysis/TrigInDetAnalysisExample/src/Analysis_Tier0.cxx +++ b/Trigger/TrigAnalysis/TrigInDetAnalysisExample/src/Analysis_Tier0.cxx @@ -226,7 +226,11 @@ void Analysis_Tier0::initialise() { addHistogram( h_nsihits_lb ); addHistogram( h_nsihits_lb_rec ); + m_h_layer_rec = new TH1F("layer_rec" , "hit layers", 32, -0.5, 31.5 ); + m_h_layer = new TH1F("layer", "hit layers", 32, -0.5, 31.5 ); + addHistogram(m_h_layer); + addHistogram(m_h_layer_rec); @@ -476,7 +480,10 @@ void Analysis_Tier0::execute(const std::vector<TIDA::Track*>& referenceTracks, h_chain->Fill(2.5); - + for ( size_t ilayer=0 ; ilayer<32 ; ilayer++ ) { + if ( (*reference)->hitPattern()&(1<<ilayer) ) m_h_layer_rec->Fill( ilayer ); + } + if (test) { h_chain->Fill(3.5); @@ -504,6 +511,10 @@ void Analysis_Tier0::execute(const std::vector<TIDA::Track*>& referenceTracks, h_trkvtx_y_lb->Fill( event()->lumi_block(), beamTesty() ); h_trkvtx_z_lb->Fill( event()->lumi_block(), beamTestz() ); + for ( size_t ilayer=0 ; ilayer<32 ; ilayer++ ) { + if ( test->hitPattern()&(1<<ilayer) ) m_h_layer_rec->Fill( ilayer ); + } + // std::cout << "SUTT beam x " << beamTestx() << " " << "\tx " << beamTesty() << " " << "\ty " << beamTestz() << std::endl; #if 0