CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
List of all members | Classes | Public Member Functions | Private Member Functions | Private Attributes
TrackingTruthAccumulator Class Reference

Replacement for TrackingTruthProducer in the new pileup mixing setup. More...

#include <TrackingTruthAccumulator.h>

Inheritance diagram for TrackingTruthAccumulator:
DigiAccumulatorMixMod

Classes

struct  OutputCollections
 

Public Member Functions

 TrackingTruthAccumulator (const edm::ParameterSet &config, edm::one::EDProducerBase &mixMod, edm::ConsumesCollector &iC)
 
- Public Member Functions inherited from DigiAccumulatorMixMod
virtual void beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup)
 
virtual void beginRun (edm::Run const &run, edm::EventSetup const &setup)
 
 DigiAccumulatorMixMod ()
 
virtual void endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup)
 
virtual void endRun (edm::Run const &run, edm::EventSetup const &setup)
 
virtual void finalizeBunchCrossing (edm::Event &event, edm::EventSetup const &setup, int bunchCrossing)
 
virtual PileupMixingContentgetEventPileupInfo ()
 
virtual void initializeBunchCrossing (edm::Event const &event, edm::EventSetup const &setup, int bunchCrossing)
 
virtual void StorePileupInformation (std::vector< int > &numInteractionList, std::vector< int > &bunchCrossingList, std::vector< float > &TrueInteractionList, std::vector< edm::EventID > &eventList, int bunchSpace)
 
virtual ~DigiAccumulatorMixMod ()
 

Private Member Functions

virtual void accumulate (const edm::Event &event, const edm::EventSetup &setup)
 
virtual void accumulate (const PileUpEventPrincipal &event, const edm::EventSetup &setup, edm::StreamID const &)
 
template<class T >
void accumulateEvent (const T &event, const edm::EventSetup &setup, const edm::Handle< edm::HepMCProduct > &hepMCproduct)
 Both forms of accumulate() delegate to this templated method. More...
 
template<class T >
void fillSimHits (std::vector< const PSimHit * > &returnValue, const T &event, const edm::EventSetup &setup)
 Fills the supplied vector with pointers to the SimHits, checking for bad modules if required. More...
 
virtual void finalizeEvent (edm::Event &event, const edm::EventSetup &setup)
 
virtual void initializeEvent (const edm::Event &event, const edm::EventSetup &setup)
 

Private Attributes

const bool addAncestors_
 Whether or not to add the full parentage of any TrackingParticle that is inserted in the collection. More...
 
bool allowDifferentProcessTypeForDifferentDetectors_
 When counting hits, allows hits in different detectors to have a different process type. More...
 
bool chargedOnly_
 Uses the same config as selector_, but can be used to drop out early since selector_ requires the TrackingParticle to be created first. More...
 
std::vector< edm::InputTagcollectionTags_
 
const bool createMergedCollection_
 
const bool createUnmergedCollection_
 If bremsstrahlung merging, whether to also add the unmerged collection to the event or not. More...
 
edm::InputTag genParticleLabel_
 
edm::InputTag hepMCproductLabel_
 Needed to add HepMC::GenVertex to SimVertex. More...
 
const bool ignoreTracksOutsideVolume_
 
const unsigned int maximumPreviousBunchCrossing_
 
const unsigned int maximumSubsequentBunchCrossing_
 
OutputCollections mergedOutput_
 
const std::string messageCategory_
 The message category used to send messages to MessageLogger. More...
 
const bool removeDeadModules_
 As of 11/Feb/2013 this option hasn't been implemented yet. More...
 
TrackingParticleSelector selector_
 
bool selectorFlag_
 
bool signalOnly_
 Uses the same config as selector_, but can be used to drop out early since selector_ requires the TrackingParticle to be created first. More...
 
const edm::InputTag simTrackLabel_
 
const edm::InputTag simVertexLabel_
 
OutputCollections unmergedOutput_
 
const double vertexDistanceCut_
 maximum distance for HepMC::GenVertex to be added to SimVertex More...
 
const double volumeRadius_
 
const double volumeZ_
 

Detailed Description

Replacement for TrackingTruthProducer in the new pileup mixing setup.

