TreeProducerCalibSimul Class Reference

#include <TreeProducerCalibSimul.h>

Inheritance diagram for TreeProducerCalibSimul:

Public Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob ()
virtual void	endJob ()
	TreeProducerCalibSimul (const edm::ParameterSet &)
	~TreeProducerCalibSimul ()
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
SerialTaskQueue *	globalLuminosityBlocksQueue ()
SerialTaskQueue *	globalRunsQueue ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
	~EDAnalyzer () 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 &&)=default
	EDConsumerBase (EDConsumerBase const &)=delete
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 &	operator= (EDConsumerBase &&)=default
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
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 Attributes
std::string	EBRecHitCollection_
std::string	eventHeaderCollection_
std::string	eventHeaderProducer_
std::string	hodoRecInfoCollection_
std::string	hodoRecInfoProducer_
TreeMatrixCalib *	myTree
int	noHeader
int	noHits
int	noHodo
int	noTdc
double	posCluster_
std::string	RecHitProducer_
std::string	rootfile_
std::string	tdcRecInfoCollection_
std::string	tdcRecInfoProducer_
int	tot_events
int	tot_events_ok
std::string	txtfile_
int	xtalInBeam

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &)
static bool	wantsGlobalLuminosityBlocks ()
static bool	wantsGlobalRuns ()
static bool	wantsStreamLuminosityBlocks ()
static bool	wantsStreamRuns ()
Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 59 of file TreeProducerCalibSimul.h.

Constructor & Destructor Documentation

TreeProducerCalibSimul()

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

explicit

Definition at line 8 of file TreeProducerCalibSimul.cc.

{
  // now do what ever initialization is needed
  xtalInBeam                 = iConfig.getUntrackedParameter<int>("xtalInBeam",-1000);
  rootfile_                  = iConfig.getUntrackedParameter<std::string>("rootfile","mySimMatrixTree.root");
  txtfile_                   = iConfig.getUntrackedParameter<std::string>("txtfile", "mySimMatrixTree.txt");
  EBRecHitCollection_        = iConfig.getParameter<std::string>("EBRecHitCollection");
  RecHitProducer_            = iConfig.getParameter<std::string>("RecHitProducer");
  hodoRecInfoCollection_     = iConfig.getParameter<std::string>("hodoRecInfoCollection");
  hodoRecInfoProducer_       = iConfig.getParameter<std::string>("hodoRecInfoProducer");
  tdcRecInfoCollection_      = iConfig.getParameter<std::string>("tdcRecInfoCollection");
  tdcRecInfoProducer_        = iConfig.getParameter<std::string>("tdcRecInfoProducer");
  eventHeaderCollection_     = iConfig.getParameter<std::string>("eventHeaderCollection");
  eventHeaderProducer_       = iConfig.getParameter<std::string>("eventHeaderProducer");
 
  // summary
  cout << endl;
  cout <<"Constructor" << endl;
  cout << endl;
  cout << "TreeProducerCalibSimul"    << endl;
  cout << "xtal in beam = " << xtalInBeam << endl;
  cout <<"Fetching hitCollection: "   << EBRecHitCollection_.c_str()     << " prod by " << RecHitProducer_.c_str()      <<endl; 
  cout <<"Fetching hodoCollection: "  << hodoRecInfoCollection_.c_str()  << " prod by " << hodoRecInfoProducer_.c_str() <<endl;
  cout <<"Fetching tdcCollection: "   << tdcRecInfoCollection_.c_str()   << " prod by " << tdcRecInfoProducer_.c_str()  <<endl;
  cout <<"Fetching evHeaCollection: " << eventHeaderCollection_.c_str()  << " prod by " << eventHeaderProducer_.c_str() <<endl;
  cout << endl;
}

References gather_cfg::cout, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

~TreeProducerCalibSimul()

TreeProducerCalibSimul::~TreeProducerCalibSimul

