#include <AlignmentPrescaler.h>

Inheritance diagram for AlignmentPrescaler:

Public Member Functions
	AlignmentPrescaler (const edm::ParameterSet &iConfig)

void	beginJob () override

void	endJob () override

void	produce (edm::Event &iEvent, const edm::EventSetup &iSetup) override

	~AlignmentPrescaler () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDProducer () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
int	layerFromId (const DetId &id, const TrackerTopology *tTopo) const

Private Attributes
unsigned int	detid_

TFile *	fpresc_

float	hitPrescFactor_

TRandom3 *	myrand_

float	overlapPrescFactor_

std::string	prescfilename_

std::string	presctreename_

edm::InputTag	src_

edm::InputTag	srcQualityMap_

int	totnhitspxl_

TTree *	tpresc_

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

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 23 of file AlignmentPrescaler.h.

Constructor & Destructor Documentation

AlignmentPrescaler::AlignmentPrescaler ( const edm::ParameterSet & iConfig )

Definition at line 28 of file AlignmentPrescaler.cc.

     : src_(iConfig.getParameter<edm::InputTag>("src")),
       srcQualityMap_(iConfig.getParameter<edm::InputTag>("assomap")),
       prescfilename_(iConfig.getParameter<std::string>("PrescFileName")),
       presctreename_(iConfig.getParameter<std::string>("PrescTreeName")) {
   // issue the produce<>
   produces<AliClusterValueMap>();
   produces<AliTrackTakenClusterValueMap>();
 }

AlignmentPrescaler::~AlignmentPrescaler ( )

override

Definition at line 38 of file AlignmentPrescaler.cc.

                                         {
   //
 }

Member Function Documentation

void AlignmentPrescaler::beginJob ( void )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 42 of file AlignmentPrescaler.cc.

References gather_cfg::cout, detid_, fpresc_, hitPrescFactor_, myrand_, overlapPrescFactor_, prescfilename_, presctreename_, and tpresc_.

                                   {
   //
   std::cout << "in AlignmentPrescaler::beginJob" << std::flush;
   fpresc_ = new TFile(prescfilename_.c_str(), "READ");
   tpresc_ = (TTree*)fpresc_->Get(presctreename_.c_str());
   tpresc_->BuildIndex("DetId");
   tpresc_->SetBranchStatus("*", false);
   tpresc_->SetBranchStatus("DetId", true);
   tpresc_->SetBranchStatus("PrescaleFactor", true);
   tpresc_->SetBranchStatus("PrescaleFactorOverlap", true);
   cout << " Branches activated " << std::flush;
   detid_ = 0;
   hitPrescFactor_ = 99.0;
   overlapPrescFactor_ = 88.0;
 
   tpresc_->SetBranchAddress("DetId", &detid_);
   tpresc_->SetBranchAddress("PrescaleFactor", &hitPrescFactor_);
   tpresc_->SetBranchAddress("PrescaleFactorOverlap", &overlapPrescFactor_);
   cout << " addressed " << std::flush;
   myrand_ = new TRandom3();
   //   myrand_->SetSeed();
   cout << " ok " << std::endl;
 }

void AlignmentPrescaler::endJob ( void )

overridevirtual

Reimplemented from edm::EDProducer.

Definition at line 66 of file AlignmentPrescaler.cc.

References fpresc_, myrand_, and tpresc_.

Referenced by o2olib.O2ORunMgr::executeJob().

                                 {
   delete tpresc_;
   fpresc_->Close();
   delete fpresc_;
   delete myrand_;
 }

int AlignmentPrescaler::layerFromId	(	const DetId &	id,
		const TrackerTopology *	tTopo
	)		const

private

Definition at line 230 of file AlignmentPrescaler.cc.

References DEFINE_FWK_MODULE, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, TrackerTopology::pxbLayer(), TrackerTopology::pxfDisk(), TrackerTopology::pxfSide(), StripSubdetector::TEC, TrackerTopology::tecWheel(), StripSubdetector::TIB, TrackerTopology::tibLayer(), StripSubdetector::TID, TrackerTopology::tidSide(), TrackerTopology::tidWheel(), StripSubdetector::TOB, and TrackerTopology::tobLayer().

                                                                                        {
   if (uint32_t(id.subdetId()) == PixelSubdetector::PixelBarrel) {
     return tTopo->pxbLayer(id);
   } else if (uint32_t(id.subdetId()) == PixelSubdetector::PixelEndcap) {
     return tTopo->pxfDisk(id) + (3 * (tTopo->pxfSide(id) - 1));
   } else if (id.subdetId() == StripSubdetector::TIB) {
     return tTopo->tibLayer(id);
   } else if (id.subdetId() == StripSubdetector::TOB) {
     return tTopo->tobLayer(id);
   } else if (id.subdetId() == StripSubdetector::TEC) {
     return tTopo->tecWheel(id) + (9 * (tTopo->pxfSide(id) - 1));
   } else if (id.subdetId() == StripSubdetector::TID) {
     return tTopo->tidWheel(id) + (3 * (tTopo->tidSide(id) - 1));
   }
   return -1;
 
 }  //end layerfromId