The configuration parameters are:

Parameter name Type Description
volumeRadius double The volume radius in cm used if ignoreTracksOutsideVolume is true.
volumeZ double The volume z in cm used if ignoreTracksOutsideVolume is true.
ignoreTracksOutsideVolume bool If true, sim tracks that have a production vertex outside the volume specified by volumeRadius and volumeZ won't be turned into TrackingParticles. Doesn't make much difference to be honest, over a huge range of volume sizes so there must be a cut earlier in the simulation.
maximumPreviousBunchCrossing unsigned int Bunch crossings before this number (inclusive; use positive integer) won't be included. Setting to zero means only in-time.
maximumSubsequentBunchCrossing unsigned int Bunch crossings after this won't create any TrackingParticles.
createUnmergedCollection bool Whether to create the TrackingParticle collection without bremsstrahlung merged.
createMergedBremsstrahlung bool Whether to create the TrackingParticle collection with bremsstrahlung merged. At least one of createUnmergedCollection or createMergedBremsstrahlung should be true otherwise nothing will be produced.
alwaysAddAncestors bool If a sim track passes selection and is turned into a TrackingParticle, all of it's parents will also be created even if they fail the selection. This was the default behaviour for the old TrackingParticleProducer.
removeDeadModules bool Hasn't been implemented yet (as of 22/May/2013).
simTrackCollection edm::InputTag The input SimTrack collection
simVertexCollection edm::InputTag The input SimVerted collection
simHitCollections edm::ParameterSet A ParameterSet of vectors of InputTags that are the input PSimHits
genParticleCollection edm::InputTag The input reco::GenParticle collection. Note that there's a difference between reco::GenParticle and HepMC::GenParticle; the old TrackingTruthProducer used to use HepMC::GenParticle.
allowDifferentSimHitProcesses bool Should be false for FullSim and true for FastSim. There's more documentation in the code if you're really interested.
select edm::ParameterSet A ParameterSet used to configure a TrackingParticleSelector. If the TrackingParticle doesn't pass this selector then it's not added to the output.
Author
Mark Grimes (mark..nosp@m.grim.nosp@m.es@br.nosp@m.isto.nosp@m.l.ac..nosp@m.uk)
Date
11/Oct/2012

Definition at line 68 of file TrackingTruthAccumulator.h.

Constructor & Destructor Documentation

TrackingTruthAccumulator::TrackingTruthAccumulator ( const edm::ParameterSet config,
edm::one::EDProducerBase mixMod,
edm::ConsumesCollector iC 
)
explicit

Definition at line 217 of file TrackingTruthAccumulator.cc.

References chargedOnly_, collectionTags_, edm::ConsumesCollector::consumes(), createMergedCollection_, createUnmergedCollection_, edm::ParameterSet::exists(), genParticleLabel_, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterNames(), edm::ParameterSet::getParameterSet(), hepMCproductLabel_, messageCategory_, selector_, selectorFlag_, signalOnly_, simTrackLabel_, simVertexLabel_, and o2o::tags.

