MCPizeroAnalyzer Class Reference

#include <MCPizeroAnalyzer.h>

Inheritance diagram for MCPizeroAnalyzer:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	beginJob () override
void	endJob () override
	MCPizeroAnalyzer (const edm::ParameterSet &)
	~MCPizeroAnalyzer () override
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () 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
float	etaTransformation (float a, float b)
float	phiNormalization (float &a)

Private Attributes
TFile *	fOutputFile_
std::string	fOutputFileName_
TH1F *	h_BremEnergy_
TH1F *	h_BremFrac_
TH2F *	h_EleEvsPhoE_
TH1F *	h_MCConvPhoE_
TH1F *	h_MCConvPhoEta_
TH1F *	h_MCConvPhoPhi_
TH1F *	h_MCConvPhoR_
TH1F *	h_MCEleE_
TH1F *	h_MCEleEta_
TH1F *	h_MCElePhi_
TH1F *	h_MCPhoE_
TH1F *	h_MCPhoEta_
TH1F *	h_MCPhoPhi_
TH1F *	h_MCPiz1ConEta_
TH1F *	h_MCPiz2ConEta_
TH1F *	h_MCPizE_
TH1F *	h_MCPizEta_
TH1F *	h_MCPizMass1_
TH1F *	h_MCPizMass2_
TH1F *	h_MCPizPhi_
TH1F *	h_MCPizUnEta_
std::string	HepMCLabel
double	mcEta_
double	mcPhi_
	global variable for the MC photon More...
int	nEvt_
int	nMatched_
std::string	SimHitLabel
std::string	SimTkLabel
std::string	SimVtxLabel
PizeroMCTruthFinder *	thePizeroMCTruthFinder_

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
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 22 of file MCPizeroAnalyzer.h.

Constructor & Destructor Documentation

MCPizeroAnalyzer::MCPizeroAnalyzer ( const edm::ParameterSet &  pset )

explicit

Definition at line 45 of file MCPizeroAnalyzer.cc.

   : fOutputFileName_( pset.getUntrackedParameter<string>("HistOutFile",std::string("TestConversions.root")) ),
     fOutputFile_(nullptr)
{


}

MCPizeroAnalyzer::~MCPizeroAnalyzer ( )

override

Definition at line 55 of file MCPizeroAnalyzer.cc.

References thePizeroMCTruthFinder_.

                                    {


  delete thePizeroMCTruthFinder_;

}

Member Function Documentation

void MCPizeroAnalyzer::analyze ( const edm::Event &  e, const edm::EventSetup &    )

override

Definition at line 162 of file MCPizeroAnalyzer.cc.

References gather_cfg::cout, PizeroMCTruthFinder::find(), edm::Event::getByLabel(), h_BremEnergy_, h_BremFrac_, h_EleEvsPhoE_, h_MCConvPhoE_, h_MCConvPhoEta_, h_MCConvPhoPhi_, h_MCConvPhoR_, h_MCEleE_, h_MCEleEta_, h_MCElePhi_, h_MCPhoE_, h_MCPhoEta_, h_MCPhoPhi_, h_MCPiz1ConEta_, h_MCPiz2ConEta_, h_MCPizE_, h_MCPizEta_, h_MCPizMass1_, h_MCPizPhi_, h_MCPizUnEta_, edm::EventBase::id(), nEvt_, mathSSE::sqrt(), and thePizeroMCTruthFinder_.