(

)

Definition at line 39 of file TreeProducerCalibSimul.cc.

{
  cout << endl;
  cout << "Deleting" << endl;
  cout << endl;
 
  delete myTree;
}

References gather_cfg::cout.

Member Function Documentation

analyze()

void TreeProducerCalibSimul::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

virtual

Implements edm::EDAnalyzer.

Definition at line 95 of file TreeProducerCalibSimul.cc.

{
  using namespace edm;
  using namespace cms;
 
  // counting events
  tot_events++;
 
  if ( tot_events%5000 == 0){ cout << "event " << tot_events << endl;}
  
 
  // ---------------------------------------------------------------------
  // taking what I need: hits
  Handle< EBRecHitCollection > pEBRecHits ;
  const EBRecHitCollection*  EBRecHits = 0 ;
  //try {
  iEvent.getByLabel (RecHitProducer_, EBRecHitCollection_, pEBRecHits) ;
  EBRecHits = pEBRecHits.product(); 
  //} catch ( std::exception& ex ) {
  //std::cout<<"Error! can't get the product " << EBRecHitCollection_.c_str () << std::endl ;
  //std::cerr<<"Error! can't get the product " << EBRecHitCollection_.c_str () << std::endl ;
  //}
 
  // taking what I need: hodoscopes
  Handle<EcalTBHodoscopeRecInfo> pHodo;
  const EcalTBHodoscopeRecInfo* recHodo=0;
  //try {
  iEvent.getByLabel( hodoRecInfoProducer_, hodoRecInfoCollection_, pHodo);
  recHodo = pHodo.product();
  //} catch ( std::exception& ex ) {
  //std::cout<<"Error! can't get the product "<<hodoRecInfoCollection_.c_str() << std::endl;
  //std::cerr<<"Error! can't get the product "<< hodoRecInfoCollection_.c_str() << std::endl;
  //}
  
  // taking what I need: tdc
  Handle<EcalTBTDCRecInfo> pTDC;
  const EcalTBTDCRecInfo* recTDC=0;
  //try {
  iEvent.getByLabel( tdcRecInfoProducer_, tdcRecInfoCollection_, pTDC);
  recTDC = pTDC.product(); 
  //} catch ( std::exception& ex ) {
  //std::cout<<"Error! can't get the product " << tdcRecInfoCollection_.c_str() << std::endl;
  //std::cerr<<"Error! can't get the product " << tdcRecInfoCollection_.c_str() << std::endl;
  //}
 
  // taking what I need: event header
  Handle<EcalTBEventHeader> pEventHeader;
  const EcalTBEventHeader* evtHeader=0;
  //try {
  iEvent.getByLabel( eventHeaderProducer_ , pEventHeader );
  evtHeader = pEventHeader.product(); 
  //} catch ( std::exception& ex ) {
  //std::cout << "Error! can't get the event header " <<std::endl;
  //std::cerr << "Error! can't get the event header " <<std::endl;
  //}
   
  // checking everything is there and fine
  if ( (!EBRecHits) || (EBRecHits->size() == 0)){ noHits++;  return; }
  if (!recTDC)    { noTdc++;     return; }        
  if (!recHodo)   { noHodo++;    return; }                  
  if (!evtHeader) { noHeader++;  return; }                
  tot_events_ok++;
 
 
 
  // ---------------------------------------------------------------------
  // info on the event
  int run   = -999;
  int tbm   = -999;
  int event = evtHeader->eventNumber();
 
  // ---------------------------------------------------------------------
  // xtal-in-beam  
  int nomXtalInBeam  = -999;
  int nextXtalInBeam = -999;
 
  EBDetId xtalInBeamId(1,xtalInBeam, EBDetId::SMCRYSTALMODE); 
  if (xtalInBeamId==EBDetId(0)){ return; }
  int mySupCry = xtalInBeamId.ic();
  int mySupEta = xtalInBeamId.ieta();
  int mySupPhi = xtalInBeamId.iphi();
 
  
  // ---------------------------------------------------------------------
  // hodoscope information
  double x  = recHodo->posX();
  double y  = recHodo->posY();
  double sx = recHodo->slopeX();
  double sy = recHodo->slopeY();
  double qx = recHodo->qualX();
  double qy = recHodo->qualY();
    
 
  // ---------------------------------------------------------------------
  // tdc information
  double tdcOffset = recTDC->offset();
  
 
  // ---------------------------------------------------------------------
  // Find EBDetId in a 7x7 Matrix
  EBDetId Xtals7x7[49];
  double energy[49];
  int crystal[49];
  int allMatrix = 1; 
  for (unsigned int icry=0; icry<49; icry++){
    unsigned int row    = icry/7;
    unsigned int column = icry%7;
    //try
    //  {
    Xtals7x7[icry]=EBDetId(xtalInBeamId.ieta()+column-3, xtalInBeamId.iphi()+row-3, EBDetId::ETAPHIMODE);
    
    if ( Xtals7x7[icry].ism() == 1){ 
      energy[icry] = EBRecHits->find(Xtals7x7[icry])->energy(); 
      crystal[icry] = Xtals7x7[icry].ic();
    } else {
      energy[icry] = -100.;
      crystal[icry] = -100;
      allMatrix = 0; 
    }
    /* 
    catch (...)
      {
    // can not construct 7x7 matrix 
    energy[icry]  = -100.; 
    crystal[icry] = -100;
    allMatrix = 0;  
      }
    */
  }
 
 
  // ---------------------------------------------------------------------
  // Looking for the max energy crystal
  double maxEne      = -999.;
  int maxEneCry      = 9999;
  int maxEneInMatrix = -999;
  for (int ii=0; ii<49; ii++){ if (energy[ii] > maxEne){ 
    maxEne         = energy[ii]; 
    maxEneCry      = crystal[ii]; 
    maxEneInMatrix = ii;} 
  }  
  
 
 
  // Position reconstruction - skipped here
  double Xcal   = -999.;
  double Ycal   = -999.;
 
  // filling the tree
  myTree->fillInfo(run, event, mySupCry, maxEneCry, nomXtalInBeam, nextXtalInBeam, mySupEta, mySupPhi, tbm, x, y, Xcal, Ycal, sx, sy, qx, qy, tdcOffset, allMatrix, energy, crystal);
  myTree->store();
}

