HcalTestAnalysis Class Reference

#include <HcalTestAnalysis.h>

Inheritance diagram for HcalTestAnalysis:

Public Member Functions
	HcalTestAnalysis (const edm::ParameterSet &p)
virtual	~HcalTestAnalysis ()
Public Member Functions inherited from SimWatcher
	SimWatcher ()
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 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	fill (const EndOfEvent *ev)
void	layerAnalysis ()
std::vector< int >	layerGrouping (int)
void	qieAnalysis (CLHEP::HepRandomEngine *)
double	timeOfFlight (int det, int layer, double eta)
std::vector< int >	towersToAdd (int centre, int nadd)
void	update (const BeginOfJob *run)
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfRun *run)
	This routine will be called when the appropriate signal arrives. More...
void	update (const BeginOfEvent *evt)
	This routine will be called when the appropriate signal arrives. More...
void	update (const EndOfEvent *evt)
	This routine will be called when the appropriate signal arrives. More...
void	update (const G4Step *step)
	This routine will be called when the appropriate signal arrives. More...

Private Attributes
int	addTower
std::vector< CaloHit >	caloHitCache
int	centralTower
unsigned int	count
double	edepEB
double	edepEE
double	edepHB
double	edepHE
double	edepHO
double	edepl [20]
double	eta0
std::string	fileName
std::vector< int >	group_
double	mudist [20]
HcalQie *	myqie
std::vector< std::string >	names
int	nGroup
int	nTower
HcalNumberingFromDDD *	numberingFromDDD
HcalTestNumberingScheme *	org
double	phi0
std::vector< int >	tower_
HcalTestHistoClass *	tuples
std::auto_ptr < HcalTestHistoManager >	tuplesManager

Additional Inherited Members

Detailed Description

Definition at line 34 of file HcalTestAnalysis.h.

Constructor & Destructor Documentation

HcalTestAnalysis::HcalTestAnalysis

(

const edm::ParameterSet &

p

)

Definition at line 32 of file HcalTestAnalysis.cc.

References addTower, centralTower, count, eta0, fileName, edm::ParameterSet::getParameter(), group_, i, layerGrouping(), myqie, names, nGroup, nTower, phi0, AlCaHLTBitMon_QueryRunRegistry::string, tower_, and towersToAdd().

                                                          : 
  myqie(0), addTower(3), tuplesManager(0), tuples(0), numberingFromDDD(0), org(0) {

  edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("HcalTestAnalysis");
  eta0         = m_Anal.getParameter<double>("Eta0");
  phi0         = m_Anal.getParameter<double>("Phi0");
  int laygroup = m_Anal.getParameter<int>("LayerGrouping");
  centralTower = m_Anal.getParameter<int>("CentralTower");
  names        = m_Anal.getParameter<std::vector<std::string> >("Names");
  fileName     = m_Anal.getParameter<std::string>("FileName");

  edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Initialised as observer of "
              << "begin/end events and of G4step";

  count  = 0;
  group_ = layerGrouping(laygroup);
  nGroup = 0;
  for (unsigned int i=0; i<group_.size(); i++) 
    if (group_[i]>nGroup) nGroup = group_[i];
  tower_ = towersToAdd(centralTower, addTower);
  nTower = tower_.size()/2;

  edm::LogInfo("HcalSim") << "HcalTestAnalysis:: initialised for " << nGroup 
              << " Longitudinal groups and " << nTower 
              << " towers";

  // qie
  myqie  = new HcalQie(p);
} 

HcalTestAnalysis::~HcalTestAnalysis ( )

virtual

Definition at line 62 of file HcalTestAnalysis.cc.

References count, myqie, numberingFromDDD, org, and tuples.

                                    {
  edm::LogInfo("HcalSim") << "HcalTestAnalysis: -------->  Total number of "
              << "selected entries : " << count;
  edm::LogInfo("HcalSim") << "HcalTestAnalysis: Pointers:: HcalQie " << myqie 
              << ", HistoClass " << tuples << ", Numbering Scheme "
              << org << " and FromDDD " << numberingFromDDD;
  if (myqie)  {
    edm::LogInfo("HcalSim") << "HcalTestAnalysis: Delete HcalQie";
    delete myqie;
  }
  if (numberingFromDDD) {
    edm::LogInfo("HcalSim") << "HcalTestAnalysis: Delete HcalNumberingFromDDD";
    delete numberingFromDDD;
  }
}

Member Function Documentation

void HcalTestAnalysis::fill ( const EndOfEvent *  ev )