{


  using namespace edm;
  //UNUSED const float etaPhiDistance=0.01;
  // Fiducial region
  //UNUSED const float TRK_BARL =0.9;
  //UNUSED const float BARL = 1.4442; // DAQ TDR p.290
  //UNUSED const float END_LO = 1.566;
  //UNUSED const float END_HI = 2.5;
 // Electron mass
  //UNUSED const Float_t mElec= 0.000511;


  nEvt_++;
  LogInfo("MCPizeroAnalyzer") << "MCPizeroAnalyzer Analyzing event number: " << e.id() << " Global Counter " << nEvt_ <<"\n";
  //  LogDebug("MCPizeroAnalyzer") << "MCPizeroAnalyzer Analyzing event number: "  << e.id() << " Global Counter " << nEvt_ <<"\n";
  std::cout << "MCPizeroAnalyzer Analyzing event number: "  << e.id() << " Global Counter " << nEvt_ <<"\n";

  std::cout  << " MCPizeroAnalyzer Looking for MC truth " << "\n";

  //get simtrack info
  std::vector<SimTrack> theSimTracks;
  std::vector<SimVertex> theSimVertices;

  edm::Handle<SimTrackContainer> SimTk;
  edm::Handle<SimVertexContainer> SimVtx;
  e.getByLabel("g4SimHits",SimTk);
  e.getByLabel("g4SimHits",SimVtx);

  theSimTracks.insert(theSimTracks.end(),SimTk->begin(),SimTk->end());
  theSimVertices.insert(theSimVertices.end(),SimVtx->begin(),SimVtx->end());
  std::cout << " MCPizeroAnalyzer This Event has " <<  theSimTracks.size() << " sim tracks " << std::endl;
  std::cout << " MCPizeroAnalyzer This Event has " <<  theSimVertices.size() << " sim vertices " << std::endl;
  if (  theSimTracks.empty() ) std::cout << " Event number " << e.id() << " has NO sim tracks " << std::endl;


  std::vector<PizeroMCTruth> MCPizeroeros=thePizeroMCTruthFinder_->find (theSimTracks,  theSimVertices);
  std::cout << " MCPizeroAnalyzer MCPizeroeros size " <<  MCPizeroeros.size() << std::endl;

  for ( std::vector<PizeroMCTruth>::const_iterator iPiz=MCPizeroeros.begin(); iPiz !=MCPizeroeros.end(); ++iPiz ){

    h_MCPizE_->Fill  ( (*iPiz).fourMomentum().e() );
    h_MCPizEta_->Fill  ( (*iPiz).fourMomentum().pseudoRapidity() );
    h_MCPizPhi_->Fill  ( (*iPiz).fourMomentum().phi() );

    std::vector<PhotonMCTruth> mcPhotons=(*iPiz).photons();
    std::cout << " MCPizeroAnalyzer mcPhotons size " << mcPhotons.size() << std::endl;

    float px = mcPhotons[0].fourMomentum().x() + mcPhotons[1].fourMomentum().x();
    float py = mcPhotons[0].fourMomentum().y() + mcPhotons[1].fourMomentum().y();
    float pz = mcPhotons[0].fourMomentum().z() + mcPhotons[1].fourMomentum().z();
    float e  = mcPhotons[0].fourMomentum().e() + mcPhotons[1].fourMomentum().e();
    float invM =  sqrt( e*e - px*px -py*py - pz*pz)*1000;
    h_MCPizMass1_ ->Fill (invM);

    int converted=0;
    for ( std::vector<PhotonMCTruth>::const_iterator iPho=mcPhotons.begin(); iPho !=mcPhotons.end(); ++iPho ){
      h_MCPhoE_->Fill  ( (*iPho).fourMomentum().e() );
      h_MCPhoEta_->Fill  ( (*iPho).fourMomentum().pseudoRapidity() );
      h_MCPhoPhi_->Fill  ( (*iPho).fourMomentum().phi() );
      if (  (*iPho).isAConversion() ) {

        converted++;

    h_MCConvPhoE_->Fill  ( (*iPho).fourMomentum().e() );
    h_MCConvPhoEta_->Fill  ( (*iPho).fourMomentum().pseudoRapidity() );
    h_MCConvPhoPhi_->Fill  ( (*iPho).fourMomentum().phi() );
    h_MCConvPhoR_->Fill  ( (*iPho).vertex().perp() );

    std::vector<ElectronMCTruth> mcElectrons=(*iPho).electrons();
    std::cout << " MCPizeroAnalyzer mcElectrons size " << mcElectrons.size() << std::endl;

    for ( std::vector<ElectronMCTruth>::const_iterator iEl=mcElectrons.begin(); iEl !=mcElectrons.end(); ++iEl ){

          if ( (*iEl).fourMomentum().e()  < 30 ) continue;
      h_MCEleE_->Fill  ( (*iEl).fourMomentum().e() );
      h_MCEleEta_->Fill  ( (*iEl).fourMomentum().pseudoRapidity() );
      h_MCElePhi_->Fill  ( (*iEl).fourMomentum().phi() );


      h_EleEvsPhoE_->Fill ( (*iPho).fourMomentum().e(), (*iEl).fourMomentum().e() );

      float totBrem=0;
      for ( unsigned int iBrem=0; iBrem < (*iEl).bremVertices().size(); ++iBrem )
        totBrem +=  (*iEl).bremMomentum()[iBrem].e();

      h_BremFrac_->Fill( totBrem/(*iEl).fourMomentum().e() );
      h_BremEnergy_->Fill (  totBrem  );
    }

      }

    }


      if ( converted > 0 ) {
    h_MCPiz1ConEta_->Fill  ( (*iPiz).fourMomentum().pseudoRapidity() );
        if ( converted==2) h_MCPiz2ConEta_->Fill  ( (*iPiz).fourMomentum().pseudoRapidity() );
      } else {
    h_MCPizUnEta_->Fill  ( (*iPiz).fourMomentum().pseudoRapidity() );
      }



  }




}

