#include <CastorTestAnalysis.h>

Inheritance diagram for CastorTestAnalysis:

Public Member Functions
	CastorTestAnalysis (const edm::ParameterSet &p)
virtual	~CastorTestAnalysis ()
Private Member Functions
void	Finish ()
void	getCastorBranchData (const CaloG4HitCollection *hc)
void	update (const EndOfRun *run)
	This routine will be called when the appropriate signal arrives.
void	update (const G4Step *step)
	This routine will be called when the appropriate signal arrives.
void	update (const EndOfEvent *evt)
	This routine will be called when the appropriate signal arrives.
void	update (const BeginOfRun *run)
	This routine will be called when the appropriate signal arrives.
void	update (const BeginOfJob *run)
	This routine will be called when the appropriate signal arrives.
void	update (const BeginOfEvent *evt)
	This routine will be called when the appropriate signal arrives.
Private Attributes
Float_t	castoreventarray [11]
TNtuple *	castoreventntuple
TFile *	castorOutputEventFile
TFile *	castorOutputStepFile
Float_t	castorsteparray [14]
TNtuple *	castorstepntuple
int	doNTcastorevent
int	doNTcastorstep
int	eventGlobalHit
int	eventIndex
std::string	eventNtFileName
int	stepIndex
std::string	stepNtFileName
CastorNumberingScheme *	theCastorNumScheme
int	verbosity

Detailed Description

Definition at line 65 of file CastorTestAnalysis.h.

Constructor & Destructor Documentation

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

Definition at line 35 of file CastorTestAnalysis.cc.

References castoreventntuple, castorstepntuple, gather_cfg::cout, doNTcastorevent, doNTcastorstep, eventNtFileName, edm::ParameterSet::getParameter(), stepNtFileName, and verbosity.

                                                               {

  edm::ParameterSet m_Anal = p.getParameter<edm::ParameterSet>("CastorTestAnalysis");
  verbosity                = m_Anal.getParameter<int>("Verbosity");
  doNTcastorstep           = m_Anal.getParameter<int>("StepNtupleFlag");
  doNTcastorevent          = m_Anal.getParameter<int>("EventNtupleFlag");
  stepNtFileName           = m_Anal.getParameter<std::string>("StepNtupleFileName");
  eventNtFileName          = m_Anal.getParameter<std::string>("EventNtupleFileName");

  if (verbosity > 0)
   std::cout<<std::endl;
   std::cout<<"============================================================================"<<std::endl;
   std::cout << "CastorTestAnalysis:: Initialized as observer"<< std::endl;
   if (doNTcastorstep  > 0){
     std::cout <<" Step Ntuple will be created"<< std::endl;
     std::cout <<" Step Ntuple file: "<<stepNtFileName<<std::endl;
   }
   if (doNTcastorevent > 0){
     std::cout <<" Event Ntuple will be created"<< std::endl;
     std::cout <<" Step Ntuple file: "<<stepNtFileName<<std::endl;
   }
   std::cout<<"============================================================================"<<std::endl;
   std::cout<<std::endl;
  
  if (doNTcastorstep  > 0)  
  castorstepntuple = new TNtuple("NTcastorstep","NTcastorstep","evt:trackid:charge:pdgcode:x:y:z:stepl:stepe:eta:phi:vpx:vpy:vpz");
  
  if (doNTcastorevent  > 0)
  castoreventntuple = new TNtuple("NTcastorevent","NTcastorevent","evt:ihit:detector:sector:module:enem:totalenergy:hitenergy:x:y:z");
}

CastorTestAnalysis::~CastorTestAnalysis ( ) [virtual]

Definition at line 66 of file CastorTestAnalysis.cc.

References gather_cfg::cout, Finish(), and verbosity.

                                        {
  //destructor of CastorTestAnalysis
    
    Finish();
  if (verbosity > 0) {
    std::cout << std::endl << "End of CastorTestAnalysis"
              << std::endl; 
  }
  
  std::cout<<"CastorTestAnalysis: End of process"<<std::endl;
  
}

Member Function Documentation

void CastorTestAnalysis::Finish ( ) [private]