private

Definition at line 340 of file HcalTestAnalysis.cc.

References caloHitCache, alignCSCRings::e, eta(), CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getPosition(), CaloG4Hit::getTimeSlice(), HcalTestNumberingScheme::getUnitID(), CaloG4Hit::getUnitID(), GeV, HcalEndcap, HcalForward, j, log, LogDebug, names, numberingFromDDD, org, HcalTestNumberingScheme::packHcalIndex(), phi, AlCaHLTBitMon_QueryRunRegistry::string, funct::tan(), theta(), cond::rpcobgas::time, timeOfFlight(), HcalNumberingFromDDD::unitID(), HcalTestNumberingScheme::unpackHcalIndex(), and ecaldqm::zside().

Referenced by update().

                                                  {

  LogDebug("HcalSim") << "HcalTestAnalysis: Fill event " 
              << (*evt)()->GetEventID();
  
  // access to the G4 hit collections 
  G4HCofThisEvent* allHC = (*evt)()->GetHCofThisEvent();
  
  int nhc = 0, neb = 0, nef = 0, j = 0;    
  caloHitCache.erase (caloHitCache.begin(), caloHitCache.end());
 
  // Hcal
  int HCHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names[0]);
  CaloG4HitCollection* theHCHC = (CaloG4HitCollection*) allHC->GetHC(HCHCid);
  LogDebug("HcalSim") << "HcalTestAnalysis :: Hit Collection for " << names[0] 
              << " of ID " << HCHCid << " is obtained at " << theHCHC;
  if (HCHCid >= 0 && theHCHC > 0) {
    for (j = 0; j < theHCHC->entries(); j++) {

      CaloG4Hit* aHit = (*theHCHC)[j]; 
    
      double e        = aHit->getEnergyDeposit()/GeV;
      double time     = aHit->getTimeSlice();
      
      math::XYZPoint pos   = aHit->getPosition();
      double theta    = pos.theta();
      double   eta    = -log(tan(theta * 0.5));
      double   phi    = pos.phi();
    
      uint32_t unitID = aHit->getUnitID();
      int subdet, zside, layer, etaIndex, phiIndex, lay;
      org->unpackHcalIndex(unitID,subdet,zside,layer,etaIndex,phiIndex,lay);
      double jitter   = time-timeOfFlight(subdet,lay,eta);
      if (jitter<0) jitter = 0;
      CaloHit hit(subdet,lay,e,eta,phi,jitter,unitID);
      caloHitCache.push_back(hit);
      nhc++;

      std::string det =  "HB";
      if (subdet == static_cast<int>(HcalForward)) {
    det = "HF";
      } else if (subdet == static_cast<int>(HcalEndcap)) {
    if (etaIndex <= 20) {
      det  = "HES";
    } else {
      det = "HED";
    }
      }

      LogDebug("HcalSim") << "HcalTest: " << det << "  layer " << std::setw(2) 
              << layer  << " time " << std::setw(6) << time 
              << " theta "  << std::setw(8) << theta << " eta " 
              << std::setw(8) << eta  << " phi " << std::setw(8) 
              << phi << " e " << std::setw(8) << e;
    }
  }
  LogDebug("HcalSim") << "HcalTestAnalysis::HCAL hits : " << nhc << "\n";

  // EB
  int EBHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names[1]);
  CaloG4HitCollection* theEBHC = (CaloG4HitCollection*) allHC->GetHC(EBHCid);
  LogDebug("HcalSim") << "HcalTestAnalysis :: Hit Collection for " << names[1]
              << " of ID " << EBHCid << " is obtained at " << theEBHC;
  if (EBHCid >= 0 && theEBHC > 0) {
    for (j = 0; j < theEBHC->entries(); j++) {

      CaloG4Hit* aHit = (*theEBHC)[j]; 
    
      double e        = aHit->getEnergyDeposit()/GeV;
      double time     = aHit->getTimeSlice();
      std::string det =  "EB";
      
      math::XYZPoint pos  = aHit->getPosition();
      double theta    = pos.theta();
      double   eta    = -log(tan(theta/2.));
      double   phi    = pos.phi();

      HcalNumberingFromDDD::HcalID id = numberingFromDDD->unitID(eta,phi,1,1);
      uint32_t unitID = org->getUnitID(id);
      int subdet, zside, layer, ieta, iphi, lay;
      org->unpackHcalIndex(unitID,subdet,zside,layer,ieta,iphi,lay);
      subdet = 10;
      layer  = 0;
      unitID = org->packHcalIndex(subdet,zside,layer,ieta,iphi,lay);
      CaloHit hit(subdet,lay,e,eta,phi,time,unitID);
      caloHitCache.push_back(hit);
      neb++;
      LogDebug("HcalSim") << "HcalTest: " << det << "  layer " << std::setw(2) 
              << layer << " time " << std::setw(6) << time 
              << " theta " << std::setw(8) << theta << " eta " 
              << std::setw(8) << eta << " phi " << std::setw(8) 
              << phi << " e " << std::setw(8) << e;
    }
  }
  LogDebug("HcalSim") << "HcalTestAnalysis::EB hits : " << neb << "\n";

  // EE
  int EEHCid = G4SDManager::GetSDMpointer()->GetCollectionID(names[2]);
  CaloG4HitCollection* theEEHC = (CaloG4HitCollection*) allHC->GetHC(EEHCid);
  LogDebug("HcalSim") << "HcalTestAnalysis :: Hit Collection for " << names[2]
              << " of ID " << EEHCid << " is obtained at " << theEEHC;
  if (EEHCid >= 0 && theEEHC > 0) {
    for (j = 0; j < theEEHC->entries(); j++) {

      CaloG4Hit* aHit = (*theEEHC)[j]; 
    
      double e        = aHit->getEnergyDeposit()/GeV;
      double time     = aHit->getTimeSlice();
      std::string det = "EE";
      
      math::XYZPoint pos  = aHit->getPosition();
      double theta    = pos.theta();
      double   eta    = -log(tan(theta/2.));
      double   phi    = pos.phi();

      HcalNumberingFromDDD::HcalID id = numberingFromDDD->unitID(eta,phi,1,1);
      uint32_t unitID = org->getUnitID(id);
      int subdet, zside, layer, ieta, iphi, lay;
      org->unpackHcalIndex(unitID,subdet,zside,layer,ieta,iphi,lay);
      subdet = 11;
      layer  = 0;
      unitID = org->packHcalIndex(subdet,zside,layer,ieta,iphi,lay);
      CaloHit hit(subdet,lay,e,eta,phi,time,unitID);
      caloHitCache.push_back(hit);
      nef++;
      LogDebug("HcalSim") << "HcalTest: " << det << "  layer " << std::setw(2)
              << layer << " time " << std::setw(6) << time 
              << " theta " << std::setw(8) << theta << " eta " 
              << std::setw(8) << eta  << " phi " << std::setw(8) 
              << phi << " e " << std::setw(8) << e;
    }
  }
  LogDebug("HcalSim") << "HcalTestAnalysis::EE hits : " << nef << "\n";
}

