HaloTrigger Class Reference

#include <HaloTrigger.h>

Inheritance diagram for HaloTrigger:

Public Member Functions
	HaloTrigger (const edm::ParameterSet &ps)
virtual	~HaloTrigger ()
Public Member Functions inherited from edm::EDAnalyzer
	EDAnalyzer ()
std::string	workerType () const
virtual	~EDAnalyzer ()

Public Attributes
unsigned int	CscHalo_Gmt
bool	first
unsigned int	hltHaloOver1
unsigned int	hltHaloOver2
unsigned int	hltHaloRing23
unsigned int	hltHaloTriggers
unsigned int	majikNumber0
unsigned int	majikNumber1
unsigned int	majikNumber2
unsigned int	majikNumber3
std::vector< std::string >	Namen
std::string	SimVtxLabel

Protected Member Functions
void	analyze (const edm::Event &e, const edm::EventSetup &es)
void	beginJob ()
void	endJob (void)
Protected Member Functions inherited from edm::EDAnalyzer
CurrentProcessingContext const *	currentContext () const

Private Attributes
edm::InputTag	cscRecHitLabel
DQMStore *	dbe
float	denCountMinus [50]
float	denCountPlus [50]
edm::InputTag	GMTInputTag
int	gtHaloBit
edm::InputTag	HLTriggerTag
edm::InputTag	lctProducer
edm::ESHandle< CSCGeometry >	m_cscGeometry
MonitorElement *	MinusEff
MonitorElement *	MinusEffDen
MonitorElement *	MinusEffNum
MonitorElement *	MinusEffProj3
MonitorElement *	MinusMe1BeamHaloOcc
MonitorElement *	MinusMe1BeamHaloOccRad
MonitorElement *	MinusMe1BeamHaloOccRing1
MonitorElement *	MinusMe1BeamHaloOccRing2
MonitorElement *	MinusMe1BeamHaloOccRing2or3
MonitorElement *	MinusMe2BeamHaloOcc
MonitorElement *	MinusMe2BeamHaloOccRing1
MonitorElement *	MinusMe2BeamHaloOccRing1Unfold
MonitorElement *	MinusMe2BeamHaloOccRing2
MonitorElement *	MinusMe2BeamHaloOccRing2or3
MonitorElement *	MinusMe2BeamHaloOccRing2Unfold
MonitorElement *	MinusMe3BeamHaloOcc
MonitorElement *	MinusMe3BeamHaloOccRing1
MonitorElement *	MinusMe3BeamHaloOccRing1Unfold
MonitorElement *	MinusMe3BeamHaloOccRing2
MonitorElement *	MinusMe3BeamHaloOccRing2or3
MonitorElement *	MinusMe3BeamHaloOccRing2Unfold
MonitorElement *	MinusMe4BeamHaloOcc
MonitorElement *	MinusMe4BeamHaloOccRing1
MonitorElement *	MinusMe4BeamHaloOccRing2
MonitorElement *	MinusMe4BeamHaloOccRing2or3
float	numCountMinus [50]
float	numCountPlus [50]
std::string	outFile
MonitorElement *	PlusEff
MonitorElement *	PlusEffDen
MonitorElement *	PlusEffNum
MonitorElement *	PlusEffProj3
MonitorElement *	PlusMe1BeamHaloOcc
MonitorElement *	PlusMe1BeamHaloOccRad
MonitorElement *	PlusMe1BeamHaloOccRing1
MonitorElement *	PlusMe1BeamHaloOccRing2
MonitorElement *	PlusMe1BeamHaloOccRing2or3
MonitorElement *	PlusMe2BeamHaloOcc
MonitorElement *	PlusMe2BeamHaloOccRing1
MonitorElement *	PlusMe2BeamHaloOccRing1Unfold
MonitorElement *	PlusMe2BeamHaloOccRing2
MonitorElement *	PlusMe2BeamHaloOccRing2or3
MonitorElement *	PlusMe2BeamHaloOccRing2Unfold
MonitorElement *	PlusMe3BeamHaloOcc
MonitorElement *	PlusMe3BeamHaloOccRing1
MonitorElement *	PlusMe3BeamHaloOccRing1Unfold
MonitorElement *	PlusMe3BeamHaloOccRing2
MonitorElement *	PlusMe3BeamHaloOccRing2or3
MonitorElement *	PlusMe3BeamHaloOccRing2Unfold
MonitorElement *	PlusMe4BeamHaloOcc
MonitorElement *	PlusMe4BeamHaloOccRing1
MonitorElement *	PlusMe4BeamHaloOccRing2
MonitorElement *	PlusMe4BeamHaloOccRing2or3

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
typedef WorkerT< EDAnalyzer >	WorkerType
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 28 of file HaloTrigger.h.

Constructor & Destructor Documentation

HaloTrigger::HaloTrigger

(

const edm::ParameterSet &

ps

)

Definition at line 49 of file HaloTrigger.cc.

References first, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and cmsCodeRules.cppFunctionSkipper::operator.

