#include <Blinder.h>
Inheritance diagram for Blinder:
Public Member Functions | |
| Blinder () | |
| ~Blinder () | |
| virtual JobCResult | Ana (const MomNavigator *mom) |
| virtual JobCResult | Reco (MomNavigator *mom) |
| const Registry & | DefaultConfig () const |
| void | Config (const Registry &r) |
Private Member Functions | |
| Bool_t | GetInfo (const MomNavigator *mom, Int_t &run, Int_t &snarl, Int_t &epoch, Float_t &length, Float_t &energy) |
| const REROOT_NeuKin * | GetTruth (const MomNavigator *mom) |
Private Attributes | |
| Int_t | fSeed |
| TFile * | fFile |
| TNtuple * | fTree |
| TRandom * | fRandom |
| TRandom * | fParameterChooser |
| BlindingFunction | fBlindingFunction |
|
|
Definition at line 81 of file Blinder.cxx. References fBlindingFunction, fFile, fParameterChooser, fRandom, fSeed, fTree, and Blinder::BlindingFunction::SetParameters(). 00082 {
00083 fSeed = 0;
00084 fFile = 0;
00085 fTree = 0;
00086
00087 fRandom = new TRandom3;
00088 fParameterChooser = new TRandom3;
00089 fBlindingFunction.SetParameters(fParameterChooser,fSeed);
00090 }
|
|
|
Definition at line 94 of file Blinder.cxx. 00095 {
00096 if(fFile) {
00097 fFile->cd();
00098 fTree->Write();
00099 fFile->Write();
00100 fFile->Close();
00101 }
00102 if(fRandom) delete fRandom;
00103 if(fParameterChooser) delete fParameterChooser;
00104 }
|
|
|
This function returns kPassed if the event is allowed in the open dataset, or kFailed if the event should be in the blind data set. If the relevant info isn't available, it returns kPassed. Reimplemented from JobCModule. Definition at line 108 of file Blinder.cxx. References Blinder::BlindingFunction::EventBlindProb(), fBlindingFunction, fParameterChooser, fRandom, fSeed, GetInfo(), run(), and Blinder::BlindingFunction::SetParameters(). 00109 {
00116 Int_t run;
00117 Int_t snarl;
00118 Int_t epoch;
00119 Float_t length;
00120 Float_t energy;
00121
00122
00123 if(GetInfo(mom,run,snarl,epoch,length,energy)) {
00124 // Set the parameters from the epoch:
00125 fBlindingFunction.SetParameters(fParameterChooser,fSeed + 3 * epoch);
00126
00127 // Find the probability of blinding this event:
00128 Float_t prob = fBlindingFunction.EventBlindProb(length,energy);
00129
00130 // Get a pseudo-random number specific for this event:
00131 fRandom->SetSeed(run*1000 + snarl);
00132 Float_t r = fRandom->Rndm();
00133
00134 // Allow event to be open if r<prob
00135 if(r<=prob) return JobCResult::kPassed;
00136 else return JobCResult::kFailed;
00137 };
00138
00139 // No blindness info obtained; return default. (This is probably a monitoring block.)
00140 return JobCResult::kPassed;
00141 }
|
|
|
Return the actual configuration. If your module directly pulls its configuration from the fConfig Registry, you don't need to override this. Override if you have local config variables. Reimplemented from JobCModule. Definition at line 370 of file Blinder.cxx. References fBlindingFunction, fParameterChooser, fSeed, Registry::Get(), and Blinder::BlindingFunction::SetParameters(). 00371 {
00372 //======================================================================
00373 // Configure the module given the Registry r
00374 //======================================================================
00375
00376 r.Get("Seed",fSeed);
00377
00378 // reset the paramaters:
00379 fBlindingFunction.SetParameters(fParameterChooser,fSeed);
00380 }
|
|
|
Get the default configuration registry. This should normally be overridden. One useful idiom is to implement it like: const Registry& MyModule::DefaultConfig() const { static Registry cfg; // never is destroyed if (cfg.Size()) return cfg; // already filled it // set defaults: cfg.Set("TheAnswer",42); cfg.Set("Units","unknown"); return cfg; } Reimplemented from JobCModule. Definition at line 348 of file Blinder.cxx. References JobCModule::GetName(), Registry::LockValues(), Registry::Set(), and Registry::UnLockValues(). 00349 {
00350 //======================================================================
00351 // Supply the default configuration for the module
00352 //======================================================================
00353 static Registry r; // Default configuration for module
00354
00355 // Set name of config
00356 std::string name = this->GetName();
00357 name += ".config.default";
00358 r.SetName(name.c_str());
00359
00360 // Set values in configuration
00361 r.UnLockValues();
00362 r.Set("Seed", 10127); // Seed with Nathaniel's fnal ID number.. that's pretty random.
00363 r.LockValues();
00364
00365 return r;
00366 }
|
|
||||||||||||||||||||||||||||
|
Definition at line 144 of file Blinder.cxx. References digit(), RawDigit::GetADC(), RawDigit::GetChannel(), RawDigitDataBlock::GetDatumIter(), VldContext::GetDetector(), MomNavigator::GetFragment(), PlexSEIdAltL::GetPlane(), RawRecord::GetRawBlockIter(), RawRecord::GetRawHeader(), RawDaqHeader::GetRun(), VldContext::GetSimFlag(), RawDaqSnarlHeader::GetSnarl(), VldContext::GetTimeStamp(), RecMinosHdr::GetVldContext(), MAXMSG, and run(). Referenced by Ana(), and Reco(). 00150 {
00151 //
00152 // This function pulls the neccessary info from the current Mom to find
00153 // the run, snarl, event length (with a basic clustering algorithm) and
00154 // event energy (in raw ADC counts).
00155 //
00156 run = 0;
00157 snarl = 0;
00158 length = 0;
00159 energy = 0;
00160 epoch = 0;
00161
00162 // Find raw data.
00163 const RawRecord* rr = dynamic_cast<const RawRecord*>(mom->GetFragment("RawRecord"));
00164 if(rr==0) return false;
00165
00166 const RawDaqSnarlHeader* snarlHead= dynamic_cast<const RawDaqSnarlHeader*>(rr->GetRawHeader());
00167 if(snarlHead==0) return false;
00168
00169
00170 TIter rdbit = rr->GetRawBlockIter();
00171 TObject *tob;
00172 const RawDigitDataBlock *rddb = 0;
00173 while ((tob = rdbit())) {
00174 if (tob->InheritsFrom("RawDigitDataBlock")) {
00175 rddb = (const RawDigitDataBlock *) tob;
00176 }
00177 }
00178 if(rddb==0) return false;
00179
00180 if(snarlHead->GetVldContext().GetDetector()!=Detector::kFar) return false;
00181
00182 if(snarlHead->GetVldContext().GetSimFlag()!=SimFlag::kData) epoch = 10;
00183
00184 // Put epoch here:
00185 // These lines change the blinding function at certain
00186 // watershed dates.
00187 const VldTimeStamp kShutdown2006(2006,4,15,0,0,0); // April 15, 2006, midnight
00188 if(snarlHead->GetVldContext().GetTimeStamp() > kShutdown2006) epoch +=1;
00189
00190 MAXMSG("Blinder",Msg::kInfo,1) << "Blinding FD data. Epoch = " << epoch << endl;
00191
00192 // Get plex info.
00193 PlexHandle plex(snarlHead->GetVldContext());
00194
00195 // Get the header data.
00196 run = snarlHead->GetRun();
00197 snarl = snarlHead->GetSnarl();
00198
00199 // Map of the planes.
00200 std::vector<int> planes(500,0);
00201
00202 // Loop the digits.
00203 TIter rdit = rddb->GetDatumIter();
00204 const RawDigit *digit;
00205 Int_t rawdigitindex = -1;
00206 while ((digit = (const RawDigit *) rdit())) {
00207 ++rawdigitindex;
00208
00209 if(plex.GetReadoutType(digit->GetChannel()) != ReadoutType::kScintStrip)
00210 continue;
00211
00212 // Find plane.
00213 PlexSEIdAltL altl = plex.GetSEIdAltL(digit->GetChannel());
00214 Int_t plane = altl.GetPlane();
00215
00216 // Find energy.
00217 Int_t adc = digit->GetADC();
00218
00219 // add it up.
00220 if((plane>=0)&&(plane<490))
00221 planes[plane] ++;
00222
00223 energy+=(float)adc;
00224 }
00225
00226 // do the number-of-planes matching. Look for the largest 'cluster'
00227 int start = 0;
00228 int end = 0;
00229 int longest = 0;
00230 while(start<500) {
00231 if(planes[start]==0) {
00232 start++;
00233 continue;
00234 };
00235
00236 // Found a starting place.
00237 end = start;
00238 int gap = 0;
00239 int thislen = 0;
00240 while(end<490) {
00241 end++;
00242 if(planes[end]==0) gap++;
00243 else thislen = end-start;
00244 if(gap>2) break;
00245 }
00246
00247 if(thislen > longest) longest = thislen;
00248
00249 start=end;
00250 }
00251 length = longest;
00252
00253 return true;
00254 }
|
|
|
Definition at line 323 of file Blinder.cxx. References RecDataRecord< T >::FindComponent(), and MomNavigator::FragmentIter(). Referenced by Reco(). 00324 {
00325 // Look for the simSnarl in Mom.
00326 SimSnarlRecord* simsnarl = 0;
00327 TObject* tobj;
00328 TIter fragiter = mom->FragmentIter();
00329 while( ( tobj = fragiter.Next() ) ) {
00330 simsnarl = dynamic_cast<SimSnarlRecord*>(tobj);
00331 if(simsnarl) break;
00332 }
00333
00334 if(simsnarl==0) return 0;
00335
00336 // Get the NeuKin list.
00337 const TClonesArray* neuKins = dynamic_cast<const TClonesArray*>(simsnarl->FindComponent("TClonesArray","NeuKinList"));
00338 if(neuKins==0) return 0;
00339
00340 return dynamic_cast<const REROOT_NeuKin*>((*neuKins)[0]);
00341 }
|
|
|
This function serves to dump the blinding info to an ntuple file for doing research on how well the blinding function is doing. Reimplemented from JobCModule. Definition at line 262 of file Blinder.cxx. References Blinder::BlindingFunction::EventBlindProb(), fBlindingFunction, fFile, fRandom, fTree, GetInfo(), GetTruth(), REROOT_NeuKin::IAction(), REROOT_NeuKin::INu(), REROOT_NeuKin::P4El1(), REROOT_NeuKin::P4El2(), REROOT_NeuKin::P4Mu1(), REROOT_NeuKin::P4Mu2(), REROOT_NeuKin::P4Neu(), REROOT_NeuKin::P4Shw(), REROOT_NeuKin::P4Tau(), and run(). 00263 {
00268 Int_t run;
00269 Int_t snarl;
00270 Int_t epoch;
00271 Float_t length;
00272 Float_t energy;
00273
00274 if(GetInfo(mom,run,epoch,snarl,length,energy)) {
00275 float prob = fBlindingFunction.EventBlindProb(length,energy);
00276
00277 fRandom->SetSeed(run*1000 + snarl);
00278 Float_t r = fRandom->Rndm();
00279
00280 float blind = 0;
00281 if(r<=prob) blind = 1;
00282
00283 float enu=0;
00284 float evis=0;
00285 float cc=0;
00286 float emu=0;
00287 float esh=0;
00288 float nu=0;
00289
00290 const REROOT_NeuKin* truth = GetTruth(mom);
00291 if(truth) {
00292 enu = fabs(truth->P4Neu()[3]);
00293 cc = truth->IAction();
00294 emu = fabs(truth->P4Mu1()[3]) + fabs(truth->P4Mu2()[3]);
00295 esh = fabs(truth->P4Shw()[3]);
00296 evis = emu + esh
00297 + fabs(truth->P4El1()[3])
00298 + fabs(truth->P4El2()[3])
00299 + fabs(truth->P4Tau()[3]);
00300 nu = truth->INu();
00301 }
00302
00303 if(!fTree) {
00304 fFile = new TFile("blind.root","RECREATE");
00305 fTree = new TNtuple("blind","blind",
00306 "run:snarl:len:energy"
00307 ":enu:evis:cc:nu:emu:esh"
00308 ":prob:blind"
00309 );
00310 }
00311
00312 if(fTree) {
00313 fTree->Fill(run,snarl,length,energy
00314 ,enu,evis,cc,nu,emu,esh
00315 ,prob,blind
00316 );
00317 }
00318 };
00319
00320 return JobCResult::kPassed; // kNoDecision, kFailed, etc.
00321 }
|
|
|
|
|
|
Definition at line 71 of file Blinder.h. Referenced by Blinder(), Reco(), and ~Blinder(). |
|
|
|
|
|
|
|
|
|
|
|
Definition at line 72 of file Blinder.h. Referenced by Blinder(), Reco(), and ~Blinder(). |
1.3.9.1