void HcalTestAnalysis::layerAnalysis ( )

private

Definition at line 596 of file HcalTestAnalysis.cc.

References edepEB, edepEE, edepHB, edepHE, edepHO, edepl, HcalTestHistoClass::fillLayers(), i, LogDebug, MeV, mudist, and tuples.

Referenced by update().

                                    {

  int i = 0;
  LogDebug("HcalSim") << "\n ===>>> HcalTestAnalysis: Energy deposit in MeV " 
              << "\n at EB : " << std::setw(6) << edepEB/MeV 
              << "\n at EE : " << std::setw(6) << edepEE/MeV 
              << "\n at HB : " << std::setw(6) << edepHB/MeV
              << "\n at HE : " << std::setw(6) << edepHE/MeV
              << "\n at HO : " << std::setw(6) << edepHO/MeV
              << "\n ---- HcalTestAnalysis: Energy deposit in Layers"; 
  for (i = 0; i < 20; i++) 
    LogDebug("HcalSim") << " Layer " << std::setw(2) << i << " E " 
            << std::setw(8) << edepl[i]/MeV  << " MeV";

  tuples->fillLayers(edepl, edepHO, edepHB+edepHE, mudist);  
}

std::vector< int > HcalTestAnalysis::layerGrouping ( int  group )

private

Definition at line 78 of file HcalTestAnalysis.cc.

References i, and groupFilesInBlocks::temp.

