Skip to content
Snippets Groups Projects
Commit 55311859 authored by Will Leight's avatar Will Leight
Browse files

Added protection against trying to retrieve RPC sector logic container when... 

Added protection against trying to retrieve RPC sector logic container when running on MC for RpcLv1RawDataEfficiency and RpcLv1RawDataSectorLogic.
This requires a configurable flag to determine whether the code is running on MC or not: the configuration is updated accordingly.


Former-commit-id: 69d11d74
parent 1534c0b8
No related branches found
No related tags found
No related merge requests found
Hide whitespace changes
Inline Side-by-side
Showing
with 178 additions and 165 deletions
MuonSpectrometer/MuonValidation/MuonDQA/MuonRawDataMonitoring/RpcRawDataMonitoring/RpcRawDataMonitoring/RpcLv1RawDataEfficiency.h
+ 2
− 0
View file @ 55311859
...@@ -201,6 +201,8 @@ class RpcLv1RawDataEfficiency: public ManagedMonitorToolBase { ...@@ -201,6 +201,8 @@ class RpcLv1RawDataEfficiency: public ManagedMonitorToolBase {
int m_lumiblock; int m_lumiblock;
int m_BCID; int m_BCID;
bool m_isMC;
SG::ReadHandleKey<Muon::RpcCoinDataContainer> m_rpcCoinKey{this,"RpcCoinKey","RPC_triggerHits","RPC coincidences"}; SG::ReadHandleKey<Muon::RpcCoinDataContainer> m_rpcCoinKey{this,"RpcCoinKey","RPC_triggerHits","RPC coincidences"};
SG::ReadHandleKey<xAOD::EventInfo> m_eventInfo{this,"EventInfo","EventInfo","EventInfo"}; SG::ReadHandleKey<xAOD::EventInfo> m_eventInfo{this,"EventInfo","EventInfo","EventInfo"};
SG::ReadHandleKey<RpcSectorLogicContainer> m_sectorLogicContainerKey{this,"RPCSec","RPC_SECTORLOGIC","RPC sector logic"}; SG::ReadHandleKey<RpcSectorLogicContainer> m_sectorLogicContainerKey{this,"RPCSec","RPC_SECTORLOGIC","RPC sector logic"};
......
MuonSpectrometer/MuonValidation/MuonDQA/MuonRawDataMonitoring/RpcRawDataMonitoring/RpcRawDataMonitoring/RpcLv1RawDataSectorLogic.h
+ 1
− 0
View file @ 55311859
...@@ -78,6 +78,7 @@ class RpcLv1RawDataSectorLogic: public ManagedMonitorToolBase { ...@@ -78,6 +78,7 @@ class RpcLv1RawDataSectorLogic: public ManagedMonitorToolBase {
int m_rpc_event_inarea; int m_rpc_event_inarea;
bool m_lumiblockhist ; bool m_lumiblockhist ;
bool m_isMC;
// Declare Histograms // Declare Histograms
TH1* m_rpclv1_Hits_per_TriggerSector; TH1* m_rpclv1_Hits_per_TriggerSector;
......
MuonSpectrometer/MuonValidation/MuonDQA/MuonRawDataMonitoring/RpcRawDataMonitoring/share/RpcRawBS_MonitoringOptions.py
+ 2
− 1
View file @ 55311859
...@@ -37,7 +37,8 @@ rpcLv1SLRawMonMan = AthenaMonManager(name="RpcLv1SLRawMonManager", ...@@ -37,7 +37,8 @@ rpcLv1SLRawMonMan = AthenaMonManager(name="RpcLv1SLRawMonManager",
OutputLevel = muonOutputLevel) OutputLevel = muonOutputLevel)
from RpcRawDataMonitoring.RpcRawDataMonitoringConf import RpcLv1RawDataSectorLogic from RpcRawDataMonitoring.RpcRawDataMonitoringConf import RpcLv1RawDataSectorLogic
rpcLV1RawDataSectorLogic = RpcLv1RawDataSectorLogic(name='rpcLV1RawDataSectorLogic',OutputLevel = OutputMessageLevel) rpcLV1RawDataSectorLogic = RpcLv1RawDataSectorLogic(name='rpcLV1RawDataSectorLogic',OutputLevel = OutputMessageLevel)
if globalflags.DataSource() != 'data':
rpcLV1RawDataSectorLogic.isMC = True
ToolSvc += rpcLV1RawDataSectorLogic ToolSvc += rpcLV1RawDataSectorLogic
rpcLv1SLRawMonMan.AthenaMonTools += [ rpcLV1RawDataSectorLogic ] rpcLv1SLRawMonMan.AthenaMonTools += [ rpcLV1RawDataSectorLogic ]
topSequence += rpcLv1SLRawMonMan topSequence += rpcLv1SLRawMonMan
......
MuonSpectrometer/MuonValidation/MuonDQA/MuonRawDataMonitoring/RpcRawDataMonitoring/share/RpcRawESD_MonitoringOptions.py
+ 2
− 0
View file @ 55311859
...@@ -93,6 +93,8 @@ if not hasattr(topSequence,"RpcClusterBuilderPRD"): ...@@ -93,6 +93,8 @@ if not hasattr(topSequence,"RpcClusterBuilderPRD"):
topSequence += RpcClusterBuilderPRD(ClusterTimeSpread = 10000, topSequence += RpcClusterBuilderPRD(ClusterTimeSpread = 10000,
OutputLevel = FATAL) OutputLevel = FATAL)
rpcLv1RawDataEfficiency = RpcLv1RawDataEfficiency(name='rpcLv1RawDataEfficiency') rpcLv1RawDataEfficiency = RpcLv1RawDataEfficiency(name='rpcLv1RawDataEfficiency')
if globalflags.DataSource() != 'data':
rpcLv1RawDataEfficiency.isMC = True
ToolSvc += rpcLv1RawDataEfficiency ToolSvc += rpcLv1RawDataEfficiency
rpcLv1RawEfficiencyMonMan.AthenaMonTools += [ rpcLv1RawDataEfficiency ] rpcLv1RawEfficiencyMonMan.AthenaMonTools += [ rpcLv1RawDataEfficiency ]
topSequence += rpcLv1RawEfficiencyMonMan topSequence += rpcLv1RawEfficiencyMonMan
......
MuonSpectrometer/MuonValidation/MuonDQA/MuonRawDataMonitoring/RpcRawDataMonitoring/src/RpcLv1RawDataEfficiency.cxx
+ 31
− 28
View file @ 55311859
...@@ -39,7 +39,9 @@ RpcLv1RawDataEfficiency::RpcLv1RawDataEfficiency( const std::string & type, ...@@ -39,7 +39,9 @@ RpcLv1RawDataEfficiency::RpcLv1RawDataEfficiency( const std::string & type,
const std::string & name, const std::string & name,
const IInterface* parent ) const IInterface* parent )
:ManagedMonitorToolBase( type, name, parent ) :ManagedMonitorToolBase( type, name, parent )
{} {
declareProperty("isMC" , m_isMC = false );
}
//================================================================================================================================ //================================================================================================================================
RpcLv1RawDataEfficiency::~RpcLv1RawDataEfficiency() RpcLv1RawDataEfficiency::~RpcLv1RawDataEfficiency()
...@@ -84,7 +86,7 @@ StatusCode RpcLv1RawDataEfficiency::initialize() ...@@ -84,7 +86,7 @@ StatusCode RpcLv1RawDataEfficiency::initialize()
ATH_CHECK(m_rpcCoinKey.initialize()); ATH_CHECK(m_rpcCoinKey.initialize());
ATH_CHECK(m_eventInfo.initialize()); ATH_CHECK(m_eventInfo.initialize());
ATH_CHECK(m_sectorLogicContainerKey.initialize()); ATH_CHECK(m_sectorLogicContainerKey.initialize(m_isMC));
// Ignore the checking code // Ignore the checking code
ManagedMonitorToolBase::initialize().ignore(); ManagedMonitorToolBase::initialize().ignore();
...@@ -215,34 +217,35 @@ StatusCode RpcLv1RawDataEfficiency::fillHistograms( ) ...@@ -215,34 +217,35 @@ StatusCode RpcLv1RawDataEfficiency::fillHistograms( )
// == Filling the Histograms // == Filling the Histograms
//--------------------Sector Hits--------------------------------- //--------------------Sector Hits---------------------------------
// Retrieve the Sector Logic container if(!m_isMC){
SG::ReadHandle<RpcSectorLogicContainer> sectorLogicContainer(m_sectorLogicContainerKey); // Retrieve the Sector Logic container
SG::ReadHandle<RpcSectorLogicContainer> sectorLogicContainer(m_sectorLogicContainerKey);
RpcSectorLogicContainer::const_iterator it = sectorLogicContainer -> begin(); RpcSectorLogicContainer::const_iterator it = sectorLogicContainer -> begin();
for ( ; it != sectorLogicContainer -> end() ; ++it ) for ( ; it != sectorLogicContainer -> end() ; ++it )
{ {
int i_sectorid = (*it)->sectorId(); int i_sectorid = (*it)->sectorId();
// Loop over the trigger hits of each sector // Loop over the trigger hits of each sector
RpcSectorLogic::const_iterator ithit = (*it) -> begin(); RpcSectorLogic::const_iterator ithit = (*it) -> begin();
for ( ; ithit != (*it) -> end() ; ++ithit ) for ( ; ithit != (*it) -> end() ; ++ithit )
{ {
bool b_isInput = (*ithit) -> isInput(); bool b_isInput = (*ithit) -> isInput();
int i_ptid = (*ithit) -> ptId(); int i_ptid = (*ithit) -> ptId();
if (b_isInput == true) // fill histograms only if there was a trigger hit if (b_isInput == true) // fill histograms only if there was a trigger hit
{ {
for( int iThresholdIndex=0; iThresholdIndex < 6; iThresholdIndex++) { for( int iThresholdIndex=0; iThresholdIndex < 6; iThresholdIndex++) {
if (i_ptid == (iThresholdIndex + 1)) if (i_ptid == (iThresholdIndex + 1))
{ {
m_rpclv1_sectorhits_all[iThresholdIndex]->Fill(m_lumiblock,i_sectorid); m_rpclv1_sectorhits_all[iThresholdIndex]->Fill(m_lumiblock,i_sectorid);
m_rpclv1_sectorhits_A[iThresholdIndex]->Fill(m_lumiblock,i_sectorid); m_rpclv1_sectorhits_A[iThresholdIndex]->Fill(m_lumiblock,i_sectorid);
m_rpclv1_sectorhits_C[iThresholdIndex]->Fill(m_lumiblock,i_sectorid); m_rpclv1_sectorhits_C[iThresholdIndex]->Fill(m_lumiblock,i_sectorid);
} }
}
} }
} }
} }
} }
return StatusCode::SUCCESS; return StatusCode::SUCCESS;
} }
......
MuonSpectrometer/MuonValidation/MuonDQA/MuonRawDataMonitoring/RpcRawDataMonitoring/src/RpcLv1RawDataSectorLogic.cxx
+ 140
− 136
View file @ 55311859
...@@ -37,6 +37,7 @@ RpcLv1RawDataSectorLogic::RpcLv1RawDataSectorLogic( const std::string & type, ...@@ -37,6 +37,7 @@ RpcLv1RawDataSectorLogic::RpcLv1RawDataSectorLogic( const std::string & type,
{ {
// Declare the properties // Declare the properties
declareProperty("LumiblockHist" , m_lumiblockhist = false ); declareProperty("LumiblockHist" , m_lumiblockhist = false );
declareProperty("isMC" , m_isMC = false );
} }
//*********************************************************************************************************************** //***********************************************************************************************************************
...@@ -68,7 +69,7 @@ StatusCode RpcLv1RawDataSectorLogic::initialize() ...@@ -68,7 +69,7 @@ StatusCode RpcLv1RawDataSectorLogic::initialize()
m_Tower_in = 0 ; m_Tower_in = 0 ;
m_in_sectorid = 0 ; m_in_sectorid = 0 ;
ATH_CHECK(m_sectorLogicContainerKey.initialize()); ATH_CHECK(m_sectorLogicContainerKey.initialize(m_isMC));
ATH_CHECK(m_eventInfo.initialize()); ATH_CHECK(m_eventInfo.initialize());
// Ignore the checking code // Ignore the checking code
...@@ -175,150 +176,153 @@ StatusCode RpcLv1RawDataSectorLogic::fillHistograms( ) ...@@ -175,150 +176,153 @@ StatusCode RpcLv1RawDataSectorLogic::fillHistograms( )
m_nTriggerHits = 0; m_nTriggerHits = 0;
// Retrieve the Sector Logic container if(!m_isMC){
SG::ReadHandle<RpcSectorLogicContainer> sectorLogicContainer(m_sectorLogicContainerKey); // Retrieve the Sector Logic container
/////////////////////////////////////////// SG::ReadHandle<RpcSectorLogicContainer> sectorLogicContainer(m_sectorLogicContainerKey);
// Loop over the Sector Logic containers // ///////////////////////////////////////////
/////////////////////////////////////////// // Loop over the Sector Logic containers //
RpcSectorLogicContainer::const_iterator it = sectorLogicContainer -> begin(); ///////////////////////////////////////////
for ( ; it != sectorLogicContainer -> end() ; ++it ) RpcSectorLogicContainer::const_iterator it = sectorLogicContainer -> begin();
{ for ( ; it != sectorLogicContainer -> end() ; ++it )
int i_sectorid = (*it)->sectorId(); {
m_nTriggerHitsperSector = 0; int i_sectorid = (*it)->sectorId();
m_nTriggerHitsperSector = 0;
// Loop over the trigger hits of each sector
RpcSectorLogic::const_iterator ithit = (*it) -> begin(); // Loop over the trigger hits of each sector
for ( ; ithit != (*it) -> end() ; ++ithit ) RpcSectorLogic::const_iterator ithit = (*it) -> begin();
{ for ( ; ithit != (*it) -> end() ; ++ithit )
// from RpcSLTriggerHit {
bool b_isInput = (*ithit) -> isInput(); // from RpcSLTriggerHit
int i_rowinBcid = (*ithit) -> rowinBcid();//readout window BCid bool b_isInput = (*ithit) -> isInput();
int i_padid = (*ithit) -> padId();//tower int i_rowinBcid = (*ithit) -> rowinBcid();//readout window BCid
int i_ptid = (*ithit) -> ptId();//threshold int i_padid = (*ithit) -> padId();//tower
int i_roi = (*ithit) -> roi();//region of interest int i_ptid = (*ithit) -> ptId();//threshold
int i_triggerBcid = (*ithit) -> triggerBcid(); int i_roi = (*ithit) -> roi();//region of interest
int i_triggerBcid = (*ithit) -> triggerBcid();
m_Diff_triggerBCid = 0;
m_in_triggerBCid = 0; m_Diff_triggerBCid = 0;
m_out_triggerBCid = 0; m_in_triggerBCid = 0;
m_in_rowinBCid = 0; m_out_triggerBCid = 0;
m_out_rowinBCid = 0; m_in_rowinBCid = 0;
m_Tower_in = -1; m_out_rowinBCid = 0;
m_Tower_in = -1;
if (b_isInput == true){
m_nTriggerHits++; if (b_isInput == true){
m_nTriggerHitsperSector++; m_nTriggerHits++;
m_nTriggerHitsperSector++;
//Fill hits per trigger sector histogram //Fill hits per trigger sector histogram
m_rpclv1_Hits_per_TriggerSector -> Fill(float(i_sectorid)); m_rpclv1_Hits_per_TriggerSector -> Fill(float(i_sectorid));
//per lumi block //per lumi block
if(m_lumiblockhist)m_rpclv1_Hits_per_TriggerSector_LB -> Fill(float(i_sectorid)); if(m_lumiblockhist)m_rpclv1_Hits_per_TriggerSector_LB -> Fill(float(i_sectorid));
//Fill Pad vs. tigger sector histograms
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
if (i_ptid == 1){ //Fill Pad vs. tigger sector histograms
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt1 -> Fill(float(-(i_padid + 1)), float(i_sectorid)); if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt1 -> Fill(float((i_padid + 1)), float(i_sectorid - 32)); else m_rpclv1_TriggerSector_vs_Pad -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
}
if (i_ptid == 2){
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt2 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt2 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
}
if (i_ptid == 3){
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt3 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt3 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
}
if (i_ptid == 4){
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt4 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt4 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
}
if (i_ptid == 5){
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt5 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt5 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
}
if (i_ptid == 6){
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt6 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt6 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
}
//Fill trigger sector vs. rowin BCid histogram if (i_ptid == 1){
m_rpclv1_rowinBCid_vs_TriggerSector -> Fill(float(i_sectorid), float(i_rowinBcid)); if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt1 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt1 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
//Fill Threshold vs. Pad(Tower) histogram }
m_rpclv1_ptid_vs_Tower -> Fill(float(i_padid), float(i_ptid)); if (i_ptid == 2){
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt2 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
//Trigger BCid in else m_rpclv1_TriggerSector_vs_Pad_Pt2 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
m_in_triggerBCid = i_triggerBcid; }
//ReadOutWindow BCid in if (i_ptid == 3){
m_in_rowinBCid = i_rowinBcid; if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt3 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt3 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
//Tower in }
m_Tower_in = i_padid; if (i_ptid == 4){
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt4 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt4 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
}
if (i_ptid == 5){
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt5 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt5 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
}
if (i_ptid == 6){
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_Pt6 -> Fill(float(-(i_padid + 1)), float(i_sectorid));
else m_rpclv1_TriggerSector_vs_Pad_Pt6 -> Fill(float((i_padid + 1)), float(i_sectorid - 32));
}
m_Tower_out = -1; //Fill trigger sector vs. rowin BCid histogram
// Loop over the trigger hits of each sector m_rpclv1_rowinBCid_vs_TriggerSector -> Fill(float(i_sectorid), float(i_rowinBcid));
RpcSectorLogic::const_iterator ithiti = (*it) -> begin();
for ( ; ithiti != (*it) -> end() ; ++ithiti ) //Fill Threshold vs. Pad(Tower) histogram
{ m_rpclv1_ptid_vs_Tower -> Fill(float(i_padid), float(i_ptid));
bool b_isInput1 = (*ithiti) -> isInput();
int i_roi1 = (*ithiti) -> roi();//region of interest //Trigger BCid in
int i_triggerBcid1 = (*ithiti) -> triggerBcid(); m_in_triggerBCid = i_triggerBcid;
int i_rowinBcid1 = (*ithiti) -> rowinBcid(); //ReadOutWindow BCid in
if (b_isInput1 == false){ m_in_rowinBCid = i_rowinBcid;
//Tower out calculated from ROI out //Tower in
if ( i_roi1 < 4 ) m_Tower_out = 0; m_Tower_in = i_padid;
if ( i_roi1 >= 4 && i_roi < 8 ) m_Tower_out = 1;
if ( i_roi1 >= 8 && i_roi < 12 ) m_Tower_out = 2; m_Tower_out = -1;
if ( i_roi1 >= 12 && i_roi < 16 ) m_Tower_out = 3; // Loop over the trigger hits of each sector
if ( i_roi1 >= 16 && i_roi < 20 ) m_Tower_out = 4; RpcSectorLogic::const_iterator ithiti = (*it) -> begin();
if ( i_roi1 >= 20 && i_roi < 24 ) m_Tower_out = 5; for ( ; ithiti != (*it) -> end() ; ++ithiti )
if ( i_roi1 >= 24 ) m_Tower_out = 6; {
bool b_isInput1 = (*ithiti) -> isInput();
int i_roi1 = (*ithiti) -> roi();//region of interest
int i_triggerBcid1 = (*ithiti) -> triggerBcid();
int i_rowinBcid1 = (*ithiti) -> rowinBcid();
if (b_isInput1 == false){
//Tower out calculated from ROI out
if ( i_roi1 < 4 ) m_Tower_out = 0;
if ( i_roi1 >= 4 && i_roi < 8 ) m_Tower_out = 1;
if ( i_roi1 >= 8 && i_roi < 12 ) m_Tower_out = 2;
if ( i_roi1 >= 12 && i_roi < 16 ) m_Tower_out = 3;
if ( i_roi1 >= 16 && i_roi < 20 ) m_Tower_out = 4;
if ( i_roi1 >= 20 && i_roi < 24 ) m_Tower_out = 5;
if ( i_roi1 >= 24 ) m_Tower_out = 6;
//Trigger BCid out
m_out_triggerBCid = i_triggerBcid1;
//ReadOutWindow BCid in
m_out_rowinBCid = i_rowinBcid1;
//Trigger BCid out //Calculate Timing difference in trigger BCid
m_out_triggerBCid = i_triggerBcid1; if ((m_Tower_out == m_Tower_in) && (m_out_rowinBCid == m_in_rowinBCid)){
//ReadOutWindow BCid in if(m_in_triggerBCid - m_out_triggerBCid >= 0) m_Diff_triggerBCid = m_in_triggerBCid - m_out_triggerBCid;
m_out_rowinBCid = i_rowinBcid1; else m_Diff_triggerBCid = (m_in_triggerBCid - m_out_triggerBCid) + 8;
}
//Calculate Timing difference in trigger BCid }
if ((m_Tower_out == m_Tower_in) && (m_out_rowinBCid == m_in_rowinBCid)){
if(m_in_triggerBCid - m_out_triggerBCid >= 0) m_Diff_triggerBCid = m_in_triggerBCid - m_out_triggerBCid;
else m_Diff_triggerBCid = (m_in_triggerBCid - m_out_triggerBCid) + 8;
}
}
} }
} }
//Fill in/out dependent histograms & histograms with timing difference only for IN //Fill in/out dependent histograms & histograms with timing difference only for IN
if (b_isInput == true){ if (b_isInput == true){
m_rpclv1_triggerBCid_inout_vs_Tower -> Fill(float(m_Tower_out), float(m_Diff_triggerBCid)); m_rpclv1_triggerBCid_inout_vs_Tower -> Fill(float(m_Tower_out), float(m_Diff_triggerBCid));
// if (i_padid == 0){
m_rpclv1_triggerBCid_inout_vs_TriggerSector -> Fill(float(i_sectorid), float(m_Diff_triggerBCid));
m_rpclv1_triggerBCid_inout -> Fill(float(m_Diff_triggerBCid));
//per lumi block
if(m_lumiblockhist)m_rpclv1_triggerBCid_inout_LB -> Fill(float(m_Diff_triggerBCid));
// }
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_triggerBCid_inout -> Fill(float(-(i_padid + 1)), float(i_sectorid), float(m_Diff_triggerBCid));
else m_rpclv1_TriggerSector_vs_Pad_triggerBCid_inout -> Fill(float((i_padid + 1)), float(i_sectorid - 32), float(m_Diff_triggerBCid));
}
} // End Loop over the trigger hits of each sector // if (i_padid == 0){
m_rpclv1_triggerBCid_inout_vs_TriggerSector -> Fill(float(i_sectorid), float(m_Diff_triggerBCid));
m_rpclv1_triggerBCid_inout -> Fill(float(m_Diff_triggerBCid));
//per lumi block
if(m_lumiblockhist)m_rpclv1_triggerBCid_inout_LB -> Fill(float(m_Diff_triggerBCid));
// }
if (i_sectorid < 32) m_rpclv1_TriggerSector_vs_Pad_triggerBCid_inout -> Fill(float(-(i_padid + 1)), float(i_sectorid), float(m_Diff_triggerBCid));
else m_rpclv1_TriggerSector_vs_Pad_triggerBCid_inout -> Fill(float((i_padid + 1)), float(i_sectorid - 32), float(m_Diff_triggerBCid));
}
} // End Loop over the trigger hits of each sector
// Fill the trigger hits per event per sector histogram here // Fill the trigger hits per event per sector histogram here
if (m_nTriggerHitsperSector > 0) { if (m_nTriggerHitsperSector > 0) {
m_rpclv1_TriggerHitsperEventperTriggerSector -> Fill(float(i_sectorid), float(m_nTriggerHitsperSector)); m_rpclv1_TriggerHitsperEventperTriggerSector -> Fill(float(i_sectorid), float(m_nTriggerHitsperSector));
//per lumi block //per lumi block
if(m_lumiblockhist)m_rpclv1_TriggerHitsperEventperTriggerSector_LB -> Fill(float(i_sectorid), float(m_nTriggerHitsperSector)); if(m_lumiblockhist)m_rpclv1_TriggerHitsperEventperTriggerSector_LB -> Fill(float(i_sectorid), float(m_nTriggerHitsperSector));
} }
} // End Loop over the Sector Logic containers } // End Loop over the Sector Logic containers
}
// Fill the trigger hits per event histogram here // Fill the trigger hits per event histogram here
m_rpclv1_TriggerHitsperEvent -> Fill(float(m_nTriggerHits)); m_rpclv1_TriggerHitsperEvent -> Fill(float(m_nTriggerHits));
......
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