#include <CheckSecondary.h>

Inheritance diagram for CheckSecondary:

Public Member Functions
	CheckSecondary (const edm::ParameterSet &p)

	~CheckSecondary () override

Public Member Functions inherited from SimWatcher
	SimWatcher ()

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)

	CheckSecondary (const CheckSecondary &)=delete

void	endTree ()

const CheckSecondary &	operator= (const CheckSecondary &)=delete

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

Detailed Description

Definition at line 25 of file CheckSecondary.h.

Constructor & Destructor Documentation

CheckSecondary::CheckSecondary ( const edm::ParameterSet & p )

Definition at line 25 of file CheckSecondary.cc.

References bookTree(), deltae, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), mass, nsec, procids, procs, px, py, pz, 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::~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;
 }

CheckSecondary::CheckSecondary ( const CheckSecondary & )

privatedelete

Member Function Documentation

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

private

Definition at line 91 of file CheckSecondary.cc.

References deltae, file, mass, nsec, procids, procs, px, py, and pz.

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;
 }

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;
 }

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

privatedelete

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 LogDebug.

                                                    {
   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();
 }

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.

                                                    {
   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;
 }

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, mps_fire::i, createfilelist::int, LogDebug, funct::m, mag(), dataset::name, nHad, storeIt, AlCaHLTBitMon_QueryRunRegistry::string, l1t::tracks, TreatSecondary::tracks(), and treatSecondary.

                                                {
   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);
       }
     }
   }
 }

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.

                                                  {
   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();
 }