calhit::CalHit Class Reference

Base class for creation of calibrated hits. More...

#include <CalHit.h>

Inheritance diagram for calhit::CalHit:

List of all members.

Public Member Functions
	CalHit (const char *version)
virtual	~CalHit ()
jobc::Result	Reco (edm::EventHandle &evt)
	read/write access to event
void	Update (const cfg::Config &)
Protected Member Functions
virtual jobc::Result	ADC2PEHit (recobase::CellHit *)
	base clase convertor
virtual jobc::Result	ADC2PEACHit (recobase::CellHit *)
	base clase convertor
virtual jobc::Result	ADC2MIPHit (recobase::CellHit *)
	base clase convertor
virtual jobc::Result	TDC2NSHit (recobase::CellHit *)
	base clase convertor
virtual jobc::Result	TDC2TT0Hit (recobase::CellHit *)
	base clase convertor
virtual jobc::Result	TDC2TTransHit (recobase::CellHit *)
	base clase convertor
virtual void	AddAttCorr (recobase::CellHit *)
	add cell hit to attenuation info
Private Member Functions
jobc::Result	GetData (edm::EventHandle &)
	get input data as needed
jobc::Result	AddToData (edm::EventHandle &)
	add calibrated information
jobc::Result	SaveData (edm::EventHandle &)
	save modified data to the record
void	PrintCalibs ()
	print out calibration status
jobc::Result	Raw (edm::EventHandle &evt)
	create new CalHits
jobc::Result	Read (edm::EventHandle &evt)
	read in existing CalHits
bool	DoCalib (recobase::CellHit_t, int=-1)
	Check or set whether calib type is requested.
Private Attributes
std::string	fInputFolder
	folder input
std::string	fOutputFolder
	folder for output
int	fDoCalib
bool	fWarnedPE
bool	fWarnedPEAC
bool	fWarnedMIP
bool	fWarnedNS
bool	fWarnedTT0
bool	fWarnedTTrans
std::vector< recobase::CellHit * >	fCellVec
	Working vector of cell hits.
std::vector< const recobase::CellHit * >	fCellInput
	const buffer of vector of cell hits

---

Detailed Description

Base class for creation of calibrated hits.

Definition at line 20 of file CalHit.h.

---

Constructor & Destructor Documentation

CalHit::CalHit (  const char *  version  )

create calibrated hits. This is a base class that does most of the work, expect for higher levels of refinement in the calibration techniques reflected in the derived classes. Definition at line 12 of file CalHit.cxx. References jobc::Module::SetCfgVersion(). : jobc::Module("CalHit"), fInputFolder("."), fOutputFolder("Hits"), fDoCalib(1), fWarnedPE(false), fWarnedPEAC(false), fWarnedMIP(false), fWarnedNS(false), fWarnedTT0(false), fWarnedTTrans(false) { this->SetCfgVersion(version); }

CalHit::~CalHit (   )  [virtual]

Definition at line 28 of file CalHit.cxx. { }

---

Member Function Documentation

jobc::Result CalHit::ADC2MIPHit (  recobase::CellHit *   )  [protected, virtual]

