---

SubsystemNeutronWriter Class Reference

theSubsystemNeutronWriter stores "events" which consist of a list of SimHits, grouped by detector type. More...

#include <SimMuon/Neutron/interface/SubsystemNeutronWriter.h>

Inheritance diagram for SubsystemNeutronWriter:

List of all members.

Public Member Functions
virtual void	analyze (edm::Event const &e, edm::EventSetup const &c)
virtual int	chamberId (int globalDetId) const =0
virtual int	chamberType (int globalDetId) const =0
void	initialize (int chamberType)
	good practice to do once for each chamber type
virtual int	localDetId (int globalDetId) const =0
void	printStats ()
	SubsystemNeutronWriter (edm::ParameterSet const &pset)
virtual	~SubsystemNeutronWriter ()
	destructor prints statistics on number of events written
Protected Member Functions
void	adjust (PSimHit &h, float timeOffset)
	helper to add time offsets and local det ID
void	updateCount (int chamberType)
	updates the counter
virtual void	writeHits (int chamberType, edm::PSimHitContainer &allSimHits)
Private Attributes
bool	initialized
std::map< int, int >	theCountPerChamberType
NeutronWriter *	theHitWriter
edm::InputTag	theInputTag
double	theNeutronTimeCut
int	theNEvents
double	theTimeWindow

---

Detailed Description

theSubsystemNeutronWriter stores "events" which consist of a list of SimHits, grouped by detector type.

These can then be read back to model neutron background int muon chambers.

You can specify the cut on how long after the signal event to define something as a Neutron Event with the configurable Muon:NeutronTimeCut

Definition at line 25 of file SubsystemNeutronWriter.h.

---

Constructor & Destructor Documentation

SubsystemNeutronWriter::SubsystemNeutronWriter

(

edm::ParameterSet const &

pset

)

[explicit]

Definition at line 18 of file SubsystemNeutronWriter.cc.

References Exception, edm::ParameterSet::getParameter(), output(), and theHitWriter.

: theInputTag(pset.getParameter<edm::InputTag>("input")),
  theNeutronTimeCut(pset.getParameter<double>("neutronTimeCut")),
  theTimeWindow(pset.getParameter<double>("timeWindow")),
  theNEvents(0),
  initialized(false)
{
  string writer = pset.getParameter<string>("writer");
  string output = pset.getParameter<string>("output");
  if(writer == "ASCII")
  {
    theHitWriter = new AsciiNeutronWriter(output);
  }
  else if (writer == "ROOT")
  {
    theHitWriter = new RootNeutronWriter(output);
  }
  else 
  {
    throw cms::Exception("NeutronWriter") << "Bad writer: "
      << writer;
  } 
}

SubsystemNeutronWriter::~SubsystemNeutronWriter

(

)

[virtual]

destructor prints statistics on number of events written

Definition at line 43 of file SubsystemNeutronWriter.cc.

References printStats(), and theHitWriter.

{
  printStats();
  delete theHitWriter;
}

---

Member Function Documentation

void SubsystemNeutronWriter::adjust	(	PSimHit &	h,
		float	timeOffset
	)			[protected]

helper to add time offsets and local det ID

Definition at line 134 of file SubsystemNeutronWriter.cc.

References PSimHit::detUnitId(), PSimHit::energyLoss(), PSimHit::entryPoint(), PSimHit::exitPoint(), localDetId(), PSimHit::momentumAtEntry(), PSimHit::pabs(), PSimHit::particleType(), PV3DBase< T, PVType, FrameType >::phi(), PV3DBase< T, PVType, FrameType >::theta(), PSimHit::timeOfFlight(), and PSimHit::trackId().

Referenced by writeHits().

                                                                 {
  
  h = PSimHit( h.entryPoint(), h.exitPoint(), h.pabs(), 
               h.timeOfFlight() + timeOffset,
               h.energyLoss(), h.particleType(), 
               localDetId(h.detUnitId()), h.trackId(),
               h.momentumAtEntry().theta(),
               h.momentumAtEntry().phi() );
}

void SubsystemNeutronWriter::analyze	(	edm::Event const &	e,
		edm::EventSetup const &	c
	)			[virtual]

Implements edm::EDAnalyzer.

Definition at line 60 of file SubsystemNeutronWriter.cc.

References chamberId(), chamberType(), edm::Event::getByLabel(), theInputTag, theNEvents, and writeHits().

{
  ++theNEvents;
  edm::Handle<edm::PSimHitContainer> hits;
  e.getByLabel(theInputTag, hits);

  // sort hits by chamber
  map<int, edm::PSimHitContainer> hitsByChamber;
  for(edm::PSimHitContainer::const_iterator hitItr = hits->begin();
      hitItr != hits->end(); ++hitItr)
  {
    int chamberIndex = chamberId(hitItr->detUnitId());
    hitsByChamber[chamberIndex].push_back(*hitItr);
  }

  // now write out each chamber's contents
  for(map<int, edm::PSimHitContainer>::iterator hitsByChamberItr = hitsByChamber.begin();
      hitsByChamberItr != hitsByChamber.end(); ++hitsByChamberItr)
  {
    int chambertype = chamberType(hitsByChamberItr->first);
    writeHits(chambertype, hitsByChamberItr->second);
  }
}

virtual int SubsystemNeutronWriter::chamberId

(

int

globalDetId

)

const [pure virtual]

Implemented in CSCNeutronWriter.

Referenced by analyze().

virtual int SubsystemNeutronWriter::chamberType