Referenced by HcalTestAnalysis().

                                                        {

  std::vector<int> temp(19);
  if (group <= 1) {
    int grp[19] = {1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2};
    for (int i=0; i<19; i++)
      temp[i] = grp[i];
  } else  if (group == 2) {
    int grp[19] = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19};
    for (int i=0; i<19; i++)
      temp[i] = grp[i];
  } else if (group == 3) {
    int grp[19] = {1, 1, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 5, 5, 6, 6, 6, 7, 7};
    for (int i=0; i<19; i++)
      temp[i] = grp[i];
  } else if (group == 4) {
    int grp[19] = {1, 1, 1, 2, 2, 3, 3, 4, 4, 5, 5, 5, 6, 6, 6, 6, 6, 7, 7};
    for (int i=0; i<19; i++)
      temp[i] = grp[i];
  } else {
    int grp[19] = {1, 1, 2, 3, 3, 4, 4, 5, 5, 5, 6, 6, 6, 6, 6, 6, 6, 7, 7};
    for (int i=0; i<19; i++)
      temp[i] = grp[i];
  }

  edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Layer Grouping ";
  for (int i=0; i<19; i++)
    edm::LogInfo("HcalSim") << "HcalTestAnalysis: Group[" << i << "] = "
                  << temp[i];
  return temp;
}

void HcalTestAnalysis::qieAnalysis ( CLHEP::HepRandomEngine *  engine )

private

Definition at line 476 of file HcalTestAnalysis.cc.

References addTower, caloHitCache, centralTower, CaloHit::det(), CaloHit::e(), eta0, HcalTestHistoClass::fillHits(), HcalTestHistoClass::fillQie(), HcalQie::getCode(), HcalQie::getEnergy(), HcalTestNumberingScheme::getUnitID(), watchdog::group, group_, HcalBarrel, HcalEndcap, i, CaloHit::id(), customizeTrackingMonitorSeedNumber::idx, CaloHit::layer(), LogDebug, myqie, nGroup, nt, nTower, numberingFromDDD, org, phi0, AlCaHLTBitMon_QueryRunRegistry::string, tower_, tuples, HcalNumberingFromDDD::unitID(), HcalTestNumberingScheme::unpackHcalIndex(), and ecaldqm::zside().

Referenced by update().

                                                               {

  //Fill tuple with hit information
  int hittot = caloHitCache.size();
  tuples->fillHits(caloHitCache);

  //Get the index of the central tower
  HcalNumberingFromDDD::HcalID id = numberingFromDDD->unitID(eta0,phi0,1,1);
  uint32_t unitID = org->getUnitID(id);
  int subdet, zside, layer, ieta, iphi, lay;
  org->unpackHcalIndex(unitID,subdet,zside,layer,ieta,iphi,lay);
  int      laymax = 0;
  std::string det = "Unknown";
  if (subdet == static_cast<int>(HcalBarrel)) {
    laymax = 4; det = "HB";
  } else if (subdet == static_cast<int>(HcalEndcap)) {
    laymax = 2; det = "HES";
  }
  LogDebug("HcalSim") << "HcalTestAnalysis::Qie: " << det << " Eta " << ieta 
              << " Phi " << iphi << " Laymax " << laymax << " Hits "
              << hittot;

  if (laymax>0 && hittot>0) {
    std::vector<CaloHit>  hits(hittot);
    std::vector<double> eqielay(80,0.0),  esimlay(80,0.0),  esimtot(4,0.0);
    std::vector<double> eqietow(200,0.0), esimtow(200,0.0), eqietot(4,0.0);
    int etac = (centralTower/100)%100;
    int phic = (centralTower%100);

    for (int layr=0; layr<nGroup; layr++) {
      /*
      int layx, layy=20;
      for (int i=0; i<20; i++) 
    if (group_[i] == layr+1 && i < layy) layy = i+1;
      if (subdet == static_cast<int>(HcalBarrel)) {
    if      (layy < 2)    layx = 0;
    else if (layy < 17)   layx = 1;
    else if (layy == 17)  layx = 2;
    else                  layx = 3;
      } else {
    if      (layy < 2)    layx = 0;
    else                  layx = 1;
      }
      */
      for (int it=0; it<nTower; it++) {
    int    nhit = 0;
    double esim = 0;
    for (int k1 = 0; k1 < hittot; k1++) {
      CaloHit hit = caloHitCache[k1];
      int     subdetc = hit.det();
      int     layer   = hit.layer();
      int     group   = 0;
      if (layer > 0 && layer < 20) group = group_[layer];
      if (subdetc == subdet && group == layr+1) {
        int zsidec, ietac, iphic, idx;
        unitID = hit.id();
        org->unpackHcalIndex(unitID,subdetc,zsidec,layer,ietac,iphic,lay);
        if (etac > 0 && phic > 0) {
          idx = ietac*100 + iphic;
        } else if (etac > 0) {
          idx = ietac*100;
        } else if (phic > 0) {
          idx = iphic;
        } else {
          idx = 0;
        }
        if (zsidec==zside && idx==tower_[it]) {
          hits[nhit] = hit;
          LogDebug("HcalSim") << "HcalTest: Hit " << nhit << " " << hit;
          nhit++;
          esim += hit.e();
        }
      }
    }

    std::vector<int> cd = myqie->getCode(nhit, hits, engine);
    double         eqie = myqie->getEnergy(cd);

    LogDebug("HcalSim") << "HcalTestAnalysis::Qie: Energy in layer " 
                << layr << " Sim " << esim << " After QIE " <<eqie;
    for (int i=0; i<4; i++) {
      if (tower_[nTower+it] <= i) {
        esimtot[i]         += esim;
        eqietot[i]         += eqie;
        esimlay[20*i+layr] += esim;
        eqielay[20*i+layr] += eqie;
        esimtow[50*i+it]   += esim;
        eqietow[50*i+it]   += eqie;
      }
    }
      }
    }
    LogDebug("HcalSim") << "HcalTestAnalysis::Qie: Total energy " << esimtot[3]
            << " (SimHit) " << eqietot[3] << " (After QIE)";

    std::vector<double> latphi(10);
    int nt = 2*addTower + 1;
    for (int it=0; it<nt; it++)
      latphi[it] = it-addTower;
    for (int i=0; i<4; i++) {
      double scals=1, scalq=1;
      std::vector<double> latfs(10,0.), latfq(10,0.), longs(20), longq(20);
      if (esimtot[i]>0) scals = 1./esimtot[i];
      if (eqietot[i]>0) scalq = 1./eqietot[i];
      for (int it=0; it<nTower; it++) {
    int phib = it%nt; 
    latfs[phib] += scals*esimtow[50*i+it];
    latfq[phib] += scalq*eqietow[50*i+it];
      }
      for (int layr=0; layr<=nGroup; layr++) {
    longs[layr] = scals*esimlay[20*i+layr];
    longq[layr] = scalq*eqielay[20*i+layr];
      }
      tuples->fillQie(i,esimtot[i],eqietot[i],nGroup,longs,longq,
              nt,latphi,latfs,latfq);
    }
  }
}