217  :
218  messageCategory_("TrackingTruthAccumulator"),
219  volumeRadius_( config.getParameter<double>("volumeRadius") ),
220  volumeZ_( config.getParameter<double>("volumeZ") ),
221  vertexDistanceCut_( config.getParameter<double>("vertexDistanceCut") ),
222  ignoreTracksOutsideVolume_( config.getParameter<bool>("ignoreTracksOutsideVolume") ),
223  maximumPreviousBunchCrossing_( config.getParameter<unsigned int>("maximumPreviousBunchCrossing") ),
224  maximumSubsequentBunchCrossing_( config.getParameter<unsigned int>("maximumSubsequentBunchCrossing") ),
225  createUnmergedCollection_( config.getParameter<bool>("createUnmergedCollection") ),
226  createMergedCollection_(config.getParameter<bool>("createMergedBremsstrahlung") ),
227  addAncestors_( config.getParameter<bool>("alwaysAddAncestors") ),
228  removeDeadModules_( config.getParameter<bool>("removeDeadModules") ),
229  simTrackLabel_( config.getParameter<edm::InputTag>("simTrackCollection") ),
230  simVertexLabel_( config.getParameter<edm::InputTag>("simVertexCollection") ),
231  collectionTags_( ),
232  genParticleLabel_( config.getParameter<edm::InputTag>("genParticleCollection") ),
233  hepMCproductLabel_( config.getParameter<edm::InputTag>("HepMCProductLabel") ),
234  allowDifferentProcessTypeForDifferentDetectors_( config.getParameter<bool>("allowDifferentSimHitProcesses") )
235 {
236  //
237  // Make sure at least one of the merged and unmerged collections have been set
238  // to be created.
239  //
241  edm::LogError(messageCategory_) << "Both \"createUnmergedCollection\" and \"createMergedBremsstrahlung\" have been"
242  << "set to false, which means no collections will be created";
243 
244  // Initialize selection for building TrackingParticles
245  //
246  if( config.exists( "select" ) )
247  {
248  edm::ParameterSet param=config.getParameter<edm::ParameterSet>("select");
249  selector_=TrackingParticleSelector( param.getParameter<double>( "ptMinTP" ),
250  param.getParameter<double>( "minRapidityTP" ),
251  param.getParameter<double>( "maxRapidityTP" ),
252  param.getParameter<double>( "tipTP" ),
253  param.getParameter<double>( "lipTP" ),
254  param.getParameter<int>( "minHitTP" ),
255  param.getParameter<bool>( "signalOnlyTP" ),
256  param.getParameter<bool>( "chargedOnlyTP" ),
257  param.getParameter<bool>( "stableOnlyTP" ),
258  param.getParameter<std::vector<int> >("pdgIdTP") );
259  selectorFlag_=true;
260 
261  // Also set these two variables, which are used to drop out early if the SimTrack doesn't conform.
262  // The selector_ requires a full TrackingParticle object, but these two variables can veto things early.
263  chargedOnly_=param.getParameter<bool>( "chargedOnlyTP" );
264  signalOnly_=param.getParameter<bool>( "signalOnlyTP" );
265  }
266  else
267  {
268  selectorFlag_=false;
269  chargedOnly_=false;
270  signalOnly_=false;
271  }
272 
273  //
274  // Need to state what collections are going to be added to the event. This
275  // depends on which of the merged and unmerged collections have been configured
276  // to be created.
277  //
279  {
280  mixMod.produces<TrackingVertexCollection>();
281  mixMod.produces<TrackingParticleCollection>();
282  }
283 
285  {
286  mixMod.produces<TrackingParticleCollection>("MergedTrackTruth");
287  mixMod.produces<TrackingVertexCollection>("MergedTrackTruth");
288  }
289 
290  iC.consumes<std::vector<SimTrack> >(simTrackLabel_);
291  iC.consumes<std::vector<SimVertex> >(simVertexLabel_);
292  iC.consumes<std::vector<reco::GenParticle> >(genParticleLabel_);
293  iC.consumes<std::vector<int> >(genParticleLabel_);
294  iC.consumes<std::vector<int> >(hepMCproductLabel_);
295 
296  // Fill the collection tags
297  const edm::ParameterSet& simHitCollectionConfig=config.getParameterSet("simHitCollections");
298  std::vector<std::string> parameterNames=simHitCollectionConfig.getParameterNames();
299 
300  for( const auto& parameterName : parameterNames )
301  {
302  std::vector<edm::InputTag> tags=simHitCollectionConfig.getParameter<std::vector<edm::InputTag> >(parameterName);
303  collectionTags_.insert(collectionTags_.end(), tags.begin(), tags.end());
304  }
305 
306  for( const auto& collectionTag : collectionTags_ ) {
307  iC.consumes<std::vector<PSimHit> >(collectionTag);
308  }
309 
310 }
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
T getParameter(std::string const &) const
std::vector< TrackingParticle > TrackingParticleCollection
const bool removeDeadModules_
As of 11/Feb/2013 this option hasn&#39;t been implemented yet.
const bool addAncestors_
Whether or not to add the full parentage of any TrackingParticle that is inserted in the collection...
bool exists(std::string const &parameterName) const
checks if a parameter exists
const double vertexDistanceCut_
maximum distance for HepMC::GenVertex to be added to SimVertex
bool signalOnly_
Uses the same config as selector_, but can be used to drop out early since selector_ requires the Tra...
const unsigned int maximumSubsequentBunchCrossing_
SingleObjectSelector< TrackingParticleCollection,::TrackingParticleSelector > TrackingParticleSelector
const std::string messageCategory_
The message category used to send messages to MessageLogger.
std::vector< std::string > getParameterNames() const
bool allowDifferentProcessTypeForDifferentDetectors_
When counting hits, allows hits in different detectors to have a different process type...
tuple tags
Definition: o2o.py:248
ParameterSet const & getParameterSet(std::string const &) const
std::vector< TrackingVertex > TrackingVertexCollection
const unsigned int maximumPreviousBunchCrossing_
bool chargedOnly_
Uses the same config as selector_, but can be used to drop out early since selector_ requires the Tra...
const bool createUnmergedCollection_
If bremsstrahlung merging, whether to also add the unmerged collection to the event or not...
TrackingParticleSelector selector_
edm::InputTag hepMCproductLabel_
Needed to add HepMC::GenVertex to SimVertex.
std::vector< edm::InputTag > collectionTags_