void MCPizeroAnalyzer::beginJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 63 of file MCPizeroAnalyzer.cc.

References fOutputFile_, fOutputFileName_, h_BremEnergy_, h_BremFrac_, h_EleEvsPhoE_, h_MCConvPhoE_, h_MCConvPhoEta_, h_MCConvPhoPhi_, h_MCConvPhoR_, h_MCEleE_, h_MCEleEta_, h_MCElePhi_, h_MCPhoE_, h_MCPhoEta_, h_MCPhoPhi_, h_MCPiz1ConEta_, h_MCPiz2ConEta_, h_MCPizE_, h_MCPizEta_, h_MCPizMass1_, h_MCPizMass2_, h_MCPizPhi_, h_MCPizUnEta_, nEvt_, mathSSE::return(), and thePizeroMCTruthFinder_.

{


  nEvt_=0;

  thePizeroMCTruthFinder_ = new PizeroMCTruthFinder();

  fOutputFile_   = new TFile( fOutputFileName_.c_str(), "RECREATE" ) ;

  h_MCPizE_ = new TH1F("MCPizE","MC piz energy",100,0.,200.);
  h_MCPizPhi_ = new TH1F("MCPizPhi","MC piz phi",40,-3.14, 3.14);
  h_MCPizEta_ = new TH1F("MCPizEta","MC piz eta",40,-3., 3.);
  h_MCPizUnEta_ = new TH1F("MCPizUnEta","MC un piz eta",40,-3., 3.);
  h_MCPiz1ConEta_ = new TH1F("MCPiz1ConEta","MC con piz eta: at least one converted photon",40,-3., 3.);
  h_MCPiz2ConEta_ = new TH1F("MCPiz2ConEta","MC con piz eta: two converted photons",40,-3., 3.);
  h_MCPizMass1_= new TH1F("MCPizMass1","Piz mass unconverted ",100, 0., 200);
  h_MCPizMass2_= new TH1F("MCPizMass2","Piz mass converted ",100, 0., 200);

  // All Photons from Pizeros
  h_MCPhoE_ = new TH1F("MCPhoE","MC photon energy",100,0.,200.);
  h_MCPhoPhi_ = new TH1F("MCPhoPhi","MC photon phi",40,-3.14, 3.14);
  h_MCPhoEta_ = new TH1F("MCPhoEta","MC photon eta",40,-3., 3.);

 // Converted photons
  h_MCConvPhoE_ = new TH1F("MCConvPhoE","MC converted photon energy",100,0.,200.);
  h_MCConvPhoPhi_ = new TH1F("MCConvPhoPhi","MC converted photon phi",40,-3.14, 3.14);
  h_MCConvPhoEta_ = new TH1F("MCConvPhoEta","MC converted photon eta",40,-3., 3.);
  h_MCConvPhoR_ = new TH1F("MCConvPhoR","MC converted photon R",120,0.,120.);
  // Electrons from converted photons
  h_MCEleE_ = new TH1F("MCEleE","MC ele energy",100,0.,200.);
  h_MCElePhi_ = new TH1F("MCElePhi","MC ele phi",40,-3.14, 3.14);
  h_MCEleEta_ = new TH1F("MCEleEta","MC ele eta",40,-3., 3.);
  h_BremFrac_ = new TH1F("bremFrac","brem frac ", 50, 0., 1.);
  h_BremEnergy_ = new TH1F("bremE","Brem energy",100,0.,200.);


  h_EleEvsPhoE_ = new TH2F ("eleEvsPhoE","eleEvsPhoE",100,0.,200.,100,0.,200.);


  return ;
}

void MCPizeroAnalyzer::endJob ( void  )

overridevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 279 of file MCPizeroAnalyzer.cc.

References gather_cfg::cout, fOutputFile_, nEvt_, and mathSSE::return().

{




   fOutputFile_->Write() ;
   fOutputFile_->Close() ;

   edm::LogInfo("MCPizeroAnalyzer") << "Analyzed " << nEvt_  << "\n";
   std::cout  << "MCPizeroAnalyzer::endJob Analyzed " << nEvt_ << " events " << "\n";

   return ;
}

float MCPizeroAnalyzer::etaTransformation ( float  a, float  b  )

private

Definition at line 108 of file MCPizeroAnalyzer.cc.