base clase convertor Definition at line 293 of file CalHit.cxx. References ADC2PEACHit(), fWarnedMIP, recobase::CellHit::PEAC(), jobc::Result, and recobase::CellHit::SetMIP(). Referenced by AddToData(). { // Boneheaded base clase - represents MIP-corrected ph if (!fWarnedMIP) { std::cerr<<"\nWarning! Using baseclass CalHit::ADC2MIPHit. "<< " Do not use this for Physics!"<<std::endl; fWarnedMIP=true; } // require PE attcor'd number before converting to MIPs float pe; if (!ch->PEAC(pe)) { if (!ADC2PEACHit(ch)) return jobc::kFailed; ch->PEAC(pe); } ch->SetMIP(pe/30.); return jobc::kPassed; } // ADC2MIPHit

jobc::Result CalHit::ADC2PEACHit (  recobase::CellHit *   )  [protected, virtual]

base clase convertor Reimplemented in calhit::MyHit. Definition at line 270 of file CalHit.cxx. References ADC2PEHit(), fWarnedPEAC, recobase::CellHit::PE(), jobc::Result, and recobase::CellHit::SetPEAC(). Referenced by ADC2MIPHit(), and AddToData(). { // Boneheaded base clase - represents attenuation-corrected PEs if (!fWarnedPEAC) { std::cerr<<"\nWarning! Using baseclass CalHit::ADC2PEACHit. "<< " Do not use this for Physics!"<<std::endl; fWarnedPEAC=true; } // require PE number before applying attenuation correction float pe; if (!ch->PE(pe)) { if (!ADC2PEHit(ch)) return jobc::kFailed; ch->PE(pe); } ch->SetPEAC(pe/0.55); return jobc::kPassed; } // ADC2PEACHit

jobc::Result CalHit::ADC2PEHit (  recobase::CellHit *   )  [protected, virtual]

base clase convertor Definition at line 248 of file CalHit.cxx. References recobase::CellHit::ADC(), recobase::CellHit::Channel(), fWarnedPE, jobc::Result, and recobase::CellHit::SetPE(). Referenced by calhit::MyHit::ADC2PEACHit(), ADC2PEACHit(), and AddToData(). { // Boneheaded base clase if (!fWarnedPE) { std::cerr<<"\nWarning! Using baseclass CalHit::ADC2PEHit. "<< " Do not use this for Physics!"<<std::endl; fWarnedPE=true; } unsigned short adc; if (!ch->ADC(adc)) { std::cerr<<"CalHit::ADC2PEHit(): raw adc not set in cell hit for channel "<< ch->Channel()<<std::endl; return jobc::kFailed; } ch->SetPE(adc/4.095);// dummy value return jobc::kPassed; } // ADC2PEHit

void CalHit::AddAttCorr (  recobase::CellHit *   )  [protected, virtual]

add cell hit to attenuation info Definition at line 239 of file CalHit.cxx. Referenced by Raw(), and Read(). { // add this cell to attenuation correction // but a dummy for the time being... }

jobc::Result CalHit::AddToData (  edm::EventHandle &   )  [private]

add calibrated information Definition at line 154 of file CalHit.cxx. References ADC2MIPHit(), ADC2PEACHit(), ADC2PEHit(), DoCalib(), fCellVec, jobc::Result, TDC2NSHit(), TDC2TT0Hit(), and TDC2TTransHit(). Referenced by Reco(). { std::vector<recobase::CellHit*>::iterator iter=fCellVec.begin(); while (iter!=fCellVec.end()) { recobase::CellHit* ch=*iter; // create Photo Electron equivalent ADC if (DoCalib(recobase::kPE)) if (ADC2PEHit(ch)==jobc::kFailed) return jobc::kFailed; //... with Attenuation correction if (DoCalib(recobase::kPEAC)) if (ADC2PEACHit(ch)==jobc::kFailed) return jobc::kFailed; //... in MIP units if (DoCalib(recobase::kMIP)) if (ADC2MIPHit(ch)==jobc::kFailed) return jobc::kFailed; // tdc in ns if (DoCalib(recobase::kTNS)) if (TDC2NSHit(ch)==jobc::kFailed) return jobc::kFailed; // tdc in ns, with transit time correction if (DoCalib(recobase::kTTrans)) if (TDC2TTransHit(ch)==jobc::kFailed) return jobc::kFailed; // tdc in ns, with transit time correction, from a T0 if (DoCalib(recobase::kTT0)) if (TDC2TT0Hit(ch)==jobc::kFailed) return jobc::kFailed; iter++; } return jobc::kPassed; }

bool CalHit::DoCalib (  recobase::CellHit_t  , int  = -1 )  [private]

Check or set whether calib type is requested. Definition at line 220 of file CalHit.cxx. References recobase::CellHit_t, and fDoCalib. Referenced by AddToData(), GetData(), PrintCalibs(), Raw(), Read(), and Update(). { // if set<0, then caller is asking what calibrations are done // set=0, caller wants to turn off a calibration // set>0, caller wants to turn on a calibration if (set < 0 ) { return (fDoCalib & (1<<level) ); } else if (set == 0 ) { fDoCalib = fDoCalib&=(~(1<<level)); return false; } else { fDoCalib = fDoCalib|=(1<<level); return true; } }

jobc::Result CalHit::GetData (  edm::EventHandle &   )  [private]

get input data as needed Definition at line 59 of file CalHit.cxx. References DoCalib(), fCellInput, fCellVec, Raw(), Read(), and jobc::Result. Referenced by Reco(). { // clear the local copy fCellVec.clear(); fCellInput.clear(); // get the data from the appropriate place if (DoCalib(recobase::kRAW)) { // create cal hits, and copy over raw data if (Raw(evt)==jobc::kFailed) return jobc::kFailed; } else { // start with existing cal hits if (Read(evt)==jobc::kFailed) return jobc::kFailed; } return jobc::kPassed; }

void CalHit::PrintCalibs (   )  [private]

print out calibration status Definition at line 433 of file CalHit.cxx. References DoCalib(). Referenced by Update(). { std::cerr<<"Calibrations Requested:"<<std::endl; std::cerr<<" Starting from Raw Data: " ; if (DoCalib(recobase::kRAW)) std::cerr<<" Yes"<<std::endl; else std::cerr<<" NO"<<std::endl; std::cerr<<" PE Units: "; if (DoCalib(recobase::kPE)) std::cerr<<" Yes"<<std::endl; else std::cerr<<" NO"<<std::endl; std::cerr<<" Attenuation Correction (to PEs): "; if (DoCalib(recobase::kPEAC)) std::cerr<<" Yes"<<std::endl; else std::cerr<<" NO"<<std::endl; std::cerr<<" MIP Units: "; if (DoCalib(recobase::kMIP)) std::cerr<<" Yes"<<std::endl; else std::cerr<<" NO"<<std::endl; std::cerr<<" Time in ns units: " ; if (DoCalib(recobase::kTNS)) std::cerr<<" Yes"<<std::endl; else std::cerr<<" NO"<<std::endl; std::cerr<<" T-zero corrected time in ns: " ; if (DoCalib(recobase::kTT0)) std::cerr<<" Yes"<<std::endl; else std::cerr<<" NO"<<std::endl; std::cerr<<" Transit time corrected time in ns: " ; if (DoCalib(recobase::kTTrans)) std::cerr<<" Yes"<<std::endl; else std::cerr<<" NO"<<std::endl; }

jobc::Result CalHit::Raw (  edm::EventHandle &  evt  )  [private]

create new CalHits Definition at line 81 of file CalHit.cxx. References rawdata::RawDigit::ADC(), AddAttCorr(), rawdata::RawDigit::Channel(), DoCalib(), fCellVec, fInputFolder, rawdata::RawDigit::GetCell(), rawdata::RawDigit::GetPlane(), edm::EventHandle::Raw(), jobc::Result, recobase::CellHit::SetCell(), recobase::CellHit::SetPlane(), and rawdata::RawDigit::TDC(). Referenced by GetData(). { // get raw data std::vector<const rawdata::RawDigit*> diglist; try { evt.Raw().Get(fInputFolder.c_str(),diglist); } catch (edm::Exception e) { std::cerr << "Error retrieving digit list, while looking for raw "<< "digits in CalHit::Raw()"<<std::endl; return jobc::kFailed; } // loop over all raw digits, and create basic calibrated hits in local // record std::vector<const rawdata::RawDigit*>::iterator iter=diglist.begin(); while (iter!=diglist.end()) { const rawdata::RawDigit *rd=*iter; recobase::CellHit* ch=new recobase::CellHit(rd->Channel(), rd->ADC(0),rd->TDC(0)); //printf("%d %d\n",rd->GetCell(),rd->GetPlane()); //printf("%d ",rd->ADC(0)); ch->SetCell(rd->GetCell()); ch->SetPlane(rd->GetPlane()); fCellVec.push_back(ch); // if we'll be doing an attenuation correction, collect // necessary information here if (DoCalib(recobase::kPEAC)) AddAttCorr(ch); iter++; } return jobc::kPassed; }

jobc::Result CalHit::Read (  edm::EventHandle &  evt  )  [private]

read in existing CalHits Definition at line 122 of file CalHit.cxx. References AddAttCorr(), DoCalib(), fCellInput, fCellVec, fInputFolder, edm::EventHandle::Reco(), and jobc::Result. Referenced by GetData(). { // get existing CalHit's try { evt.Reco().Get(fInputFolder.c_str(),fCellInput); } catch (edm::Exception e) { std::cerr << "Error retrieving CellHit list in "<<fInputFolder.c_str() <<std::endl; return jobc::kFailed; } // if we will need attenuation correction, accumulate necessary info // here std::vector<const recobase::CellHit*>::iterator iter=fCellInput.begin(); std::cerr<<"Found "<<fCellInput.size()<<" cell hits in input"<<std::endl; while (iter!=fCellInput.end()) { const recobase::CellHit* oldch=*iter; recobase::CellHit* ch=new recobase::CellHit(*oldch); fCellVec.push_back(ch); if (DoCalib(recobase::kPEAC)) { AddAttCorr(ch); } iter++; } return jobc::kPassed; }

jobc::Result CalHit::Reco (  edm::EventHandle &  evt  )  [virtual]

read/write access to event Reimplemented from jobc::Module. Definition at line 33 of file CalHit.cxx. References AddToData(), GetData(), jobc::Result, and SaveData(). { // base class Reco method. if (jobc::kFailed==GetData(evt)) { std::cerr <<"CalHit::GetData failed. Returning"<<std::endl; return jobc::kFailed; } if (jobc::kFailed==AddToData(evt)) { std::cerr <<"CalHit::AddToData failed. Returning"<<std::endl; return jobc::kFailed; } if (jobc::kFailed==SaveData(evt)) { std::cerr <<"CalHit::SaveData failed. Returning"<<std::endl; return jobc::kFailed; } return jobc::kPassed; }

jobc::Result CalHit::SaveData (  edm::EventHandle &   )  [private]

save modified data to the record Definition at line 193 of file CalHit.cxx. References fCellVec, fOutputFolder, edm::EventHandle::Reco(), and jobc::Result. Referenced by Reco(). { if (evt.Reco().GetFolder(fOutputFolder.c_str())) { std::cerr << "Output Folder reco/"<<fOutputFolder.c_str()<< " already exists"<<std::endl; return jobc::kFailed; } TFolder* folder=evt.Reco().MakeFolder(fOutputFolder.c_str()); int nsave=evt.Reco().PutVector(fCellVec,fOutputFolder.c_str()); if (nsave!=fCellVec.size()) { std::cerr <<"CalHit::SaveData: wrong number CellHits saved: "<< nsave<<" instead of expected "<<fCellVec.size()<<std::endl; return jobc::kFailed; } float pe,peac; fCellVec[0]->PE(pe); fCellVec[0]->PEAC(peac); std::cerr << "saved "<<nsave<<" CellHits. first pe="<<pe<<", atcor="<<peac<< std::endl; return jobc::kPassed; }

jobc::Result CalHit::TDC2NSHit (  recobase::CellHit *  ch  )  [protected, virtual]

base clase convertor < CALIB_HOOK Definition at line 316 of file CalHit.cxx. References recobase::CellHit::Channel(), fWarnedNS, jobc::Result, recobase::CellHit::SetTNS(), and recobase::CellHit::TDC(). Referenced by AddToData(), and TDC2TT0Hit(). { // Boneheaded base clase - represents time in ns if (!fWarnedNS) { std::cerr<<"\nWarning! Using baseclass CalHit::TDC2NSHit. "<< " Do not use this for Physics!"<<std::endl; fWarnedNS=true; } unsigned short tdc; if (!ch->TDC(tdc)) { std::cerr<<"CalHit::TDC2NSHit(): raw tdc not set in cell hit for channel "<< ch->Channel()<<std::endl; return jobc::kFailed; } ch->SetTNS(tdc*1.5e3); return jobc::kPassed; } // TDC2NSHit

jobc::Result CalHit::TDC2TT0Hit (  recobase::CellHit *  ch  )  [protected, virtual]

base clase convertor < CALIB_HOOK Definition at line 338 of file CalHit.cxx. References fWarnedTT0, jobc::Result, recobase::CellHit::SetTT0(), TDC2NSHit(), and recobase::CellHit::TNS(). Referenced by AddToData(), and TDC2TTransHit(). { // Boneheaded base clase - represents hit times corrected for global event // t0 if (!fWarnedTT0) { std::cerr<<"\nWarning! Using baseclass CalHit::TDC2TT0Hit. "<< " Do not use this for Physics!"<<std::endl; fWarnedTT0=true; } // require NS number before applying t0 double ns; if (!ch->TNS(ns)) { if (!TDC2NSHit(ch)) return jobc::kFailed; ch->TNS(ns); // should be there now... } ch->SetTT0(ns+24.6); return jobc::kPassed; } // TDC2TT0Hit

jobc::Result CalHit::TDC2TTransHit (  recobase::CellHit *  ch  )  [protected, virtual]

base clase convertor < CALIB_HOOK Definition at line 362 of file CalHit.cxx. References fWarnedTTrans, jobc::Result, recobase::CellHit::SetTTrans(), TDC2TT0Hit(), and recobase::CellHit::TT0(). Referenced by AddToData(). { // Boneheaded base clase - represents transit-time corrected hits if (!fWarnedTTrans) { std::cerr<<"\nWarning! Using baseclass CalHit::TDC2TTransHit. "<< " Do not use this for Physics!"<<std::endl; fWarnedTTrans=true; } // require t0-corrected time before handling transit time double tt0; if (!ch->TT0(tt0)) { if (!TDC2TT0Hit(ch)) return jobc::kFailed; // try setting ns ch->TT0(tt0); // should be there now... } ch->SetTTrans(tt0-7.5*5); return jobc::kPassed; } // TDC2TTransHit

void CalHit::Update (  const cfg::Config &   )  [virtual]

Implements cfg::Observer. Definition at line 384 of file CalHit.cxx. References DoCalib(), fDoCalib, fInputFolder, fOutputFolder, and PrintCalibs(). { c("Input"). Get(fInputFolder); c("Output"). Get(fOutputFolder); if (fInputFolder==fOutputFolder) { std::cerr<<"CalHit::CalHit() - output folder ("<<fOutputFolder.c_str()<< ") is not allowed to equal input folder ("<<fInputFolder.c_str() <<")"<< std::endl; assert(0); } std::cerr<<"Set Input Folder ="<<fInputFolder.c_str()<<std::endl; std::cerr<<"Set Output Folder ="<<fOutputFolder.c_str()<<std::endl; // zero out calibrations fDoCalib=0; std::string calibs; c("Calibs"). Get(calibs); // set calibrations according to the string std::istringstream icalibs(calibs); while (icalibs.good()) { std::string tmp; icalibs>>tmp; if (tmp=="RAW") DoCalib(recobase::kRAW,1); else if (tmp=="PE") DoCalib(recobase::kPE,1); else if (tmp=="PEAC") DoCalib(recobase::kPEAC,1); else if (tmp=="MIP") DoCalib(recobase::kMIP,1); else if (tmp=="TNS") DoCalib(recobase::kTNS,1); else if (tmp=="TT0") DoCalib(recobase::kTT0,1); else if (tmp=="TTRANS") DoCalib(recobase::kTTrans,1); else { std::cerr<<"Unknown calibration requested: "<<tmp.c_str()<<std::endl; abort(); } } // parse calibration string PrintCalibs(); }

---

Member Data Documentation

std::vector<const recobase::CellHit*> calhit::CalHit::fCellInput [private]

const buffer of vector of cell hits Definition at line 64 of file CalHit.h. Referenced by GetData(), and Read().

std::vector<recobase::CellHit*> calhit::CalHit::fCellVec [private]

Working vector of cell hits. Definition at line 63 of file CalHit.h. Referenced by AddToData(), GetData(), Raw(), Read(), and SaveData().

int calhit::CalHit::fDoCalib [private]

Definition at line 55 of file CalHit.h. Referenced by DoCalib(), and Update().

std::string calhit::CalHit::fInputFolder [private]

folder input Definition at line 51 of file CalHit.h. Referenced by Raw(), Read(), and Update().

std::string calhit::CalHit::fOutputFolder [private]

folder for output Definition at line 52 of file CalHit.h. Referenced by SaveData(), and Update().

bool calhit::CalHit::fWarnedMIP [private]

Definition at line 58 of file CalHit.h. Referenced by ADC2MIPHit().

bool calhit::CalHit::fWarnedNS [private]

Definition at line 59 of file CalHit.h. Referenced by TDC2NSHit().

bool calhit::CalHit::fWarnedPE [private]

Definition at line 56 of file CalHit.h. Referenced by ADC2PEHit().

bool calhit::CalHit::fWarnedPEAC [private]

Definition at line 57 of file CalHit.h. Referenced by ADC2PEACHit().

bool calhit::CalHit::fWarnedTT0 [private]

Definition at line 60 of file CalHit.h. Referenced by TDC2TT0Hit().

bool calhit::CalHit::fWarnedTTrans [private]

Definition at line 61 of file CalHit.h. Referenced by TDC2TTransHit().

---

The documentation for this class was generated from the following files:

• CalHit.h
• CalHit.cxx

---