loon: NuFluxChain Class Reference

00143 {
00144   TString fFirstFileName;
00145   //check if there is actually a wildcard
00146   if (!files.MaybeWildcard()) {
00147     fFirstFileName = files;
00148   }
00149   else {
00150     //this was ripped off from ROOT's TChain::AddFile
00151     //Ripped off again from NuBase
00152     
00153     vector<string> fileList;
00154     
00155     //wildcarding used in name
00156     TString basename = files;
00157     
00158     Int_t dotslashpos = basename.Index(".root/");
00159     TString behind_dot_root;
00160     if (dotslashpos>=0) {
00161       // Copy the tree name specification
00162       behind_dot_root = basename(dotslashpos+6,basename.Length()-dotslashpos+6);
00163       // and remove it from basename
00164       basename.Remove(dotslashpos+5);
00165     }
00166     
00167     Int_t slashpos = basename.Last('/');
00168     TString directory;
00169     if (slashpos>=0) {
00170       directory = basename(0,slashpos); // Copy the directory name
00171       basename.Remove(0,slashpos+1);      // and remove it from basename
00172     } else {
00173       directory = gSystem->WorkingDirectory();
00174     }
00175     
00176     const char *file;
00177     void *dir = gSystem->OpenDirectory(gSystem->ExpandPathName(directory.Data()));
00178     
00179     if (dir) {
00180       //create a TList to store the file names (not yet sorted)
00181       TList l;
00182       TRegexp re(basename,kTRUE);
00183       while ((file = gSystem->GetDirEntry(dir))) {
00184         if (!strcmp(file,".") || !strcmp(file,"..")) continue;
00185         TString s = file;
00186         if ( (basename!=file) && s.Index(re) == kNPOS) continue;
00187         l.Add(new TObjString(file));
00188       }
00189       gSystem->FreeDirectory(dir);
00190       //sort the files in alphanumeric order
00191       l.Sort();
00192       TIter next(&l);
00193       TObjString *obj;
00194       while ((obj = (TObjString*)next())) {
00195         file = obj->GetName();
00196         if (behind_dot_root.Length()!=0){
00197           string fileName=Form("%s/%s/%s",directory.Data(),
00198                                file,behind_dot_root.Data());
00199           fileList.push_back(fileName);
00200           //nf += AddFile(Form("%s/%s/%s",directory.Data(),file,behind_dot_root.Data()),kBigNumber);
00201         }
00202         else {
00203           string fileName=Form("%s/%s",directory.Data(),file);
00204           fileList.push_back(fileName);
00205           //nf += AddFile(Form("%s/%s",directory.Data(),file),kBigNumber);
00206         }
00207       }
00208       l.Delete();
00209     }
00210     
00211     //check if any files were found
00212     if (fileList.begin()!=fileList.end()){
00213       cout<<"Used wildcard expansion to find first file name="
00214       <<*fileList.begin()<<endl;
00215       fFirstFileName = *fileList.begin();
00216     }
00217     else {
00218       return;//null string
00219     }
00220   }
00221   
00222   if (fFirstFileName.IsNull()) return;
00223   
00224   //filename: fluka05_le010z185i_extra flugg_le010z185i_run1_extra
00225     
00226   if (fFirstFileName.Contains("le010z-185i"))      fBeamType = BeamType::kL010z185i_rev;
00227   else if (fFirstFileName.Contains("le010z000i"))  fBeamType = BeamType::kL010z000i;
00228   else if (fFirstFileName.Contains("le010z170i"))  fBeamType = BeamType::kL010z170i;
00229   else if (fFirstFileName.Contains("le010z185i"))  fBeamType = BeamType::kL010z185i;
00230   else if (fFirstFileName.Contains("le010z200i"))  fBeamType = BeamType::kL010z200i;
00231   else if (fFirstFileName.Contains("le100z200i"))  fBeamType = BeamType::kL100z200i;
00232   else if (fFirstFileName.Contains("le150z200i"))  fBeamType = BeamType::kL150z200i;
00233   else if (fFirstFileName.Contains("le250z200i"))  fBeamType = BeamType::kL250z200i;
00234 
00235   if (fFirstFileName.Contains("run1")) fRunPeriod = SKZPWeightCalculator::kRunI;
00236   if (fFirstFileName.Contains("run2")) fRunPeriod = SKZPWeightCalculator::kRunII;
00237   if (fFirstFileName.Contains("run3")) fRunPeriod = SKZPWeightCalculator::kRunIII;
00238   if (fFirstFileName.Contains("run4")) fRunPeriod = SKZPWeightCalculator::kRunIV;
00239   
00240   // Correct Run 4 beam type
00241   if (fBeamType == BeamType::kL010z185i && fRunPeriod == SKZPWeightCalculator::kRunIV) {
00242     fBeamType = BeamType::kL010z185i_rev;
00243   }
00244 }


