#include <CleanAndMergeProducer.h>

Inheritance diagram for CleanAndMergeProducer:

Public Member Functions
	CleanAndMergeProducer (const edm::ParameterSet &ps)

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

	~CleanAndMergeProducer () override

Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default

bool	hasAbilityToProduceInLumis () const final

bool	hasAbilityToProduceInRuns () const final

Private Attributes
std::string	bcCollection_

edm::EDGetTokenT< reco::SuperClusterCollection >	cleanScToken_

std::string	refScCollection_

std::string	scCollection_

edm::EDGetTokenT< reco::SuperClusterCollection >	uncleanScToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T... >	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 18 of file CleanAndMergeProducer.h.

Constructor & Destructor Documentation

CleanAndMergeProducer::CleanAndMergeProducer ( const edm::ParameterSet & ps )

Definition at line 53 of file CleanAndMergeProducer.cc.

References bcCollection_, cleanScToken_, edm::ParameterSet::getParameter(), refScCollection_, scCollection_, AlCaHLTBitMon_QueryRunRegistry::string, and uncleanScToken_.

 {
 
         // get the parameters
         // the cleaned collection:
         cleanScToken_   = 
             consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("cleanScInputTag"));
     uncleanScToken_ = 
         consumes<reco::SuperClusterCollection>(ps.getParameter<edm::InputTag>("uncleanScInputTag"));          
 
         // the names of the products to be produced:
         bcCollection_ = ps.getParameter<std::string>("bcCollection");
         scCollection_ = ps.getParameter<std::string>("scCollection");
         refScCollection_ = ps.getParameter<std::string>("refScCollection");
 
         std::map<std::string,double> providedParameters;  
         providedParameters.insert(std::make_pair("LogWeighted",ps.getParameter<bool>("posCalc_logweight")));
         providedParameters.insert(std::make_pair("T0_barl",ps.getParameter<double>("posCalc_t0")));
         providedParameters.insert(std::make_pair("W0",ps.getParameter<double>("posCalc_w0")));
         providedParameters.insert(std::make_pair("X0",ps.getParameter<double>("posCalc_x0")));
 
 
         // the products:
         produces< reco::BasicClusterCollection >(bcCollection_);
         produces< reco::SuperClusterCollection >(scCollection_);
         produces< reco::SuperClusterRefVector >(refScCollection_);
 
 }

CleanAndMergeProducer::~CleanAndMergeProducer ( )

override

Definition at line 82 of file CleanAndMergeProducer.cc.

82 {;}

Member Function Documentation

void CleanAndMergeProducer::produce	(	edm::Event &	evt,
		const edm::EventSetup &	es
	)

override

Definition at line 84 of file CleanAndMergeProducer.cc.

