MDTMonitorAlgorithm Namespace Reference

Functions
def	MdtMonitoringConfig (inputFlags)

Variables
	flags = initConfigFlags()
	Files
	HISTFileName
	GeometryMuon
	GeometryCSC
	GeometryRPC
	GeometryTGC
	GeometrysTGC
	GeometryMM
	GeometryMDT
	UseILines
	UseALines
	UseBLines
	useTrigger
	cfg = MainServicesCfg(flags)
def	mdtMonitorAcc = MdtMonitoringConfig(flags)

Function Documentation

MdtMonitoringConfig()

def MDTMonitorAlgorithm.MdtMonitoringConfig

(

inputFlags

)

Function to configures some algorithms in the monitoring system.

Definition at line 12 of file MDTMonitorAlgorithm.py.

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
 
 

Variable Documentation

cfg

MDTMonitorAlgorithm.cfg = MainServicesCfg(flags)

Definition at line 869 of file MDTMonitorAlgorithm.py.

Files

MDTMonitorAlgorithm.Files

Definition at line 834 of file MDTMonitorAlgorithm.py.

flags

MDTMonitorAlgorithm.flags = initConfigFlags()

Definition at line 833 of file MDTMonitorAlgorithm.py.

GeometryCSC

MDTMonitorAlgorithm.GeometryCSC

Definition at line 848 of file MDTMonitorAlgorithm.py.

GeometryMDT

MDTMonitorAlgorithm.GeometryMDT

Definition at line 854 of file MDTMonitorAlgorithm.py.

GeometryMM

MDTMonitorAlgorithm.GeometryMM

Definition at line 852 of file MDTMonitorAlgorithm.py.

GeometryMuon

MDTMonitorAlgorithm.GeometryMuon

Definition at line 847 of file MDTMonitorAlgorithm.py.

GeometryRPC

MDTMonitorAlgorithm.GeometryRPC

Definition at line 849 of file MDTMonitorAlgorithm.py.

GeometrysTGC

MDTMonitorAlgorithm.GeometrysTGC

Definition at line 851 of file MDTMonitorAlgorithm.py.

GeometryTGC

MDTMonitorAlgorithm.GeometryTGC

Definition at line 850 of file MDTMonitorAlgorithm.py.

HISTFileName

MDTMonitorAlgorithm.HISTFileName

Definition at line 838 of file MDTMonitorAlgorithm.py.

mdtMonitorAcc

def MDTMonitorAlgorithm.mdtMonitorAcc = MdtMonitoringConfig(flags)

Definition at line 875 of file MDTMonitorAlgorithm.py.

UseALines

MDTMonitorAlgorithm.UseALines

Definition at line 857 of file MDTMonitorAlgorithm.py.

UseBLines

MDTMonitorAlgorithm.UseBLines

Definition at line 858 of file MDTMonitorAlgorithm.py.

UseILines

MDTMonitorAlgorithm.UseILines

Definition at line 856 of file MDTMonitorAlgorithm.py.

useTrigger

MDTMonitorAlgorithm.useTrigger

Definition at line 860 of file MDTMonitorAlgorithm.py.