Member Function Documentation

void TrackingTruthAccumulator::accumulate ( const edm::Event event,
const edm::EventSetup setup 
)
privatevirtual

create handle to edm::HepMCProduct here because event.getByLabel with edm::HepMCProduct only works for edm::Event but not for PileUpEventPrincipal; PileUpEventPrincipal::getByLabel tries to call T::value_type and T::iterator (where T is the type of the object one wants to get a handle to) which is only implemented for container-like objects like std::vector but not for edm::HepMCProduct!

Implements DigiAccumulatorMixMod.

Definition at line 336 of file TrackingTruthAccumulator.cc.

References accumulateEvent(), and hepMCproductLabel_.

337 {
338  // Call the templated version that does the same for both signal and pileup events
339 
341  event.getByLabel(hepMCproductLabel_, hepmc);
342 
343  accumulateEvent( event, setup, hepmc );
344 }
void accumulateEvent(const T &event, const edm::EventSetup &setup, const edm::Handle< edm::HepMCProduct > &hepMCproduct)
Both forms of accumulate() delegate to this templated method.
edm::InputTag hepMCproductLabel_
Needed to add HepMC::GenVertex to SimVertex.
void TrackingTruthAccumulator::accumulate ( const PileUpEventPrincipal event,
const edm::EventSetup setup,
edm::StreamID const &   
)
privatevirtual

Implements DigiAccumulatorMixMod.

Definition at line 346 of file TrackingTruthAccumulator.cc.

References accumulateEvent(), PileUpEventPrincipal::bunchCrossing(), maximumPreviousBunchCrossing_, maximumSubsequentBunchCrossing_, and messageCategory_.

347 {
348  // If this bunch crossing is outside the user configured limit, don't do anything.
349  if( event.bunchCrossing()>=-static_cast<int>(maximumPreviousBunchCrossing_) && event.bunchCrossing()<=static_cast<int>(maximumSubsequentBunchCrossing_) )
350  {
351  //edm::LogInfo(messageCategory_) << "Analysing pileup event for bunch crossing " << event.bunchCrossing();
352 
353  //simply create empty handle as we do not have a HepMCProduct in PU anyway
355  accumulateEvent( event, setup, hepmc );
356  }
357  else edm::LogInfo(messageCategory_) << "Skipping pileup event for bunch crossing " << event.bunchCrossing();
358 }
const unsigned int maximumSubsequentBunchCrossing_
void accumulateEvent(const T &event, const edm::EventSetup &setup, const edm::Handle< edm::HepMCProduct > &hepMCproduct)
Both forms of accumulate() delegate to this templated method.
const std::string messageCategory_
The message category used to send messages to MessageLogger.
const unsigned int maximumPreviousBunchCrossing_
template<class T >
void TrackingTruthAccumulator::accumulateEvent ( const T event,
const edm::EventSetup setup,
const edm::Handle< edm::HepMCProduct > &  hepMCproduct 
)
private

