CheckSecondary Class Reference

#include <CheckSecondary.h>

Inheritance diagram for CheckSecondary:

Public Member Functions
	CheckSecondary (const edm::ParameterSet &p)
	CheckSecondary (const CheckSecondary &)=delete
const CheckSecondary &	operator= (const CheckSecondary &)=delete
	~CheckSecondary () override
Public Member Functions inherited from SimWatcher
virtual void	beginRun (edm::EventSetup const &)
bool	isMT () const
const SimWatcher &	operator= (const SimWatcher &)=delete
virtual void	registerConsumes (edm::ConsumesCollector)
	SimWatcher ()
	SimWatcher (const SimWatcher &)=delete
virtual	~SimWatcher ()
Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()
void	slotForUpdate (const BeginOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfTrack *>
	Observer ()
void	slotForUpdate (const BeginOfTrack * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const EndOfEvent *>
	Observer ()
void	slotForUpdate (const EndOfEvent * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const G4Step *>
	Observer ()
void	slotForUpdate (const G4Step * iT)
virtual	~Observer ()

Private Member Functions
TTree *	bookTree (std::string)
void	endTree ()
void	update (const BeginOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfTrack *trk) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const G4Step *step) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const EndOfEvent *evt) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
std::vector< double > *	deltae
TFile *	file
std::vector< double > *	mass
int	nHad
std::vector< int > *	nsec
std::vector< int > *	procids
std::vector< std::string > *	procs
std::vector< double > *	px
std::vector< double > *	py
std::vector< double > *	pz
bool	saveToTree
bool	storeIt
TreatSecondary *	treatSecondary
TTree *	tree
G4ProcessTypeEnumerator *	typeEnumerator

Additional Inherited Members
Protected Member Functions inherited from SimWatcher
void	setMT (bool val)

Detailed Description

Definition at line 25 of file CheckSecondary.h.

Constructor & Destructor Documentation

CheckSecondary() [1/2]

CheckSecondary::CheckSecondary

(

const edm::ParameterSet &

p

)

Definition at line 25 of file CheckSecondary.cc.

References bookTree(), deltae, edm::ParameterSet::getUntrackedParameter(), mass, nsec, AlCaHLTBitMon_ParallelJobs::p, procids, procs, px, py, pz, fed_dqm_sourceclient-live_cfg::saveFile, saveToTree, AlCaHLTBitMon_QueryRunRegistry::string, treatSecondary, and typeEnumerator.

    : treatSecondary(nullptr),
      typeEnumerator(nullptr),
      nsec(nullptr),
      procids(nullptr),
      px(nullptr),
      py(nullptr),
      pz(nullptr),
      mass(nullptr),
      deltae(nullptr),
      procs(nullptr),
      file(nullptr),
      tree(nullptr) {
  edm::ParameterSet m_p = p.getParameter<edm::ParameterSet>("CheckSecondary");
  std::string saveFile = m_p.getUntrackedParameter<std::string>("SaveInFile", "None");
  treatSecondary = new TreatSecondary(m_p);
  typeEnumerator = new G4ProcessTypeEnumerator();

  nsec = new std::vector<int>();
  px = new std::vector<double>();
  py = new std::vector<double>();
  pz = new std::vector<double>();
  mass = new std::vector<double>();
  deltae = new std::vector<double>();
  procids = new std::vector<int>();
  procs = new std::vector<std::string>();

  if (saveFile != "None") {
    saveToTree = true;
    tree = bookTree(saveFile);
    edm::LogInfo("CheckSecondary") << "Instantiate CheckSecondary with first"
                                   << " hadronic interaction information"
                                   << " to be saved in file " << saveFile;
  } else {
    saveToTree = false;
    edm::LogInfo("CheckSecondary") << "Instantiate CheckSecondary with first"
                                   << " hadronic interaction information"
                                   << " not saved";
  }
}

CheckSecondary() [2/2]

CheckSecondary::CheckSecondary ( const CheckSecondary &  )

delete

~CheckSecondary()

CheckSecondary::~CheckSecondary ( )

override

Definition at line 66 of file CheckSecondary.cc.

References deltae, endTree(), mass, nsec, procids, procs, px, py, pz, saveToTree, treatSecondary, and typeEnumerator.

                                {
  if (saveToTree)
    endTree();
  if (nsec)
    delete nsec;
  if (px)
    delete px;
  if (py)
    delete py;
  if (pz)
    delete pz;
  if (mass)
    delete mass;
  if (deltae)
    delete deltae;
  if (procs)
    delete procs;
  if (procids)
    delete procids;
  if (typeEnumerator)
    delete typeEnumerator;
  if (treatSecondary)
    delete treatSecondary;
}

Member Function Documentation

bookTree()

TTree * CheckSecondary::bookTree ( std::string  fileName )

private

Definition at line 91 of file CheckSecondary.cc.

References deltae, file, MillePedeFileConverter_cfg::fileName, mass, nsec, procids, procs, px, py, pz, and RandomServiceHelper::t1.

Referenced by CheckSecondary().

                                                  {
  file = new TFile(fileName.c_str(), "RECREATE");
  file->cd();

  TTree *t1 = new TTree("T1", "Secondary Particle Information");
  t1->Branch("SecondaryPx", "std::vector<double>", &px);
  t1->Branch("SecondaryPy", "std::vector<double>", &py);
  t1->Branch("SecondaryPz", "std::vector<double>", &pz);
  t1->Branch("SecondaryMass", "std::vector<double>", &mass);
  t1->Branch("NumberSecondaries", "std::vector<int>", &nsec);
  t1->Branch("DeltaEinInteract", "std::vector<double>", &deltae);
  t1->Branch("ProcessID", "std::vector<int>", &procids);
  t1->Branch("ProcessNames", "std::vector<std::string>", &procs);
  return t1;
}

endTree()

void CheckSecondary::endTree ( )

private

Definition at line 107 of file CheckSecondary.cc.

References file.

Referenced by ~CheckSecondary().

                             {
  edm::LogInfo("CheckSecondary") << "Save the Secondary Tree " << tree->GetName() << " (" << tree << ") in file "
                                 << file->GetName() << " (" << file << ")";
  file->cd();
  tree->Write();
  file->Close();
  delete file;
}

operator=()

const CheckSecondary& CheckSecondary::operator= ( const CheckSecondary &  )

delete

update() [1/4]

void CheckSecondary::update ( const BeginOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 116 of file CheckSecondary.cc.

References caHitNtupletGeneratorKernels::iev, and LogDebug.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                   {
  int iev = (*evt)()->GetEventID();
  LogDebug("CheckSecondary") << "CheckSecondary::=====> Begin of event = " << iev;

  (*nsec).clear();
  (*procs).clear();
  (*procids).clear();
  (*deltae).clear();
  (*px).clear();
  (*py).clear();
  (*pz).clear();
  (*mass).clear();
}

update() [2/4]

void CheckSecondary::update ( const BeginOfTrack *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfTrack *>.

Definition at line 130 of file CheckSecondary.cc.

References TreatSecondary::initTrack(), LogDebug, nHad, storeIt, and treatSecondary.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                   {
  const G4Track *thTk = (*trk)();
  treatSecondary->initTrack(thTk);
  if (thTk->GetParentID() <= 0)
    storeIt = true;
  else
    storeIt = false;
  nHad = 0;
  LogDebug("CheckSecondary") << "CheckSecondary:: Track " << thTk->GetTrackID() << " Parent " << thTk->GetParentID()
                             << " Flag " << storeIt;
}

update() [3/4]

void CheckSecondary::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 142 of file CheckSecondary.cc.

References ALCARECOTkAlJpsiMuMu_cff::charge, EgHLTOffHistBins_cfi::deltaE, mps_fire::i, createfilelist::int, LogDebug, visualization-live-secondInstance_cfg::m, mag(), Skims_PA_cff::name, nHad, fileinputsource_cfi::sec, storeIt, AlCaHLTBitMon_QueryRunRegistry::string, TreatSecondary::tracks(), DiMuonV_cfg::tracks, and treatSecondary.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                               {
  std::string name;
  int procID;
  bool hadrInt;
  double deltaE;
  std::vector<int> charge;
  std::vector<math::XYZTLorentzVector> tracks = treatSecondary->tracks(aStep, name, procID, hadrInt, deltaE, charge);
  if (storeIt && hadrInt) {
    double pInit = (aStep->GetPreStepPoint()->GetMomentum()).mag();
    double pEnd = (aStep->GetPostStepPoint()->GetMomentum()).mag();
    nHad++;
    int sec = (int)(tracks.size());
    LogDebug("CheckSecondary") << "CheckSecondary:: Hadronic Interaction " << nHad << " of type " << name << " ID "
                               << procID << " with " << sec << " secondaries "
                               << " and Momentum (MeV/c) at start " << pInit << " and at end " << pEnd;
    (*nsec).push_back(sec);
    (*procs).push_back(name);
    (*procids).push_back(procID);
    (*deltae).push_back(deltaE);
    if (nHad == 1) {
      for (int i = 0; i < sec; i++) {
        double m = tracks[i].M();
        if (charge[i] < 0)
          m = -m;
        (*px).push_back(tracks[i].Px());
        (*py).push_back(tracks[i].Py());
        (*pz).push_back(tracks[i].Pz());
        (*mass).push_back(m);
      }
    }
  }
}

update() [4/4]

void CheckSecondary::update ( const EndOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent *>.

Definition at line 175 of file CheckSecondary.cc.

References mps_fire::i, LogDebug, G4ProcessTypeEnumerator::processG4Name(), procids, saveToTree, and typeEnumerator.

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), progressbar.ProgressBar::finish(), and MatrixUtil.Steps::overwrite().

                                                 {
  LogDebug("CheckSecondary") << "CheckSecondary::EndofEvent =====> Event " << (*evt)()->GetEventID() << " with "
                             << (*nsec).size() << " hadronic collisions with"
                             << " secondaries produced in each step";
  for (unsigned int i = 0; i < (*nsec).size(); i++)
    LogDebug("CheckSecondary") << " " << (*nsec)[i] << " from " << (*procs)[i] << " ID " << (*procids)[i] << " ("
                               << typeEnumerator->processG4Name((*procids)[i]) << ") deltaE = " << (*deltae)[i]
                               << " MeV";
  LogDebug("CheckSecondary") << "And " << (*mass).size() << " secondaries "
                             << "produced in the first interactions";
  for (unsigned int i = 0; i < (*mass).size(); i++)
    LogDebug("CheckSecondary") << "Secondary " << i << " (" << (*px)[i] << ", " << (*py)[i] << ", " << (*pz)[i] << ", "
                               << (*mass)[i] << ")";

  if (saveToTree)
    tree->Fill();
}

Member Data Documentation

deltae

std::vector<double> * CheckSecondary::deltae

private

Definition at line 51 of file CheckSecondary.h.

Referenced by bookTree(), CheckSecondary(), and ~CheckSecondary().

file

TFile* CheckSecondary::file

private

Definition at line 53 of file CheckSecondary.h.

Referenced by bookTree(), and endTree().

mass

std::vector<double> * CheckSecondary::mass

private

Definition at line 51 of file CheckSecondary.h.

Referenced by Particle.Particle::__str__(), ZMuMuRochCorAnalyzer.DiMuon::__str__(), DiObject.DiMuon::__str__(), bookTree(), CheckSecondary(), datamodel.Object::p4(), and ~CheckSecondary().

nHad

int CheckSecondary::nHad

private

Definition at line 49 of file CheckSecondary.h.

Referenced by update().

nsec

std::vector<int>* CheckSecondary::nsec

private

Definition at line 50 of file CheckSecondary.h.

Referenced by bookTree(), CheckSecondary(), and ~CheckSecondary().

procids

std::vector<int> * CheckSecondary::procids

private

Definition at line 50 of file CheckSecondary.h.

Referenced by bookTree(), CheckSecondary(), update(), and ~CheckSecondary().

procs

std::vector<std::string>* CheckSecondary::procs

private

Definition at line 52 of file CheckSecondary.h.

Referenced by bookTree(), CheckSecondary(), and ~CheckSecondary().

px

std::vector<double>* CheckSecondary::px

private

Definition at line 51 of file CheckSecondary.h.

Referenced by bookTree(), CheckSecondary(), and ~CheckSecondary().

py

std::vector<double> * CheckSecondary::py

private

Definition at line 51 of file CheckSecondary.h.

Referenced by bookTree(), CheckSecondary(), and ~CheckSecondary().

pz

std::vector<double> * CheckSecondary::pz

private

Definition at line 51 of file CheckSecondary.h.

Referenced by bookTree(), CheckSecondary(), and ~CheckSecondary().

saveToTree

bool CheckSecondary::saveToTree

private

Definition at line 48 of file CheckSecondary.h.

Referenced by CheckSecondary(), update(), and ~CheckSecondary().

storeIt

bool CheckSecondary::storeIt

private

Definition at line 48 of file CheckSecondary.h.

Referenced by update().

treatSecondary

TreatSecondary* CheckSecondary::treatSecondary

private

Definition at line 46 of file CheckSecondary.h.

Referenced by CheckSecondary(), update(), and ~CheckSecondary().

tree

TTree* CheckSecondary::tree

private

Definition at line 54 of file CheckSecondary.h.

Referenced by svgfig.SVG::__str__(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::beginLoop(), MetTreeProducer.MetTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::declareVariables(), core.AutoFillTreeProducer.AutoFillTreeProducer::fillTree(), ZJetsTreeAnalyzer.ZJetsTreeAnalyzer::process(), MetTreeProducer.MetTreeProducer::process(), and python.cmstools.EventTree::SetAlias().

typeEnumerator

G4ProcessTypeEnumerator* CheckSecondary::typeEnumerator

private

Definition at line 47 of file CheckSecondary.h.

Referenced by CheckSecondary(), update(), and ~CheckSecondary().