ClusterizerFP420.cc

Go to the documentation of this file.

// File: ClusterizerFP420.cc
// Date: 12.2006
// Description: ClusterizerFP420 for FP420
// Modifications: 
#include <memory>
#include <string>

#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/EDProducer.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Framework/interface/ESHandle.h"

#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/CommonDetUnit/interface/GeomDetType.h"
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/TrackerNumberingBuilder/interface/GeometricDet.h"

#include "DataFormats/Common/interface/DetSetVector.h"
//#include "SimGeneral/HepPDTRecord/interface/ParticleDataTable.h"


#include "RecoRomanPot/RecoFP420/interface/ClusterizerFP420.h"
#include "DataFormats/FP420Digi/interface/DigiCollectionFP420.h"
#include "DataFormats/FP420Cluster/interface/ClusterCollectionFP420.h"
#include "RecoRomanPot/RecoFP420/interface/ClusterNoiseFP420.h"

//Random Number
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "FWCore/Utilities/interface/RandomNumberGenerator.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "CLHEP/Random/RandomEngine.h"

#include <cstdlib> 
#include <iostream> 
using namespace std;

namespace cms
{
  ClusterizerFP420::ClusterizerFP420(const edm::ParameterSet& conf):
    conf_(conf) {
    
    
    std::string alias ( conf.getParameter<std::string>("@module_label") );
    
    produces<ClusterCollectionFP420>().setBranchAlias( alias );
    
    trackerContainers.clear();
    trackerContainers = conf.getParameter<std::vector<std::string> >("ROUList");
    
    verbosity = conf_.getUntrackedParameter<int>("VerbosityLevel");
    dn0   = conf_.getParameter<int>("NumberFP420Detectors");
    sn0   = conf_.getParameter<int>("NumberFP420Stations");
    pn0 = conf_.getParameter<int>("NumberFP420SPlanes");
    rn0 = 7;
    if (verbosity > 0) {
      std::cout << "Creating a ClusterizerFP420" << std::endl;
      std::cout << "ClusterizerFP420: dn0=" << dn0 << " sn0=" << sn0 << " pn0=" << pn0 <<  " rn0=" << rn0 << std::endl;
    }
    
    sClusterizerFP420_ = new FP420ClusterMain(conf_,dn0,sn0,pn0,rn0);
  }
  
  // Virtual destructor needed.
  ClusterizerFP420::~ClusterizerFP420() { 
    delete sClusterizerFP420_;
  }  
  
  //Get at the beginning
  void ClusterizerFP420::beginJob() {
    if (verbosity > 0) {
      std::cout << "BeginJob method " << std::endl;
    }
    //Getting Calibration data (Noises and BadElectrodes Flag)
    //    bool UseNoiseBadElectrodeFlagFromDB_=conf_.getParameter<bool>("UseNoiseBadElectrodeFlagFromDB");  
    //    if (UseNoiseBadElectrodeFlagFromDB_==true){
    //      iSetup.get<ClusterNoiseFP420Rcd>().get(noise);// AZ: do corrections for noise here
    //=========================================================
    // 
    // Debug: show noise for DetIDs
    //       ElectrodNoiseMapIterator mapit = noise->m_noises.begin();
    //       for (;mapit!=noise->m_noises.end();mapit++)
    //  {
    //    unsigned int detid = (*mapit).first;
    //    std::cout << "detid " <<  detid << " # Electrode " << (*mapit).second.size()<<std::endl;
    //    //ElectrodNoiseVector theElectrodVector =  (*mapit).second;     
    //    const ElectrodNoiseVector theElectrodVector =  noise->getElectrodNoiseVector(detid);
    
    
    //    int electrode=0;
    //    ElectrodNoiseVectorIterator iter=theElectrodVector.begin();
    //    //for(; iter!=theElectrodVector.end(); iter++)
    //    {
    //      std::cout << " electrode " << electrode++ << " =\t"
    //            << iter->getNoise()     << " \t" 
    //            << iter->getDisable()   << " \t" 
    //            << std::endl;         
    //    } 
    //       }
    //===========================================================
    //    }
  }
  