Both forms of accumulate() delegate to this templated method.

Definition at line 383 of file TrackingTruthAccumulator.cc.

References addAncestors_, allowDifferentProcessTypeForDifferentDetectors_, EncodedEventId::bunchCrossing(), CoreSimTrack::charge(), chargedOnly_, createMergedCollection_, createUnmergedCollection_, EncodedEventId::event(), CoreSimTrack::eventId(), cppFunctionSkipper::exception, fillSimHits(), genParticleLabel_, edm::EventSetup::get(), ignoreTracksOutsideVolume_, cmsHarvester::index, mergedOutput_, NULL, CoreSimVertex::position(), edm::ESHandle< class >::product(), selector_, selectorFlag_, signalOnly_, simTrackLabel_, simVertexLabel_, unmergedOutput_, vertexDistanceCut_, volumeRadius_, and volumeZ_.

Referenced by accumulate().

384 {
385  //
386  // Get the collections
387  //
389  edm::Handle<std::vector<SimVertex> > hSimVertices;
391  edm::Handle< std::vector<int> > hGenParticleIndices;
392 
393  event.getByLabel( simTrackLabel_, hSimTracks );
394  event.getByLabel( simVertexLabel_, hSimVertices );
395 
396  try
397  {
398  event.getByLabel( genParticleLabel_, hGenParticles );
399  event.getByLabel( genParticleLabel_, hGenParticleIndices );
400  }
401  catch( cms::Exception& exception )
402  {
403  //
404  // The Monte Carlo is not always available, e.g. for pileup events. The information
405  // is only used if it's available, but for some reason the PileUpEventPrincipal
406  // wrapper throws an exception here rather than waiting to see if the handle is
407  // used (as is the case for edm::Event). So I just want to catch this exception
408  // and use the normal handle checking later on.
409  //
410  }
411 
412  //Retrieve tracker topology from geometry
413  edm::ESHandle<TrackerTopology> tTopoHandle;
414  setup.get<IdealGeometryRecord>().get(tTopoHandle);
415  const TrackerTopology* const tTopo = tTopoHandle.product();
416 
417 
418  // Run through the collections and work out the decay chain of each track/vertex. The
419  // information in SimTrack and SimVertex only allows traversing upwards, but this will
420  // allow traversal in both directions. This is required for things like grouping electrons
421  // that bremsstrahlung as one TrackingParticle if "mergedBremsstrahlung" is set in the
422  // config file.
423  DecayChain decayChain( *hSimTracks, *hSimVertices );
424 
425  // I only want to create these collections if they're actually required
426  std::auto_ptr< ::OutputCollectionWrapper> pUnmergedCollectionWrapper;
427  std::auto_ptr< ::OutputCollectionWrapper> pMergedCollectionWrapper;
428  if( createUnmergedCollection_ ) pUnmergedCollectionWrapper.reset( new ::OutputCollectionWrapper( decayChain, unmergedOutput_ ) );
429  if( createMergedCollection_ ) pMergedCollectionWrapper.reset( new ::OutputCollectionWrapper( decayChain, mergedOutput_ ) );
430 
431  std::vector<const PSimHit*> simHitPointers;
432  fillSimHits( simHitPointers, event, setup );
433  TrackingParticleFactory objectFactory( decayChain, hGenParticles, hepMCproduct, hGenParticleIndices, simHitPointers, volumeRadius_, volumeZ_, vertexDistanceCut_, allowDifferentProcessTypeForDifferentDetectors_ );
434 
435  // While I'm testing, perform some checks.
436  // TODO - drop this call once I'm happy it works in all situations.
437  //decayChain.integrityCheck();
438 
439  TrackingParticleSelector* pSelector=NULL;
440  if( selectorFlag_ ) pSelector=&selector_;
441 
442  // Run over all of the SimTracks, but because I'm interested in the decay hierarchy
443  // do it through the DecayChainTrack objects. These are looped over in sequence here
444  // but they have the hierarchy information for the functions called to traverse the
445  // decay chain.
446 
447  for( size_t index=0; index<decayChain.decayTracksSize; ++index )
448  {
449  ::DecayChainTrack* pDecayTrack=&decayChain.decayTracks[index];
450  const SimTrack& simTrack=hSimTracks->at(pDecayTrack->simTrackIndex);
451 
452 
453  // Perform some quick checks to see if we can drop out early. Note that these are
454  // a subset of the cuts in the selector_ so the created TrackingParticle could still
455  // fail. The selector_ requires the full TrackingParticle to be made however, which
456  // can be computationally expensive.
457  if( chargedOnly_ && simTrack.charge()==0 ) continue;
458  if( signalOnly_ && (simTrack.eventId().bunchCrossing()!=0 || simTrack.eventId().event()!=0) ) continue;
459 
460  // Also perform a check to see if the production vertex is inside the tracker volume (if required).
462  {
463  const SimVertex& simVertex=hSimVertices->at( pDecayTrack->pParentVertex->simVertexIndex );
464  if( !objectFactory.vectorIsInsideVolume( simVertex.position() ) ) continue;
465  }
466 
467 
468  // This function creates the TrackinParticle and adds it to the collection if it
469  // passes the selection criteria specified in the configuration. If the config
470  // specifies adding ancestors, the function is called recursively to do that.
471  ::addTrack( pDecayTrack, pSelector, pUnmergedCollectionWrapper.get(), pMergedCollectionWrapper.get(), objectFactory, addAncestors_, tTopo );
472  }
473 }
EncodedEventId eventId() const
Definition: CoreSimTrack.h:31
int event() const
get the contents of the subdetector field (should be protected?)
const bool addAncestors_
Whether or not to add the full parentage of any TrackingParticle that is inserted in the collection...
const double vertexDistanceCut_
maximum distance for HepMC::GenVertex to be added to SimVertex
#define NULL
Definition: scimark2.h:8
float charge() const
charge
Definition: CoreSimTrack.cc:18
bool signalOnly_
Uses the same config as selector_, but can be used to drop out early since selector_ requires the Tra...
void fillSimHits(std::vector< const PSimHit * > &returnValue, const T &event, const edm::EventSetup &setup)
Fills the supplied vector with pointers to the SimHits, checking for bad modules if required...
int bunchCrossing() const
get the detector field from this detid
const math::XYZTLorentzVectorD & position() const
Definition: CoreSimVertex.h:26
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
bool allowDifferentProcessTypeForDifferentDetectors_
When counting hits, allows hits in different detectors to have a different process type...
const T & get() const
Definition: EventSetup.h:55
T const * product() const
Definition: ESHandle.h:86
bool chargedOnly_
Uses the same config as selector_, but can be used to drop out early since selector_ requires the Tra...
const bool createUnmergedCollection_
If bremsstrahlung merging, whether to also add the unmerged collection to the event or not...
TrackingParticleSelector selector_
template<class T >
void TrackingTruthAccumulator::fillSimHits ( std::vector< const PSimHit * > &  returnValue,
const T event,
const edm::EventSetup setup 
)
private