References gather_cfg::cout, filterRecHits_cfi::EBRecHits, HCALHighEnergyHPDFilter_cfi::energy, EBDetId::ETAPHIMODE, EcalTBEventHeader::eventNumber(), EBDetId::ic(), EBDetId::ieta(), iEvent, cuy::ii, EBDetId::iphi(), ecalpyutils::ism(), HcalResponse_cfi::maxEne, EcalTBTDCRecInfo::offset(), EcalTBHodoscopeRecInfo::posX(), EcalTBHodoscopeRecInfo::posY(), edm::Handle< T >::product(), EcalTBHodoscopeRecInfo::qualX(), EcalTBHodoscopeRecInfo::qualY(), writedatasetfile::run, EcalTBHodoscopeRecInfo::slopeX(), EcalTBHodoscopeRecInfo::slopeY(), EBDetId::SMCRYSTALMODE, fftjetcommon_cfi::sx, and fftjetcommon_cfi::sy.

beginJob()

void TreeProducerCalibSimul::beginJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 52 of file TreeProducerCalibSimul.cc.

{
  cout << endl;
  cout << "BeginJob" << endl;
  cout << endl;
 
  // tree
  myTree = new TreeMatrixCalib(rootfile_.c_str());
 
  // counters
  tot_events      = 0;
  tot_events_ok   = 0;
  noHits          = 0;
  noHodo          = 0;
  noTdc           = 0;
  noHeader        = 0;
} 