double HcalTestAnalysis::timeOfFlight ( int  det, int  layer, double  eta  )

private

Definition at line 615 of file HcalTestAnalysis.cc.

References funct::cos(), create_public_lumi_plots::exp, HcalBarrel, LogDebug, funct::sin(), theta(), and tmp.

Referenced by fill().

                                                                    {

  double theta = 2.0*atan(exp(-eta));
  double dist  = 0.;
  if (det == static_cast<int>(HcalBarrel)) {
    const double rLay[19] = {
      1836.0, 1902.0, 1962.0, 2022.0, 2082.0, 2142.0, 2202.0, 2262.0, 2322.0, 
      2382.0, 2448.0, 2514.0, 2580.0, 2646.0, 2712.0, 2776.0, 2862.5, 3847.0,
      4052.0};
    if (layer>0 && layer<20) dist += rLay[layer-1]*mm/sin(theta);
  } else {
    const double zLay[19] = {
      4034.0, 4032.0, 4123.0, 4210.0, 4297.0, 4384.0, 4471.0, 4558.0, 4645.0, 
      4732.0, 4819.0, 4906.0, 4993.0, 5080.0, 5167.0, 5254.0, 5341.0, 5428.0,
      5515.0};
    if (layer>0 && layer<20) dist += zLay[layer-1]*mm/cos(theta);
  }
  double tmp = dist/c_light/ns;
  LogDebug("HcalSim") << "HcalTestAnalysis::timeOfFlight " << tmp 
              << " for det/lay " << det << " " << layer 
              << " eta/theta " << eta << " " << theta/deg << " dist " 
              << dist;
  return tmp;
}

std::vector< int > HcalTestAnalysis::towersToAdd ( int  centre, int  nadd  )

private

Definition at line 110 of file HcalTestAnalysis.cc.

References funct::abs(), eta(), i, bookConverter::max, phi, phimax, phimin, and groupFilesInBlocks::temp.

