#include <HLTPixelAsymmetryFilter.h>

Inheritance diagram for HLTPixelAsymmetryFilter:

Public Member Functions
	HLTPixelAsymmetryFilter (const edm::ParameterSet &)

	~HLTPixelAsymmetryFilter ()

Public Member Functions inherited from HLTFilter
	HLTFilter (const edm::ParameterSet &config)

int	module () const

const std::string *	moduleLabel () const

int	path () const

const std::string *	pathName () const

std::pair< int, int >	pmid () const

bool	saveTags () const

virtual	~HLTFilter ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Private Member Functions
virtual bool	hltFilter (edm::Event &, const edm::EventSetup &, trigger::TriggerFilterObjectWithRefs &filterproduct)

Private Attributes
double	bmincharge_

double	clus_thresh_

edm::InputTag	inputTag_

double	max_asym_

double	min_asym_

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

typedef WorkerT< EDFilter >	WorkerType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Static Public Member Functions inherited from HLTFilter
static void	makeHLTFilterDescription (edm::ParameterSetDescription &desc)

Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDFilter
CurrentProcessingContext const *	currentContext () const

Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Definition at line 50 of file HLTPixelAsymmetryFilter.h.

Constructor & Destructor Documentation

HLTPixelAsymmetryFilter::HLTPixelAsymmetryFilter ( const edm::ParameterSet & config )

explicit

Definition at line 17 of file HLTPixelAsymmetryFilter.cc.

References bmincharge_, clus_thresh_, inputTag_, LogDebug, max_asym_, and min_asym_.

                                                                               : HLTFilter(config),
   inputTag_ (config.getParameter<edm::InputTag>("inputTag")),
   min_asym_ (config.getParameter<double>("MinAsym")),
   max_asym_ (config.getParameter<double>("MaxAsym")),
   clus_thresh_ (config.getParameter<double>("MinCharge")),
   bmincharge_ (config.getParameter<double>("MinBarrel"))
 {
   LogDebug("") << "Using the " << inputTag_ << " input collection";
   LogDebug("") << "Requesting events with a charge repartition asymmetry between " << min_asym_ << " and " << max_asym_;
   LogDebug("") << "Mean cluster charge in the barrel should be higher than" << bmincharge_ << " electrons ";
   LogDebug("") << "Only clusters with a charge larger than " << clus_thresh_ << " electrons will be used ";
 }

HLTPixelAsymmetryFilter::~HLTPixelAsymmetryFilter ( )

Definition at line 30 of file HLTPixelAsymmetryFilter.cc.

31 {

32 }

Member Function Documentation

bool HLTPixelAsymmetryFilter::hltFilter	(	edm::Event &	event,
		const edm::EventSetup &	iSetup,
		trigger::TriggerFilterObjectWithRefs &	filterproduct
	)

privatevirtual

Implements HLTFilter.

Definition at line 39 of file HLTPixelAsymmetryFilter.cc.

References accept(), trigger::TriggerFilterObjectWithRefs::addCollectionTag(), Reference_intrackfit_cff::barrel, begin, bmincharge_, clus_thresh_, cond::rpcobgas::detid, end, Reference_intrackfit_cff::endcap, PixelEndcapName::halfCylinder(), i, inputTag_, LogDebug, PixelEndcapName::mI, min_asym_, PixelEndcapName::mO, PixelEndcapName::pI, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, PixelEndcapName::pO, position, HLTFilter::saveTags(), and DetId::subdetId().

 {
   // All HLT filters must create and fill an HLT filter object,
   // recording any reconstructed physics objects satisfying (or not)
   // this HLT filter, and place it in the Event.
 
   // The filter object
   if (saveTags()) filterproduct.addCollectionTag(inputTag_);
 
   // get hold of products from Event
   edm::Handle<edmNew::DetSetVector<SiPixelCluster> > clusterColl;
   event.getByLabel(inputTag_, clusterColl);
 
   unsigned int clusterSize = clusterColl->dataSize();
   LogDebug("") << "Number of clusters accepted: " << clusterSize;
 
   bool accept = (clusterSize >= 2); // Not necessary to go further in this case
 
   if (!accept) return accept;
 
   double asym_pix_1  = -1;
   double asym_pix_2  = -1;
 
   int n_clus[3]   = {0,0,0};
   double e_clus[3] = {0.,0.,0.};
 
   for (edmNew::DetSetVector<SiPixelCluster>::const_iterator DSViter=clusterColl->begin(); DSViter!=clusterColl->end();DSViter++ ) 
   {
     edmNew::DetSet<SiPixelCluster>::const_iterator begin=DSViter->begin();
     edmNew::DetSet<SiPixelCluster>::const_iterator end  =DSViter->end();
     uint32_t detid = DSViter->id();
 
     bool barrel = DetId(detid).subdetId() == static_cast<int>(PixelSubdetector::PixelBarrel);
     bool endcap = DetId(detid).subdetId() == static_cast<int>(PixelSubdetector::PixelEndcap);
 
     int detpos = -1;
 
     // First we look if we are in the endcap or in the barrel PIXELS
     // Asymmetry is computed with endcap pixels only
 
     if (endcap)
     {
       PixelEndcapName::HalfCylinder position = PixelEndcapName(DetId(detid)).halfCylinder();
 
       if (position == PixelEndcapName::mI || position == PixelEndcapName::mO) // FPIX-
     detpos = 0;
 
       if (position == PixelEndcapName::pI || position == PixelEndcapName::pO) // FPIX+
     detpos = 2;
     }
 
     if (barrel)
       detpos = 1;
 
     if (detpos<0) continue;
 
     for(edmNew::DetSet<SiPixelCluster>::const_iterator iter=begin;iter!=end;++iter) 
     {
       if (iter->charge()>clus_thresh_ )
       {
     ++n_clus[detpos];
     e_clus[detpos]+=iter->charge();
       }    
     }
   } // End of cluster loop
 
   for (int i=0;i<3;++i) 
   {
     if (n_clus[i])
       e_clus[i] /= n_clus[i];
   }
 
   if (e_clus[1] < bmincharge_) return false; // Reject event if the Barrel mean charge is too low
 
 
   if (e_clus[0]+e_clus[2] != 0) // Compute asyms, if applicable
   {
     asym_pix_1 = e_clus[0]/(e_clus[0]+e_clus[2]);
     asym_pix_2 = e_clus[2]/(e_clus[0]+e_clus[2]);
   }
   else  // Otherwise reject event
   {
     return false;
   }
 
   bool pkam_1 = (asym_pix_1 <= max_asym_ && asym_pix_1 >= min_asym_);
   bool pkam_2 = (asym_pix_2 <= max_asym_ && asym_pix_2 >= min_asym_);
 
   if (pkam_1 || pkam_2) return accept; // Final selection
 
   return false;
 }