References gather_cfg::cout.

endJob()

void TreeProducerCalibSimul::endJob ( void  )

virtual

Reimplemented from edm::EDAnalyzer.

Definition at line 74 of file TreeProducerCalibSimul.cc.

{
  cout << endl;
  cout << "EndJob" << endl;
  cout << endl;
 
  ofstream *MyOut = new ofstream(txtfile_.c_str());   
  *MyOut << "total events: "                                   << tot_events      << endl;
  *MyOut << "events skipped because of no hits: "              << noHits          << endl;
  *MyOut << "events skipped because of no hodos: "             << noHodo          << endl;
  *MyOut << "events skipped because of no tdc: "               << noTdc           << endl;
  *MyOut << "events skipped because of no header: "            << noHeader        << endl;
  *MyOut << "total OK events (passing the basic selection): "  << tot_events_ok   << endl;
  MyOut->close();
  delete MyOut;
}

References gather_cfg::cout.

Member Data Documentation

EBRecHitCollection_

std::string TreeProducerCalibSimul::EBRecHitCollection_

private

Definition at line 72 of file TreeProducerCalibSimul.h.

eventHeaderCollection_

std::string TreeProducerCalibSimul::eventHeaderCollection_

private

Definition at line 78 of file TreeProducerCalibSimul.h.

eventHeaderProducer_

std::string TreeProducerCalibSimul::eventHeaderProducer_

private

Definition at line 79 of file TreeProducerCalibSimul.h.

hodoRecInfoCollection_

std::string TreeProducerCalibSimul::hodoRecInfoCollection_

private

Definition at line 74 of file TreeProducerCalibSimul.h.

hodoRecInfoProducer_

std::string TreeProducerCalibSimul::hodoRecInfoProducer_

private

Definition at line 75 of file TreeProducerCalibSimul.h.

myTree

TreeMatrixCalib* TreeProducerCalibSimul::myTree

private

Definition at line 82 of file TreeProducerCalibSimul.h.

noHeader

int TreeProducerCalibSimul::noHeader

private

Definition at line 90 of file TreeProducerCalibSimul.h.

noHits

int TreeProducerCalibSimul::noHits

private

Definition at line 87 of file TreeProducerCalibSimul.h.

noHodo

int TreeProducerCalibSimul::noHodo

private

Definition at line 88 of file TreeProducerCalibSimul.h.

noTdc

int TreeProducerCalibSimul::noTdc

private

Definition at line 89 of file TreeProducerCalibSimul.h.

posCluster_

double TreeProducerCalibSimul::posCluster_

private

Definition at line 80 of file TreeProducerCalibSimul.h.

RecHitProducer_

std::string TreeProducerCalibSimul::RecHitProducer_

private

Definition at line 73 of file TreeProducerCalibSimul.h.

rootfile_

std::string TreeProducerCalibSimul::rootfile_

private

Definition at line 70 of file TreeProducerCalibSimul.h.

tdcRecInfoCollection_

std::string TreeProducerCalibSimul::tdcRecInfoCollection_

private

Definition at line 76 of file TreeProducerCalibSimul.h.

tdcRecInfoProducer_

std::string TreeProducerCalibSimul::tdcRecInfoProducer_

private

Definition at line 77 of file TreeProducerCalibSimul.h.

tot_events

int TreeProducerCalibSimul::tot_events

private

Definition at line 85 of file TreeProducerCalibSimul.h.

tot_events_ok

int TreeProducerCalibSimul::tot_events_ok

private

Definition at line 86 of file TreeProducerCalibSimul.h.

txtfile_

std::string TreeProducerCalibSimul::txtfile_

private

Definition at line 71 of file TreeProducerCalibSimul.h.

xtalInBeam

int TreeProducerCalibSimul::xtalInBeam

private

Definition at line 84 of file TreeProducerCalibSimul.h.