Referenced by HcalTestAnalysis().

                                                                 {

  int etac = (centre/100)%100;
  int phic = (centre%100);
  int etamin, etamax, phimin, phimax;
  if (etac>0) {
    etamin = etac-nadd;
    etamax = etac+nadd;
  } else {
    etamin = etac;
    etamax = etac;
  }
  if (phic>0) {
    phimin = phic-nadd;
    phimax = phic+nadd;
  } else {
    phimin = phic;
    phimax = phic;
  }

  int nbuf, kount=0;
  nbuf = (etamax-etamin+1)*(phimax-phimin+1);
  std::vector<int> temp(2*nbuf);
  for (int eta=etamin; eta<=etamax; eta++) {
    for (int phi=phimin; phi<=phimax; phi++) {
      temp[kount] = (eta*100 + phi);
      temp[kount+nbuf] = std::max(abs(eta-etac),abs(phi-phic));
      kount++;
    }
  }

  edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Towers to be considered for"
              << " Central " << centre << " and " << nadd 
              << " on either side";
  for (int i=0; i<nbuf; i++)
    edm::LogInfo("HcalSim") << "HcalTestAnalysis: Tower[" << std::setw(3) << i 
                << "] " << temp[i] << " " << temp[nbuf+i];
  return temp;
}

void HcalTestAnalysis::update ( const BeginOfJob *  )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfJob * >.

Definition at line 151 of file HcalTestAnalysis.cc.

References fileName, names, numberingFromDDD, org, and tuplesManager.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                    {

  // Numbering From DDD
  edm::ESTransientHandle<DDCompactView> pDD;
  (*job)()->get<IdealGeometryRecord>().get(pDD);
  edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Initialise "
              << "HcalNumberingFromDDD for " << names[0];
  numberingFromDDD = new HcalNumberingFromDDD(names[0], (*pDD));

  // Ntuples
  tuplesManager.reset(new HcalTestHistoManager(fileName));

  // Numbering scheme
  org    = new HcalTestNumberingScheme(false);

}

void HcalTestAnalysis::update ( const BeginOfRun *  )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfRun * >.

Definition at line 169 of file HcalTestAnalysis.cc.

References HcalNumberingFromDDD::cell(), centralTower, HcalCellType::HcalCell::eta, eta(), eta0, HcalBarrel, HcalEndcap, HcalForward, cmsHarvester::loop, names, numberingFromDDD, HcalCellType::HcalCell::ok, org, phi, HcalCellType::HcalCell::phi, phi0, sd, HCalSD::setNumberingScheme(), AlCaHLTBitMon_QueryRunRegistry::string, and tmp.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                    {

  int irun = (*run)()->GetRunID();
  edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Begin of Run = " << irun;

  bool loop = true, eta = true, phi = true;
  int  etac = (centralTower/100)%100;
  if (etac == 0) {
    etac = 1;
    eta  = false;
  }
  int  phic = (centralTower%100);
  if (phic == 0) {
    phic = 1;
    phi  = false;
  }
  int  idet = static_cast<int>(HcalBarrel);
  while (loop) {
    HcalCellType::HcalCell tmp = numberingFromDDD->cell(idet,1,1,etac,phic);
    if (tmp.ok) {
      if (eta) eta0 = tmp.eta;
      if (phi) phi0 = tmp.phi;
      loop = false;
    } else if (idet == static_cast<int>(HcalBarrel)) {
      idet = static_cast<int>(HcalEndcap);
    } else if (idet == static_cast<int>(HcalEndcap)) {
      idet = static_cast<int>(HcalForward);
    } else {
      loop = false;
    }
  }

  edm::LogInfo("HcalSim") << "HcalTestAnalysis:: Central Tower " 
              << centralTower << " corresponds to eta0 = " << eta0 
              << " phi0 = " << phi0;
 
  std::string sdname = names[0];
  G4SDManager* sd = G4SDManager::GetSDMpointerIfExist();
  if (sd != 0) {
    G4VSensitiveDetector* aSD = sd->FindSensitiveDetector(sdname);
    if (aSD==0) {
      edm::LogWarning("HcalSim") << "HcalTestAnalysis::beginOfRun: No SD with "
                 << "name " << sdname << " in this Setup";
    } else {
      HCalSD* theCaloSD = dynamic_cast<HCalSD*>(aSD);
      edm::LogInfo("HcalSim") << "HcalTestAnalysis::beginOfRun: Finds SD with "
                  << "name " << theCaloSD->GetName() 
                  << " in this Setup";
      if (org) {
        theCaloSD->setNumberingScheme(org);
    edm::LogInfo("HcalSim") << "HcalTestAnalysis::beginOfRun: set a new "
                << "numbering scheme";
      }
    }
  } else {
    edm::LogWarning("HcalSim") << "HcalTestAnalysis::beginOfRun: Could not get"
                   << " SD Manager!";
  }

}