  void ClusterizerFP420::produce(edm::Event& iEvent, const edm::EventSetup& iSetup)
  {
    //  beginJob;
    // be lazy and include the appropriate namespaces
    using namespace edm; 
    using namespace std;   
    if (verbosity > 0) {
      std::cout << "ClusterizerFP420: produce" << std::endl;
    }
    
    // Get input
    //A
    //   edm::Handle<DigiCollectionFP420> icf_simhit;
    /*    
    Handle<DigiCollectionFP420> cf_simhit;
    std::vector<const DigiCollectionFP420 *> cf_simhitvec;
    for(uint32_t i = 0; i< trackerContainers.size();i++){
      iEvent.getByLabel( trackerContainers[i], cf_simhit);
      cf_simhitvec.push_back(cf_simhit.product());   }
    std::auto_ptr<DigiCollectionFP420 > digis(new DigiCollectionFP420(cf_simhitvec));

    std::vector<HDigiFP420> input;
    DigiCollectionFP420::iterator isim;
    for (isim=digis->begin(); isim!= digis->end();isim++) {
      input.push_back(*isim);
     }
*/    
    //B
       
    Handle<DigiCollectionFP420> input;
    try{
      //      iEvent.getByLabel( "FP420Digi" , digis);
      iEvent.getByLabel( trackerContainers[0] , input);
    } catch(const Exception&)
      {
    // in principal, should never happen, as it's taken care of by Framework
    throw cms::Exception("InvalidReference")
      << "Invalid reference to DigiCollectionFP420 \n";
      }
    
    if (verbosity > 0) {
      std::cout << "ClusterizerFP420: OK1" << std::endl;
    }

    // Step C: create empty output collection
    std::auto_ptr<ClusterCollectionFP420> soutput(new ClusterCollectionFP420);
    /*    
   std::vector<SimVertex> input;
   Handle<DigiCollectionFP420> digis;
   iEvent.getByLabel("FP420Digi",digis);
   input.insert(input.end(),digis->begin(),digis->end());



    std::vector<HDigiFP420> input;
    for(std::vector<HDigiFP420>::const_iterator vsim=digis->begin();
    vsim!=digis->end(); ++vsim){
      input.push_back(*vsim);
    }
   theSimTracks.insert(theSimTracks.end(),SimTk->begin(),SimTk->end());
*/
    //     std::vector<HDigiFP420> input;
    //   DigiCollectionFP420  input;
       //input.push_back(digis);
    //   input.insert(input.end(), digis->begin(), digis->end());

    /*
    std::vector<HDigiFP420> input;
    input.clear();
    DigiCollectionFP420::ContainerIterator sort_begin = digis->begin();
    DigiCollectionFP420::ContainerIterator sort_end = digis->end();
    for ( ;sort_begin != sort_end; ++sort_begin ) {
      input.push_back(*sort_begin);
    } // for
*/
    
    
    //    put zero to container info from the beginning (important! because not any detID is updated with coming of new event     !!!!!!   
    // clean info of container from previous event
      for (int det=1; det<dn0; det++) {
    for (int sector=1; sector<sn0; sector++) {
      for (int zmodule=1; zmodule<pn0; zmodule++) {
        for (int zside=1; zside<rn0; zside++) {
          // intindex is a continues numbering of FP420
          unsigned int detID = theFP420NumberingScheme->FP420NumberingScheme::packMYIndex(rn0, pn0, sn0, det, zside, sector, zmodule);
          std::vector<ClusterFP420> collector;
          collector.clear();
          ClusterCollectionFP420::Range inputRange;
          inputRange.first = collector.begin();
          inputRange.second = collector.end();
          
          soutput->putclear(inputRange,detID);
          
        }//for
      }//for
    }//for
      }//for
      
    
    //                                                                                                                      !!!!!!   
    // if we want to keep Cluster container/Collection for one event --->   uncomment the line below and vice versa
    soutput->clear();   //container_.clear() --> start from the beginning of the container
    
    //                                RUN now:                                                                                 !!!!!!     
    //  sClusterizerFP420_.run(input, soutput, noise);
    if (verbosity > 0) {
      std::cout << "ClusterizerFP420: OK2" << std::endl;
    }
      sClusterizerFP420_->run(input, soutput, noise);

    if (verbosity > 0) {
      std::cout << "ClusterizerFP420: OK3" << std::endl;
    }

    //  if(collectorZS.data.size()>0){

    //  std::cout <<"=======           ClusterizerFP420:                    end of produce     " <<  std::endl;
    
    // Step D: write output to file
    iEvent.put(soutput);
    if (verbosity > 0) {
      std::cout << "ClusterizerFP420: OK4" << std::endl;
    }
  }//produce
  
} // namespace cms