(

int

globalDetId

)

const [pure virtual]

Implemented in CSCNeutronWriter.

Referenced by analyze().

void SubsystemNeutronWriter::initialize

(

int

chamberType

)

good practice to do once for each chamber type

Definition at line 85 of file SubsystemNeutronWriter.cc.

References NeutronWriter::initialize(), and theHitWriter.

Referenced by CSCNeutronWriter::CSCNeutronWriter().

{
  // should instantiate one of every chamber type, just so
  // ROOT knows what file to associate them with
  theHitWriter->initialize(chamberType);
}

virtual int SubsystemNeutronWriter::localDetId

(

int

globalDetId

)

const [pure virtual]

Implemented in CSCNeutronWriter.

Referenced by adjust().

void SubsystemNeutronWriter::printStats

(

)

Definition at line 50 of file SubsystemNeutronWriter.cc.

References theCountPerChamberType, and theNEvents.

Referenced by ~SubsystemNeutronWriter().

                                        {
  edm::LogInfo("SubsystemNeutronWriter") << "SubsystemNeutronWriter Statistics:\n";
  for(map<int,int>::iterator mapItr = theCountPerChamberType.begin();
      mapItr != theCountPerChamberType.end();  ++mapItr) {
     edm::LogInfo("SubsystemNeutronWriter") << "   theEventOccupancy[" << mapItr->first << "] = "
         << mapItr->second << " / " << theNEvents << " / NCT \n";
  }
}

void SubsystemNeutronWriter::updateCount

(

int

chamberType

)

[protected]

updates the counter

Definition at line 145 of file SubsystemNeutronWriter.cc.

References theCountPerChamberType.

Referenced by writeHits().

                                                        {
  map<int,int>::iterator entry = theCountPerChamberType.find(chamberType);
  if(entry == theCountPerChamberType.end()) {
    theCountPerChamberType.insert( pair<int,int>(chamberType, 1) );
  } else {
    ++(entry->second);
  }
}

void SubsystemNeutronWriter::writeHits	(	int	chamberType,
		edm::PSimHitContainer &	allSimHits
	)			[protected, virtual]

Definition at line 93 of file SubsystemNeutronWriter.cc.

References adjust(), GenMuonPlsPt100GeV_cfg::cout, PSimHit::detUnitId(), lat::endl(), i, LogDebug, PSimHit::pabs(), PSimHit::particleType(), python::multivaluedict::sort(), theHitWriter, theNeutronTimeCut, theTimeWindow, PSimHit::tof(), updateCount(), and NeutronWriter::writeEvent().

Referenced by analyze().

{

  sort(input.begin(), input.end(), SortByTime());
  edm::PSimHitContainer current;
  float startTime = -1000.;
  for(size_t i = 0; i < input.size(); ++i) {
    PSimHit hit = input[i];
std::cout << hit << std::endl;
    float tof = hit.tof();
    LogDebug("SubsystemNeutronWriter") << "found hit from part type " << hit.particleType()
                   << " at tof " << tof << " p " << hit.pabs() 
                   << " on det " << hit.detUnitId() 
                   << " chamber type " << chamberType;
    if(tof > theNeutronTimeCut) {
      if(tof > (startTime + theTimeWindow) ) { // 1st in cluster
        startTime = tof;
        if(!current.empty()) {
          LogDebug("SubsystemNeutronWriter") << "filling old cluster";
          theHitWriter->writeEvent(chamberType, current);
          updateCount(chamberType);
          current.clear();
        }
        LogDebug("SubsystemNeutronWriter") << "starting neutron cluster at time " << startTime 
          << " on detType " << chamberType;
      }
      // set the time to be 0 at start of event
std::cout << "ADJUST " << startTime << std::endl;
      adjust(hit, -1.*startTime);
      current.push_back( hit );
std::cout << "NEXT HIT" << std::endl;
    }
  }
std::cout << "LOOPED OVER HITS " << theHitWriter << std::endl;
  if(!current.empty()) {
    theHitWriter->writeEvent(chamberType, current);
    updateCount(chamberType);
  }
}

---

Member Data Documentation

bool SubsystemNeutronWriter::initialized [private]

Definition at line 64 of file SubsystemNeutronWriter.h.

std::map<int, int> SubsystemNeutronWriter::theCountPerChamberType [private]

Definition at line 65 of file SubsystemNeutronWriter.h.

Referenced by printStats(), and updateCount().

NeutronWriter* SubsystemNeutronWriter::theHitWriter [private]

Definition at line 59 of file SubsystemNeutronWriter.h.

Referenced by initialize(), SubsystemNeutronWriter(), writeHits(), and ~SubsystemNeutronWriter().

edm::InputTag SubsystemNeutronWriter::theInputTag [private]

Definition at line 60 of file SubsystemNeutronWriter.h.

Referenced by analyze().

double SubsystemNeutronWriter::theNeutronTimeCut [private]

Definition at line 61 of file SubsystemNeutronWriter.h.

Referenced by writeHits().

int SubsystemNeutronWriter::theNEvents [private]

Definition at line 63 of file SubsystemNeutronWriter.h.

Referenced by analyze(), and printStats().

double SubsystemNeutronWriter::theTimeWindow [private]

Definition at line 62 of file SubsystemNeutronWriter.h.

Referenced by writeHits().

---

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

• SimMuon/Neutron/interface/SubsystemNeutronWriter.h
• SimMuon/Neutron/src/SubsystemNeutronWriter.cc

---