Fills the supplied vector with pointers to the SimHits, checking for bad modules if required.

Definition at line 475 of file TrackingTruthAccumulator.cc.

References collectionTags_.

Referenced by accumulateEvent().

476 {
477  // loop over the collections
478  for( const auto& collectionTag : collectionTags_ )
479  {
481  event.getByLabel( collectionTag, hSimHits );
482 
483  // TODO - implement removing the dead modules
484  for( const auto& simHit : *hSimHits )
485  {
486  returnValue.push_back( &simHit );
487  }
488 
489  } // end of loop over InputTags
490 }
std::vector< edm::InputTag > collectionTags_
void TrackingTruthAccumulator::finalizeEvent ( edm::Event event,
const edm::EventSetup setup 
)
privatevirtual

Implements DigiAccumulatorMixMod.

Definition at line 360 of file TrackingTruthAccumulator.cc.

References createMergedCollection_, createUnmergedCollection_, mergedOutput_, TrackingTruthAccumulator::OutputCollections::pTrackingParticles, TrackingTruthAccumulator::OutputCollections::pTrackingVertices, and unmergedOutput_.

361 {
362 
364  {
365  edm::LogInfo("TrackingTruthAccumulator") << "Adding " << unmergedOutput_.pTrackingParticles->size() << " TrackingParticles and " << unmergedOutput_.pTrackingVertices->size()
366  << " TrackingVertexs to the event.";
367 
368  event.put( unmergedOutput_.pTrackingParticles );
369  event.put( unmergedOutput_.pTrackingVertices );
370  }
371 
373  {
374  edm::LogInfo("TrackingTruthAccumulator") << "Adding " << mergedOutput_.pTrackingParticles->size() << " merged TrackingParticles and " << mergedOutput_.pTrackingVertices->size()
375  << " merged TrackingVertexs to the event.";
376 
377  event.put( mergedOutput_.pTrackingParticles, "MergedTrackTruth" );
378  event.put( mergedOutput_.pTrackingVertices, "MergedTrackTruth" );
379  }
380 
381 }
std::auto_ptr< TrackingParticleCollection > pTrackingParticles
std::auto_ptr< TrackingVertexCollection > pTrackingVertices
const bool createUnmergedCollection_
If bremsstrahlung merging, whether to also add the unmerged collection to the event or not...
void TrackingTruthAccumulator::initializeEvent ( const edm::Event event,
const edm::EventSetup setup 
)
privatevirtual

