MDTMonitorAlgorithm.py

Go to the documentation of this file.

#
#Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
#
 
from AthenaConfiguration.ComponentFactory import CompFactory
from AthenaConfiguration.Enums import LHCPeriod
from MdtRawDataMonitoring.MdtMonUtils import getMDTLabel, getMDTLabelx
from MdtRawDataMonitoring.MDTTubeMax import tubeMax, tubeMax_smdt
from MdtRawDataMonitoring.MDTChambers import mdtBA,mdtBC,mdtEA,mdtEC
 
 
def MdtMonitoringConfig(inputFlags):
    '''Function to configures some algorithms in the monitoring system.'''
 
    
    from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
    result = ComponentAccumulator()
 
    # Make sure muon geometry is configured
    from MuonConfig.MuonGeometryConfig import MuonGeoModelCfg
    result.merge(MuonGeoModelCfg(inputFlags))
 
    # Fixup for R21 ESD
    from MuonConfig.MuonSegmentNameFixConfig import MuonSegmentNameFixCfg
    result.merge(MuonSegmentNameFixCfg(inputFlags))
 
    # The following class will make a sequence, configure algorithms, and link
    # them to GenericMonitoringTools
    from AthenaMonitoring import AthMonitorCfgHelper
    helper = AthMonitorCfgHelper(inputFlags,'MdtAthMonitorCfg')
    
    
    mdtMonAlg = helper.addAlgorithm(CompFactory.MdtRawDataMonAlg,'MdtMonAlg')
    mdtMonAlg.DoMdtEsd = True
    mdtMonAlg.DoChamberHist=True
    mdtMonAlg.maskNoisyTubes = True
    mdtMonAlg.ADCCut = 80.
    mdtMonAlg.do_mdtChamberHits=True
    mdtMonAlg.do_mdttdccut_sector=True
    mdtMonAlg.do_mdtchamberstatphislice=True
    if (inputFlags.GeoModel.Run == LHCPeriod.Run3):
        mdtMonAlg.do_Run3Geometry=True
    
    
    if inputFlags.Input.Release and inputFlags.Input.Release.find("21.") != -1:
        mdtMonAlg.SegmentKey = [ "TrackMuonSegments"]
        
    if not inputFlags.DQ.triggerDataAvailable:
        mdtMonAlg.L1RoiKey=''
    # Add a generic monitoring tool (a "group" in old language). The returned 
    # object here is the standard GenericMonitoringTool.
    mdtGroup = helper.addGroup(mdtMonAlg,'MdtMonitor','Muon/MuonRawDataMonitoring/MDT/')
 
    # Add a GMT for the other example monitor algorithm
    
 
    mdtGroup.defineHistogram('run3geo;GeometryFlag',  type='TH1F', cutmask='firstEvent',
                             title='GeometryFlag;flag;',
                             path='Overview',   xbins=5, xmin=0., xmax=5)
 
    mdtGroup.defineHistogram('mdt_tube_z_barrel,mdt_tube_perp_barrel;Number_of_BarrelMDTHits_inRZView_Global_ADCCut', type='TH2F',
                            title='Number_of_BarrelMDTHits_inRZView_Global_ADCCut;MDT-GlobalZ(mm);MDT-GlobalR(mm)',
                            path='Overview/Expert/Geometry',  xbins=250,xmin=-25000.,xmax=25000., ybins=120, ymin=0., ymax=12000.   
                            )
 
    mdtGroup.defineHistogram('mdt_tube_x_barrel,mdt_tube_y_barrel;Number_of_BarrelMDTHits_inXYView_Global_ADCCut',   type='TH2F',
                            title='Number_of_BarrelMDTHits_inXYView_Global_ADCCut;MDT-GlobalX(mm);MDT-GlobalY(mm)',
                            path='Overview/Expert/Geometry',   xbins=150, xmin=-15000., xmax=15000., ybins=150, ymin=-15000., ymax=15000.
                            )
 
    mdtGroup.defineHistogram('mdt_tube_z_ovl,mdt_tube_perp_ovl;Number_of_OverlapMDTHits_inRZView_Global_ADCCut',  type='TH2F',
                            title='Number_of_OverlapMDTHits_inRZView_Global_ADCCut;MDT-GlobalZ(mm);MDT-GlobalR(mm)',
                            path='Overview/Expert/Geometry',  xbins=250,xmin=-25000.,xmax=25000., ybins=120, ymin=0., ymax=12000.   
                            )
 
    mdtGroup.defineHistogram('mdt_tube_x_ovl,mdt_tube_y_ovl;Number_of_OverlapMDTHits_inXYView_Global_ADCCut',   type='TH2F',
                            title='Number_of_OverlapMDTHits_inXYView_Global_ADCCut;MDT-GlobalX(mm);MDT-GlobalY(mm)',
                            path='Overview/Expert/Geometry',   xbins=150, xmin=-15000., xmax=15000., ybins=150, ymin=-15000., ymax=15000.
                            )
 
    mdtGroup.defineHistogram('mdt_tube_z_endcap,mdt_tube_perp_endcap;Number_of_EndcapMDTHits_inRZView_Global_ADCCut', type='TH2F',
                            title='Number_of_EndcapMDTHits_inRZView_Global_ADCCut;MDT-GlobalZ(mm);MDT-GlobalR(mm)',
                            path='Overview/Expert/Geometry',  xbins=250,xmin=-25000.,xmax=25000., ybins=120, ymin=0., ymax=12000.   
                            )
 
    mdtGroup.defineHistogram('mdt_tube_x_endcap,mdt_tube_y_endcap;Number_of_EndcapMDTHits_inXYView_Global_ADCCut',  type='TH2F',
                            title='Number_of_EndcapMDTHits_inXYView_Global_ADCCut;MDT-GlobalX(mm);MDT-GlobalY(mm)',
                            path='Overview/Expert/Geometry',   xbins=150, xmin=-15000., xmax=15000., ybins=150, ymin=-15000., ymax=15000.
                            )
 
    mdtGroup.defineHistogram('tdc_mon,adc_mon;Overall_TDCADC_spectrum',    type='TH2F',
                            title='Overall_TDCADC_spectrum;[nsec];[adc counts]',
                            path='Overview',   xbins=50, xmin=0., xmax=2000., ybins=20, ymin=0., ymax=400.
                            )
 
    labels_allSec=list()
    for phi in range(0, 16):
        sectPhi=str(phi+1)
        if phi+1 < 10: 
            sectPhi="0"+sectPhi
 
        if mdtMonAlg.do_Run3Geometry :
            thisLabelx=getMDTLabelx("run3_labels_sectorPhi"+sectPhi)
        else:
            thisLabelx=getMDTLabelx("labels_sectorPhi"+sectPhi)            
        labels_allSec=labels_allSec+thisLabelx
        title="MDTHitsOccup_ADCCut_Sector"+sectPhi
        maxbin=len(thisLabelx)
        nbins=maxbin
        mdtGroup.defineHistogram('hits_phi_'+sectPhi+';'+title,  type='TH1F',
                                 title=title+';;Number of hits',
                                 path='Overview/HitOccupancies_Sectors',   xbins=nbins, xmin=0., xmax=maxbin, xlabels=thisLabelx
                             )
    maxbin=len(labels_allSec)
    nbins=maxbin
    mdtGroup.defineHistogram('hits_allphi;Number_of_MDT_hits_per_chamber_ADCCut',  type='TH1F',
                             title='Number_of_MDT_hits_per_chamber_ADCCut;;Number of hits',
                            path='Overview',   xbins=nbins, xmin=0., xmax=maxbin, xlabels=labels_allSec
                             )
 
    mdtGroup.defineHistogram('adc_mon_nosel;Overall_ADC_spectrum',  type='TH1F',
                            title='Overall_ADC_spectrum;[adc counts];Number of Entries',
                            path='Overview',   xbins=100, xmin=0., xmax=400.
                            )
 
    mdtGroup.defineHistogram('adc_mon_nosel;Overall_ADC_spectrum_NoiseBursts',  type='TH1F',
                            cutmask='noiseBurst',
                            title='Overall_ADC_spectrum_NoiseBursts;[adc counts];Number of Entries',
                            path='Overview',   xbins=100, xmin=0., xmax=400.
                            )
 
    
    mdtGroup.defineHistogram('adc_segs_overall_mon;Overall_ADC_OnSegm_spectrum',  type='TH1F',
                             title='Overall_ADC_spectrum;[adc counts];Number of Entries',
                             path='Overview',   xbins=100, xmin=0., xmax=400.)
 
    
    mdtGroup.defineHistogram('tdc_mon_nosel,adc_mon_nosel;Overall_TDCADC_spectrum_NoiseBursts',    type='TH2F',
                            cutmask='noiseBurst',
                            title='Overall_TDCADC_spectrum_NoiseBursts;[nsec];[adc counts]',
                            path='Overview',   xbins=50, xmin=0., xmax=2000., ybins=20, ymin=0., ymax=400.
                            )
 
    mdtGroup.defineHistogram('tdc_mon_nosel;Overall_TDC_spectrum_NoiseBursts',  type='TH1F',
                            cutmask='noiseBurst',
                            title='Overall_TDC_spectrum_NoiseBursts;[nsec];Number of Entries',
                            path='Overview',   xbins=120, xmin=0., xmax=2000.
                            )
 
    
    mdtGroup.defineHistogram('tdc_mon_noiseBurst_adcCut;Overall_TDC_spectrum_NoiseBursts_ADCCut',  type='TH1F',
                            title='Overall_TDC_spectrum_NoiseBursts_ADCCut;[nsec];Number of Entries',
                            path='Overview',   xbins=120, xmin=0., xmax=2000.
                            )
 
    mdtGroup.defineHistogram('tdc_mon_adcCut;Overall_TDC_ADCCUT_spectrum',  type='TH1F',
                            title='Overall_TDC_spectrum_ADCCut;[nsec];Number of Entries',
                            path='Overview',   xbins=120, xmin=0., xmax=2000.
                            )
 
    mdtGroup.defineHistogram('tdc_segs_overall_mon;Overall_TDC_onSegm_ADCCut_spectrum',  type='TH1F',
                            title='Overall_TDC_onSegm_ADCCut_spectrum;[nsec];Number of Entries',
                            path='Overview',   xbins=120, xmin=0., xmax=2000.
                            )
 
    mdtGroup.defineHistogram('nHighOccChambers_mon;number_of_Chambers_with_high_occupancy_per_event',  type='TH1F',
                             title='number_of_Chambers_with_high_occupancy_per_event;[Number_of_MDT_chambers_with_high_occupancy];Number of Entries',
                             path='Overview',   xbins=200, xmin=0., xmax=800.
                            )
    mdtGroup.defineHistogram('nPrd_mon,Nhitsrpc_mon;MdtNHitsvsRpcNHits',  type='TH2F',
                             title='MdtNHitsvsRpcNHits;[Number of MDT hits];[Number of RPC hits]',
                             path='Overview',  xbins=1000, xmin=0., xmax=100000., ybins=100, ymin=0., ymax=10000.
                            )
    mdtGroup.defineHistogram('nPrdcut_mon;TotalNumber_of_MDT_hits_per_event_ADCCut',  type='TH1F',
                             title='TotalNumber_of_MDT_hits_per_event_ADCCut;[counts];Number of Event',
                             path='Overview',  xbins=400, xmin=0., xmax=10000.
                            )
 
    mdtGroup.defineHistogram('nPrdcut_mon;TotalNumber_of_MDT_hits_per_event_ADCCut_big',  type='TH1F',
                             title='TotalNumber_of_MDT_hits_per_event_ADCCut_big;[counts];Number of Event',
                             path='Overview/Expert',  xbins=200, xmin=0., xmax=100000.
                            )
 
    mdtGroup.defineHistogram('nPrd_mon;TotalNumber_of_MDT_hits_per_event',  type='TH1F',
                             title='TotalNumber_of_MDT_hits_per_event;[counts];Number of Event',
                             path='Overview',  xbins=500, xmin=0., xmax=20000.
                            )
 
    mdtGroup.defineHistogram('nPrd_mon;TotalNumber_of_MDT_hits_per_event_big',  type='TH1F',
                             title='TotalNumber_of_MDT_hits_per_event_big;[counts];Number of Event',
                             path='Overview/Expert',  xbins=200, xmin=0., xmax=100000.
                            )
 
    mdtGroup.defineHistogram('nColl_mon;number_of_Chambers_with_hits_per_event',  type='TH1F',
                             title='number_of_Chambers_with_hits_per_event;[Number_of_MDT_chambers_with_hits];Number of Entries',
                             path='Overview/Expert',  xbins=400, xmin=0., xmax=1600.
                            )
 
    mdtGroup.defineHistogram('nColl_ADCCut_mon;number_of_Chambers_with_hits_per_event_ADCCut',  type='TH1F',
                             title='number_of_Chambers_with_hits_per_event_ADCCut;[Number_of_MDT_chambers_with_hits];Number of Entries',
                             path='Overview/Expert',  xbins=400, xmin=0., xmax=1600.
                            )
 
    bigRegions = ["Barrel","EndCap"]
    regions = ["BA","BC","EA","EC"]
    layers = ["Inner","Middle","Outer","Extra"]
    crates  = ["01","02","03","04"]
  
    for iregion in bigRegions:
        theTitle = "NumberOfHitsIn"+iregion+"PerChamber_ADCCut_forpp"
        theTitle_noise = "NumberOfHitsIn"+iregion+"PerChamber_ADCCutt_NoiseBurst"
        theTitle_segs = "NumberOfHitsIn"+iregion+"PerChamber_onSegms_ADCCut_forpp"
        theTitle_eff = "effsIn"+iregion+"PerChamber_ADCCut"
        thisLabelx,thisLabely=getMDTLabel(theTitle+"_labelx",theTitle+"_labely")
        if iregion=="Barrel":
            etamin=-8.
            etamax=9.
            etabins=17
            phimin=0.
            phimax=58. #as defined in binMdtGlobal() c++ function
            phibins=58
            hits="x_mon_barrel,y_mon_barrel;"
            hits_noise="x_mon_barrel_noise,y_mon_barrel_noise;"
            hits_segs="x_segs_mon_barrel,y_segs_mon_barrel;"
            hits_eff="x_eff_mon_barrel,y_eff_mon_barrel;"
        if iregion=="EndCap":
            etamin=-6.
            etamax=6.
            etabins=12
            phimin=0.
            phimax=64. #as defined in binMdtGlobal() c++ function       
            phibins=64
            hits="x_mon_endcap,y_mon_endcap;"
            hits_noise="x_mon_endcap_noise,y_mon_endcap_noise;"
            hits_segs="x_segs_mon_endcap,y_segs_mon_endcap;"
            hits_eff="x_eff_mon_endcap,y_eff_mon_endcap;"
 
        mdtGroup.defineHistogram(hits+theTitle, type='TH2F', title=theTitle+";Eta;LayerPhi", path='Overview',
                                 xbins=etabins, xmin=etamin, xmax=etamax,  ybins=phibins, ymin=phimin, ymax=phimax, xlabels=thisLabelx, ylabels=thisLabely)  
 
        mdtGroup.defineHistogram(hits_noise+theTitle_noise, type='TH2F', title=theTitle_noise+";Eta;LayerPhi", path='Overview',
                                 xbins=etabins, xmin=etamin, xmax=etamax,  ybins=phibins, ymin=phimin, ymax=phimax, xlabels=thisLabelx, ylabels=thisLabely)  
 
        mdtGroup.defineHistogram(hits_segs+theTitle_segs, type='TH2F', title=theTitle_segs+";Eta;LayerPhi", path='Overview/Expert/RecoMonitoring',
                                 xbins=etabins, xmin=etamin, xmax=etamax,  ybins=phibins, ymin=phimin, ymax=phimax, xlabels=thisLabelx, ylabels=thisLabely)  
 
        mdtGroup.defineHistogram(hits_eff+theTitle_eff, type='TH2F', title=theTitle_eff+";Eta;LayerPhi", path='Overview/Expert/RecoMonitoring',
                                 xbins=etabins, xmin=etamin, xmax=etamax,  ybins=phibins, ymin=phimin, ymax=phimax, xlabels=thisLabelx, ylabels=thisLabely)  
        
 
 
    for ilayer in layers:
        if(ilayer=="Extra"):
            continue
        hits="x_mon_"+ilayer+",y_mon_"+ilayer+";"
        theTitle = "NumberOfHitsInMDT"+ilayer+"_ADCCut"
        thisLabelx,thisLabely=getMDTLabel(theTitle+"_labelx",theTitle+"_labely")
        if(ilayer=="Inner"):
            etabins=38
            phibins=36
        if(ilayer=="Middle"):
            etabins=24
            phibins=32
        if(ilayer=="Outer"):
            etabins=31
            phibins=32
        mdtGroup.defineHistogram(hits+theTitle, type='TH2F', title=theTitle+";Eta station;Phi station", path='Overview',
                                 xbins=etabins, xmin=0, xmax=etabins,  ybins=phibins, ymin=0, ymax=phibins, xlabels=thisLabelx, ylabels=thisLabely)  
        
 
 
    phimax = 16
    for iregion in regions:
        MDT_regionGroup="MDT_regionGroup{0}".format(iregion)
        mdtRegionGroup=helper.addGroup(mdtMonAlg, MDT_regionGroup, "Muon/MuonRawDataMonitoring/MDT/MDT"+iregion)
        MDT_regionGroup_bycrate="MDT_regionGroup_bycrate{0}".format(iregion)
        mdtRegionGroup_bycrate=helper.addGroup(mdtMonAlg, MDT_regionGroup_bycrate, "Muon/MuonRawDataMonitoring/MDT/MDT"+iregion)
 
        #histograms from segments 
 
        title_ADC_segs_spectrum="MDTADC_segm_Summary_"+iregion
        mdtRegionGroup.defineHistogram("adc_segs_mon;"+title_ADC_segs_spectrum, title=title_ADC_segs_spectrum+";[adc counts];Number of Entries", type='TH1F',
                                             path='Overview', xbins=100, xmin=0.5, xmax=400.5)
 
        title_TDCADC_spectrum="Overall_TDCADC_spectrum_"+iregion
        title_TDCADC_spectrum_noiseBurst="Overall_TDCADC_spectrum_"+iregion+"_NoiseBurst"
        mdtRegionGroup.defineHistogram("tdc_mon,adc_mon;"+title_TDCADC_spectrum, title=title_TDCADC_spectrum+";[nsec];[adc counts]", type='TH2F',
                                             path='Overview', xbins=50, xmin=0., xmax=2000., ybins=20, ymin=0., ymax=400.)
        mdtRegionGroup.defineHistogram("tdc_mon_nb2,adc_mon_nb2;"+title_TDCADC_spectrum_noiseBurst, title=title_TDCADC_spectrum_noiseBurst+";[nsec];[adc counts]",
                                                       type='TH2F', path='Overview', xbins=50, xmin=0., xmax=2000., ybins=20, ymin=0., ymax=400.)
        title_ADCSummaryRegion="MDTADC_Summary_"+iregion
        mdtRegionGroup.defineHistogram("adc_mon;"+title_ADCSummaryRegion,  type='TH1F',
                                       title=title_ADCSummaryRegion+";[adc counts];Number of Entries",
                                       path='Overview',   xbins=100, xmin=0., xmax=400.)
 
        title_ADCSummaryRegion="MDTADC_Summary_ADCCUT_"+iregion
        mdtRegionGroup.defineHistogram("adc_mon_adccut;"+title_ADCSummaryRegion,  type='TH1F',
                                       title=title_ADCSummaryRegion+";[adc counts];Number of Entries",
                                       path='Overview',   xbins=100, xmin=0., xmax=400.)
 
        title_TDCSummaryRegion="MDTTDC_Summary_ADCCut_"+iregion
        mdtRegionGroup.defineHistogram("tdc_mon_adccut;"+title_TDCSummaryRegion,  type='TH1F',
                                       title=title_TDCSummaryRegion+";[nsec];Number of Entries",
                                       path='Overview',   xbins=120, xmin=0., xmax=2000.)
 
        title_TDCSummaryRegion_segs="MDTTDC_segm_Summary_ADCCut_"+iregion
        mdtRegionGroup.defineHistogram("tdc_segs_region_mon;"+title_TDCSummaryRegion_segs,  type='TH1F',
                                       title=title_TDCSummaryRegion_segs+";[nsec];Number of Entries",
                                       path='Overview',   xbins=120, xmin=0., xmax=2000.)
        
        title_TDC_NoiseBursts = "MDTTDC_NoiseBursts_"+iregion
        mdtRegionGroup.defineHistogram("tdc_mon_nb1;"+title_TDC_NoiseBursts,  type='TH1F',
                                       title=title_TDC_NoiseBursts+";[nsec];Number of Entries",
                                       path='Overview',   xbins=120, xmin=0., xmax=2000.)
 
        title_TDC_NoiseBursts = "MDTTDC_NoiseBursts_ADCCUT_"+iregion
        mdtRegionGroup.defineHistogram("tdc_mon_nb3;"+title_TDC_NoiseBursts,  type='TH1F',
                                       title=title_TDC_NoiseBursts+";[nsec];Number of Entries",
                                       path='Overview',   xbins=120, xmin=0., xmax=2000.)
 
        title_ADC_NoiseBursts = "MDTADC_NoiseBursts_"+iregion
        mdtRegionGroup.defineHistogram("adc_mon_nb1;"+title_ADC_NoiseBursts,  type='TH1F',
                                       title=title_ADC_NoiseBursts+";[adc counts];Number of Entries",
                                       path='Overview',   xbins=100, xmin=0., xmax=400.)
 
        title_TDC_RPCtrig = "MDTTDC_Summary_ADCCut_"+iregion+"_RPC"
        mdtRegionGroup.defineHistogram("tdc_mon_rpc;"+title_TDC_RPCtrig,  type='TH1F',
                                       title=title_TDC_RPCtrig+";[nsec];Number of Entries",
                                       path='Overview',   xbins=120, xmin=0., xmax=2000.)
 
        title_TDC_TGCtrig = "MDTTDC_Summary_ADCCut_"+iregion+"_TGC"
        mdtRegionGroup.defineHistogram("tdc_mon_tgc;"+title_TDC_TGCtrig,  type='TH1F',
                                       title=title_TDC_TGCtrig+";[nsec];Number of Entries",
                                       path='Overview',   xbins=120, xmin=0., xmax=2000.)
 
        for ilayer in layers:
            title_MDTNumberOfHitsPerML="NumberOfHitsIn"+iregion+ilayer+"PerMultiLayer_ADCCut_forpp"
            thisLabelx,thisLabely=getMDTLabel(title_MDTNumberOfHitsPerML+"_labelx",title_MDTNumberOfHitsPerML+"_labely")
            phimaxML=phimax*2
            if (iregion=="BA" or iregion=="BC") and ilayer=="Extra": phimaxML=phimax//2
            if (iregion=="BA" or iregion=="BC") and ilayer=="Inner": phimaxML=phimax*2+4
            etamax=2 #extra layer
            if (iregion=="BA" or iregion=="BC") and ilayer=="Inner": etamax=8
            if (iregion=="BA" or iregion=="BC") and ilayer=="Middle": etamax=6
            if iregion=="BA" and  ilayer=="Outer": etamax=9
            if iregion=="BC" and  ilayer=="Outer": etamax=8
            if (iregion=="EA" or iregion=="EC") and ilayer=="Inner": etamax=5
            if (iregion=="EA" or iregion=="EC") and ilayer=="Middle": etamax=5
            if (iregion=="EA" or iregion=="EC") and ilayer=="Outer": etamax=6
            var="x_mon_"+iregion+"_"+ilayer+",y_mon_"+iregion+"_"+ilayer+";"
            mdtRegionGroup.defineHistogram(var+title_MDTNumberOfHitsPerML, title=title_MDTNumberOfHitsPerML+";[Eta];[Phi - Multilayer]", type='TH2F',
                                           path='Overview',  xbins=etamax, xmin=0, xmax=etamax, ybins=phimaxML, ymin=0., ymax=phimaxML, xlabels=thisLabelx, ylabels=thisLabely)
 
            for phi in range(1, phimax+1):
                if not (mdtMonAlg.do_Run3Geometry and (iregion=="EA" or iregion=="EC") and (phi%2==0)):
                    if not ((iregion=="BA" or iregion=="BC") and ilayer=="Extra" and (phi%2!=0)):
                        var="stEta_"+iregion+"_"+ilayer+"_phi"+str(phi)+";"
                        title_MDTHitSummary="MDTHits_ADCCut_"+iregion+"_"+ilayer+"_StPhi"+str(phi)
                    
                        max=7
                        if ilayer=="Inner" and (iregion=="BA" or iregion=="BC"):
                            max=11
                        
                        mdtRegionGroup.defineHistogram(var+title_MDTHitSummary, title=title_MDTHitSummary+";StationEta;Number of Entries", type='TH1F',
                                                    path='Overview/Hits',  xbins=max-1, xmin=1., xmax=max)
    
                    
                        vartdc="tdc_segs_"+iregion+"_"+ilayer+"_phi"+str(phi)+";"
                        title_MDTTDCSummary="MDTTDC_ADCCut_"+iregion+"_"+ilayer+"_StPhi"+str(phi)
                        mdtRegionGroup.defineHistogram(vartdc+title_MDTTDCSummary, title=title_MDTTDCSummary+";[nsec];Number of Entries", type='TH1F',
                                                       path='Overview/TDC',  xbins=100, xmin=0, xmax=2000)
                    
 
            if(ilayer=="Outer"):
                labelsY=getMDTLabelx("labelY_OccupancyVsLB_"+iregion+"_OuterExtra")
                titleOccvsLbPerRegionPerLayer = "OccupancyVsLB_"+iregion+"OuterPlusExtra"
                var="lb_mon,y_mon_bin_"+iregion+"_OuterPlusExtra;"+titleOccvsLbPerRegionPerLayer
                if(iregion=="BA"):
                     maxy=118 # outer sideA
                if(iregion=="BC"):
                    maxy=116 #  outer sideC
                if(iregion=="EA"):
                     maxy=103 # outer sideA
                if(iregion=="EC"):
                     maxy=127 # outer sideA
            elif(ilayer=="Middle" or ilayer=="Inner"):
                labelsY=getMDTLabelx("labelY_OccupancyVsLB_"+iregion+"_"+ilayer)
                titleOccvsLbPerRegionPerLayer = "OccupancyVsLB_"+iregion+ilayer
                var="lb_mon,y_mon_bin_"+iregion+"_"+ilayer+";"+titleOccvsLbPerRegionPerLayer
                if(ilayer=="Inner"):
                    if(iregion=="BA" or iregion=="BC"):
                        maxy=122
                        if (mdtMonAlg.do_Run3Geometry and iregion=="BA"):
                            labelsY=getMDTLabelx("run3_labelY_OccupancyVsLB_"+iregion+"_"+ilayer)
                            maxy=114
                    if(iregion=="EA" or iregion=="EC"):
                        maxy=50
                        if (mdtMonAlg.do_Run3Geometry):
                            maxy=10
                            labelsY=getMDTLabelx("run3_labelY_OccupancyVsLB_"+iregion+"_"+ilayer)
                elif (ilayer=="Middle"):
                    if(iregion=="BA" or iregion=="BC"):
                        maxy=95
                    if(iregion=="EA" or iregion=="EC"):
                        maxy=80
 
            if ilayer != 'Extra':
                mdtRegionGroup.defineHistogram(var, title=titleOccvsLbPerRegionPerLayer+";LB;[Eta - Phi]", type='TH2F',
                                               path='Overview', xbins=834, xmin=1, xmax=2502, ybins=maxy, ymin=0, ymax=maxy, ylabels=labelsY, opt='kAddBinsDynamically', merge='merge')
 
        for icrate in crates:
            maxy=122
            labelsY=getMDTLabelx("labelY_OccupancyVsLB_"+iregion+icrate)
            if(iregion=="BA" or iregion=="BC"):
                if(icrate=="01"):
                    maxy=73
                elif(icrate=="02"):
                    maxy=72
                elif(icrate=="03"):
                    maxy= 80
                elif(icrate=="04"):
                    maxy= 79
                    if(iregion=="BA"):
                        maxy= 80
            elif(iregion=="EA" or iregion=="EC"):
                if mdtMonAlg.do_Run3Geometry:
                    labelsY=getMDTLabelx("run3_labelY_OccupancyVsLB_"+iregion+icrate)
                if(icrate=="01"):
                    maxy=73
                    if(mdtMonAlg.do_Run3Geometry and iregion=="EA"):
                        maxy=61
                    if(mdtMonAlg.do_Run3Geometry  and iregion=="EC"):
                        maxy=63
                elif(icrate=="02"):
                    maxy=71
                    if(mdtMonAlg.do_Run3Geometry and iregion=="EA"):
                        maxy=59
                    if(mdtMonAlg.do_Run3Geometry  and iregion=="EC"):
                        maxy=61
                elif(icrate=="03"):
                    maxy= 73
                    if(mdtMonAlg.do_Run3Geometry and iregion=="EA"):
                        maxy=61
                    if(mdtMonAlg.do_Run3Geometry  and iregion=="EC"):
                        maxy=63
                elif(icrate=="04"):
                    maxy= 72
                    if(mdtMonAlg.do_Run3Geometry and iregion=="EA"):
                        maxy=60
                    if(mdtMonAlg.do_Run3Geometry  and iregion=="EC"):
                        maxy=62
 
            titleOccvsLbPerRegionPerCrate = "OccupancyVsLB_"+iregion+icrate 
            var="lb_mon,y_mon_bin_bycrate_"+iregion+"_"+icrate+";"+titleOccvsLbPerRegionPerCrate
            mdtRegionGroup_bycrate.defineHistogram(var, type='TH2F', title=titleOccvsLbPerRegionPerCrate+";LB;[Eta-Phi]", 
                                                   path='Overview', xbins=33, xmin=-0.5, xmax=98.5, ybins=maxy, ymin=0, ymax=maxy, ylabels=labelsY, opt='kAddBinsDynamically', merge='merge')
                                           
            titleOccvsLbPerRegionPerCrate_ontrack = "OccupancyVsLB_ontrack_"+iregion+icrate
            var="lb_mon,y_mon_bin_bycrate_ontrack_"+iregion+"_"+icrate+";"+titleOccvsLbPerRegionPerCrate_ontrack
            mdtRegionGroup_bycrate.defineHistogram(var, title=titleOccvsLbPerRegionPerCrate_ontrack+";LB;[Eta - Phi]", type='TH2F',
                                           path='Overview', xbins=33, xmin=-0.5, xmax=98.5, ybins=maxy, ymin=0, ymax=maxy, ylabels=labelsY, opt='kAddBinsDynamically', merge='merge') 
 
    labelsY=getMDTLabelx("labelY_OccupancyVsLB")
    mdtGroup.defineHistogram('lb_mon,sector;OccupancyPerSectorVsLB', type='TH2F',
                             title='OccupancyPerSectorVsLB;LB;[Phi]',
                             path='Overview', xbins=33, xmin= -0.5, xmax=98.5, ybins=64, ymin=0., ymax=64., ylabels=labelsY, opt='kAddBinsDynamically', merge='merge')
    
 
    #histograms per chambers
    mdtPerChamberBAGroup=helper.addGroup(mdtMonAlg,"MdtMonPerChamberBA", "Muon/MuonRawDataMonitoring/MDT/MDTBA/Chambers")    
    for ch in mdtBA:
        title_mdttdc= ch+"_MDT_Station_TDC"
        var="tdc_perch_"+ch+";"+title_mdttdc
        mdtPerChamberBAGroup.defineHistogram(var, title=title_mdttdc+";[nsec];Number of Entries", 
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdttdc_ml1_adccut=ch+"_MDT_Station_TDC_ML1_ADCCut"
        var="tdc_perch_"+ch+";"+title_mdttdc_ml1_adccut
        mdtPerChamberBAGroup.defineHistogram(var, title=title_mdttdc_ml1_adccut+";[nsec];Number of Entries",
                                             cutmask='ml1_adccut',
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdttdc_ml2_adccut=ch+"_MDT_Station_TDC_ML2_ADCCut"
        var="tdc_perch_"+ch+";"+title_mdttdc_ml2_adccut
        mdtPerChamberBAGroup.defineHistogram(var, title=title_mdttdc_ml2_adccut+";[nsec];Number of Entries",
                                             cutmask='ml2_adccut',
                                           type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdtadc= ch+"_MDT_Station_ADC"
        var="adc_perch_"+ch+";"+title_mdtadc
        mdtPerChamberBAGroup.defineHistogram(var,  type='TH1F',
                                             title=title_mdtadc+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
 
        title_mdtadc_segs_ml1= ch+"_MDT_Station_ADC_onSegm_ML1"
        var="adc_segs_perch_"+ch+";"+title_mdtadc_segs_ml1
        mdtPerChamberBAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adc_ml1',
                                             title=title_mdtadc_segs_ml1+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
        title_mdtadc_segs_ml2= ch+"_MDT_Station_ADC_onSegm_ML2"
        var="adc_segs_perch_"+ch+";"+title_mdtadc_segs_ml2
        mdtPerChamberBAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adc_ml2',
                                             title=title_mdtadc_segs_ml2+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
        
        title_mdtdcadc= ch+"_MDT_Station_TDCADC"
        var="tdc_perch_"+ch+",adc_perch_"+ch+";"+title_mdtdcadc
        mdtPerChamberBAGroup.defineHistogram(var, type='TH2F',
                                             title=title_mdtdcadc+';[nsec];[adc counts]',
                                             path=ch,   xbins=50, xmin=0., xmax=2000., ybins=20, ymin=0., ymax=400.)
        
 
        title_mdtlayer= ch+"_MDT_Station_LAYER_ADCCut"
        var="layer_perch_"+ch+";"+title_mdtlayer
        mdtPerChamberBAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adccut_nonoise',
                                             title=title_mdtlayer+";layerID;Number of Entries",
                                             path=ch,   xbins=10, xmin=0., xmax=10.)
        
        title_mdttube= ch+"_MDT_Station_TUBE_ADCCut"
        var="tube_perch_"+ch+";"+title_mdttube
        binmax=tubeMax[ch]
        if mdtMonAlg.do_Run3Geometry and ch[0:5]=="BIS7A":
            binmax=tubeMax_smdt[ch]
        mdtPerChamberBAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adccut',
                                             title=title_mdttube+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
        
        title_effentries=ch+"_MDT_Station_EFFENTRIES"
        var="tube_perch_segs_"+ch+";"+title_effentries
        mdtPerChamberBAGroup.defineHistogram(var,  type='TH1F',
                                             title=title_effentries+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_effcounts=ch+"_MDT_Station_EFFCOUNTS"
        var="tube_perch_segs_"+ch+";"+title_effcounts
        mdtPerChamberBAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='hitcut',
                                             title=title_effcounts+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_eff=ch+"_MDT_Station_EFFPERTUBE"
        var_eff="hitcut,tube_perch_segs_"+ch+";"+title_eff
        mdtPerChamberBAGroup.defineHistogram(var_eff,
                                             title=title_eff+";tubeID;Efficiency",
                                             type='TEfficiency',
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
                                             
 
 
        title_mdtmezz= ch+"_MDT_Station_MEZZ_ADCCut"
        mdtPerChamberBAGroup.defineHistogram("mezz_perch_"+ch+";"+title_mdtmezz,  type='TH1F',
                                             title=title_mdtmezz+";mezzID;Number of Entries",
                                             path=ch,   xbins=20, xmin=0., xmax=20.)
 
    mdtPerChamberBCGroup=helper.addGroup(mdtMonAlg, "MdtMonPerChamberBC", "Muon/MuonRawDataMonitoring/MDT/MDTBC/Chambers")    
    for ch in mdtBC:
        title_mdttdc=ch+"_MDT_Station_TDC"
        var="tdc_perch_"+ch+";"+title_mdttdc
        mdtPerChamberBCGroup.defineHistogram(var, title=title_mdttdc+";[nsec];Number of Entries",
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdttdc_ml1_adccut=ch+"_MDT_Station_TDC_ML1_ADCCut"
        var="tdc_perch_"+ch+";"+title_mdttdc_ml1_adccut
        mdtPerChamberBCGroup.defineHistogram(var, title=title_mdttdc_ml1_adccut+";[nsec];Number of Entries",
                                             cutmask='ml1_adccut',
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdttdc_ml2_adccut=ch+"_MDT_Station_TDC_ML2_ADCCut"
        var="tdc_perch_"+ch+";"+title_mdttdc_ml2_adccut
        mdtPerChamberBCGroup.defineHistogram(var, title=title_mdttdc_ml2_adccut+";[nsec];Number of Entries",
                                             cutmask='ml2_adccut',
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdtadc= ch+"_MDT_Station_ADC"
        var="adc_perch_"+ch+";"+title_mdtadc
        mdtPerChamberBCGroup.defineHistogram(var, title=title_mdtadc+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
 
        title_mdtadc_segs_ml1= ch+"_MDT_Station_ADC_onSegm_ML1"
        var="adc_segs_perch_"+ch+";"+title_mdtadc_segs_ml1
        mdtPerChamberBCGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adc_ml1',
                                             title=title_mdtadc_segs_ml1+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
        title_mdtadc_segs_ml2= ch+"_MDT_Station_ADC_onSegm_ML2"
        var="adc_segs_perch_"+ch+";"+title_mdtadc_segs_ml2
        mdtPerChamberBCGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adc_ml2',
                                             title=title_mdtadc_segs_ml2+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
 
 
        title_mdtdcadc= ch+"_MDT_Station_TDCADC"
        var="tdc_perch_"+ch+",adc_perch_"+ch+";"+title_mdtdcadc
        mdtPerChamberBCGroup.defineHistogram(var,    type='TH2F',
                                             title=title_mdtdcadc+';[nsec];[adc counts]',
                                             path=ch,   xbins=50, xmin=0., xmax=2000., ybins=20, ymin=0., ymax=400.)
        title_mdtlayer= ch+"_MDT_Station_LAYER_DCCut"
        var="layer_perch_"+ch+";"+title_mdtlayer
        mdtPerChamberBCGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adccut_nonoise',
                                             title=title_mdtlayer+";layerID;Number of Entries",
                                             path=ch,   xbins=10, xmin=0., xmax=10.)
        title_mdttube= ch+"_MDT_Station_TUBE_ADCCut"
        var="tube_perch_"+ch+";"+title_mdttube
        binmax=tubeMax[ch]
        mdtPerChamberBCGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adccut',
                                             title=title_mdttube+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_effentries=ch+"_MDT_Station_EFFENTRIES"
        var="tube_perch_segs_"+ch+";"+title_effentries
        mdtPerChamberBCGroup.defineHistogram(var,  type='TH1F',
                                             title=title_effentries+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_effcounts=ch+"_MDT_Station_EFFCOUNTS"
        var="tube_perch_segs_"+ch+";"+title_effcounts
        mdtPerChamberBCGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='hitcut',
                                             title=title_effcounts+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
        title_eff=ch+"_MDT_Station_EFFPERTUBE"
        var_eff="hitcut,tube_perch_segs_"+ch+";"+title_eff
        mdtPerChamberBCGroup.defineHistogram(var_eff,
                                             title=title_eff+";tubeID;Efficiency",
                                             type='TEfficiency',
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_mdtmezz= ch+"_MDT_Station_MEZZ_ADCCut"
        mdtPerChamberBCGroup.defineHistogram("mezz_perch_"+ch+";"+title_mdtmezz,  type='TH1F',
                                             title=title_mdtmezz+";mezzID;Number of Entries",
                                             path=ch,   xbins=20, xmin=0., xmax=20.)
 
    mdtPerChamberEAGroup=helper.addGroup(mdtMonAlg, "MdtMonPerChamberEA", "Muon/MuonRawDataMonitoring/MDT/MDTEA/Chambers")           
    for ch in mdtEA:
        title_mdttdc=ch+"_MDT_Station_TDC"
        var="tdc_perch_"+ch+";"+title_mdttdc
        mdtPerChamberEAGroup.defineHistogram(var, title=title_mdttdc+";[nsec];Number of Entries", 
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdttdc_ml1_adccut=ch+"_MDT_Station_TDC_ML1_ADCCut"
        var="tdc_perch_"+ch+";"+title_mdttdc_ml1_adccut
        mdtPerChamberEAGroup.defineHistogram(var, title=title_mdttdc_ml1_adccut+";[nsec];Number of Entries",
                                             cutmask='ml1_adccut',
                                            type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdttdc_ml2_adccut=ch+"_MDT_Station_TDC_ML2_ADCCut"
        var="tdc_perch_"+ch+";"+title_mdttdc_ml2_adccut
        mdtPerChamberEAGroup.defineHistogram(var, title=title_mdttdc_ml2_adccut+";[nsec];Number of Entries",
                                             cutmask='ml2_adccut',
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdtadc= ch+"_MDT_Station_ADC"
        var="adc_perch_"+ch+";"+title_mdtadc
        mdtPerChamberEAGroup.defineHistogram(var,  type='TH1F',
                                             title=title_mdtadc+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
      
        title_mdtadc_segs_ml1= ch+"_MDT_Station_ADC_onSegm_ML1"
        var="adc_segs_perch_"+ch+";"+title_mdtadc_segs_ml1
        mdtPerChamberEAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adc_ml1',
                                             title=title_mdtadc_segs_ml1+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
        title_mdtadc_segs_ml2= ch+"_MDT_Station_ADC_onSegm_ML2"
        var="adc_segs_perch_"+ch+";"+title_mdtadc_segs_ml2
        mdtPerChamberEAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adc_ml2',
                                             title=title_mdtadc_segs_ml2+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
       
 
 
        title_mdtdcadc= ch+"_MDT_Station_TDCADC"
        var="tdc_perch_"+ch+",adc_perch_"+ch+";"+title_mdtdcadc
        mdtPerChamberEAGroup.defineHistogram(var,    type='TH2F',
                                             title=title_mdtdcadc+";[nsec];[adc counts]",
                                             path=ch,   xbins=50, xmin=0., xmax=2000., ybins=20, ymin=0., ymax=400.)
        
        title_mdtlayer= ch+"_MDT_Station_LAYER_DCCut"
        var="layer_perch_"+ch+";"+title_mdtlayer
        mdtPerChamberEAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adccut_nonoise',
                                             title=title_mdtlayer+";layerID;Number of Entries",
                                             path=ch,   xbins=10, xmin=0., xmax=10.)
 
        title_mdttube= ch+"_MDT_Station_TUBE_ADCCut"
        var="tube_perch_"+ch+";"+title_mdttube
        binmax=tubeMax[ch]
        mdtPerChamberEAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adccut',
                                             title=title_mdttube+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_effentries=ch+"_MDT_Station_EFFENTRIES"
        var="tube_perch_segs_"+ch+";"+title_effentries
        mdtPerChamberEAGroup.defineHistogram(var,  type='TH1F',
                                             title=title_effentries+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_effcounts=ch+"_MDT_Station_EFFCOUNTS"
        var="tube_perch_segs_"+ch+";"+title_effcounts
        mdtPerChamberEAGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='hitcut',
                                             title=title_effcounts+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_eff=ch+"_MDT_Station_EFFPERTUBE"
        var_eff="hitcut,tube_perch_segs_"+ch+";"+title_eff
        mdtPerChamberEAGroup.defineHistogram(var_eff,
                                             title=title_eff+";tubeID;Efficiency",
                                             type='TEfficiency',
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_mdtmezz= ch+"_MDT_Station_MEZZ_ADCCut"
        mdtPerChamberEAGroup.defineHistogram("mezz_perch_"+ch+";"+title_mdtmezz,  type='TH1F',
                                             cutmask='adccut',
                                             title=title_mdtmezz+";mezzID;Number of Entries",
                                             path=ch,   xbins=20, xmin=0., xmax=20.)
    
    mdtPerChamberECGroup=helper.addGroup(mdtMonAlg, "MdtMonPerChamberEC", "Muon/MuonRawDataMonitoring/MDT/MDTEC/Chambers")    
    for ch in mdtEC:
        title_mdttdc=ch+"_MDT_Station_TDC"
        var="tdc_perch_"+ch+";"+title_mdttdc
        mdtPerChamberECGroup.defineHistogram(var, title=title_mdttdc+";[nsec];Number of Entries", 
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdttdc_ml1_adccut=ch+"_MDT_Station_TDC_ML1_ADCCut"
        var="tdc_perch_"+ch+";"+title_mdttdc_ml1_adccut
        mdtPerChamberECGroup.defineHistogram(var, title=title_mdttdc_ml1_adccut+";[nsec];Number of Entries",
                                             cutmask='ml1_adccut',
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdttdc_ml2_adccut=ch+"_MDT_Station_TDC_ML2_ADCCut"
        var="tdc_perch_"+ch+";"+title_mdttdc_ml2_adccut
        mdtPerChamberECGroup.defineHistogram(var, title=title_mdttdc_ml2_adccut+";[nsec];Number of Entries",
                                             cutmask='ml2_adccut',
                                             type='TH1F', path=ch, xbins=100, xmin=0., xmax=2000.)
        title_mdtadc= ch+"_MDT_Station_ADC"
        var="adc_perch_"+ch+";"+title_mdtadc
        mdtPerChamberECGroup.defineHistogram(var,  type='TH1F',
                                             title=title_mdtadc+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
        title_mdtadc_segs_ml1= ch+"_MDT_Station_ADC_onSegm_ML1"
        var="adc_segs_perch_"+ch+";"+title_mdtadc_segs_ml1
        mdtPerChamberECGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adc_ml1',
                                             title=title_mdtadc_segs_ml1+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
        title_mdtadc_segs_ml2= ch+"_MDT_Station_ADC_onSegm_ML2"
        var="adc_segs_perch_"+ch+";"+title_mdtadc_segs_ml2
        mdtPerChamberECGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adc_ml2',
                                             title=title_mdtadc_segs_ml2+";[adc counts];Number of Entries",
                                             path=ch,   xbins=100, xmin=0., xmax=400.)
 
 
 
        title_mdtdcadc= ch+"_MDT_Station_TDCADC"
        var="tdc_perch_"+ch+",adc_perch_"+ch+";"+title_mdtdcadc
        mdtPerChamberECGroup.defineHistogram(var,    type='TH2F',
                                             title=title_mdtdcadc+';[nsec];[adc counts]',
                                             path=ch,   xbins=50, xmin=0., xmax=2000., ybins=20, ymin=0., ymax=400.)
        #to do
        title_mdtlayer= ch+"_MDT_Station_LAYER_ADCCut"
        var="layer_perch_"+ch+";"+title_mdtlayer
        mdtPerChamberECGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adccut_nonoise',
                                             title=title_mdtlayer+";layerID;Number of Entries",
                                             path=ch,   xbins=10, xmin=0., xmax=10.)
        title_mdttube= ch+"_MDT_Station_TUBE_ADCCut"
        var="tube_perch_"+ch+";"+title_mdttube
        binmax=tubeMax[ch]
        mdtPerChamberECGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='adccut',
                                             title=title_mdttube+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=0., xmax=binmax)
 
        title_effentries=ch+"_MDT_Station_EFFENTRIES"
        var="tube_perch_segs_"+ch+";"+title_effentries
        mdtPerChamberECGroup.defineHistogram(var,  type='TH1F',
                                             title=title_effentries+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
        title_effcounts=ch+"_MDT_Station_EFFCOUNTS"
        var="tube_perch_segs_"+ch+";"+title_effcounts
        mdtPerChamberECGroup.defineHistogram(var,  type='TH1F',
                                             cutmask='hitcut',
                                             title=title_effcounts+";tubeID;Number of Entries",
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
        
        title_eff=ch+"_MDT_Station_EFFPERTUBE"
        var_eff="hitcut,tube_perch_segs_"+ch+";"+title_eff
        mdtPerChamberECGroup.defineHistogram(var_eff,
                                             title=title_eff+";tubeID;Efficiency",
                                             type='TEfficiency',
                                             path=ch,   xbins=binmax, xmin=1., xmax=binmax+1)
 
 
        title_mdtmezz= ch+"_MDT_Station_MEZZ_ADCCut"
        mdtPerChamberECGroup.defineHistogram("mezz_perch_"+ch+";"+title_mdtmezz,  type='TH1F',
                                             title=title_mdtmezz+";mezzID;Number of Entries",
                                             path=ch,   xbins=20, xmin=0., xmax=20.)
 
 
 
    
        
    
 
 
    
    # # Otherwise, merge with result object and return
    #acc, seq = helper.result()
    acc = helper.result()
    result.merge(acc)
    return result
 
 
if __name__=='__main__':
 
    # Set the Athena configuration flags
    from AthenaConfiguration.AllConfigFlags import initConfigFlags
    from AthenaConfiguration.TestDefaults import defaultTestFiles
    flags = initConfigFlags()
    flags.Input.Files = defaultTestFiles.ESD
    
    #flags.Input.isMC = True
    #flags.Common.isOnline = True
    flags.Output.HISTFileName = 'MdtMonitorOutput.root'
    #flags.Output.AODFileName='aod.root'
 
    #flags.IOVDb.GlobalTag = "CONDBR2-BLKPA-2018-12"
    #flags.IOVDb.GlobalTag =  "CONDBR2-BLKPA-2018-10" 
    #flags.GeoModel.AtlasVersion = "ATLAS-R2-2016-01-00-01"
    #flags.IOVDb.DatabaseInstance=""
    #flags.GeoModel.AtlasVersion = "ATLAS-R2-2016-01-00-01  "
 
    flags.Detector.GeometryMuon = False
    flags.Detector.GeometryCSC = False
    flags.Detector.GeometryRPC = False
    flags.Detector.GeometryTGC = False
    flags.Detector.GeometrysTGC = False
    flags.Detector.GeometryMM = False
 
    flags.Detector.GeometryMDT = True
 
    flags.Muon.Align.UseILines = False
    flags.Muon.Align.UseALines = False
    flags.Muon.Align.UseBLines = False
 
    flags.DQ.useTrigger=False
 
 
    flags.lock()
    flags.dump()
    
    # Initialize configuration object, add accumulator, merge, and run.
    from AthenaConfiguration.MainServicesConfig import MainServicesCfg 
    from AthenaPoolCnvSvc.PoolReadConfig import PoolReadCfg
    cfg = MainServicesCfg(flags)
 
    cfg.merge(PoolReadCfg(flags))
    
 
    #mdtMonitorAcc,mdtMonitorAlg =  MdtMonitoringConfig(flags)
    mdtMonitorAcc =  MdtMonitoringConfig(flags)
    
    #mdtMonitorAcc.OutputLevel=DEBUG
    
    cfg.merge(mdtMonitorAcc)
 
    #cfg.printConfig(withDetails=True, summariseProps = True)
    
    cfg.run(5000)