{
//  dbe = NULL;
//  if( ps.getUntrackedParameter<bool>("DQMStore", false) )
//  {
        dbe = Service<DQMStore>().operator->();
        dbe->setVerbose(0);
//  }

    first = false;
    outFile = ps.getUntrackedParameter<string>("outFile");
    
    HLTriggerTag = ps.getParameter< edm::InputTag >("HLTriggerTag");
    GMTInputTag = ps.getParameter< edm::InputTag >("GMTInputTag");
    lctProducer = ps.getParameter< edm::InputTag >("LCTInputTag");
    cscRecHitLabel = ps.getParameter< edm::InputTag >("RecInputTag");

}

HaloTrigger::~HaloTrigger ( )

virtual

Definition at line 68 of file HaloTrigger.cc.

{}

Member Function Documentation

void HaloTrigger::analyze ( const edm::Event &  e, const edm::EventSetup &  es  )

protectedvirtual

Implements edm::EDAnalyzer.

Definition at line 210 of file HaloTrigger.cc.

References PV3DBase< T, PVType, FrameType >::eta(), first, edm::EventSetup::get(), edm::Event::getByLabel(), i, j, n, L1TEmulatorMonitor_cff::p, PV3DBase< T, PVType, FrameType >::phi(), edm::Handle< T >::product(), CosmicsPD_Skims::radius, mathSSE::sqrt(), edm::TriggerNames::triggerNames(), edm::Event::triggerNames(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

{


    es.get<MuonGeometryRecord>().get(m_cscGeometry);
    
    edm::Handle<edm::SimVertexContainer> simvertices_handle;
    e.getByLabel("g4SimHits",simvertices_handle);
    edm::SimVertexContainer const* simvertices = simvertices_handle.product();
 
    edm::Handle<edm::SimTrackContainer> simtracks_handle;
    e.getByLabel("g4SimHits",simtracks_handle);
    edm::SimTrackContainer const* simtracks = simtracks_handle.product();

    float vtxArray[7];
    for(int jj=0; jj<7; jj++)   vtxArray[jj] = 0;
  edm::SimTrackContainer::const_iterator isimtr;
    for(isimtr=simtracks->begin(); isimtr!=simtracks->end(); isimtr++) {
        if( (*simvertices)[(*isimtr).vertIndex()].position().z() > 2299 ){
            vtxArray[1] = 1;
            vtxArray[2] = (*simvertices)[(*isimtr).vertIndex()].position().x();
            vtxArray[3] = (*simvertices)[(*isimtr).vertIndex()].position().y();
            vtxArray[4] = (*isimtr).momentum().px();
            vtxArray[5] = (*isimtr).momentum().py();
            vtxArray[6] = (*isimtr).momentum().pz();
        }
        
        if( (*simvertices)[(*isimtr).vertIndex()].position().z() < -2299 ){
            vtxArray[1] = -1;
            vtxArray[2] = (*simvertices)[(*isimtr).vertIndex()].position().x();
            vtxArray[3] = (*simvertices)[(*isimtr).vertIndex()].position().y();
            vtxArray[4] = (*isimtr).momentum().px();
            vtxArray[5] = (*isimtr).momentum().py();
            vtxArray[6] = (*isimtr).momentum().pz();
        }
    }

    if( HLTriggerTag.label() != "null" )
    {
        edm::Handle<TriggerResults> trh;
        e.getByLabel(HLTriggerTag,trh);

        // Initialize HLT Paths //
        if (!first)
        {
            first = true;
            edm::TriggerNames const& triggerNames = e.triggerNames(*trh);
            Namen = triggerNames.triggerNames();
            const unsigned int n(Namen.size());
            std::string searchName_Halo0 ("HLT_CSCBeamHalo");
            std::string searchName_Halo1 ("HLT_CSCBeamHaloOverlapRing1");
            std::string searchName_Halo2 ("HLT_CSCBeamHaloOverlapRing2");
            std::string searchName_Halo3 ("HLT_CSCBeamHaloRing2or3");
            for( unsigned int i=0; i!=n; ++i){
                int compy = Namen[i].compare(searchName_Halo0);
                if( compy == 0) majikNumber0 = i;
                
                compy = Namen[i].compare(searchName_Halo1);
                if( compy == 0) majikNumber1 = i;
                
                compy = Namen[i].compare(searchName_Halo2);
                if( compy == 0) majikNumber2 = i;
                
                compy = Namen[i].compare(searchName_Halo3);
                if( compy == 0) majikNumber3 = i;
            }// loop over all trigger indices searching for HLT of interest
        }
        
        // ID Halo Triggers //
        bool haloTriggerAny = false;
        bool halo0 = false, halo1 = false, halo2 = false, halo3 = false;
        const unsigned int n(Namen.size());
        for( unsigned int j=0; j!=n; ++j)
        {
            if( (j == majikNumber0) && (trh->accept(j)) )
            {
                halo0 = true;
                haloTriggerAny = true;
            }
            
            if( (j == majikNumber1) && (trh->accept(j)) )
            {
                halo1 = true;
                haloTriggerAny = true;
            }
            
            if( (j == majikNumber2) && (trh->accept(j)) )
            {
                halo2 = true;
                haloTriggerAny = true;
            }
            
            if( (j == majikNumber3) && (trh->accept(j)) )
            {
                halo3 = true;
                haloTriggerAny = true;
            }
        }// loop over all triggers in every event, check if path was run
        
        float radius = sqrt( vtxArray[2]*vtxArray[2] + vtxArray[3]*vtxArray[3] );
        
        if( vtxArray[1] == 1){
            float projX = ((vtxArray[4]/vtxArray[6])*(-1375)) + vtxArray[2];
            float projY = ((vtxArray[5]/vtxArray[6])*(-1375)) + vtxArray[3];
            float projRad = sqrt(projX*projX + projY*projY);
            int radInd = 0;
            for(int j=0; j<50; j++)
            { 
                if((projRad >= (j)*15.0) && (projRad < (j+1)*15)) radInd = j;
            }
            PlusEffDen->Fill(radius);
            if(haloTriggerAny == true){
                PlusEff->Fill(radius,1);
                PlusEffNum->Fill(radius);
                if((vtxArray[6] < 10) && (radInd != 0))
                {
                    numCountPlus[radInd]+=1.0;
                    denCountPlus[radInd]+=1.0;
                }
            } else {
                PlusEff->Fill(radius,0);
                if((vtxArray[6] < 50) && (radInd !=0))
                {
                    denCountPlus[radInd]+=1.0;
                }
            }
        }
        
        if( vtxArray[1] == -1){
            float projX = ((vtxArray[4]/vtxArray[6])*(1375)) + vtxArray[2];
            float projY = ((vtxArray[5]/vtxArray[6])*(1375)) + vtxArray[3];
            float projRad = sqrt(projX*projX + projY*projY);
            int radInd = 0;
            for(int j=0; j<50; j++) 
                if((projRad >= (j)*15.0) && (projRad < (j+1)*15)) radInd = j;
            MinusEffDen->Fill(radius);
            if(haloTriggerAny == true){
                MinusEff->Fill(radius,1);
                if((vtxArray[6] > 50) && (radInd != 0)){
                    numCountMinus[radInd]+=1.0;
                    denCountMinus[radInd]+=1.0;
                }
                MinusEffNum->Fill(radius);
            } else {
                MinusEff->Fill(radius,0);
                if((vtxArray[6] > 10) && (radInd != 0))
                {
                    denCountMinus[radInd]+=1.0;
                }
            }
        }
        
        // Rec Hit Info //
        if( haloTriggerAny == true )
        {
            float cscHit[9][8];
            for(int j = 1; j < 9; j++) cscHit[j][1] = 0;
            
            edm::Handle<CSCRecHit2DCollection> cscRecHits;
            e.getByLabel(cscRecHitLabel, cscRecHits);
            CSCRecHit2DCollection::const_iterator hit = cscRecHits->begin();
            for(; hit != cscRecHits->end(); hit++)
            {
                LocalPoint p = hit->localPosition();
                CSCDetId id((hit)->geographicalId().rawId());
                GlobalPoint gP = m_cscGeometry->idToDet(id)->toGlobal(hit->localPosition());
                int hitIndex = id.station() + 5*(id.endcap() - 1);
                if( cscHit[hitIndex][1] == 0 )
                {
                    cscHit[hitIndex][1] = 1;
                    cscHit[hitIndex][2] = gP.x();
                    cscHit[hitIndex][3] = gP.y();
                    cscHit[hitIndex][4] = gP.z();
                    cscHit[hitIndex][5] = sqrt( gP.x()*gP.x() + gP.y()*gP.y() );
                    float joe1 = ((gP.phi()+3.2)*180*100/(3.14));
                    int joe2 = int(joe1) % 1000;
                    float joe3 = (1.0*joe2)/100;
                    cscHit[hitIndex][6] = joe3;
                    cscHit[hitIndex][7] = gP.eta();
                }//cscHit[x][1] == 0
            }//2d rec hit loop
            
            // Fill Histograms //
            if(halo0 == true)
            {
                if(cscHit[1][1] == 1)
                { 
                    PlusMe1BeamHaloOcc->Fill(cscHit[1][2],cscHit[1][3]);
                    PlusMe1BeamHaloOccRad->Fill(cscHit[1][5]);
                }
                if(cscHit[2][1] == 1) PlusMe2BeamHaloOcc->Fill(cscHit[2][2],cscHit[2][3]);
                if(cscHit[3][1] == 1) PlusMe3BeamHaloOcc->Fill(cscHit[3][2],cscHit[3][3]);
                if(cscHit[4][1] == 1) PlusMe4BeamHaloOcc->Fill(cscHit[4][2],cscHit[4][3]);
                if(cscHit[5][1] == 1)
                {
                    MinusMe1BeamHaloOcc->Fill(cscHit[5][2],cscHit[5][3]);
                    MinusMe1BeamHaloOccRad->Fill(cscHit[5][5]);
                }
                if(cscHit[6][1] == 1) MinusMe2BeamHaloOcc->Fill(cscHit[6][2],cscHit[6][3]);
                if(cscHit[7][1] == 1) MinusMe3BeamHaloOcc->Fill(cscHit[7][2],cscHit[7][3]);
                if(cscHit[8][1] == 1) MinusMe4BeamHaloOcc->Fill(cscHit[8][2],cscHit[8][3]); 
        
            }//HLT_CSCBeamHalo
            
            if(halo1 == true)
            {
                if(cscHit[1][1] == 1) PlusMe1BeamHaloOccRing1->Fill(cscHit[1][2],cscHit[1][3]);
                if(cscHit[2][1] == 1)
                {
                    PlusMe2BeamHaloOccRing1->Fill(cscHit[2][2],cscHit[2][3]);
                    PlusMe2BeamHaloOccRing1Unfold->Fill(cscHit[2][6],cscHit[2][7]);
                }
                if(cscHit[3][1] == 1)
                {
                    PlusMe3BeamHaloOccRing1->Fill(cscHit[3][2],cscHit[3][3]);
                    PlusMe3BeamHaloOccRing1Unfold->Fill(cscHit[3][6],cscHit[3][7]);
                }
                if(cscHit[4][1] == 1) PlusMe4BeamHaloOccRing1->Fill(cscHit[4][2],cscHit[4][3]);
                if(cscHit[5][1] == 1) MinusMe1BeamHaloOccRing1->Fill(cscHit[5][2],cscHit[5][3]);
                if(cscHit[6][1] == 1)
                {
                    MinusMe2BeamHaloOccRing1->Fill(cscHit[6][2],cscHit[6][3]);
                    MinusMe2BeamHaloOccRing1Unfold->Fill(cscHit[6][6],cscHit[6][7]);
                }
                if(cscHit[7][1] == 1)
                {
                    MinusMe3BeamHaloOccRing1->Fill(cscHit[7][2],cscHit[7][3]);
                    MinusMe3BeamHaloOccRing1Unfold->Fill(cscHit[7][6],cscHit[7][7]);
                }
                if(cscHit[8][1] == 1) MinusMe4BeamHaloOccRing1->Fill(cscHit[8][2],cscHit[8][3]);
            }//HLT_CSCBeamHalo
            
            if(halo2 == true)
            {
                if(cscHit[1][1] == 1) PlusMe1BeamHaloOccRing2->Fill(cscHit[1][2],cscHit[1][3]);
                if(cscHit[2][1] == 1) 
                {
                    PlusMe2BeamHaloOccRing2->Fill(cscHit[2][2],cscHit[2][3]);
                    PlusMe2BeamHaloOccRing2Unfold->Fill(cscHit[2][6],cscHit[2][7]);
                }
                if(cscHit[3][1] == 1)
                {
                    PlusMe3BeamHaloOccRing2->Fill(cscHit[3][2],cscHit[3][3]);
                    PlusMe3BeamHaloOccRing2Unfold->Fill(cscHit[2][6],cscHit[2][7]);
                }
                if(cscHit[4][1] == 1) PlusMe4BeamHaloOccRing2->Fill(cscHit[4][2],cscHit[4][3]);
                if(cscHit[5][1] == 1) MinusMe1BeamHaloOccRing2->Fill(cscHit[5][2],cscHit[5][3]);
                if(cscHit[6][1] == 1)
                {
                    MinusMe2BeamHaloOccRing2->Fill(cscHit[6][2],cscHit[6][3]);
                    MinusMe2BeamHaloOccRing2Unfold->Fill(cscHit[6][6],cscHit[6][7]);
                }
                if(cscHit[7][1] == 1)
                {
                    MinusMe3BeamHaloOccRing2->Fill(cscHit[7][2],cscHit[7][3]);
                    MinusMe3BeamHaloOccRing2Unfold->Fill(cscHit[7][6],cscHit[7][7]);
                }
                if(cscHit[8][1] == 1) MinusMe4BeamHaloOccRing2->Fill(cscHit[8][2],cscHit[8][3]);
            }//HLT_CSCBeamHalo
            
            if(halo3 == true)
            {
                if(cscHit[1][1] == 1) PlusMe1BeamHaloOccRing2or3->Fill(cscHit[1][2],cscHit[1][3]);
                if(cscHit[2][1] == 1) PlusMe2BeamHaloOccRing2or3->Fill(cscHit[2][2],cscHit[2][3]);
                if(cscHit[3][1] == 1) PlusMe3BeamHaloOccRing2or3->Fill(cscHit[3][2],cscHit[3][3]);
                if(cscHit[4][1] == 1) PlusMe4BeamHaloOccRing2or3->Fill(cscHit[4][2],cscHit[4][3]);
                if(cscHit[5][1] == 1) MinusMe1BeamHaloOccRing2or3->Fill(cscHit[5][2],cscHit[5][3]);
                if(cscHit[6][1] == 1) MinusMe2BeamHaloOccRing2or3->Fill(cscHit[6][2],cscHit[6][3]);
                if(cscHit[7][1] == 1) MinusMe3BeamHaloOccRing2or3->Fill(cscHit[7][2],cscHit[7][3]);
                if(cscHit[8][1] == 1) MinusMe4BeamHaloOccRing2or3->Fill(cscHit[8][2],cscHit[8][3]);
            }//HLT_CSCBeamHalo
        }//haloTrigger true
    }//HLTriggerTag != null
}

void HaloTrigger::beginJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 71 of file HaloTrigger.cc.

References DQMStore::book1D(), DQMStore::book2D(), DQMStore::bookProfile(), i, cmsCodeRules.cppFunctionSkipper::operator, MonitorElement::setAxisTitle(), MonitorElement::setBinLabel(), and DQMStore::setCurrentFolder().

{
    // Initialize Counters //
    hltHaloTriggers = 0;
    hltHaloOver1 = 0;
    hltHaloOver2 = 0;
    hltHaloRing23 = 0;
    CscHalo_Gmt = 0;    
    
    // DQM Directory Structure //
    DQMStore * dbe = 0;
    dbe = Service<DQMStore>().operator->();
    if( dbe ){
        dbe->setCurrentFolder("HLTriggerOffline/special");
        
        // Define Monitor Elements a.k.a. Plots //
        PlusMe1BeamHaloOcc = dbe->book2D("PlusMe1BeamHaloOcc","Positive Endcap, Me1 Beam Halo Occupancy",1000,-500,500,1000,-500,500);
        PlusMe1BeamHaloOccRing1 = dbe->book2D("PlusMe1BeamHaloOccRing1","Positive Endcap, Me1 Beam Halo Occupancy, HLT Trig Ring1",1000,-500,500,1000,-500,500);
        PlusMe1BeamHaloOccRing2 = dbe->book2D("PlusMe1BeamHaloOccRing2","Positive Endcap, Me1 Beam Halo Occupancy, HLT Trig Ring2",1000,-500,500,1000,-500,500);
        PlusMe1BeamHaloOccRing2or3 = dbe->book2D("PlusMe1BeamHaloOccRing2or3","Positive Endcap, Me1 Beam Halo Occupancy, HLT Trig Ring 2 or 3",1000,-500,500,1000,-500,500);
        PlusMe1BeamHaloOccRad = dbe->book1D("PlusMe1BeamHaloOccRad","Positive Endcap, Me1 Beam Halo Radial Occupancy",750,0,750);
    
        PlusMe2BeamHaloOcc = dbe->book2D("PlusMe2BeamHaloOcc","Positive Endcap, Me2 Beam Halo Occupancy",1000,-500,500,1000,-500,500);
        PlusMe2BeamHaloOccRing1 = dbe->book2D("PlusMe2BeamHaloOccRing1","Positive Endcap, Me2 Beam Halo Occupancy, HLT Trig Ring1",1000,-500,500,1000,-500,500);
        PlusMe2BeamHaloOccRing2 = dbe->book2D("PlusMe2BeamHaloOccRing2","Positive Endcap, Me2 Beam Halo Occupancy, HLT Trig Ring2",1000,-500,500,1000,-500,500);
        PlusMe2BeamHaloOccRing2or3 = dbe->book2D("PlusMe2BeamHaloOccRing2or3","Positive Endcap, Me2 Beam Halo Occupancy, HLT Trig Ring 2 or 3",1000,-500,500,1000,-500,500);
    
        PlusMe3BeamHaloOcc = dbe->book2D("PlusMe3BeamHaloOcc","Positive Endcap, Me3 Beam Halo Occupancy",1000,-500,500,1000,-500,500);
        PlusMe3BeamHaloOccRing1 = dbe->book2D("PlusMe3BeamHaloOccRing1","Positive Endcap, Me3 Beam Halo Occupancy, HLT Trig Ring1",1000,-500,500,1000,-500,500);
        PlusMe3BeamHaloOccRing2 = dbe->book2D("PlusMe3BeamHaloOccRing2","Positive Endcap, Me3 Beam Halo Occupancy, HLT Trig Ring2",1000,-500,500,1000,-500,500);
        PlusMe3BeamHaloOccRing2or3 = dbe->book2D("PlusMe3BeamHaloOccRing2or3","Positive Endcap, Me3 Beam Halo Occupancy, HLT Trig Ring 2 or 3",1000,-500,500,1000,-500,500);
    
        PlusMe4BeamHaloOcc = dbe->book2D("PlusMe4BeamHaloOcc","Positive Endcap, Me4 Beam Halo Occupancy",1000,-500,500,1000,-500,500);
        PlusMe4BeamHaloOccRing1 = dbe->book2D("PlusMe4BeamHaloOccRing1","Positive Endcap, Me4 Beam Halo Occupancy, HLT Trig Ring1",1000,-500,500,1000,-500,500);
        PlusMe4BeamHaloOccRing2 = dbe->book2D("PlusMe4BeamHaloOccRing2","Positive Endcap, Me4 Beam Halo Occupancy, HLT Trig Ring2",1000,-500,500,1000,-500,500);
        PlusMe4BeamHaloOccRing2or3 = dbe->book2D("PlusMe4BeamHaloOccRing2or3","Positive Endcap, Me4 Beam Halo Occupancy, HLT Trig Ring 2 or 3",1000,-500,500,1000,-500,500);
    
        MinusMe1BeamHaloOcc = dbe->book2D("MinusMe1Me1BeamHaloOcc","Negative Endcap, Me1 Beam Halo Occupancy",1000,-500,500,1000,-500,500);
        MinusMe1BeamHaloOccRing1 = dbe->book2D("MinusBeamHaloOccRing1","Negative Endcap, Me1 Beam Halo Occupancy, HLT Trig Ring1",1000,-500,500,1000,-500,500);
        MinusMe1BeamHaloOccRing2 = dbe->book2D("MinusMe1BeamHaloOccRing2","Negative Endcap, Me1 Beam Halo Occupancy, HLT Trig Ring2",1000,-500,500,1000,-500,500);
        MinusMe1BeamHaloOccRing2or3 = dbe->book2D("MinusMe1BeamHaloOccRing2or3","Negative Endcap, Me1 Beam Halo Occupancy, HLT Trig Ring 2 or 3",1000,-500,500,1000,-500,500);
        MinusMe1BeamHaloOccRad = dbe->book1D("MinusMe1BeamHaloOccRad","Minus Endcap, Me1 Beam Halo Radial Occupancy",750,0,750);
        
        MinusMe2BeamHaloOcc = dbe->book2D("MinusMe2BeamHaloOcc","Negative Endcap, Me2 Beam Halo Occupancy",1000,-500,500,1000,-500,500);
        MinusMe2BeamHaloOccRing1 = dbe->book2D("MinusMe2BeamHaloOccRing1","Negative Endcap, Me2 Beam Halo Occupancy, HLT Trig Ring1",1000,-500,500,1000,-500,500);
        MinusMe2BeamHaloOccRing2 = dbe->book2D("MinusMe2BeamHaloOccRing2","Negative Endcap, Me2 Beam Halo Occupancy, HLT Trig Ring2",1000,-500,500,1000,-500,500);
        MinusMe2BeamHaloOccRing2or3 = dbe->book2D("MinusMe2BeamHaloOccRing2or3","Negative Endcap, Me2 Beam Halo Occupancy, HLT Trig Ring 2 or 3",1000,-500,500,1000,-500,500);
        
        MinusMe3BeamHaloOcc = dbe->book2D("MinusMe3BeamHaloOcc","Negative Endcap, Me3 Beam Halo Occupancy",1000,-500,500,1000,-500,500);
        MinusMe3BeamHaloOccRing1 = dbe->book2D("MinusMe3BeamHaloOccRing1","Negative Endcap, Me3 Beam Halo Occupancy, HLT Trig Ring1",1000,-500,500,1000,-500,500);
        MinusMe3BeamHaloOccRing2 = dbe->book2D("MinusMe3BeamHaloOccRing2","Negative Endcap, Me3 Beam Halo Occupancy, HLT Trig Ring2",1000,-500,500,1000,-500,500);
        MinusMe3BeamHaloOccRing2or3 = dbe->book2D("MinusMe3BeamHaloOccRing2or3","Negative Endcap, Me3 Beam Halo Occupancy, HLT Trig Ring 2 or 3",1000,-500,500,1000,-500,500);
        
        MinusMe4BeamHaloOcc = dbe->book2D("MinusMe4BeamHaloOcc","Negative Endcap, Me4 Beam Halo Occupancy",1000,-500,500,1000,-500,500);
        MinusMe4BeamHaloOccRing1 = dbe->book2D("MinusMe4BeamHaloOccRing1","Negative Endcap, Me4 Beam Halo Occupancy, HLT Trig Ring1",1000,-500.,500.,1000,-500.,500.);
        MinusMe4BeamHaloOccRing2 = dbe->book2D("MinusMe4BeamHaloOccRing2","Negative Endcap, Me4 Beam Halo Occupancy, HLT Trig Ring2",1000,-500,500,1000,-500,500);
        MinusMe4BeamHaloOccRing2or3 = dbe->book2D("MinusMe4BeamHaloOccRing2or3","Negative Endcap, Me4 Beam Halo Occupancy, HLT Trig Ring 2 or 3",1000,-500,500,1000,-500,500);
    
        MinusMe3BeamHaloOccRing2Unfold = dbe->book2D("MinusMe3BeamHaloOccRing2Unfold","Negative Endcap Me3, Ring 2 Unfolded",120,-1,11,500,-2.4,-0.9);
        MinusMe2BeamHaloOccRing2Unfold = dbe->book2D("MinusMe2BeamHaloOccRing2Unfold","Negative Endcap Me2, Ring 2 Unfolded",120,-1,11,500,-2.4,-0.9);
        PlusMe3BeamHaloOccRing2Unfold = dbe->book2D("PlusMe3BeamHaloOccRing2Unfold","Negative Endcap Me3, Ring 2 Unfolded",120,-1,11,500,0.9,2.4);
        PlusMe2BeamHaloOccRing2Unfold = dbe->book2D("PlusMe2BeamHaloOccRing2Unfold","Negative Endcap Me2, Ring 2 Unfolded",120,-1,11,500,0.9,2.4);
        
        MinusMe3BeamHaloOccRing1Unfold = dbe->book2D("MinusMe3BeamHaloOccRing1Unfold","Negative Endcap Me3, Ring 1 Unfolded",120,-1,11,500,-2.5,-0.9);
        MinusMe2BeamHaloOccRing1Unfold = dbe->book2D("MinusMe2BeamHaloOccRing1Unfold","Negative Endcap Me2, Ring 1 Unfolded",120,-1,11,500,-2.5,-0.9);
        PlusMe3BeamHaloOccRing1Unfold = dbe->book2D("PlusMe3BeamHaloOccRing1Unfold","Negative Endcap Me3, Ring 1 Unfolded",120,-1,11,500,0.9,2.5);
        PlusMe2BeamHaloOccRing1Unfold = dbe->book2D("PlusMe2BeamHaloOccRing1Unfold","Negative Endcap Me2, Ring 1 Unfolded",120,-1,11,500,0.9,2.5);
    
        MinusEff = dbe->bookProfile("MinusEff", "Minus Endcap Efficiency by Vertex",50,0,750,50,0,1);
        MinusEff->setAxisTitle("Radius",1);
        MinusEff->setAxisTitle("Efficiency",2);
        //MinusEffProj3 = dbe->bookProfile("MinusEffProj3", "Minus Endcap Efficiency Projected to statoin 3",50,0,750,50,0,1);
        MinusEffProj3 = dbe->book1D("MinusEffProj3","Minus Endcap Efficiency by projection",50,0,750);
        MinusEffProj3->setBinLabel(1,"0",1);
        MinusEffProj3->setBinLabel(10,"150",1);
        MinusEffProj3->setBinLabel(20,"300",1);
        MinusEffProj3->setBinLabel(30,"450",1);
        MinusEffProj3->setBinLabel(40,"600",1);
        MinusEffProj3->setBinLabel(50,"750",1);
        MinusEffProj3->setAxisTitle("Radius",1);
        MinusEffProj3->setAxisTitle("Efficiency",2);
        MinusEffNum = dbe->book1D("MinusEffNum", "Minus Endcap numerator",50,0,750);
        MinusEffDen = dbe->book1D("MinusEffDen", "Minus Endcap denomenator", 50, 0, 750);
        PlusEff  = dbe->bookProfile("PlusEff",  "Plus Endcap Efficiency by Vertex",50,0,750,50,0,1);
        PlusEff->setAxisTitle("Radius",1);
        PlusEff->setAxisTitle("Efficiency",2);
        //PlusEffProj3  = dbe->bookProfile("PlusEffProj3",  "Plus Endcap Efficiency Projected to station 3",50,0,750,50,0,1);
        PlusEffProj3 = dbe->book1D("PlusEffProj3","Plus Endcap Efficiency by projection",50,0,750);
        PlusEffProj3->setAxisTitle("Radius",1);
        PlusEffProj3->setAxisTitle("Efficiency",2);
        PlusEffProj3->setBinLabel(1,"0",1);
        PlusEffProj3->setBinLabel(10,"150",1);
        PlusEffProj3->setBinLabel(20,"300",1);
        PlusEffProj3->setBinLabel(30,"450",1);
        PlusEffProj3->setBinLabel(40,"600",1);
        PlusEffProj3->setBinLabel(50,"750",1);
        PlusEffNum = dbe->book1D("PlusEffNum", "Plus Endcap numerator",50,0,750);
        PlusEffDen = dbe->book1D("PlusEffDen", "Plus Endcap denomenator", 50, 0, 750);
    }
    
    //es.get<MuonGeometryRecord>().get(m_cscGeometry);
    for( int i=0; i<50; ++i)
    {
        numCountPlus[i] = 0;
        denCountPlus[i] = 0;
        numCountMinus[i] = 0;
        denCountMinus[i] = 0;
    }
}

void HaloTrigger::endJob ( void  )

protectedvirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 187 of file HaloTrigger.cc.

References j.

{
    for(int j=1;j<51;j++)
    {
        if(denCountPlus[j] != 0)
        {
            float fillValuePlus = numCountPlus[j]/denCountPlus[j];
            PlusEffProj3->setBinContent(j,fillValuePlus);
        }
        
        if(denCountMinus[j] != 0)
        {
            float fillValueMinus = numCountMinus[j]/denCountMinus[j];
            MinusEffProj3->setBinContent(j,fillValueMinus);
        }
    }
        


    dbe->save(outFile); 
    return;
}

Member Data Documentation

unsigned int HaloTrigger::CscHalo_Gmt

Definition at line 38 of file HaloTrigger.h.

edm::InputTag HaloTrigger::cscRecHitLabel

private

Definition at line 109 of file HaloTrigger.h.

DQMStore* HaloTrigger::dbe

private

Definition at line 47 of file HaloTrigger.h.

float HaloTrigger::denCountMinus[50]

private

Definition at line 105 of file HaloTrigger.h.

float HaloTrigger::denCountPlus[50]

private

Definition at line 103 of file HaloTrigger.h.

bool HaloTrigger::first

Definition at line 36 of file HaloTrigger.h.

edm::InputTag HaloTrigger::GMTInputTag

private

Definition at line 109 of file HaloTrigger.h.

int HaloTrigger::gtHaloBit

private

Definition at line 112 of file HaloTrigger.h.

unsigned int HaloTrigger::hltHaloOver1

Definition at line 38 of file HaloTrigger.h.

unsigned int HaloTrigger::hltHaloOver2

Definition at line 38 of file HaloTrigger.h.

unsigned int HaloTrigger::hltHaloRing23

Definition at line 38 of file HaloTrigger.h.

unsigned int HaloTrigger::hltHaloTriggers

Definition at line 38 of file HaloTrigger.h.

edm::InputTag HaloTrigger::HLTriggerTag

private

Definition at line 109 of file HaloTrigger.h.

edm::InputTag HaloTrigger::lctProducer

private

Definition at line 109 of file HaloTrigger.h.

edm::ESHandle<CSCGeometry> HaloTrigger::m_cscGeometry

private

Definition at line 108 of file HaloTrigger.h.

unsigned int HaloTrigger::majikNumber0

Definition at line 39 of file HaloTrigger.h.

unsigned int HaloTrigger::majikNumber1

Definition at line 39 of file HaloTrigger.h.

unsigned int HaloTrigger::majikNumber2

Definition at line 39 of file HaloTrigger.h.

unsigned int HaloTrigger::majikNumber3

Definition at line 39 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusEff

private

Definition at line 94 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusEffDen

private

Definition at line 100 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusEffNum

private

Definition at line 99 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusEffProj3

private

Definition at line 96 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe1BeamHaloOcc

private

Definition at line 66 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe1BeamHaloOccRad

private

Definition at line 82 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe1BeamHaloOccRing1

private

Definition at line 67 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe1BeamHaloOccRing2

private

Definition at line 68 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe1BeamHaloOccRing2or3

private

Definition at line 69 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe2BeamHaloOcc

private

Definition at line 70 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe2BeamHaloOccRing1

private

Definition at line 71 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe2BeamHaloOccRing1Unfold

private

Definition at line 89 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe2BeamHaloOccRing2

private

Definition at line 72 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe2BeamHaloOccRing2or3

private

Definition at line 73 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe2BeamHaloOccRing2Unfold

private

Definition at line 85 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe3BeamHaloOcc

private

Definition at line 74 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe3BeamHaloOccRing1

private

Definition at line 75 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe3BeamHaloOccRing1Unfold

private

Definition at line 88 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe3BeamHaloOccRing2

private

Definition at line 76 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe3BeamHaloOccRing2or3

private

Definition at line 77 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe3BeamHaloOccRing2Unfold

private

Definition at line 84 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe4BeamHaloOcc

private

Definition at line 78 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe4BeamHaloOccRing1

private

Definition at line 79 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe4BeamHaloOccRing2

private

Definition at line 80 of file HaloTrigger.h.

MonitorElement* HaloTrigger::MinusMe4BeamHaloOccRing2or3

private

Definition at line 81 of file HaloTrigger.h.

std::vector<std::string> HaloTrigger::Namen

Definition at line 37 of file HaloTrigger.h.

float HaloTrigger::numCountMinus[50]

private

Definition at line 104 of file HaloTrigger.h.

float HaloTrigger::numCountPlus[50]

private

Definition at line 102 of file HaloTrigger.h.

std::string HaloTrigger::outFile

private

Definition at line 111 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusEff

private

Definition at line 93 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusEffDen

private

Definition at line 98 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusEffNum

private

Definition at line 97 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusEffProj3

private

Definition at line 95 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe1BeamHaloOcc

private

Definition at line 48 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe1BeamHaloOccRad

private

Definition at line 64 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe1BeamHaloOccRing1

private

Definition at line 49 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe1BeamHaloOccRing2

private

Definition at line 50 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe1BeamHaloOccRing2or3

private

Definition at line 51 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe2BeamHaloOcc

private

Definition at line 52 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe2BeamHaloOccRing1

private

Definition at line 53 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe2BeamHaloOccRing1Unfold

private

Definition at line 91 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe2BeamHaloOccRing2

private

Definition at line 54 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe2BeamHaloOccRing2or3

private

Definition at line 55 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe2BeamHaloOccRing2Unfold

private

Definition at line 87 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe3BeamHaloOcc

private

Definition at line 56 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe3BeamHaloOccRing1

private

Definition at line 57 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe3BeamHaloOccRing1Unfold

private

Definition at line 90 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe3BeamHaloOccRing2

private

Definition at line 58 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe3BeamHaloOccRing2or3

private

Definition at line 59 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe3BeamHaloOccRing2Unfold

private

Definition at line 86 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe4BeamHaloOcc

private

Definition at line 60 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe4BeamHaloOccRing1

private

Definition at line 61 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe4BeamHaloOccRing2

private

Definition at line 62 of file HaloTrigger.h.

MonitorElement* HaloTrigger::PlusMe4BeamHaloOccRing2or3

private

Definition at line 63 of file HaloTrigger.h.

std::string HaloTrigger::SimVtxLabel

Definition at line 34 of file HaloTrigger.h.