References ETA, etaBarrelEndcap, cmsBatch::log, PI, R_ECAL, funct::tan(), and Z_Endcap.

                                                                              {

//---Definitions
    const float PI    = 3.1415927;
    //UNUSED const float TWOPI = 2.0*PI;

//---Definitions for ECAL
    const float R_ECAL           = 136.5;
    const float Z_Endcap         = 328.0;
    const float etaBarrelEndcap  = 1.479;

//---ETA correction

    float Theta = 0.0  ;
        float ZEcal = R_ECAL*sinh(EtaParticle)+Zvertex;

    if(ZEcal != 0.0) Theta = atan(R_ECAL/ZEcal);
    if(Theta<0.0) Theta = Theta+PI ;
    float ETA = - log(tan(0.5*Theta));

    if( fabs(ETA) > etaBarrelEndcap )
      {
       float Zend = Z_Endcap ;
       if(EtaParticle<0.0 )  Zend = -Zend ;
       float Zlen = Zend - Zvertex ;
       float RR = Zlen/sinh(EtaParticle);
       Theta = atan(RR/Zend);
       if(Theta<0.0) Theta = Theta+PI ;
       ETA = - log(tan(0.5*Theta));
      }
//---Return the result
        return ETA;
//---end
}

float MCPizeroAnalyzer::phiNormalization ( float &  a )

private

Definition at line 143 of file MCPizeroAnalyzer.cc.

References phi, PI, and TWOPI.

{
//---Definitions
 const float PI    = 3.1415927;
 const float TWOPI = 2.0*PI;


 if(phi >  PI) {phi = phi - TWOPI;}
 if(phi < -PI) {phi = phi + TWOPI;}

 //  cout << " Float_t PHInormalization out " << PHI << endl;
 return phi;

}

Member Data Documentation

TFile* MCPizeroAnalyzer::fOutputFile_

private

Definition at line 49 of file MCPizeroAnalyzer.h.

Referenced by beginJob(), and endJob().

std::string MCPizeroAnalyzer::fOutputFileName_

private

Definition at line 48 of file MCPizeroAnalyzer.h.

Referenced by beginJob().

TH1F* MCPizeroAnalyzer::h_BremEnergy_

private

Definition at line 80 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_BremFrac_

private

Definition at line 79 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH2F* MCPizeroAnalyzer::h_EleEvsPhoE_

private

Definition at line 82 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCConvPhoE_

private

Definition at line 87 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCConvPhoEta_

private

Definition at line 88 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCConvPhoPhi_

private

Definition at line 89 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCConvPhoR_

private

Definition at line 90 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCEleE_

private

Definition at line 76 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCEleEta_

private

Definition at line 77 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCElePhi_

private

Definition at line 78 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPhoE_

private

Definition at line 84 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPhoEta_

private

Definition at line 85 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPhoPhi_

private

Definition at line 86 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPiz1ConEta_

private

Definition at line 70 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPiz2ConEta_

private

Definition at line 71 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPizE_

private

Definition at line 67 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPizEta_

private

Definition at line 68 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPizMass1_

private

Definition at line 73 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPizMass2_

private

Definition at line 74 of file MCPizeroAnalyzer.h.

Referenced by beginJob().

TH1F* MCPizeroAnalyzer::h_MCPizPhi_

private

Definition at line 72 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

TH1F* MCPizeroAnalyzer::h_MCPizUnEta_

private

Definition at line 69 of file MCPizeroAnalyzer.h.

Referenced by analyze(), and beginJob().

std::string MCPizeroAnalyzer::HepMCLabel

private

Definition at line 61 of file MCPizeroAnalyzer.h.

double MCPizeroAnalyzer::mcEta_

private

Definition at line 59 of file MCPizeroAnalyzer.h.

double MCPizeroAnalyzer::mcPhi_

private

global variable for the MC photon

Definition at line 58 of file MCPizeroAnalyzer.h.

int MCPizeroAnalyzer::nEvt_

private

Definition at line 54 of file MCPizeroAnalyzer.h.

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

int MCPizeroAnalyzer::nMatched_

private

Definition at line 55 of file MCPizeroAnalyzer.h.

std::string MCPizeroAnalyzer::SimHitLabel

private

Definition at line 64 of file MCPizeroAnalyzer.h.

std::string MCPizeroAnalyzer::SimTkLabel

private

Definition at line 62 of file MCPizeroAnalyzer.h.

std::string MCPizeroAnalyzer::SimVtxLabel

private

Definition at line 63 of file MCPizeroAnalyzer.h.

PizeroMCTruthFinder* MCPizeroAnalyzer::thePizeroMCTruthFinder_

private

Definition at line 44 of file MCPizeroAnalyzer.h.

Referenced by analyze(), beginJob(), and ~MCPizeroAnalyzer().