Definition at line 354 of file CastorTestAnalysis.cc.

References castoreventntuple, castorOutputEventFile, castorOutputStepFile, castorstepntuple, gather_cfg::cout, doNTcastorevent, and doNTcastorstep.

Referenced by ~CastorTestAnalysis().

                                {
  if (doNTcastorstep) {
    castorOutputStepFile->cd();
    castorstepntuple->Write();
    std::cout << "CastorTestAnalysis: Ntuple step  written" <<std::endl;
    castorOutputStepFile->Close();
    std::cout << "CastorTestAnalysis: Step file closed" << std::endl;
  }
  
   if (doNTcastorevent) {
   castorOutputEventFile->cd();
   castoreventntuple->Write("",TObject::kOverwrite);
   std::cout << "CastorTestAnalysis: Ntuple event written" << std::endl;   
   castorOutputEventFile->Close();
   std::cout << "CastorTestAnalysis: Event file closed" << std::endl;
 }

}

void CastorTestAnalysis::getCastorBranchData ( const CaloG4HitCollection * hc ) [private]

Definition at line 297 of file CastorTestAnalysis.cc.

References castoreventarray, castoreventntuple, eventGlobalHit, eventIndex, CaloG4Hit::getEnergyDeposit(), CaloG4Hit::getPosition(), CaloG4Hit::getUnitID(), ntcastore_detector, ntcastore_enem, ntcastore_enhad, ntcastore_evt, ntcastore_hitenergy, ntcastore_ihit, ntcastore_module, ntcastore_sector, ntcastore_x, ntcastore_y, ntcastore_z, theCastorNumScheme, and CastorNumberingScheme::unpackIndex().

Referenced by update().

                                                                           {

    int nentries = hc->entries();
    
    if (nentries > 0) {
      
      unsigned int volumeID=0;
      int det=0, zside, sector, zmodule;
      std::map<int,float,std::less<int> > themap;
      double totalEnergy = 0;
      double hitEnergy = 0;
      double en_in_sd = 0.;

      for (int ihit = 0; ihit < nentries; ihit++) {
        CaloG4Hit* aHit = (*hc)[ihit];
        totalEnergy += aHit->getEnergyDeposit();
      }
    
      for (int ihit = 0; ihit < nentries; ihit++) {
        CaloG4Hit* aHit = (*hc)[ihit];
        volumeID = aHit->getUnitID();
        hitEnergy = aHit->getEnergyDeposit();
        en_in_sd += aHit->getEnergyDeposit();
//      double enEm = aHit->getEM();
//      double enHad = aHit->getHadr();
        
        themap[volumeID] += aHit->getEnergyDeposit();
        // int det, zside, sector, zmodule;
        theCastorNumScheme->unpackIndex(volumeID, zside, sector,zmodule);

        // det = 2 ;  //  det=2/3 for CAFI/CAPL
        
        castoreventarray[ntcastore_evt]       = (float)eventIndex;
//      castoreventarray[ntcastore_ihit]      = (float)ihit;
        castoreventarray[ntcastore_ihit]      = (float)eventGlobalHit;
        castoreventarray[ntcastore_detector]  = (float)det;
        castoreventarray[ntcastore_sector]    = (float)sector;
        castoreventarray[ntcastore_module]    = (float)zmodule;
        castoreventarray[ntcastore_enem]      = en_in_sd;
        castoreventarray[ntcastore_enhad]     = totalEnergy;
        castoreventarray[ntcastore_hitenergy] = hitEnergy;
        castoreventarray[ntcastore_x]         = aHit->getPosition().x();
        castoreventarray[ntcastore_y]         = aHit->getPosition().y();
        castoreventarray[ntcastore_z]         = aHit->getPosition().z();
//      castoreventarray[ntcastore_x]         = aHit->getEntry().x();
//      castoreventarray[ntcastore_y]         = aHit->getEntry().y();
//      castoreventarray[ntcastore_z]         = aHit->getEntry().z();
        
        castoreventntuple->Fill(castoreventarray);
        
        eventGlobalHit++ ;
      }
    } // nentries > 0
}

void CastorTestAnalysis::update ( const EndOfRun * ) [private, virtual]