void HcalTestAnalysis::update ( const BeginOfEvent *  )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const BeginOfEvent * >.

Definition at line 231 of file HcalTestAnalysis.cc.

References edepEB, edepEE, edepHB, edepHE, edepHO, edepl, i, LogDebug, mudist, HcalTestHistoClass::setCounters(), and tuples.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                      {
 
  // create tuple object
  tuples = new HcalTestHistoClass();
  // Reset counters
  tuples->setCounters();
 
  int i = 0;
  edepEB = edepEE = edepHB = edepHE = edepHO = 0.;
  for (i = 0; i < 20; i++) edepl[i] = 0.;
  for (i = 0; i < 20; i++) mudist[i] = -1.;

  int iev = (*evt)()->GetEventID();
  LogDebug("HcalSim") <<"HcalTestAnalysis: Begin of event = " << iev;
}

void HcalTestAnalysis::update ( const EndOfEvent *  )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const EndOfEvent * >.

Definition at line 316 of file HcalTestAnalysis.cc.

References count, fill(), layerAnalysis(), LogDebug, qieAnalysis(), tuples, and tuplesManager.

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                    {

  count++;
  // Fill event input 
  fill(evt);
  LogDebug("HcalSim") << "HcalTestAnalysis:: ---  after Fill";

  // Qie analysis
  CLHEP::HepRandomEngine* engine = CLHEP::HepRandom::getTheEngine();
  qieAnalysis(engine);
  LogDebug("HcalSim") << "HcalTestAnalysis:: ---  after QieAnalysis";

  // Layers tuples filling
  layerAnalysis();
  LogDebug("HcalSim") << "HcalTestAnalysis:: ---  after LayerAnalysis";

  // Writing the data to the Tree
  tuplesManager->fillTree(tuples); // (no need to delete it...)
  tuples = 0; // but avoid to reuse it...
  LogDebug("HcalSim") << "HcalTestAnalysis:: --- after fillTree";

}

void HcalTestAnalysis::update ( const G4Step *  )

privatevirtual

This routine will be called when the appropriate signal arrives.

Implements Observer< const G4Step * >.

Definition at line 248 of file HcalTestAnalysis.cc.

References edepEB, edepEE, edepHB, edepHE, edepHO, edepl, eta(), eta0, log, LogDebug, MeV, mudist, mergeVDriftHistosByStation::name, NULL, phi, phi0, mathSSE::sqrt(), funct::tan(), and theta().

Referenced by progressbar.ProgressBar::__next__(), relval_steps.Matrix::__setitem__(), relval_steps.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), relval_steps.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                                                  {

  if (aStep != NULL) {
    G4VPhysicalVolume* curPV  = aStep->GetPreStepPoint()->GetPhysicalVolume();
    G4String name = curPV->GetName();
    name.assign(name,0,3);
    double edeposit = aStep->GetTotalEnergyDeposit();
    int    layer=-1;
    if (name == "EBR") {
      edepEB += edeposit;
    } else if (name == "EFR") {
      edepEE += edeposit;
    } else if (name == "HBS") {
      layer = (curPV->GetCopyNo()/10)%100;
      if (layer >= 0 && layer < 17) {
    edepHB += edeposit;
      } else {
    edm::LogWarning("HcalSim") << "HcalTestAnalysis::Error in HB "
                   << curPV->GetName() << curPV->GetCopyNo();
    layer = -1;
      }
    } else if (name == "HES") {
      layer = (curPV->GetCopyNo()/10)%100;
      if (layer >= 0 && layer < 19) {
    edepHE += edeposit;
      } else {
    edm::LogWarning("HcalSim") << "HcalTestAnalysis::Error in HE " 
                   << curPV->GetName() << curPV->GetCopyNo();
    layer = -1;
      }
    } else if (name == "HTS") {
      layer = (curPV->GetCopyNo()/10)%100;
      if (layer >= 17 && layer < 20) {
    edepHO += edeposit;
       } else {
     edm::LogWarning("HcalSim") << "HcalTestAnalysis::Error in HO " 
                    << curPV->GetName() << curPV->GetCopyNo();
    layer = -1;
       }
    }
    if (layer >= 0 && layer < 20) {
      edepl[layer] += edeposit;

      // Calculate the distance if it is a muon
      G4String part = aStep->GetTrack()->GetDefinition()->GetParticleName();
      if ((part == "mu-" || part == "mu+") && mudist[layer] < 0) {
        math::XYZPoint pos(aStep->GetPreStepPoint()->GetPosition().x(),
                           aStep->GetPreStepPoint()->GetPosition().y(),
                           aStep->GetPreStepPoint()->GetPosition().z());
        double theta   = pos.theta();
        double   eta   = -log(tan(theta * 0.5));
        double   phi   = pos.phi();
        double  dist   = sqrt ((eta-eta0)*(eta-eta0) + (phi-phi0)*(phi-phi0));
        mudist[layer]  = dist*std::sqrt(pos.perp2());
      }
    }

    if (layer >= 0 && layer < 20) {
      LogDebug("HcalSim") << "HcalTestAnalysis:: G4Step: " << name << " Layer "
              << std::setw(3) << layer << " Edep " << std::setw(6) 
              << edeposit/MeV << " MeV";
    }
  } else {
    edm::LogInfo("HcalSim") << "HcalTestAnalysis:: G4Step: Null Step";
  }
}