References bcCollection_, cleanScToken_, plotBeamSpotDB::first, edm::Event::getByToken(), mps_fire::i, edm::OrphanHandleBase::isValid(), LogDebug, LogTrace, eostools::move(), edm::RefVector< C, T, F >::push_back(), edm::PtrVector< T >::push_back(), edm::Event::put(), refScCollection_, scCollection_, SurveyInfoScenario_cff::seed, edm::RefVector< C, T, F >::size(), and uncleanScToken_.

 {
         // get the input collections
         // ______________________________________________________________________________________
         //
         // cluster collections:
 
         edm::Handle<reco::SuperClusterCollection> pCleanSC;
         edm::Handle<reco::SuperClusterCollection> pUncleanSC;
 
         evt.getByToken(cleanScToken_, pCleanSC);
         evt.getByToken(uncleanScToken_, pUncleanSC);
     
         //
         // collections are all taken now _____________________________________________________
         //
         //
         // the collections to be produced:
         reco::BasicClusterCollection basicClusters;
         reco::SuperClusterCollection superClusters;
         reco::SuperClusterRefVector * scRefs = new reco::SuperClusterRefVector;
 
         //
         // run over the uncleaned SC and check how many of them are matched to the cleaned ones
         // if you find a matched one, create a reference to the cleaned collection and store it
         // if you find an unmatched one, then keep all its basic clusters in the basic Clusters 
         // vector
         int uncleanSize =  pUncleanSC->size();
         int cleanSize =    pCleanSC->size();
 
         LogTrace("EcalCleaning") << "Size of Clean Collection: " << cleanSize 
                 << ", uncleanSize: " << uncleanSize  << std::endl;
         //
         // keep whether the SC in unique in the uncleaned collection
         std::vector<int> isUncleanOnly;     // 1 if unique in the uncleaned
         std::vector<int> basicClusterOwner; // contains the index of the SC that owns that BS
         std::vector<int> isSeed; // if this basic cluster is a seed it is 1
         for (int isc =0; isc< uncleanSize; ++isc) {
       reco::SuperClusterRef unscRef( pUncleanSC, isc );    
       const std::vector< std::pair<DetId, float> > & uhits = unscRef->hitsAndFractions();
       int uhitsSize = uhits.size();
       bool foundTheSame = false;
       for (int jsc=0; jsc < cleanSize; ++jsc) { // loop over the cleaned SC
         reco::SuperClusterRef cscRef( pCleanSC, jsc );
         const std::vector< std::pair<DetId, float> > & chits = cscRef->hitsAndFractions();
         int chitsSize =  chits.size();
         foundTheSame = true;
         if (unscRef->seed()->seed() == cscRef->seed()->seed() && chitsSize == uhitsSize) { 
           // if the clusters are exactly the same then because the clustering
           // algorithm works in a deterministic way, the order of the rechits
           // will be the same
           for (int i=0; i< chitsSize; ++i) {
         if (uhits[i].first != chits[i].first ) { foundTheSame=false;  break;}
           }
           if (foundTheSame) { // ok you have found it!
         // make the reference
         //scRefs->push_back( edm::Ref<reco::SuperClusterCollection>(pCleanSC, jsc) );
         scRefs->push_back( cscRef );
         isUncleanOnly.push_back(0);
         break;
           }
         }
       }
       if (not foundTheSame) {
         // mark it as unique in the uncleaned
         isUncleanOnly.push_back(1);
         // keep all its basic clusters
         for (reco::CaloCluster_iterator bciter = unscRef->clustersBegin(); bciter != unscRef->clustersEnd(); ++bciter) {
           // the basic clusters
           basicClusters.push_back(**bciter);
           basicClusterOwner.push_back(isc);
         }
       }
         }
         int bcSize =  basicClusters.size();
     
         LogDebug("EcalCleaning") << "Found cleaned SC: " << cleanSize <<  " uncleaned SC: " 
                 << uncleanSize << " from which " << scRefs->size() 
                 << " will become refs to the cleaned collection" ;
      
 
         // now you have the collection of basic clusters of the SC to be remain in the
         // in the clean collection, export them to the event
         // you will need to reread them later in order to make correctly the refs to the SC
         auto basicClusters_p = std::make_unique<reco::BasicClusterCollection>();
         basicClusters_p->assign(basicClusters.begin(), basicClusters.end());
         edm::OrphanHandle<reco::BasicClusterCollection> bccHandle =  
                 evt.put(std::move(basicClusters_p), bcCollection_);
         if (!(bccHandle.isValid())) {
       edm::LogWarning("MissingInput")<<"could not get a handle on the BasicClusterCollection!" << std::endl;
       return;
         }
 
         LogDebug("EcalCleaning")<< "Got the BasicClusters from the event again";
         // now you have to create again your superclusters
         // you run over the uncleaned SC, but now you know which of them are
         // the ones that are needed and which are their basic clusters
         for (int isc=0; isc< uncleanSize; ++isc) {
       if (isUncleanOnly[isc]==1) { // look for sc that are unique in the unclean collection
         // make the basic cluster collection
         reco::CaloClusterPtrVector clusterPtrVector;
         reco::CaloClusterPtr seed; // the seed is the basic cluster with the highest energy
         double energy = -1;
         for (int jbc=0; jbc< bcSize; ++jbc) {
           if (basicClusterOwner[jbc]==isc) {
         reco::CaloClusterPtr currentClu = reco::CaloClusterPtr(bccHandle, jbc);
         clusterPtrVector.push_back(currentClu);
         if (energy< currentClu->energy()) {
           energy = currentClu->energy(); seed = currentClu;
         }
           }
         }
         reco::SuperClusterRef unscRef( pUncleanSC, isc ); 
         reco::SuperCluster newSC(unscRef->energy(), unscRef->position(), seed, clusterPtrVector );
         superClusters.push_back(newSC);
       }
       
         }
         // export the collection of references to the clean collection
         std::unique_ptr<reco::SuperClusterRefVector> scRefs_p(scRefs);
         evt.put(std::move(scRefs_p), refScCollection_);
 
         // the collection of basic clusters is already in the event
         // the collection of uncleaned SC
         auto superClusters_p = std::make_unique<reco::SuperClusterCollection>();
         superClusters_p->assign(superClusters.begin(), superClusters.end());
         evt.put(std::move(superClusters_p), scCollection_);
         LogDebug("EcalCleaning")<< "Hybrid Clusters (Basic/Super) added to the Event! :-)";
 }