---

DTDigiReader Class Reference

Analyse the the muon-drift-tubes digitizer. More...

Inheritance diagram for DTDigiReader:

List of all members.

Public Member Functions
void	analyze (const Event &event, const EventSetup &eventSetup)
	DTDigiReader (const ParameterSet &pset)
virtual	~DTDigiReader ()
Private Attributes
TH1F *	DigiTimeBox
TH1F *	DigiTimeBoxW0
TH1F *	DigiTimeBoxW1
TH1F *	DigiTimeBoxW2
TFile *	file
string	label

---

Detailed Description

Analyse the the muon-drift-tubes digitizer.

Date

2006/08/01 15:24:06

Revision

1.5

Authors:

: R. Bellan

Definition at line 30 of file DTDigiReader.cc.

---

Constructor & Destructor Documentation

DTDigiReader::DTDigiReader

(

const ParameterSet &

pset

)

[inline, explicit]

Definition at line 33 of file DTDigiReader.cc.

References GenMuonPlsPt100GeV_cfg::cout, DigiTimeBox, DigiTimeBoxW0, DigiTimeBoxW1, DigiTimeBoxW2, lat::endl(), file, edm::ParameterSet::getUntrackedParameter(), and label.

                                                 {
    file = new TFile("DTDigiPlots.root","RECREATE");
    file->cd();
    DigiTimeBox = new TH1F("DigiTimeBox","Digi Time Box",2048,0,1600);
    DigiTimeBoxW0 = new TH1F("DigiTimeBoxW0","Digi Time Box W0",2000,0,1600);
    DigiTimeBoxW1 = new TH1F("DigiTimeBoxW1","Digi Time Box W1",2000,0,1600);
    DigiTimeBoxW2 = new TH1F("DigiTimeBoxW2","Digi Time Box W2",2000,0,1600);
    if(file->IsOpen()) cout<<"file open!"<<endl;
    else cout<<"*** Error in opening file ***"<<endl;
    label = pset.getUntrackedParameter<string>("label");
  }

virtual DTDigiReader::~DTDigiReader

(

)

[inline, virtual]

Definition at line 45 of file DTDigiReader.cc.

References DigiTimeBox, DigiTimeBoxW0, DigiTimeBoxW1, DigiTimeBoxW2, and file.

                         {
    file->cd();
    DigiTimeBox->Write();
    DigiTimeBoxW0->Write();
    DigiTimeBoxW1->Write();
    DigiTimeBoxW2->Write();
    file->Close();
    //    delete file;
    // delete DigiTimeBox;
  }

---

Member Function Documentation

void DTDigiReader::analyze	(	const Event &	event,
		const EventSetup &	eventSetup
	)			[inline, virtual]

Implements edm::EDAnalyzer.

Definition at line 56 of file DTDigiReader.cc.

References funct::abs(), GenMuonPlsPt100GeV_cfg::cout, DigiTimeBox, DigiTimeBoxW0, DigiTimeBoxW1, DigiTimeBoxW2, lat::endl(), label, range, and trackerHits::simHits.

                                                                 {
    cout << "--- Run: " << event.id().run()
         << " Event: " << event.id().event() << endl;
    
    Handle<DTDigiCollection> dtDigis;
    event.getByLabel(label, dtDigis);
     // event.getByLabel("MuonDTDigis", dtDigis);
    Handle<PSimHitContainer> simHits; 
    event.getByLabel("g4SimHits","MuonDTHits",simHits);    


    DTDigiCollection::DigiRangeIterator detUnitIt;
    for (detUnitIt=dtDigis->begin();
         detUnitIt!=dtDigis->end();
         ++detUnitIt){
      
      const DTLayerId& id = (*detUnitIt).first;
      const DTDigiCollection::Range& range = (*detUnitIt).second;
      
      // DTLayerId print-out
      cout<<"--------------"<<endl;
      cout<<"id: "<<id;
      
      // Loop over the digis of this DetUnit
      for (DTDigiCollection::const_iterator digiIt = range.first;
           digiIt!=range.second;
           ++digiIt){
        //      if((*digiIt).time()<703 &&(*digiIt).time()>699) {
          cout<<" Wire: "<<(*digiIt).wire()<<endl
              <<" digi time (ns): "<<(*digiIt).time()<<endl;
          
          for(vector<PSimHit>::const_iterator simHit = simHits->begin();
              simHit != simHits->end(); simHit++){
            DTWireId wireId((*simHit).detUnitId());
            if (wireId.layerId()==id && abs((*simHit).particleType())==13){
              cout<<"entry: "<<(*simHit).entryPoint()<<endl
                  <<"exit: "<<(*simHit).exitPoint()<<endl
                  <<"TOF: "<<(*simHit).timeOfFlight()<<endl;
            }
          }

          //    }

        if(id.layer()==3)
          DigiTimeBoxW0->Fill((*digiIt).time());
        else if(abs(id.superlayer())==1)
          DigiTimeBoxW1->Fill((*digiIt).time());
        else if(abs(id.superlayer())==2)
          DigiTimeBoxW2->Fill((*digiIt).time());
        else cout<<"Error"<<endl;
        DigiTimeBox->Fill((*digiIt).time());
        
      }// for digis in layer
    }// for layers
    cout<<"--------------"<<endl;
  }

---

Member Data Documentation

TH1F* DTDigiReader::DigiTimeBox [private]

Definition at line 115 of file DTDigiReader.cc.

Referenced by analyze(), DTDigiReader(), and ~DTDigiReader().

TH1F* DTDigiReader::DigiTimeBoxW0 [private]

Definition at line 116 of file DTDigiReader.cc.

Referenced by analyze(), DTDigiReader(), and ~DTDigiReader().

TH1F* DTDigiReader::DigiTimeBoxW1 [private]

Definition at line 117 of file DTDigiReader.cc.

Referenced by analyze(), DTDigiReader(), and ~DTDigiReader().

TH1F* DTDigiReader::DigiTimeBoxW2 [private]

Definition at line 118 of file DTDigiReader.cc.

Referenced by analyze(), DTDigiReader(), and ~DTDigiReader().

TFile* DTDigiReader::file [private]

Definition at line 119 of file DTDigiReader.cc.

Referenced by DTDigiReader(), and ~DTDigiReader().

string DTDigiReader::label [private]

Definition at line 114 of file DTDigiReader.cc.

Referenced by analyze(), and DTDigiReader().

---

The documentation for this class was generated from the following file:

• SimMuon/DTDigitizer/src/DTDigiReader.cc

---