82 double perp = std::sqrt(HitPoint.
y()*HitPoint.
y()+HitPoint.
z()*HitPoint.
z());
83 int radialPosition =
std::min(69,
int(perp/
cm));
92 G4double totalFilledProfileHcal = 0.0;
94 static G4int IDlateralProfile = -1;
95 if (IDlateralProfile < 0)
96 IDlateralProfile = man->
GetH1Id(
"h500");
98 for (
G4int i=0; i<70; i++) {
107 G4cout <<
"CCalAnalysis::InsertLateralProfile: Total filled Profile Hcal"
108 <<
" histo " << totalFilledProfileHcal <<
G4endl;
111 static G4int IDTimeHist = -1;
113 IDTimeHist = man->
GetH1Id(
"h300");
114 G4double totalFilledTimeProfile = 0.0;
123 static G4int IDTimeProfile = -1;
124 if (IDTimeProfile < 0)
125 IDTimeProfile = man->
GetH1Id(
"h901");
129 G4cout <<
"Fill Time profile histo 1 with " << t <<
" " <<
x <<
G4endl;
133 G4cout <<
"CCalAnalysis::InsertTime: Total filled Time profile histo "
134 << totalFilledTimeProfile <<
G4endl;
static G4RunManager * GetRunManager()
const G4UserRunAction * GetUserRunAction() const
G4StepPoint * GetPreStepPoint() const
G4int GetH1Id(const G4String &name, G4bool warn=true) const
const G4ThreeVector & GetPosition() const
G4double GetTotalEnergyDeposit() const
G4double GetGlobalTime() const
G4StepPoint * GetPostStepPoint() const
static constexpr double nanosecond
static constexpr double cm
G4GLOB_DLL std::ostream G4cout
const CCalRunAction * runAct
virtual void UserSteppingAction(const G4Step *aStep)
G4bool FillH1(G4int id, G4double value, G4double weight=1.0)
static constexpr double GeV
T min(const T t1, const T t2)
brief Return the smallest of the two arguments