Member Data Documentation

int HcalTestAnalysis::addTower

private

Definition at line 68 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), and qieAnalysis().

std::vector<CaloHit> HcalTestAnalysis::caloHitCache

private

Definition at line 79 of file HcalTestAnalysis.h.

Referenced by fill(), and qieAnalysis().

int HcalTestAnalysis::centralTower

private

Definition at line 86 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), qieAnalysis(), and update().

unsigned int HcalTestAnalysis::count

private

Definition at line 89 of file HcalTestAnalysis.h.

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

double HcalTestAnalysis::edepEB

private

Definition at line 90 of file HcalTestAnalysis.h.

Referenced by layerAnalysis(), and update().

double HcalTestAnalysis::edepEE

private

Definition at line 90 of file HcalTestAnalysis.h.

Referenced by layerAnalysis(), and update().

double HcalTestAnalysis::edepHB

private

Definition at line 90 of file HcalTestAnalysis.h.

Referenced by layerAnalysis(), and update().

double HcalTestAnalysis::edepHE

private

Definition at line 90 of file HcalTestAnalysis.h.

Referenced by layerAnalysis(), and update().

double HcalTestAnalysis::edepHO

private

Definition at line 90 of file HcalTestAnalysis.h.

Referenced by layerAnalysis(), and update().

double HcalTestAnalysis::edepl[20]

private

Definition at line 91 of file HcalTestAnalysis.h.

Referenced by layerAnalysis(), and update().

double HcalTestAnalysis::eta0

private

Definition at line 85 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), qieAnalysis(), and update().

std::string HcalTestAnalysis::fileName

private

Definition at line 64 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), and update().

std::vector<int> HcalTestAnalysis::group_

private

Definition at line 80 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), and qieAnalysis().

double HcalTestAnalysis::mudist[20]

private

Definition at line 92 of file HcalTestAnalysis.h.

Referenced by layerAnalysis(), and update().

HcalQie* HcalTestAnalysis::myqie

private

Definition at line 67 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), qieAnalysis(), and ~HcalTestAnalysis().

std::vector<std::string> HcalTestAnalysis::names

private

Definition at line 84 of file HcalTestAnalysis.h.

Referenced by fill(), HcalTestAnalysis(), and update().

int HcalTestAnalysis::nGroup

private

Definition at line 81 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), and qieAnalysis().

int HcalTestAnalysis::nTower

private

Definition at line 81 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), and qieAnalysis().

HcalNumberingFromDDD* HcalTestAnalysis::numberingFromDDD

private

Definition at line 75 of file HcalTestAnalysis.h.

Referenced by fill(), qieAnalysis(), update(), and ~HcalTestAnalysis().

HcalTestNumberingScheme* HcalTestAnalysis::org

private

Definition at line 76 of file HcalTestAnalysis.h.

Referenced by fill(), qieAnalysis(), update(), and ~HcalTestAnalysis().

double HcalTestAnalysis::phi0

private

Definition at line 85 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), qieAnalysis(), and update().

std::vector<int> HcalTestAnalysis::tower_

private

Definition at line 80 of file HcalTestAnalysis.h.

Referenced by HcalTestAnalysis(), and qieAnalysis().

HcalTestHistoClass* HcalTestAnalysis::tuples

private

Definition at line 72 of file HcalTestAnalysis.h.

Referenced by layerAnalysis(), qieAnalysis(), update(), and ~HcalTestAnalysis().

std::auto_ptr<HcalTestHistoManager> HcalTestAnalysis::tuplesManager

private

Definition at line 71 of file HcalTestAnalysis.h.

Referenced by update().