void AlignmentPrescaler::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

override

Definition at line 73 of file AlignmentPrescaler.cc.

References SiStripRecHit1D::cluster(), SiStripRecHit2D::cluster(), SiPixelRecHit::cluster(), gather_cfg::cout, dqmMemoryStats::float, edm::Event::getByLabel(), hitPrescFactor_, eostools::move(), myrand_, overlapPrescFactor_, pixelClusterTagInfos_cfi::pixelhit, edm::Event::put(), findQualityFiles::rr, src_, srcQualityMap_, tpresc_, and gather_cfg::Tracks.

Referenced by JSONExport.JsonExport::export(), HTMLExport.HTMLExport::export(), and HTMLExport.HTMLExportStatic::export().

                                                                               {
   // std::cout<<"\n\n#################\n### Starting the AlignmentPrescaler::produce ; Event: "<<iEvent.id().run() <<", "<<iEvent.id().event()<<std::endl;
   edm::Handle<reco::TrackCollection> Tracks;
   iEvent.getByLabel(src_, Tracks);
 
   //take  HitAssomap
   edm::Handle<AliClusterValueMap> hMap;
   iEvent.getByLabel(srcQualityMap_, hMap);
   AliClusterValueMap InValMap = *hMap;
 
   //prepare the output of the ValueMap flagging tracks
   std::vector<int> trackflags(Tracks->size(), 0);
 
   //int npxlhits=0;
 
   //loop on tracks
   for (std::vector<reco::Track>::const_iterator ittrk = Tracks->begin(), edtrk = Tracks->end(); ittrk != edtrk;
        ++ittrk) {
     //loop on tracking rechits
     // std::cout << "Loop on hits of track #" << (ittrk - Tracks->begin()) << std::endl;
     int nhit = 0;
     int ntakenhits = 0;
     bool firstTakenHit = false;
 
     for (auto const& hit : ittrk->recHits()) {
       if (!hit->isValid()) {
         nhit++;
         continue;
       }
       uint32_t tmpdetid = hit->geographicalId().rawId();
       tpresc_->GetEntryWithIndex(tmpdetid);
 
       //-------------
       //decide whether to take this hit or not
       bool takeit = false;
       int subdetId = hit->geographicalId().subdetId();
 
       //check first if the cluster is also in the overlap asso map
       bool isOverlapHit = false;
       //  bool first=true;
       //ugly...
       const SiPixelRecHit* pixelhit = dynamic_cast<const SiPixelRecHit*>(hit);
       const SiStripRecHit1D* stripHit1D = dynamic_cast<const SiStripRecHit1D*>(hit);
       const SiStripRecHit2D* stripHit2D = dynamic_cast<const SiStripRecHit2D*>(hit);
 
       AlignmentClusterFlag tmpflag(hit->geographicalId());
       int stripType = 0;
       if (subdetId > 2) {  // SST case
         const std::type_info& type = typeid(*hit);
         if (type == typeid(SiStripRecHit1D))
           stripType = 1;
         else if (type == typeid(SiStripRecHit2D))
           stripType = 2;
         else
           stripType = 3;
 
         if (stripType == 1) {
           //      const SiStripRecHit1D* stripHit1D = dynamic_cast<const SiStripRecHit1D*>(hit);
 
           if (stripHit1D != nullptr) {
             SiStripRecHit1D::ClusterRef stripclust(stripHit1D->cluster());
             tmpflag = InValMap[stripclust];
             tmpflag.SetDetId(hit->geographicalId());
             if (tmpflag.isOverlap())
               isOverlapHit = true;
             // std::cout<<"~*~*~* Prescale (1D) for module "<<tmpflag.detId().rawId()<<"("<<InValMap[stripclust].detId().rawId() <<") is "<<hitPrescFactor_<<std::flush;
             //  if(tmpflag.isOverlap())cout<<" (it is Overlap)"<<endl;
             // else cout<<endl;
 
           }  //end if striphit1D!=0
         } else if (stripType == 2) {
           //const SiStripRecHit2D* stripHit2D = dynamic_cast<const SiStripRecHit2D*>(hit);
           if (stripHit2D != nullptr) {
             SiStripRecHit2D::ClusterRef stripclust(stripHit2D->cluster());
             tmpflag = InValMap[stripclust];
             tmpflag.SetDetId(hit->geographicalId());
             if (tmpflag.isOverlap())
               isOverlapHit = true;
             // std::cout<<"~*~*~* Prescale (2D) for module "<<tmpflag.detId().rawId()<<"("<<InValMap[stripclust].detId().rawId() <<") is "<<hitPrescFactor_<<std::flush;
             //  if(tmpflag.isOverlap())cout<<" (it is Overlap)"<<endl;
             // else cout<<endl;
 
           }  //end if striphit2D!=0
         }
       }  //end if is a strip hit
       else {
         //  const SiPixelRecHit*   pixelhit= dynamic_cast<const SiPixelRecHit*>(hit);
         if (pixelhit != nullptr) {
           //npxlhits++;
           SiPixelClusterRefNew pixclust(pixelhit->cluster());
           tmpflag = InValMap[pixclust];
           tmpflag.SetDetId(hit->geographicalId());
           if (tmpflag.isOverlap())
             isOverlapHit = true;
         }
       }  //end else is a pixel hit
       //      tmpflag.SetDetId(hit->geographicalId());
 
       if (isOverlapHit) {
         //cout<<"  DetId="<<tmpdetid<<" is Overlap! "<<flush;
         takeit = (float(myrand_->Rndm()) <= overlapPrescFactor_);
       }
       if (!takeit) {
         float rr = float(myrand_->Rndm());
         takeit = (rr <= hitPrescFactor_);
       }
       if (takeit) {  //HIT TAKEN !
         //cout<<"  DetId="<<tmpdetid<<" taken!"<<flush;
         tmpflag.SetTakenFlag();
 
         if (subdetId > 2) {
           if (stripType == 1) {
             SiStripRecHit1D::ClusterRef stripclust(stripHit1D->cluster());
             InValMap[stripclust] = tmpflag;  //.SetTakenFlag();
           } else if (stripType == 2) {
             SiStripRecHit1D::ClusterRef stripclust(stripHit2D->cluster());
             InValMap[stripclust] = tmpflag;  //.SetTakenFlag();
           } else
             std::cout << "Unknown type of strip hit" << std::endl;
         } else {
           SiPixelClusterRefNew pixclust(pixelhit->cluster());
           InValMap[pixclust] = tmpflag;  //.SetTakenFlag();
         }
 
         if (!firstTakenHit) {
           firstTakenHit = true;
           //std::cout<<"Index of the track iterator is "<< ittrk-Tracks->begin() <<endl;
         }
         ntakenhits++;
       }  //end if take this hit
          //cout<<endl;
 
       nhit++;
       //cout<<endl;
     }  //end loop on RecHits
     trackflags[ittrk - Tracks->begin()] = ntakenhits;
 
   }  //end loop on tracks
 
   // totnhitspxl_+=ntakenhits;
   //cout<<"AlignmentPrescaler::produce says that in this event "<<ntakenhits<<" pixel clusters were taken (out of "<<npxlhits<<" total pixel hits."<<endl;
 
   //save the asso map, tracks...
   // prepare output
   auto OutVM = std::make_unique<AliClusterValueMap>();
   *OutVM = InValMap;
 
   iEvent.put(std::move(OutVM));
 
   auto trkVM = std::make_unique<AliTrackTakenClusterValueMap>();
   AliTrackTakenClusterValueMap::Filler trkmapfiller(*trkVM);
   trkmapfiller.insert(Tracks, trackflags.begin(), trackflags.end());
   trkmapfiller.fill();
   iEvent.put(std::move(trkVM));
 
 }  //end produce