Public Member Functions
	NuFluxChain ()
virtual	~NuFluxChain ()
virtual void	LoadEvent (const Int_t event) const
virtual const Int_t	GetEntries () const
virtual const Double_t	MuonParentPX (const Int_t event) const =0
virtual const Double_t	MuonParentPY (const Int_t event) const =0
virtual const Double_t	MuonParentPZ (const Int_t event) const =0
virtual const Double_t	MuonParentEnergy (const Int_t event) const =0
virtual const Double_t	NeutrinoFDEnergy (const Int_t event) const =0
virtual const Double_t	NeutrinoNDEnergy (const Int_t event) const =0
virtual const Double_t	NuECofM (const Int_t event) const =0
virtual const NuParticle::NuParticleType_t	NuType (const Int_t event) const =0
virtual const Double_t	NuImportanceWeight (const Int_t event) const =0
virtual const Double_t	NuFDWeight (const Int_t event) const =0
virtual const Int_t	NuParticleAsPDGCode (const NuParticle::NuParticleType_t nuParticle) const
virtual const Double_t	NuNDWeight (const Int_t event) const =0
virtual const Double_t	ParentDecayVtxX (const Int_t event) const =0
virtual const Double_t	ParentDecayVtxY (const Int_t event) const =0
virtual const Double_t	ParentDecayVtxZ (const Int_t event) const =0
virtual const Double_t	ParentMass (const Int_t event) const =0
virtual const NuParticle::NuParticleType_t	ParentParticleType (const Int_t event) const =0
virtual const Double_t	ParentPX (const Int_t event) const =0
virtual const Double_t	ParentPY (const Int_t event) const =0
virtual const Double_t	ParentPZ (const Int_t event) const =0
virtual const Double_t	ParentProdDXDZ (const Int_t event) const =0
virtual const Double_t	ParentProdDYDZ (const Int_t event) const =0
virtual const Double_t	ParentProdEnergy (const Int_t event) const =0
virtual const Double_t	ParentProdPZ (const Int_t event) const =0
virtual const Double_t	ParentProdVtxX (const Int_t event) const =0
virtual const Double_t	ParentProdVtxY (const Int_t event) const =0
virtual const Double_t	ParentProdVtxZ (const Int_t event) const =0
virtual const NuParticle::NuParticleType_t	ParentType (const Int_t event) const =0
virtual const Double_t	TargetParentPX (const Int_t event) const =0
virtual const Double_t	TargetParentPY (const Int_t event) const =0
virtual const Double_t	TargetParentPZ (const Int_t event) const =0
virtual const BeamType::BeamType_t	GetBeamType () const
virtual const SKZPWeightCalculator::RunPeriod_t	GetRunPeriod () const
Protected Member Functions
virtual void	ParseFileName (TString files)
	ClassDef (NuFluxChain, 0)
Protected Attributes
TChain	fchain
Int_t	fcurrEntry
Int_t	fNuType
Int_t	fparPartType
Int_t	fparType
BeamType::BeamType_t	fBeamType
SKZPWeightCalculator::RunPeriod_t	fRunPeriod

NuFluxChain Class Reference

Public Member Functions

Protected Member Functions

Protected Attributes

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation