ZdcTestAnalysis Class Reference

Inheritance diagram for ZdcTestAnalysis:

Public Member Functions
	ZdcTestAnalysis (const edm::ParameterSet &p)
	~ZdcTestAnalysis () 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 BeginOfJob *>
	Observer ()
void	slotForUpdate (const BeginOfJob * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfRun *>
	Observer ()
void	slotForUpdate (const BeginOfRun * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const EndOfRun *>
	Observer ()
void	slotForUpdate (const EndOfRun * iT)
virtual	~Observer ()
Public Member Functions inherited from Observer< const BeginOfEvent *>
	Observer ()
void	slotForUpdate (const BeginOfEvent * 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
void	finish ()
void	update (const BeginOfJob *run) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const EndOfRun *run) override
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfEvent *evt) 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...
void	update (const G4Step *step) override
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
int	doNTzdcevent
int	doNTzdcstep
int	eventIndex
std::string	eventNtFileName
int	stepIndex
std::string	stepNtFileName
ZdcNumberingScheme *	theZdcNumScheme
int	verbosity
Float_t	zdceventarray [16]
TNtuple *	zdceventntuple
TFile *	zdcOutputEventFile
TFile *	zdcOutputStepFile
Float_t	zdcsteparray [18]
TNtuple *	zdcstepntuple

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

Detailed Description

Definition at line 57 of file ZdcTestAnalysis.cc.

Constructor & Destructor Documentation

ZdcTestAnalysis()

ZdcTestAnalysis::ZdcTestAnalysis

(

const edm::ParameterSet &

p

)

Definition at line 140 of file ZdcTestAnalysis.cc.

References doNTzdcevent, doNTzdcstep, eventNtFileName, edm::ParameterSet::getParameter(), AlCaHLTBitMon_ParallelJobs::p, stepNtFileName, AlCaHLTBitMon_QueryRunRegistry::string, theZdcNumScheme, verbosity, zdceventntuple, and zdcstepntuple.

                                                         {
  //constructor
  edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("ZdcTestAnalysis");
  verbosity = m_Anal.getParameter<int>("Verbosity");
  doNTzdcstep = m_Anal.getParameter<int>("StepNtupleFlag");
  doNTzdcevent = m_Anal.getParameter<int>("EventNtupleFlag");
  stepNtFileName = m_Anal.getParameter<std::string>("StepNtupleFileName");
  eventNtFileName = m_Anal.getParameter<std::string>("EventNtupleFileName");

  if (verbosity > 0)
    edm::LogVerbatim("ZdcAnalysis") << "\n============================================================================";
  edm::LogVerbatim("ZdcAnalysis") << "ZdcTestAnalysis:: Initialized as observer";
  if (doNTzdcstep > 0) {
    edm::LogVerbatim("ZdcAnalysis") << " Step Ntuple will be created";
    edm::LogVerbatim("ZdcAnalysis") << " Step Ntuple file: " << stepNtFileName;
  }
  if (doNTzdcevent > 0) {
    edm::LogVerbatim("ZdcAnalysis") << " Event Ntuple will be created";
    edm::LogVerbatim("ZdcAnalysis") << " Step Ntuple file: " << stepNtFileName;
  }
  edm::LogVerbatim("ZdcAnalysis") << "============================================================================"
                                  << std::endl;

  if (doNTzdcstep > 0)
    zdcstepntuple =
        new TNtuple("NTzdcstep",
                    "NTzdcstep",
                    "evt:trackid:charge:pdgcode:x:y:z:stepl:stepe:eta:phi:vpx:vpy:vpz:idx:idl:pvtype:ncherphot");

  if (doNTzdcevent > 0)
    zdceventntuple =
        new TNtuple("NTzdcevent",
                    "NTzdcevent",
                    "evt:ihit:fiberid:zside:subdet:layer:fiber:channel:enem:enhad:hitenergy:x:y:z:time:etot");

  theZdcNumScheme = nullptr;
  //theZdcSD = new ZdcSD("ZDCHITSB", new ZdcNumberingScheme());
}

~ZdcTestAnalysis()

ZdcTestAnalysis::~ZdcTestAnalysis ( )

override

Definition at line 179 of file ZdcTestAnalysis.cc.

References finish(), and theZdcNumScheme.

                                  {
  // destructor
  finish();
  delete theZdcNumScheme;
}

Member Function Documentation

finish()

void ZdcTestAnalysis::finish ( )

private

Definition at line 484 of file ZdcTestAnalysis.cc.

References doNTzdcevent, doNTzdcstep, eventIndex, zdceventntuple, zdcOutputEventFile, zdcOutputStepFile, and zdcstepntuple.

Referenced by progressbar.ProgressBar::__next__(), and ~ZdcTestAnalysis().

                             {
  if (doNTzdcstep) {
    zdcOutputStepFile->cd();
    zdcstepntuple->Write();
    edm::LogVerbatim("ZdcAnalysis") << "ZdcTestAnalysis: Ntuple step  written for event: " << eventIndex;
    zdcOutputStepFile->Close();
    edm::LogVerbatim("ZdcAnalysis") << "ZdcTestAnalysis: Step file closed";
  }

  if (doNTzdcevent) {
    zdcOutputEventFile->cd();
    zdceventntuple->Write("", TObject::kOverwrite);
    edm::LogVerbatim("ZdcAnalysis") << "ZdcTestAnalysis: Ntuple event written for event: " << eventIndex;
    zdcOutputEventFile->Close();
    edm::LogVerbatim("ZdcAnalysis") << "ZdcTestAnalysis: Event file closed";
  }
}

update() [1/6]

void ZdcTestAnalysis::update ( const BeginOfJob *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob *>.

Definition at line 185 of file ZdcTestAnalysis.cc.

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

                                                  {
  //job
  edm::LogVerbatim("ZdcAnalysis") << "beggining of job";
}

update() [2/6]

void ZdcTestAnalysis::update ( const BeginOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun *>.

Definition at line 191 of file ZdcTestAnalysis.cc.

References doNTzdcevent, doNTzdcstep, eventIndex, eventNtFileName, stepNtFileName, zdcOutputEventFile, and zdcOutputStepFile.

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

                                                  {
  //run

  edm::LogVerbatim("ZdcAnalysis") << "\nZdcTestAnalysis: Begining of Run";
  if (doNTzdcstep) {
    edm::LogVerbatim("ZdcAnalysis") << "ZDCTestAnalysis: output step file created";
    TString stepfilename = stepNtFileName;
    zdcOutputStepFile = new TFile(stepfilename, "RECREATE");
  }

  if (doNTzdcevent) {
    edm::LogVerbatim("ZdcAnalysis") << "ZDCTestAnalysis: output event file created";
    TString stepfilename = eventNtFileName;
    zdcOutputEventFile = new TFile(stepfilename, "RECREATE");
  }

  eventIndex = 0;
}

update() [3/6]

void ZdcTestAnalysis::update ( const EndOfRun *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfRun *>.

Definition at line 482 of file ZdcTestAnalysis.cc.

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

{}

update() [4/6]

void ZdcTestAnalysis::update ( const BeginOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent *>.

Definition at line 210 of file ZdcTestAnalysis.cc.

References eventIndex, and stepIndex.

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

                                                    {
  //event
  edm::LogVerbatim("ZdcAnalysis") << "ZdcTest: Processing Event Number: " << eventIndex;
  eventIndex++;
  stepIndex = 0;
}

update() [5/6]

void ZdcTestAnalysis::update ( const EndOfEvent *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent *>.

Definition at line 338 of file ZdcTestAnalysis.cc.

References TauDecayModes::dec, doNTzdcevent, eventIndex, dqmMemoryStats::float, CaloG4Hit::getEM(), CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getHadr(), CaloG4Hit::getPosition(), CaloG4Hit::getTimeSliceID(), CaloG4Hit::getTrackID(), CaloG4Hit::getUnitID(), mps_fire::i, createfilelist::int, genParticles_cff::map, HLT_2022v15_cff::maxTime, npart, ntzdce_channel, ntzdce_enem, ntzdce_enhad, ntzdce_etot, ntzdce_evt, ntzdce_fiber, ntzdce_fiberid, ntzdce_hitenergy, ntzdce_ihit, ntzdce_layer, ntzdce_subdet, ntzdce_time, ntzdce_x, ntzdce_y, ntzdce_z, ntzdce_zside, funct::pow(), multPhiCorr_741_25nsDY_cfi::px, multPhiCorr_741_25nsDY_cfi::py, mathSSE::sqrt(), stepIndex, theZdcNumScheme, protons_cff::time, ZdcNumberingScheme::unpackZdcIndex(), verbosity, zdceventarray, and zdceventntuple.

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

                                                  {
  //end of event

  // Look for the Hit Collection
  edm::LogVerbatim("ZdcAnalysis") << "\nZdcTest::upDate(const EndOfEvent * evt) - event #" << (*evt)()->GetEventID()
                                  << "\n  # of aSteps followed in event = " << stepIndex;

  // access to the G4 hit collections
  G4HCofThisEvent* allHC = (*evt)()->GetHCofThisEvent();
  edm::LogVerbatim("ZdcAnalysis") << "  accessed all HC";

  int theZDCHCid = G4SDManager::GetSDMpointer()->GetCollectionID("ZDCHITS");
  edm::LogVerbatim("ZdcAnalysis") << " - theZDCHCid = " << theZDCHCid;

  CaloG4HitCollection* theZDCHC = (CaloG4HitCollection*)allHC->GetHC(theZDCHCid);
  edm::LogVerbatim("ZdcAnalysis") << " - theZDCHC = " << theZDCHC;

  if (!theZdcNumScheme) {
    theZdcNumScheme = new ZdcNumberingScheme(1);
  }

  float ETot = 0., SEnergy = 0.;
  int maxTime = 0;
  int fiberID = 0;
  unsigned int unsignedfiberID = 0;
  std::map<int, float, std::less<int> > energyInFibers;
  std::map<int, float, std::less<int> > primaries;
  float totalEnergy = 0;
  int nentries = theZDCHC->entries();
  edm::LogVerbatim("ZdcAnalysis") << "  theZDCHC has " << nentries << " entries";

  if (doNTzdcevent) {
    if (nentries > 0) {
      for (int ihit = 0; ihit < nentries; ihit++) {
        CaloG4Hit* caloHit = (*theZDCHC)[ihit];
        totalEnergy += caloHit->getEnergyDeposit();
      }

      for (int ihit = 0; ihit < nentries; ihit++) {
        CaloG4Hit* aHit = (*theZDCHC)[ihit];
        fiberID = aHit->getUnitID();
        unsignedfiberID = aHit->getUnitID();
        double enEm = aHit->getEM();
        double enHad = aHit->getHadr();
        math::XYZPoint hitPoint = aHit->getPosition();
        double hitEnergy = aHit->getEnergyDeposit();
        if (verbosity >= 1)
          edm::LogVerbatim("ZdcAnalysis")
              << " entry #" << ihit << ": fiberID=0x" << std::hex << fiberID << std::dec << "; enEm=" << enEm
              << "; enHad=" << enHad << "; hitEnergy=" << hitEnergy << "z=" << hitPoint.z();
        energyInFibers[fiberID] += enEm + enHad;
        primaries[aHit->getTrackID()] += enEm + enHad;
        float time = aHit->getTimeSliceID();
        if (time > maxTime)
          maxTime = (int)time;

        int thesubdet, thelayer, thefiber, thechannel, thez;
        theZdcNumScheme->unpackZdcIndex(fiberID, thesubdet, thelayer, thefiber, thechannel, thez);
        int unsignedsubdet, unsignedlayer, unsignedfiber, unsignedchannel, unsignedz;
        theZdcNumScheme->unpackZdcIndex(
            unsignedfiberID, unsignedsubdet, unsignedlayer, unsignedfiber, unsignedchannel, unsignedz);

        // unsigned int packidx1 = packZdcIndex(thesubdet, thelayer, thefiber, thechannel, thez);
        // unsigned int packidx1 = packZdcIndex(thesubdet, thelayer, thefiber, thechannel, thez);
        // unsigned int packidx1 = packZdcIndex(thesubdet, thelayer, thefiber, thechannel, thez);
        // unsigned int packidx1 = packZdcIndex(thesubdet, thelayer, thefiber, thechannel, thez);

        zdceventarray[ntzdce_evt] = (float)eventIndex;
        zdceventarray[ntzdce_ihit] = (float)ihit;
        zdceventarray[ntzdce_fiberid] = (float)fiberID;
        zdceventarray[ntzdce_zside] = (float)thez;
        zdceventarray[ntzdce_subdet] = (float)thesubdet;
        zdceventarray[ntzdce_layer] = (float)thelayer;
        zdceventarray[ntzdce_fiber] = (float)thefiber;
        zdceventarray[ntzdce_channel] = (float)thechannel;
        zdceventarray[ntzdce_enem] = enEm;
        zdceventarray[ntzdce_enhad] = enHad;
        zdceventarray[ntzdce_hitenergy] = hitEnergy;
        zdceventarray[ntzdce_x] = hitPoint.x();
        zdceventarray[ntzdce_y] = hitPoint.y();
        zdceventarray[ntzdce_z] = hitPoint.z();
        zdceventarray[ntzdce_time] = time;
        zdceventarray[ntzdce_etot] = totalEnergy;
        zdceventntuple->Fill(zdceventarray);
      }

      for (std::map<int, float, std::less<int> >::iterator is = energyInFibers.begin(); is != energyInFibers.end();
           is++) {
        ETot = (*is).second;
        SEnergy += ETot;
      }

      // Find Primary info:
      int trackID = 0;
      G4PrimaryParticle* thePrim = nullptr;
      G4int nvertex = (*evt)()->GetNumberOfPrimaryVertex();
      edm::LogVerbatim("ZdcAnalysis") << "Event has " << nvertex << " vertex";
      if (nvertex == 0)
        edm::LogVerbatim("ZdcAnalysis") << "ZdcTest End Of Event  ERROR: no vertex";

      for (int i = 0; i < nvertex; i++) {
        G4PrimaryVertex* avertex = (*evt)()->GetPrimaryVertex(i);
        if (avertex == nullptr) {
          edm::LogVerbatim("ZdcAnalysis") << "ZdcTest End Of Event ERR: pointer to vertex = 0";
        } else {
          edm::LogVerbatim("ZdcAnalysis") << "Vertex number :" << i;
          int npart = avertex->GetNumberOfParticle();
          if (npart == 0)
            edm::LogVerbatim("ZdcAnalysis") << "ZdcTest End Of Event ERR: no primary!";
          if (thePrim == nullptr)
            thePrim = avertex->GetPrimary(trackID);
        }
      }

      double px = 0., py = 0., pz = 0.;
      double pInit = 0.;

      if (thePrim != nullptr) {
        px = thePrim->GetPx();
        py = thePrim->GetPy();
        pz = thePrim->GetPz();
        pInit = sqrt(pow(px, 2.) + pow(py, 2.) + pow(pz, 2.));
        if (pInit == 0) {
          edm::LogVerbatim("ZdcAnalysis") << "ZdcTest End Of Event  ERR: primary has p=0 ";
        }
      } else {
        edm::LogVerbatim("ZdcAnalysis") << "ZdcTest End Of Event ERR: could not find primary ";
      }

    }  // nentries > 0

  }  // createNTzdcevent

  int iEvt = (*evt)()->GetEventID();
  if (iEvt < 10)
    edm::LogVerbatim("ZdcAnalysis") << " ZdcTest Event " << iEvt;
  else if ((iEvt < 100) && (iEvt % 10 == 0))
    edm::LogVerbatim("ZdcAnalysis") << " ZdcTest Event " << iEvt;
  else if ((iEvt < 1000) && (iEvt % 100 == 0))
    edm::LogVerbatim("ZdcAnalysis") << " ZdcTest Event " << iEvt;
  else if ((iEvt < 10000) && (iEvt % 1000 == 0))
    edm::LogVerbatim("ZdcAnalysis") << " ZdcTest Event " << iEvt;
}

update() [6/6]

void ZdcTestAnalysis::update ( const G4Step *  )

overrideprivatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step *>.

Definition at line 217 of file ZdcTestAnalysis.cc.

References doNTzdcstep, PVValHelper::eta, eventIndex, dqmMemoryStats::float, eventwithhistoryproducer_cfi::historyDepth, heavyIonCSV_trainingSettings::idx, findQualityFiles::jj, dqm-mbProfile::log, ntzdcs_charge, ntzdcs_eta, ntzdcs_evt, ntzdcs_idl, ntzdcs_idx, ntzdcs_ncherphot, ntzdcs_pdgcode, ntzdcs_phi, ntzdcs_pvtype, ntzdcs_stepe, ntzdcs_stepl, ntzdcs_trackid, ntzdcs_vpx, ntzdcs_vpy, ntzdcs_vpz, ntzdcs_x, ntzdcs_y, ntzdcs_z, PbPb_ZMuSkimMuonDPG_cff::particleType, phi, stepIndex, AlignmentTrackSelector_cfi::theCharge, verbosity, zdcsteparray, and zdcstepntuple.

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

                                                {
  //step;
  stepIndex++;

  if (doNTzdcstep) {
    G4StepPoint* preStepPoint = aStep->GetPreStepPoint();
    // G4StepPoint * postStepPoint= aStep->GetPostStepPoint();
    G4double stepL = aStep->GetStepLength();
    G4double stepE = aStep->GetTotalEnergyDeposit();

    if (verbosity >= 2)
      edm::LogVerbatim("ZdcAnalysis") << "Step " << stepL << "," << stepE;

    G4Track* theTrack = aStep->GetTrack();
    G4int theTrackID = theTrack->GetTrackID();
    G4double theCharge = theTrack->GetDynamicParticle()->GetCharge();
    G4String particleType = theTrack->GetDefinition()->GetParticleName();
    G4int pdgcode = theTrack->GetDefinition()->GetPDGEncoding();

    const G4ThreeVector& vert_mom = theTrack->GetVertexMomentumDirection();
    G4double vpx = vert_mom.x();
    G4double vpy = vert_mom.y();
    G4double vpz = vert_mom.z();
    double eta = 0.5 * log((1. + vpz) / (1. - vpz));
    double phi = atan2(vpy, vpx);

    const G4ThreeVector& hitPoint = preStepPoint->GetPosition();
    G4ThreeVector localPoint = theTrack->GetTouchable()->GetHistory()->GetTopTransform().TransformPoint(hitPoint);

    const G4VTouchable* touch = aStep->GetPreStepPoint()->GetTouchable();
    int idx = touch->GetReplicaNumber(0);
    int idLayer = -1;
    int thePVtype = -1;

    int historyDepth = touch->GetHistoryDepth();

    if (historyDepth > 0) {
      std::vector<int> theReplicaNumbers(historyDepth);
      std::vector<G4VPhysicalVolume*> thePhysicalVolumes(historyDepth);
      std::vector<G4String> thePVnames(historyDepth);
      std::vector<G4LogicalVolume*> theLogicalVolumes(historyDepth);
      std::vector<G4String> theLVnames(historyDepth);
      std::vector<G4Material*> theMaterials(historyDepth);
      std::vector<G4String> theMaterialNames(historyDepth);

      for (int jj = 0; jj < historyDepth; jj++) {
        theReplicaNumbers[jj] = touch->GetReplicaNumber(jj);
        thePhysicalVolumes[jj] = touch->GetVolume(jj);
        thePVnames[jj] = thePhysicalVolumes[jj]->GetName();
        theLogicalVolumes[jj] = thePhysicalVolumes[jj]->GetLogicalVolume();
        theLVnames[jj] = theLogicalVolumes[jj]->GetName();
        theMaterials[jj] = theLogicalVolumes[jj]->GetMaterial();
        theMaterialNames[jj] = theMaterials[jj]->GetName();
        if (verbosity >= 2)
          edm::LogVerbatim("ZdcAnalysis") << " GHD " << jj << ": " << theReplicaNumbers[jj] << "," << thePVnames[jj]
                                          << "," << theLVnames[jj] << "," << theMaterialNames[jj];
      }

      idLayer = theReplicaNumbers[1];
      if (thePVnames[0] == "ZDC_EMLayer")
        thePVtype = 1;
      else if (thePVnames[0] == "ZDC_EMAbsorber")
        thePVtype = 2;
      else if (thePVnames[0] == "ZDC_EMFiber")
        thePVtype = 3;
      else if (thePVnames[0] == "ZDC_HadLayer")
        thePVtype = 7;
      else if (thePVnames[0] == "ZDC_HadAbsorber")
        thePVtype = 8;
      else if (thePVnames[0] == "ZDC_HadFiber")
        thePVtype = 9;
      else if (thePVnames[0] == "ZDC_PhobosLayer")
        thePVtype = 11;
      else if (thePVnames[0] == "ZDC_PhobosAbsorber")
        thePVtype = 12;
      else if (thePVnames[0] == "ZDC_PhobosFiber")
        thePVtype = 13;
      else {
        thePVtype = 0;
        if (verbosity >= 2)
          edm::LogVerbatim("ZdcAnalysis") << " pvtype=0 hd=" << historyDepth << " " << theReplicaNumbers[0] << ","
                                          << thePVnames[0] << "," << theLVnames[0] << "," << theMaterialNames[0];
      }
    } else if (historyDepth == 0) {
      int theReplicaNumber = touch->GetReplicaNumber(0);
      G4VPhysicalVolume* thePhysicalVolume = touch->GetVolume(0);
      const G4String& thePVname = thePhysicalVolume->GetName();
      G4LogicalVolume* theLogicalVolume = thePhysicalVolume->GetLogicalVolume();
      const G4String& theLVname = theLogicalVolume->GetName();
      G4Material* theMaterial = theLogicalVolume->GetMaterial();
      const G4String& theMaterialName = theMaterial->GetName();
      if (verbosity >= 2)
        edm::LogVerbatim("ZdcAnalysis") << " hd=0 " << theReplicaNumber << "," << thePVname << "," << theLVname << ","
                                        << theMaterialName;
    } else {
      edm::LogVerbatim("ZdcAnalysis") << " hd<0:  hd=" << historyDepth;
    }

    double NCherPhot = -1;
    zdcsteparray[ntzdcs_evt] = (float)eventIndex;
    zdcsteparray[ntzdcs_trackid] = (float)theTrackID;
    zdcsteparray[ntzdcs_charge] = theCharge;
    zdcsteparray[ntzdcs_pdgcode] = (float)pdgcode;
    zdcsteparray[ntzdcs_x] = hitPoint.x();
    zdcsteparray[ntzdcs_y] = hitPoint.y();
    zdcsteparray[ntzdcs_z] = hitPoint.z();
    zdcsteparray[ntzdcs_stepl] = stepL;
    zdcsteparray[ntzdcs_stepe] = stepE;
    zdcsteparray[ntzdcs_eta] = eta;
    zdcsteparray[ntzdcs_phi] = phi;
    zdcsteparray[ntzdcs_vpx] = vpx;
    zdcsteparray[ntzdcs_vpy] = vpy;
    zdcsteparray[ntzdcs_vpz] = vpz;
    zdcsteparray[ntzdcs_idx] = (float)idx;
    zdcsteparray[ntzdcs_idl] = (float)idLayer;
    zdcsteparray[ntzdcs_pvtype] = thePVtype;
    zdcsteparray[ntzdcs_ncherphot] = NCherPhot;
    zdcstepntuple->Fill(zdcsteparray);
  }
}

Member Data Documentation

doNTzdcevent

int ZdcTestAnalysis::doNTzdcevent

private

Definition at line 81 of file ZdcTestAnalysis.cc.

Referenced by finish(), update(), and ZdcTestAnalysis().

doNTzdcstep

int ZdcTestAnalysis::doNTzdcstep

private

Definition at line 80 of file ZdcTestAnalysis.cc.

Referenced by finish(), update(), and ZdcTestAnalysis().

eventIndex

int ZdcTestAnalysis::eventIndex

private

Definition at line 91 of file ZdcTestAnalysis.cc.

Referenced by finish(), and update().

eventNtFileName

std::string ZdcTestAnalysis::eventNtFileName

private

Definition at line 83 of file ZdcTestAnalysis.cc.

Referenced by update(), and ZdcTestAnalysis().

stepIndex

int ZdcTestAnalysis::stepIndex

private

Definition at line 92 of file ZdcTestAnalysis.cc.

Referenced by update().

stepNtFileName

std::string ZdcTestAnalysis::stepNtFileName

private

Definition at line 82 of file ZdcTestAnalysis.cc.

Referenced by update(), and ZdcTestAnalysis().

theZdcNumScheme

ZdcNumberingScheme* ZdcTestAnalysis::theZdcNumScheme

private

Definition at line 97 of file ZdcTestAnalysis.cc.

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

verbosity

int ZdcTestAnalysis::verbosity

private

Definition at line 79 of file ZdcTestAnalysis.cc.

Referenced by update(), and ZdcTestAnalysis().

zdceventarray

Float_t ZdcTestAnalysis::zdceventarray[16]

private

Definition at line 95 of file ZdcTestAnalysis.cc.

Referenced by update().

zdceventntuple

TNtuple* ZdcTestAnalysis::zdceventntuple

private

Definition at line 89 of file ZdcTestAnalysis.cc.

Referenced by finish(), update(), and ZdcTestAnalysis().

zdcOutputEventFile

TFile* ZdcTestAnalysis::zdcOutputEventFile

private

Definition at line 85 of file ZdcTestAnalysis.cc.

Referenced by finish(), and update().

zdcOutputStepFile

TFile* ZdcTestAnalysis::zdcOutputStepFile

private

Definition at line 86 of file ZdcTestAnalysis.cc.

Referenced by finish(), and update().

zdcsteparray

Float_t ZdcTestAnalysis::zdcsteparray[18]

private

Definition at line 94 of file ZdcTestAnalysis.cc.

Referenced by update().

zdcstepntuple

TNtuple* ZdcTestAnalysis::zdcstepntuple

private

Definition at line 88 of file ZdcTestAnalysis.cc.

Referenced by finish(), update(), and ZdcTestAnalysis().