Implements DigiAccumulatorMixMod.

Definition at line 312 of file TrackingTruthAccumulator.cc.

References createMergedCollection_, createUnmergedCollection_, event(), mergedOutput_, TrackingTruthAccumulator::OutputCollections::pTrackingParticles, TrackingTruthAccumulator::OutputCollections::pTrackingVertices, TrackingTruthAccumulator::OutputCollections::refTrackingParticles, TrackingTruthAccumulator::OutputCollections::refTrackingVertexes, and unmergedOutput_.

313 {
315  {
318  unmergedOutput_.refTrackingParticles=const_cast<edm::Event&>( event ).getRefBeforePut<TrackingParticleCollection>();
319  unmergedOutput_.refTrackingVertexes=const_cast<edm::Event&>( event ).getRefBeforePut<TrackingVertexCollection>();
320  }
321 
323  {
326  mergedOutput_.refTrackingParticles=const_cast<edm::Event&>( event ).getRefBeforePut<TrackingParticleCollection>("MergedTrackTruth");
327  mergedOutput_.refTrackingVertexes=const_cast<edm::Event&>( event ).getRefBeforePut<TrackingVertexCollection>("MergedTrackTruth");
328  }
329 }
std::vector< TrackingParticle > TrackingParticleCollection
How EventSelector::AcceptEvent() decides whether to accept an event for output otherwise it is excluding the probing of A single or multiple positive and the trigger will pass if any such matching triggers are PASS or EXCEPTION[A criterion thatmatches no triggers at all is detected and causes a throw.] A single negative with an expectation of appropriate bit checking in the decision and the trigger will pass if any such matching triggers are FAIL or EXCEPTION A wildcarded negative criterion that matches more than one trigger in the trigger but the state exists so we define the behavior If all triggers are the negative crieriion will lead to accepting the event(this again matches the behavior of"!*"before the partial wildcard feature was incorporated).The per-event"cost"of each negative criterion with multiple relevant triggers is about the same as!*was in the past
std::auto_ptr< TrackingParticleCollection > pTrackingParticles
std::auto_ptr< TrackingVertexCollection > pTrackingVertices
std::vector< TrackingVertex > TrackingVertexCollection
const bool createUnmergedCollection_
If bremsstrahlung merging, whether to also add the unmerged collection to the event or not...

Member Data Documentation

const bool TrackingTruthAccumulator::addAncestors_
private

Whether or not to add the full parentage of any TrackingParticle that is inserted in the collection.

Definition at line 102 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent().

bool TrackingTruthAccumulator::allowDifferentProcessTypeForDifferentDetectors_
private

When counting hits, allows hits in different detectors to have a different process type.

Fast sim PSimHits seem to have a peculiarity where the process type (as reported by PSimHit::processType()) is different for the tracker than the muons. When counting how many hits there are, the code usually only counts the number of hits that have the same process type as the first hit. Setting this to true will also count hits that have the same process type as the first hit in the second detector.