Member Data Documentation

unsigned int AlignmentPrescaler::detid_

private

Definition at line 44 of file AlignmentPrescaler.h.

Referenced by beginJob().

TFile* AlignmentPrescaler::fpresc_

private

Definition at line 38 of file AlignmentPrescaler.h.

Referenced by beginJob(), and endJob().

float AlignmentPrescaler::hitPrescFactor_

private

Definition at line 45 of file AlignmentPrescaler.h.

Referenced by beginJob(), and produce().

TRandom3* AlignmentPrescaler::myrand_

private

Definition at line 40 of file AlignmentPrescaler.h.

Referenced by beginJob(), endJob(), and produce().

float AlignmentPrescaler::overlapPrescFactor_

private

Definition at line 45 of file AlignmentPrescaler.h.

Referenced by beginJob(), and produce().

std::string AlignmentPrescaler::prescfilename_

private

Definition at line 35 of file AlignmentPrescaler.h.

Referenced by beginJob().

std::string AlignmentPrescaler::presctreename_

private

Definition at line 36 of file AlignmentPrescaler.h.

Referenced by beginJob().

edm::InputTag AlignmentPrescaler::src_

private

Definition at line 32 of file AlignmentPrescaler.h.

Referenced by produce().

edm::InputTag AlignmentPrescaler::srcQualityMap_

private

Definition at line 33 of file AlignmentPrescaler.h.

Referenced by produce().

int AlignmentPrescaler::totnhitspxl_

private

Definition at line 46 of file AlignmentPrescaler.h.

TTree* AlignmentPrescaler::tpresc_

private

Definition at line 39 of file AlignmentPrescaler.h.

Referenced by beginJob(), endJob(), and produce().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation