d9/daa/GflashShowino_8h_source.html

#ifndef GflashShowino_H

#define GflashShowino_H


#include "SimGeneral/GFlash/interface/Gflash3Vector.h"

#include "SimGeneral/GFlash/interface/GflashNameSpace.h"

#include "SimGeneral/GFlash/interface/GflashTrajectory.h"


class GflashShowino {

public:

  //-------------------------

  // Constructor, destructor

  //-------------------------

  GflashShowino();

  ~GflashShowino();


  void initialize(int showerType,

                  double energy,

                  double globalTime,

                  double charge,

                  Gflash3Vector &position,

                  Gflash3Vector &momentum,

                  double magneticField);


  void updateShowino(double deltaStep);


  int getShowerType() { return theShowerType; }

  double getEnergy() { return theEnergy; }

  double getPathLengthOnEcal() { return thePathLengthOnEcal; }

  double getPathLengthAtShower() { return thePathLengthAtShower; }

  Gflash3Vector &getPositionAtShower() { return thePositionAtShower; }

  double getStepLengthToHcal() { return theStepLengthToHcal; }

  double getStepLengthToOut() { return theStepLengthToOut; }

  GflashTrajectory *getHelix() { return theHelix; }


  double getGlobalTime() { return theGlobalTime; }

  double getPathLength() { return thePathLength; }

  Gflash3Vector &getPosition() { return thePosition; }

  double getEnergyDeposited() { return theEnergyDeposited; }

  double getDepth() { return (thePathLength - thePathLengthAtShower); }


  void setGlobalTime(double globalTime) { theGlobalTime = globalTime; }

  void setPathLength(double pathLength) { thePathLength = pathLength; }

  void setPosition(const Gflash3Vector &position) { thePosition = position; }

  void addEnergyDeposited(double energy) { theEnergyDeposited += energy; }


private:

  Gflash3Vector &simulateFirstInteractionPoint(int showType, Gflash3Vector &pos);

  int convertShowerType(int fastSimShowerType, const Gflash3Vector &pos);

  void evaluateLengths();


private:

  // fixed at the shower starting point

  int theShowerType;

  double theEnergy;

  Gflash3Vector thePositionAtShower;

  double thePathLengthAtShower;

  double thePathLengthOnEcal;

  double theStepLengthToHcal;

  double theStepLengthToOut;


  // updated along the showino trajectory line

  double thePathLength;

  double theGlobalTime;

  Gflash3Vector thePosition;

  double theEnergyDeposited;


  GflashTrajectory *theHelix;

  GflashTrajectoryPoint theTrajectoryPoint;

};


#endif