Definition at line 127 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent().

bool TrackingTruthAccumulator::chargedOnly_
private

Uses the same config as selector_, but can be used to drop out early since selector_ requires the TrackingParticle to be created first.

Definition at line 116 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), and TrackingTruthAccumulator().

std::vector<edm::InputTag> TrackingTruthAccumulator::collectionTags_
private

Definition at line 108 of file TrackingTruthAccumulator.h.

Referenced by fillSimHits(), and TrackingTruthAccumulator().

const bool TrackingTruthAccumulator::createMergedCollection_
private
const bool TrackingTruthAccumulator::createUnmergedCollection_
private

If bremsstrahlung merging, whether to also add the unmerged collection to the event or not.

Definition at line 99 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), finalizeEvent(), initializeEvent(), and TrackingTruthAccumulator().

edm::InputTag TrackingTruthAccumulator::genParticleLabel_
private

Definition at line 109 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), and TrackingTruthAccumulator().

edm::InputTag TrackingTruthAccumulator::hepMCproductLabel_
private

Needed to add HepMC::GenVertex to SimVertex.

Definition at line 111 of file TrackingTruthAccumulator.h.

Referenced by accumulate(), and TrackingTruthAccumulator().

const bool TrackingTruthAccumulator::ignoreTracksOutsideVolume_
private

Definition at line 90 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent().

const unsigned int TrackingTruthAccumulator::maximumPreviousBunchCrossing_
private

The maximum bunch crossing BEFORE the signal crossing to create TrackinParticles for. Use positive values. If set to zero no previous bunches are added and only in-time, signal and after bunches (defined by maximumSubsequentBunchCrossing_) are used.

Definition at line 94 of file TrackingTruthAccumulator.h.

Referenced by accumulate().

const unsigned int TrackingTruthAccumulator::maximumSubsequentBunchCrossing_
private

The maximum bunch crossing AFTER the signal crossing to create TrackinParticles for. E.g. if set to zero only uses the signal and in time pileup (and previous bunches defined by the maximumPreviousBunchCrossing_ parameter).

Definition at line 97 of file TrackingTruthAccumulator.h.

Referenced by accumulate().

OutputCollections TrackingTruthAccumulator::mergedOutput_
private

Definition at line 140 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), finalizeEvent(), and initializeEvent().

const std::string TrackingTruthAccumulator::messageCategory_
private

The message category used to send messages to MessageLogger.

Definition at line 84 of file TrackingTruthAccumulator.h.

Referenced by accumulate(), and TrackingTruthAccumulator().

const bool TrackingTruthAccumulator::removeDeadModules_
private

As of 11/Feb/2013 this option hasn't been implemented yet.

Definition at line 105 of file TrackingTruthAccumulator.h.

TrackingParticleSelector TrackingTruthAccumulator::selector_
private

Definition at line 114 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), and TrackingTruthAccumulator().

bool TrackingTruthAccumulator::selectorFlag_
private

Definition at line 113 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), and TrackingTruthAccumulator().

bool TrackingTruthAccumulator::signalOnly_
private

Uses the same config as selector_, but can be used to drop out early since selector_ requires the TrackingParticle to be created first.

Definition at line 118 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), and TrackingTruthAccumulator().

const edm::InputTag TrackingTruthAccumulator::simTrackLabel_
private

Definition at line 106 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), and TrackingTruthAccumulator().

const edm::InputTag TrackingTruthAccumulator::simVertexLabel_
private

Definition at line 107 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), and TrackingTruthAccumulator().

OutputCollections TrackingTruthAccumulator::unmergedOutput_
private

Definition at line 139 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent(), finalizeEvent(), and initializeEvent().

const double TrackingTruthAccumulator::vertexDistanceCut_
private

maximum distance for HepMC::GenVertex to be added to SimVertex

Definition at line 89 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent().

const double TrackingTruthAccumulator::volumeRadius_
private

Definition at line 86 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent().

const double TrackingTruthAccumulator::volumeZ_
private

Definition at line 87 of file TrackingTruthAccumulator.h.

Referenced by accumulateEvent().